27 #define fabs_( a ) ((( a ) < 0 ) ? -( a ):( a )) 35 #define fmax_( a,b ) ( ( a ) < ( b ) ? ( b ) : ( a ) ) 43 #define fmin_( a,b ) ( ( a ) > ( b ) ? ( b ) : ( a ) ) 51 #define maximize_( maxval, val ) { if (( maxval ) < ( val )) ( maxval )= ( val ); } 59 #define minimize_( minval, val ) { if (( minval ) > ( val )) ( minval )= ( val ); } 69 #define det2_( a1,a2,b1,b2 ) (( a1 )*( b2 ) - ( a2 )*( b1 )) 80 #define det3_( a1,a2,a3,b1,b2,b3,c1,c2,c3 ) ( ( a1 )*det2_( b2,b3,c2,c3 ) \ 81 - ( b1 )*det2_( a2,a3,c2,c3 ) + ( c1 )*det2_( a2,a3,b2,b3 ) ) 94 #define dX( p1,p2 ) ( *( rows[p1] ) - *( rows[p2] )) 95 #define dY( p1,p2 ) ( *( rows[p1]+1 ) - *( rows[p2]+1 )) 96 #define dZ( p1,p2 ) ( *( rows[p1]+2 ) - *( rows[p2]+2 )) 97 #define dW( p1,p2 ) ( *( rows[p1]+3 ) - *( rows[p2]+3 )) 109 boolT bestoutside,
boolT *isoutside,
int *numpart);
160 int numpoints,
int dim,
realT *newpoints,
int newdim);
boolT qh_sethalfspace(qhT *qh, int dim, coordT *coords, coordT **nextp, coordT *normal, coordT *offset, coordT *feasible)
void qh_detroundoff(qhT *qh)
facetT * qh_findbesthorizon(qhT *qh, boolT ischeckmax, pointT *point, facetT *startfacet, boolT noupper, realT *bestdist, int *numpart)
void qh_projectinput(qhT *qh)
void qh_rotateinput(qhT *qh, realT **rows)
facetT * qh_findgooddist(qhT *qh, pointT *point, facetT *facetA, realT *distp, facetT **facetlist)
boolT qh_gram_schmidt(qhT *qh, int dim, realT **rows)
pointT * qh_facetcenter(qhT *qh, setT *vertices)
void qh_projectpoints(qhT *qh, signed char *project, int n, realT *points, int numpoints, int dim, realT *newpoints, int newdim)
void qh_joggleinput(qhT *qh)
facetT * qh_findbest(qhT *qh, pointT *point, facetT *startfacet, boolT bestoutside, boolT isnewfacets, boolT noupper, realT *dist, boolT *isoutside, int *numpart)
pointT * qh_getcenter(qhT *qh, setT *vertices)
realT qh_determinant(qhT *qh, realT **rows, int dim, boolT *nearzero)
void qh_scaleinput(qhT *qh)
realT qh_minabsval(realT *normal, int dim)
realT qh_getangle(qhT *qh, pointT *vect1, pointT *vect2)
facetT * qh_findbestnew(qhT *qh, pointT *point, facetT *startfacet, realT *dist, boolT bestoutside, boolT *isoutside, int *numpart)
void qh_sethyperplane_gauss(qhT *qh, int dim, coordT **rows, pointT *point0, boolT toporient, coordT *normal, coordT *offset, boolT *nearzero)
void qh_sethyperplane_det(qhT *qh, int dim, coordT **rows, coordT *point0, boolT toporient, coordT *normal, realT *offset, boolT *nearzero)
void qh_normalize(qhT *qh, coordT *normal, int dim, boolT toporient)
void qh_printmatrix(qhT *qh, FILE *fp, const char *string, realT **rows, int numrow, int numcol)
boolT qh_inthresholds(qhT *qh, coordT *normal, realT *angle)
realT * qh_maxabsval(realT *normal, int dim)
void qh_printpoints(qhT *qh, FILE *fp, const char *string, setT *points)
void qh_setfacetplane(qhT *qh, facetT *newfacets)
realT qh_distround(qhT *qh, int dimension, realT maxabs, realT maxsumabs)
boolT qh_orientoutside(qhT *qh, facetT *facet)
coordT * qh_sethalfspace_all(qhT *qh, int dim, int count, coordT *halfspaces, pointT *feasible)
realT qh_getdistance(qhT *qh, facetT *facet, facetT *neighbor, realT *mindist, realT *maxdist)
void qh_normalize2(qhT *qh, coordT *normal, int dim, boolT toporient, realT *minnorm, boolT *ismin)
realT qh_facetarea(qhT *qh, facetT *facet)
realT qh_detsimplex(qhT *qh, pointT *apex, setT *points, int dim, boolT *nearzero)
void qh_scalepoints(qhT *qh, pointT *points, int numpoints, int dim, realT *newlows, realT *newhighs)
pointT * qh_projectpoint(qhT *qh, pointT *point, facetT *facet, realT dist)
boolT qh_sharpnewfacets(qhT *qh)
void qh_outerinner(qhT *qh, facetT *facet, realT *outerplane, realT *innerplane)
coordT * qh_copypoints(qhT *qh, coordT *points, int numpoints, int dimension)
void qh_maxsimplex(qhT *qh, int dim, setT *maxpoints, pointT *points, int numpoints, setT **simplex)
void qh_rotatepoints(qhT *qh, realT *points, int numpoints, int dim, realT **rows)
pointT * qh_voronoi_center(qhT *qh, int dim, setT *points)
int qh_mindiff(realT *vecA, realT *vecB, int dim)
pointT * qh_getcentrum(qhT *qh, facetT *facet)
void qh_gausselim(qhT *qh, realT **rows, int numrow, int numcol, boolT *sign, boolT *nearzero)
void qh_crossproduct(int dim, realT vecA[3], realT vecB[3], realT vecC[3])
void qh_backnormal(qhT *qh, realT **rows, int numrow, int numcol, boolT sign, coordT *normal, boolT *nearzero)
realT qh_distnorm(int dim, pointT *point, pointT *normal, realT *offsetp)
void qh_scalelast(qhT *qh, coordT *points, int numpoints, int dim, coordT low, coordT high, coordT newhigh)
coordT qh_pointdist(pointT *point1, pointT *point2, int dim)
realT qh_detjoggle(qhT *qh, pointT *points, int numpoints, int dimension)
setT * qh_maxmin(qhT *qh, pointT *points, int numpoints, int dimension)
realT qh_divzero(realT numer, realT denom, realT mindenom1, boolT *zerodiv)
realT qh_facetarea_simplex(qhT *qh, int dim, coordT *apex, setT *vertices, vertexT *notvertex, boolT toporient, coordT *normal, realT *offset)
void qh_getarea(qhT *qh, facetT *facetlist)
void qh_distplane(qhT *qh, pointT *point, facetT *facet, realT *dist)
realT qh_maxouter(qhT *qh)