runQhull 求取二维数据包络 代码

下面是使用runQhull库求取二维数据包络的C代码示例：

#include <stdio.h>
#include <stdlib.h>
#include "runqhull.h"

#define DIM 2   // 二维数据
#define MAX_POINTS 1000   // 最大点数

int main() {
    int i, j, numpoints, exitcode;
    coordT *points;
    facetT *facet;
    vertexT *vertex;
    setT *vertices;

    // 读取点数据
    scanf("%d", &numpoints);
    points = (coordT*)malloc(numpoints * DIM * sizeof(coordT));
    for (i = 0; i < numpoints; i++) {
        for (j = 0; j < DIM; j++) {
            scanf("%lf", &points[i*DIM+j]);
        }
    }

    // 调用runQhull库求取凸包
    exitcode = qh_new_qhull(DIM, numpoints, points, False, "qhull d Qt", NULL, stderr);
    if (!exitcode) {
        // 输出凸包顶点
        vertices = qh_vertexintersect_new(NULL, qh facet_list, False);
        FORALLvertices {
            printf("%lf %lf\n", vertex->point[0], vertex->point[1]);
        }
        qh_freeqhull(True);
    } else {
        fprintf(stderr, "runqhull error: %d\n", exitcode);
    }

    free(points);
    return 0;
}

其中，需要注意的是：

- 二维数据点的坐标存储在`points`数组中；
- 调用`qh_new_qhull`函数以指定选项字符串来求取凸包；
- 通过`qh_vertexintersect_new`和`FORALLvertices`宏可以遍历凸包的所有顶点；
- 使用`qh_freeqhull`函数释放qhull库分配的内存。