/misc/nistmeshlab/meshlab/src/external/qhull-2003.1/src/global.c
C | 1621 lines | 1328 code | 33 blank | 260 comment | 414 complexity | bfe3ed74deac8fa26dcfd78f082e649e MD5 | raw file
- /*<html><pre> -<a href="qh-globa.htm"
- >-------------------------------</a><a name="TOP">-</a>
- global.c
- initializes all the globals of the qhull application
- see README
- see qhull.h for qh.globals and function prototypes
- see qhull_a.h for internal functions
- copyright (c) 1993-2003, The Geometry Center
- */
- #include "qhull_a.h"
- /*========= qh definition (see qhull.h) =======================*/
- #if qh_QHpointer
- qhT *qh_qh= NULL; /* pointer to all global variables */
- #else
- qhT qh_qh; /* all global variables.
- Add "= {0}" if this causes a compiler error.
- Also qh_qhstat in stat.c and qhmem in mem.c. */
- #endif
- /*-<a href ="qh-globa.htm#TOC"
- >--------------------------------</a><a name="qh_version">-</a>
- qh_version
- version string by year and date
- the revision increases on code changes only
- notes:
- change date: Changes.txt, Announce.txt, README.txt,
- qhull.man, qhull.txt, qhull-news.html, Eudora signatures,
- change version: README.txt, qh-get.htm, File_id.diz, Makefile.txt
- change year: Copying.txt
- check download size
- recompile user_eg.c, rbox.c, qhull.c, qconvex.c, qdelaun.c qvoronoi.c, qhalf.c
- */
- char *qh_version = "2003.1 2003/12/30";
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="appendprint">-</a>
- qh_appendprint( printFormat )
- append printFormat to qh.PRINTout unless already defined
- */
- void qh_appendprint (qh_PRINT format) {
- int i;
- for (i=0; i < qh_PRINTEND; i++) {
- if (qh PRINTout[i] == format && format != qh_PRINTqhull)
- break;
- if (!qh PRINTout[i]) {
- qh PRINTout[i]= format;
- break;
- }
- }
- } /* appendprint */
-
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="checkflags">-</a>
-
- qh_checkflags( commandStr, hiddenFlags )
- errors if commandStr contains hiddenFlags
- hiddenFlags starts and ends with a space and is space deliminated (checked)
- notes:
- ignores first word (e.g., "qconvex i")
- use qh_strtol/strtod since strtol/strtod may or may not skip trailing spaces
-
- see:
- qh_initflags() initializes Qhull according to commandStr
- */
- void qh_checkflags(char *command, char *hiddenflags) {
- char *s= command, *t, *chkerr, key, opt, prevopt;
- char chkkey[]= " ";
- char chkopt[]= " ";
- char chkopt2[]= " ";
- if (*hiddenflags != ' ' || hiddenflags[strlen(hiddenflags)-1] != ' ') {
- fprintf(qh ferr, "qhull error (qh_checkflags): hiddenflags must start and end with a space: \"%s\"", hiddenflags);
- qh_errexit(qh_ERRinput, NULL, NULL);
- }
- if (strpbrk(hiddenflags, ",\n\r\t")) {
- fprintf(qh ferr, "qhull error (qh_checkflags): hiddenflags contains commas, newlines, or tabs: \"%s\"", hiddenflags);
- qh_errexit(qh_ERRinput, NULL, NULL);
- }
- while (*s && !isspace(*s)) /* skip program name */
- s++;
- while (*s) {
- while (*s && isspace(*s))
- s++;
- if (*s == '-')
- s++;
- if (!*s)
- break;
- key = *s++;
- chkerr = NULL;
- if (key == '\'') { /* TO 'file name' */
- t= strchr(s, '\'');
- if (!t) {
- fprintf(qh ferr, "qhull error (qh_checkflags): missing the 2nd single-quote for:\n%s\n", s-1);
- qh_errexit(qh_ERRinput, NULL, NULL);
- }
- s= t+1;
- continue;
- }
- chkkey[1]= key;
- if (strstr(hiddenflags, chkkey)) {
- chkerr= chkkey;
- }else if (isupper(key)) {
- opt= ' ';
- prevopt= ' ';
- chkopt[1]= key;
- chkopt2[1]= key;
- while (!chkerr && *s && !isspace(*s)) {
- opt= *s++;
- if (isalpha(opt)) {
- chkopt[2]= opt;
- if (strstr(hiddenflags, chkopt))
- chkerr= chkopt;
- if (prevopt != ' ') {
- chkopt2[2]= prevopt;
- chkopt2[3]= opt;
- if (strstr(hiddenflags, chkopt2))
- chkerr= chkopt2;
- }
- }else if (key == 'Q' && isdigit(opt) && prevopt != 'b'
- && (prevopt == ' ' || islower(prevopt))) {
- chkopt[2]= opt;
- if (strstr(hiddenflags, chkopt))
- chkerr= chkopt;
- }else {
- qh_strtod (s-1, &t);
- if (s < t)
- s= t;
- }
- prevopt= opt;
- }
- }
- if (chkerr) {
- *chkerr= '\'';
- chkerr[strlen(chkerr)-1]= '\'';
- fprintf(qh ferr, "qhull error: option %s is not used with this program.\n It may be used with qhull.\n", chkerr);
- qh_errexit(qh_ERRinput, NULL, NULL);
- }
- }
- } /* checkflags */
-
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="clock">-</a>
-
- qh_clock()
- return user CPU time in 100ths (qh_SECtick)
- only defined for qh_CLOCKtype == 2
- notes:
- use first value to determine time 0
- from Stevens '92 8.15
- */
- unsigned long qh_clock (void) {
- #if (qh_CLOCKtype == 2)
- struct tms time;
- static long clktck; /* initialized first call */
- double ratio, cpu;
- unsigned long ticks;
- if (!clktck) {
- if ((clktck= sysconf (_SC_CLK_TCK)) < 0) {
- fprintf (qh ferr, "qhull internal error (qh_clock): sysconf() failed. Use qh_CLOCKtype 1 in user.h\n");
- qh_errexit (qh_ERRqhull, NULL, NULL);
- }
- }
- if (times (&time) == -1) {
- fprintf (qh ferr, "qhull internal error (qh_clock): times() failed. Use qh_CLOCKtype 1 in user.h\n");
- qh_errexit (qh_ERRqhull, NULL, NULL);
- }
- ratio= qh_SECticks / (double)clktck;
- ticks= time.tms_utime * ratio;
- return ticks;
- #else
- fprintf (qh ferr, "qhull internal error (qh_clock): use qh_CLOCKtype 2 in user.h\n");
- qh_errexit (qh_ERRqhull, NULL, NULL); /* never returns */
- return 0;
- #endif
- } /* clock */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="freebuffers">-</a>
- qh_freebuffers()
- free up global memory buffers
- notes:
- must match qh_initbuffers()
- */
- void qh_freebuffers (void) {
- trace5((qh ferr, "qh_freebuffers: freeing up global memory buffers\n"));
- /* allocated by qh_initqhull_buffers */
- qh_memfree (qh NEARzero, qh hull_dim * sizeof(realT));
- qh_memfree (qh lower_threshold, (qh input_dim+1) * sizeof(realT));
- qh_memfree (qh upper_threshold, (qh input_dim+1) * sizeof(realT));
- qh_memfree (qh lower_bound, (qh input_dim+1) * sizeof(realT));
- qh_memfree (qh upper_bound, (qh input_dim+1) * sizeof(realT));
- qh_memfree (qh gm_matrix, (qh hull_dim+1) * qh hull_dim * sizeof(coordT));
- qh_memfree (qh gm_row, (qh hull_dim+1) * sizeof(coordT *));
- qh NEARzero= qh lower_threshold= qh upper_threshold= NULL;
- qh lower_bound= qh upper_bound= NULL;
- qh gm_matrix= NULL;
- qh gm_row= NULL;
- qh_setfree (&qh other_points);
- qh_setfree (&qh del_vertices);
- qh_setfree (&qh coplanarset);
- if (qh line) /* allocated by qh_readinput, freed if no error */
- free (qh line);
- if (qh half_space)
- free (qh half_space);
- if (qh temp_malloc)
- free (qh temp_malloc);
- if (qh feasible_point) /* allocated by qh_readfeasible */
- free (qh feasible_point);
- if (qh feasible_string) /* allocated by qh_initflags */
- free (qh feasible_string);
- qh line= qh feasible_string= NULL;
- qh half_space= qh feasible_point= qh temp_malloc= NULL;
- /* usually allocated by qh_readinput */
- if (qh first_point && qh POINTSmalloc) {
- free(qh first_point);
- qh first_point= NULL;
- }
- if (qh input_points && qh input_malloc) { /* set by qh_joggleinput */
- free (qh input_points);
- qh input_points= NULL;
- }
- trace5((qh ferr, "qh_freebuffers: finished\n"));
- } /* freebuffers */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="freebuild">-</a>
- qh_freebuild( allmem )
- free global memory used by qh_initbuild and qh_buildhull
- if !allmem,
- does not free short memory (freed by qh_memfreeshort)
- design:
- free centrums
- free each vertex
- mark unattached ridges
- for each facet
- free ridges
- free outside set, coplanar set, neighbor set, ridge set, vertex set
- free facet
- free hash table
- free interior point
- free merge set
- free temporary sets
- */
- void qh_freebuild (boolT allmem) {
- facetT *facet;
- vertexT *vertex;
- ridgeT *ridge, **ridgep;
- mergeT *merge, **mergep;
- trace1((qh ferr, "qh_freebuild: free memory from qh_inithull and qh_buildhull\n"));
- if (qh del_vertices)
- qh_settruncate (qh del_vertices, 0);
- if (allmem) {
- qh_clearcenters (qh_ASnone);
- while ((vertex= qh vertex_list)) {
- if (vertex->next)
- qh_delvertex (vertex);
- else {
- qh_memfree (vertex, sizeof(vertexT));
- qh newvertex_list= qh vertex_list= NULL;
- }
- }
- }else if (qh VERTEXneighbors) {
- FORALLvertices
- qh_setfreelong (&(vertex->neighbors));
- }
- qh VERTEXneighbors= False;
- qh GOODclosest= NULL;
- if (allmem) {
- FORALLfacets {
- FOREACHridge_(facet->ridges)
- ridge->seen= False;
- }
- FORALLfacets {
- if (facet->visible) {
- FOREACHridge_(facet->ridges) {
- if (!otherfacet_(ridge, facet)->visible)
- ridge->seen= True; /* an unattached ridge */
- }
- }
- }
- while ((facet= qh facet_list)) {
- FOREACHridge_(facet->ridges) {
- if (ridge->seen) {
- qh_setfree(&(ridge->vertices));
- qh_memfree(ridge, sizeof(ridgeT));
- }else
- ridge->seen= True;
- }
- qh_setfree (&(facet->outsideset));
- qh_setfree (&(facet->coplanarset));
- qh_setfree (&(facet->neighbors));
- qh_setfree (&(facet->ridges));
- qh_setfree (&(facet->vertices));
- if (facet->next)
- qh_delfacet (facet);
- else {
- qh_memfree (facet, sizeof(facetT));
- qh visible_list= qh newfacet_list= qh facet_list= NULL;
- }
- }
- }else {
- FORALLfacets {
- qh_setfreelong (&(facet->outsideset));
- qh_setfreelong (&(facet->coplanarset));
- if (!facet->simplicial) {
- qh_setfreelong (&(facet->neighbors));
- qh_setfreelong (&(facet->ridges));
- qh_setfreelong (&(facet->vertices));
- }
- }
- }
- qh_setfree (&(qh hash_table));
- qh_memfree (qh interior_point, qh normal_size);
- qh interior_point= NULL;
- FOREACHmerge_(qh facet_mergeset) /* usually empty */
- qh_memfree (merge, sizeof(mergeT));
- qh facet_mergeset= NULL; /* temp set */
- qh degen_mergeset= NULL; /* temp set */
- qh_settempfree_all();
- } /* freebuild */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="freeqhull">-</a>
- qh_freeqhull( allmem )
- free global memory
- if !allmem,
- does not free short memory (freed by qh_memfreeshort)
- notes:
- sets qh.NOerrexit in case caller forgets to
- design:
- free global and temporary memory from qh_initbuild and qh_buildhull
- free buffers
- free statistics
- */
- void qh_freeqhull (boolT allmem) {
- trace1((qh ferr, "qh_freeqhull: free global memory\n"));
- qh NOerrexit= True; /* no more setjmp since called at exit */
- qh_freebuild (allmem);
- qh_freebuffers();
- qh_freestatistics();
- #if qh_QHpointer
- free (qh_qh);
- qh_qh= NULL;
- #else
- memset((char *)&qh_qh, 0, sizeof(qhT));
- qh NOerrexit= True;
- #endif
- } /* freeqhull */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="init_A">-</a>
- qh_init_A( infile, outfile, errfile, argc, argv )
- initialize memory and stdio files
- convert input options to option string (qh.qhull_command)
- notes:
- infile may be NULL if qh_readpoints() is not called
- errfile should always be defined. It is used for reporting
- errors. outfile is used for output and format options.
- argc/argv may be 0/NULL
- called before error handling initialized
- qh_errexit() may not be used
- */
- void qh_init_A (FILE *infile, FILE *outfile, FILE *errfile, int argc, char *argv[]) {
- qh_meminit (errfile);
- qh_initqhull_start (infile, outfile, errfile);
- qh_init_qhull_command (argc, argv);
- } /* init_A */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="init_B">-</a>
- qh_init_B( points, numpoints, dim, ismalloc )
- initialize globals for points array
- points has numpoints dim-dimensional points
- points[0] is the first coordinate of the first point
- points[1] is the second coordinate of the first point
- points[dim] is the first coordinate of the second point
- ismalloc=True
- Qhull will call free(points) on exit or input transformation
- ismalloc=False
- Qhull will allocate a new point array if needed for input transformation
- qh.qhull_command
- is the option string.
- It is defined by qh_init_B(), qh_qhull_command(), or qh_initflags
- returns:
- if qh.PROJECTinput or (qh.DELAUNAY and qh.PROJECTdelaunay)
- projects the input to a new point array
- if qh.DELAUNAY,
- qh.hull_dim is increased by one
- if qh.ATinfinity,
- qh_projectinput adds point-at-infinity for Delaunay tri.
- if qh.SCALEinput
- changes the upper and lower bounds of the input, see qh_scaleinput()
- if qh.ROTATEinput
- rotates the input by a random rotation, see qh_rotateinput()
- if qh.DELAUNAY
- rotates about the last coordinate
- notes:
- called after points are defined
- qh_errexit() may be used
- */
- void qh_init_B (coordT *points, int numpoints, int dim, boolT ismalloc) {
- qh_initqhull_globals (points, numpoints, dim, ismalloc);
- if (qhmem.LASTsize == 0)
- qh_initqhull_mem();
- /* mem.c and qset.c are initialized */
- qh_initqhull_buffers();
- qh_initthresholds (qh qhull_command);
- if (qh PROJECTinput || (qh DELAUNAY && qh PROJECTdelaunay))
- qh_projectinput();
- if (qh SCALEinput)
- qh_scaleinput();
- if (qh ROTATErandom >= 0) {
- qh_randommatrix (qh gm_matrix, qh hull_dim, qh gm_row);
- if (qh DELAUNAY) {
- int k, lastk= qh hull_dim-1;
- for (k= 0; k < lastk; k++) {
- qh gm_row[k][lastk]= 0.0;
- qh gm_row[lastk][k]= 0.0;
- }
- qh gm_row[lastk][lastk]= 1.0;
- }
- qh_gram_schmidt (qh hull_dim, qh gm_row);
- qh_rotateinput (qh gm_row);
- }
- } /* init_B */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="init_qhull_command">-</a>
- qh_init_qhull_command( argc, argv )
- build qh.qhull_command from argc/argv
- returns:
- a space-deliminated string of options (just as typed)
- notes:
- makes option string easy to input and output
- argc/argv may be 0/NULL
- */
- void qh_init_qhull_command(int argc, char *argv[]) {
- int i;
- char *s;
- if (argc) {
- if ((s= strrchr( argv[0], '\\'))) /* Borland gives full path */
- strcpy (qh qhull_command, s+1);
- else
- strcpy (qh qhull_command, argv[0]);
- if ((s= strstr (qh qhull_command, ".EXE"))
- || (s= strstr (qh qhull_command, ".exe")))
- *s= '\0';
- }
- for (i=1; i < argc; i++) {
- if (strlen (qh qhull_command) + strlen(argv[i]) + 1 < sizeof(qh qhull_command)) {
- strcat (qh qhull_command, " ");
- strcat (qh qhull_command, argv[i]);
- }else {
- fprintf (qh ferr, "qhull input error: more than %d characters in command line\n",
- (int)sizeof(qh qhull_command));
- exit (1); /* can not use qh_errexit */
- }
- }
- } /* init_qhull_command */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="initflags">-</a>
- qh_initflags( commandStr )
- set flags and initialized constants from commandStr
- returns:
- sets qh.qhull_command to command if needed
- notes:
- ignores first word (e.g., "qhull d")
- use qh_strtol/strtod since strtol/strtod may or may not skip trailing spaces
- see:
- qh_initthresholds() continues processing of 'Pdn' and 'PDn'
- 'prompt' in unix.c for documentation
- design:
- for each space-deliminated option group
- if top-level option
- check syntax
- append approriate option to option string
- set appropriate global variable or append printFormat to print options
- else
- for each sub-option
- check syntax
- append approriate option to option string
- set appropriate global variable or append printFormat to print options
- */
- void qh_initflags(char *command) {
- int k, i, lastproject;
- char *s= command, *t, *prev_s, *start, key;
- boolT isgeom= False, wasproject;
- realT r;
- if (command != &qh qhull_command[0]) {
- *qh qhull_command= '\0';
- strncat( qh qhull_command, command, sizeof( qh qhull_command));
- }
- while (*s && !isspace(*s)) /* skip program name */
- s++;
- while (*s) {
- while (*s && isspace(*s))
- s++;
- if (*s == '-')
- s++;
- if (!*s)
- break;
- prev_s= s;
- switch (*s++) {
- case 'd':
- qh_option ("delaunay", NULL, NULL);
- qh DELAUNAY= True;
- break;
- case 'f':
- qh_option ("facets", NULL, NULL);
- qh_appendprint (qh_PRINTfacets);
- break;
- case 'i':
- qh_option ("incidence", NULL, NULL);
- qh_appendprint (qh_PRINTincidences);
- break;
- case 'm':
- qh_option ("mathematica", NULL, NULL);
- qh_appendprint (qh_PRINTmathematica);
- break;
- case 'n':
- qh_option ("normals", NULL, NULL);
- qh_appendprint (qh_PRINTnormals);
- break;
- case 'o':
- qh_option ("offFile", NULL, NULL);
- qh_appendprint (qh_PRINToff);
- break;
- case 'p':
- qh_option ("points", NULL, NULL);
- qh_appendprint (qh_PRINTpoints);
- break;
- case 's':
- qh_option ("summary", NULL, NULL);
- qh PRINTsummary= True;
- break;
- case 'v':
- qh_option ("voronoi", NULL, NULL);
- qh VORONOI= True;
- qh DELAUNAY= True;
- break;
- case 'A':
- if (!isdigit(*s) && *s != '.' && *s != '-')
- fprintf(qh ferr, "qhull warning: no maximum cosine angle given for option 'An'. Ignored.\n");
- else {
- if (*s == '-') {
- qh premerge_cos= -qh_strtod (s, &s);
- qh_option ("Angle-premerge-", NULL, &qh premerge_cos);
- qh PREmerge= True;
- }else {
- qh postmerge_cos= qh_strtod (s, &s);
- qh_option ("Angle-postmerge", NULL, &qh postmerge_cos);
- qh POSTmerge= True;
- }
- qh MERGING= True;
- }
- break;
- case 'C':
- if (!isdigit(*s) && *s != '.' && *s != '-')
- fprintf(qh ferr, "qhull warning: no centrum radius given for option 'Cn'. Ignored.\n");
- else {
- if (*s == '-') {
- qh premerge_centrum= -qh_strtod (s, &s);
- qh_option ("Centrum-premerge-", NULL, &qh premerge_centrum);
- qh PREmerge= True;
- }else {
- qh postmerge_centrum= qh_strtod (s, &s);
- qh_option ("Centrum-postmerge", NULL, &qh postmerge_centrum);
- qh POSTmerge= True;
- }
- qh MERGING= True;
- }
- break;
- case 'E':
- if (*s == '-')
- fprintf(qh ferr, "qhull warning: negative maximum roundoff given for option 'An'. Ignored.\n");
- else if (!isdigit(*s))
- fprintf(qh ferr, "qhull warning: no maximum roundoff given for option 'En'. Ignored.\n");
- else {
- qh DISTround= qh_strtod (s, &s);
- qh_option ("Distance-roundoff", NULL, &qh DISTround);
- qh SETroundoff= True;
- }
- break;
- case 'H':
- start= s;
- qh HALFspace= True;
- qh_strtod (s, &t);
- while (t > s) {
- if (*t && !isspace (*t)) {
- if (*t == ',')
- t++;
- else
- fprintf (qh ferr, "qhull warning: origin for Halfspace intersection should be 'Hn,n,n,...'\n");
- }
- s= t;
- qh_strtod (s, &t);
- }
- if (start < t) {
- if (!(qh feasible_string= (char*)calloc (t-start+1, 1))) {
- fprintf(qh ferr, "qhull error: insufficient memory for 'Hn,n,n'\n");
- qh_errexit(qh_ERRmem, NULL, NULL);
- }
- strncpy (qh feasible_string, start, t-start);
- qh_option ("Halfspace-about", NULL, NULL);
- qh_option (qh feasible_string, NULL, NULL);
- }else
- qh_option ("Halfspace", NULL, NULL);
- break;
- case 'R':
- if (!isdigit(*s))
- fprintf(qh ferr, "qhull warning: missing random perturbation for option 'Rn'. Ignored\n");
- else {
- qh RANDOMfactor= qh_strtod (s, &s);
- qh_option ("Random_perturb", NULL, &qh RANDOMfactor);
- qh RANDOMdist= True;
- }
- break;
- case 'V':
- if (!isdigit(*s) && *s != '-')
- fprintf(qh ferr, "qhull warning: missing visible distance for option 'Vn'. Ignored\n");
- else {
- qh MINvisible= qh_strtod (s, &s);
- qh_option ("Visible", NULL, &qh MINvisible);
- }
- break;
- case 'U':
- if (!isdigit(*s) && *s != '-')
- fprintf(qh ferr, "qhull warning: missing coplanar distance for option 'Un'. Ignored\n");
- else {
- qh MAXcoplanar= qh_strtod (s, &s);
- qh_option ("U-coplanar", NULL, &qh MAXcoplanar);
- }
- break;
- case 'W':
- if (*s == '-')
- fprintf(qh ferr, "qhull warning: negative outside width for option 'Wn'. Ignored.\n");
- else if (!isdigit(*s))
- fprintf(qh ferr, "qhull warning: missing outside width for option 'Wn'. Ignored\n");
- else {
- qh MINoutside= qh_strtod (s, &s);
- qh_option ("W-outside", NULL, &qh MINoutside);
- qh APPROXhull= True;
- }
- break;
- /************ sub menus ***************/
- case 'F':
- while (*s && !isspace(*s)) {
- switch(*s++) {
- case 'a':
- qh_option ("Farea", NULL, NULL);
- qh_appendprint (qh_PRINTarea);
- qh GETarea= True;
- break;
- case 'A':
- qh_option ("FArea-total", NULL, NULL);
- qh GETarea= True;
- break;
- case 'c':
- qh_option ("Fcoplanars", NULL, NULL);
- qh_appendprint (qh_PRINTcoplanars);
- break;
- case 'C':
- qh_option ("FCentrums", NULL, NULL);
- qh_appendprint (qh_PRINTcentrums);
- break;
- case 'd':
- qh_option ("Fd-cdd-in", NULL, NULL);
- qh CDDinput= True;
- break;
- case 'D':
- qh_option ("FD-cdd-out", NULL, NULL);
- qh CDDoutput= True;
- break;
- case 'F':
- qh_option ("FFacets-xridge", NULL, NULL);
- qh_appendprint (qh_PRINTfacets_xridge);
- break;
- case 'i':
- qh_option ("Finner", NULL, NULL);
- qh_appendprint (qh_PRINTinner);
- break;
- case 'I':
- qh_option ("FIDs", NULL, NULL);
- qh_appendprint (qh_PRINTids);
- break;
- case 'm':
- qh_option ("Fmerges", NULL, NULL);
- qh_appendprint (qh_PRINTmerges);
- break;
- case 'M':
- qh_option ("FMaple", NULL, NULL);
- qh_appendprint (qh_PRINTmaple);
- break;
- case 'n':
- qh_option ("Fneighbors", NULL, NULL);
- qh_appendprint (qh_PRINTneighbors);
- break;
- case 'N':
- qh_option ("FNeighbors-vertex", NULL, NULL);
- qh_appendprint (qh_PRINTvneighbors);
- break;
- case 'o':
- qh_option ("Fouter", NULL, NULL);
- qh_appendprint (qh_PRINTouter);
- break;
- case 'O':
- if (qh PRINToptions1st) {
- qh_option ("FOptions", NULL, NULL);
- qh_appendprint (qh_PRINToptions);
- }else
- qh PRINToptions1st= True;
- break;
- case 'p':
- qh_option ("Fpoint-intersect", NULL, NULL);
- qh_appendprint (qh_PRINTpointintersect);
- break;
- case 'P':
- qh_option ("FPoint-nearest", NULL, NULL);
- qh_appendprint (qh_PRINTpointnearest);
- break;
- case 'Q':
- qh_option ("FQhull", NULL, NULL);
- qh_appendprint (qh_PRINTqhull);
- break;
- case 's':
- qh_option ("Fsummary", NULL, NULL);
- qh_appendprint (qh_PRINTsummary);
- break;
- case 'S':
- qh_option ("FSize", NULL, NULL);
- qh_appendprint (qh_PRINTsize);
- qh GETarea= True;
- break;
- case 't':
- qh_option ("Ftriangles", NULL, NULL);
- qh_appendprint (qh_PRINTtriangles);
- break;
- case 'v':
- /* option set in qh_initqhull_globals */
- qh_appendprint (qh_PRINTvertices);
- break;
- case 'V':
- qh_option ("FVertex-average", NULL, NULL);
- qh_appendprint (qh_PRINTaverage);
- break;
- case 'x':
- qh_option ("Fxtremes", NULL, NULL);
- qh_appendprint (qh_PRINTextremes);
- break;
- default:
- s--;
- fprintf (qh ferr, "qhull warning: unknown 'F' output option %c, rest ignored\n", (int)s[0]);
- while (*++s && !isspace(*s));
- break;
- }
- }
- break;
- case 'G':
- isgeom= True;
- qh_appendprint (qh_PRINTgeom);
- while (*s && !isspace(*s)) {
- switch(*s++) {
- case 'a':
- qh_option ("Gall-points", NULL, NULL);
- qh PRINTdots= True;
- break;
- case 'c':
- qh_option ("Gcentrums", NULL, NULL);
- qh PRINTcentrums= True;
- break;
- case 'h':
- qh_option ("Gintersections", NULL, NULL);
- qh DOintersections= True;
- break;
- case 'i':
- qh_option ("Ginner", NULL, NULL);
- qh PRINTinner= True;
- break;
- case 'n':
- qh_option ("Gno-planes", NULL, NULL);
- qh PRINTnoplanes= True;
- break;
- case 'o':
- qh_option ("Gouter", NULL, NULL);
- qh PRINTouter= True;
- break;
- case 'p':
- qh_option ("Gpoints", NULL, NULL);
- qh PRINTcoplanar= True;
- break;
- case 'r':
- qh_option ("Gridges", NULL, NULL);
- qh PRINTridges= True;
- break;
- case 't':
- qh_option ("Gtransparent", NULL, NULL);
- qh PRINTtransparent= True;
- break;
- case 'v':
- qh_option ("Gvertices", NULL, NULL);
- qh PRINTspheres= True;
- break;
- case 'D':
- if (!isdigit (*s))
- fprintf (qh ferr, "qhull input error: missing dimension for option 'GDn'\n");
- else {
- if (qh DROPdim >= 0)
- fprintf (qh ferr, "qhull warning: can only drop one dimension. Previous 'GD%d' ignored\n",
- qh DROPdim);
- qh DROPdim= qh_strtol (s, &s);
- qh_option ("GDrop-dim", &qh DROPdim, NULL);
- }
- break;
- default:
- s--;
- fprintf (qh ferr, "qhull warning: unknown 'G' print option %c, rest ignored\n", (int)s[0]);
- while (*++s && !isspace(*s));
- break;
- }
- }
- break;
- case 'P':
- while (*s && !isspace(*s)) {
- switch(*s++) {
- case 'd': case 'D': /* see qh_initthresholds() */
- key= s[-1];
- i= qh_strtol (s, &s);
- r= 0;
- if (*s == ':') {
- s++;
- r= qh_strtod (s, &s);
- }
- if (key == 'd')
- qh_option ("Pdrop-facets-dim-less", &i, &r);
- else
- qh_option ("PDrop-facets-dim-more", &i, &r);
- break;
- case 'g':
- qh_option ("Pgood-facets", NULL, NULL);
- qh PRINTgood= True;
- break;
- case 'G':
- qh_option ("PGood-facet-neighbors", NULL, NULL);
- qh PRINTneighbors= True;
- break;
- case 'o':
- qh_option ("Poutput-forced", NULL, NULL);
- qh FORCEoutput= True;
- break;
- case 'p':
- qh_option ("Pprecision-ignore", NULL, NULL);
- qh PRINTprecision= False;
- break;
- case 'A':
- if (!isdigit (*s))
- fprintf (qh ferr, "qhull input error: missing facet count for keep area option 'PAn'\n");
- else {
- qh KEEParea= qh_strtol (s, &s);
- qh_option ("PArea-keep", &qh KEEParea, NULL);
- qh GETarea= True;
- }
- break;
- case 'F':
- if (!isdigit (*s))
- fprintf (qh ferr, "qhull input error: missing facet area for option 'PFn'\n");
- else {
- qh KEEPminArea= qh_strtod (s, &s);
- qh_option ("PFacet-area-keep", NULL, &qh KEEPminArea);
- qh GETarea= True;
- }
- break;
- case 'M':
- if (!isdigit (*s))
- fprintf (qh ferr, "qhull input error: missing merge count for option 'PMn'\n");
- else {
- qh KEEPmerge= qh_strtol (s, &s);
- qh_option ("PMerge-keep", &qh KEEPmerge, NULL);
- }
- break;
- default:
- s--;
- fprintf (qh ferr, "qhull warning: unknown 'P' print option %c, rest ignored\n", (int)s[0]);
- while (*++s && !isspace(*s));
- break;
- }
- }
- break;
- case 'Q':
- lastproject= -1;
- while (*s && !isspace(*s)) {
- switch(*s++) {
- case 'b': case 'B': /* handled by qh_initthresholds */
- key= s[-1];
- if (key == 'b' && *s == 'B') {
- s++;
- r= qh_DEFAULTbox;
- qh SCALEinput= True;
- qh_option ("QbBound-unit-box", NULL, &r);
- break;
- }
- if (key == 'b' && *s == 'b') {
- s++;
- qh SCALElast= True;
- qh_option ("Qbbound-last", NULL, NULL);
- break;
- }
- k= qh_strtol (s, &s);
- r= 0.0;
- wasproject= False;
- if (*s == ':') {
- s++;
- if ((r= qh_strtod(s, &s)) == 0.0) {
- t= s; /* need true dimension for memory allocation */
- while (*t && !isspace(*t)) {
- if (toupper(*t++) == 'B'
- && k == qh_strtol (t, &t)
- && *t++ == ':'
- && qh_strtod(t, &t) == 0.0) {
- qh PROJECTinput++;
- trace2((qh ferr, "qh_initflags: project dimension %d\n", k));
- qh_option ("Qb-project-dim", &k, NULL);
- wasproject= True;
- lastproject= k;
- break;
- }
- }
- }
- }
- if (!wasproject) {
- if (lastproject == k && r == 0.0)
- lastproject= -1; /* doesn't catch all possible sequences */
- else if (key == 'b') {
- qh SCALEinput= True;
- if (r == 0.0)
- r= -qh_DEFAULTbox;
- qh_option ("Qbound-dim-low", &k, &r);
- }else {
- qh SCALEinput= True;
- if (r == 0.0)
- r= qh_DEFAULTbox;
- qh_option ("QBound-dim-high", &k, &r);
- }
- }
- break;
- case 'c':
- qh_option ("Qcoplanar-keep", NULL, NULL);
- qh KEEPcoplanar= True;
- break;
- case 'f':
- qh_option ("Qfurthest-outside", NULL, NULL);
- qh BESToutside= True;
- break;
- case 'g':
- qh_option ("Qgood-facets-only", NULL, NULL);
- qh ONLYgood= True;
- break;
- case 'i':
- qh_option ("Qinterior-keep", NULL, NULL);
- qh KEEPinside= True;
- break;
- case 'm':
- qh_option ("Qmax-outside-only", NULL, NULL);
- qh ONLYmax= True;
- break;
- case 'r':
- qh_option ("Qrandom-outside", NULL, NULL);
- qh RANDOMoutside= True;
- break;
- case 's':
- qh_option ("Qsearch-initial-simplex", NULL, NULL);
- qh ALLpoints= True;
- break;
- case 't':
- qh_option ("Qtriangulate", NULL, NULL);
- qh TRIangulate= True;
- break;
- case 'T':
- qh_option ("QTestPoints", NULL, NULL);
- if (!isdigit (*s))
- fprintf (qh ferr, "qhull input error: missing number of test points for option 'QTn'\n");
- else {
- qh TESTpoints= qh_strtol (s, &s);
- qh_option ("QTestPoints", &qh TESTpoints, NULL);
- }
- break;
- case 'u':
- qh_option ("QupperDelaunay", NULL, NULL);
- qh UPPERdelaunay= True;
- break;
- case 'v':
- qh_option ("Qvertex-neighbors-convex", NULL, NULL);
- qh TESTvneighbors= True;
- break;
- case 'x':
- qh_option ("Qxact-merge", NULL, NULL);
- qh MERGEexact= True;
- break;
- case 'z':
- qh_option ("Qz-infinity-point", NULL, NULL);
- qh ATinfinity= True;
- break;
- case '0':
- qh_option ("Q0-no-premerge", NULL, NULL);
- qh NOpremerge= True;
- break;
- case '1':
- if (!isdigit(*s)) {
- qh_option ("Q1-no-angle-sort", NULL, NULL);
- qh ANGLEmerge= False;
- break;
- }
- switch(*s++) {
- case '0':
- qh_option ("Q10-no-narrow", NULL, NULL);
- qh NOnarrow= True;
- break;
- case '1':
- qh_option ("Q11-trinormals Qtriangulate", NULL, NULL);
- qh TRInormals= True;
- qh TRIangulate= True;
- break;
- default:
- s--;
- fprintf (qh ferr, "qhull warning: unknown 'Q' qhull option 1%c, rest ignored\n", (int)s[0]);
- while (*++s && !isspace(*s));
- break;
- }
- break;
- case '2':
- qh_option ("Q2-no-merge-independent", NULL, NULL);
- qh MERGEindependent= False;
- goto LABELcheckdigit;
- break; /* no warnings */
- case '3':
- qh_option ("Q3-no-merge-vertices", NULL, NULL);
- qh MERGEvertices= False;
- LABELcheckdigit:
- if (isdigit(*s))
- fprintf (qh ferr, "qhull warning: can not follow '1', '2', or '3' with a digit. '%c' skipped.\n",
- *s++);
- break;
- case '4':
- qh_option ("Q4-avoid-old-into-new", NULL, NULL);
- qh AVOIDold= True;
- break;
- case '5':
- qh_option ("Q5-no-check-outer", NULL, NULL);
- qh SKIPcheckmax= True;
- break;
- case '6':
- qh_option ("Q6-no-concave-merge", NULL, NULL);
- qh SKIPconvex= True;
- break;
- case '7':
- qh_option ("Q7-no-breadth-first", NULL, NULL);
- qh VIRTUALmemory= True;
- break;
- case '8':
- qh_option ("Q8-no-near-inside", NULL, NULL);
- qh NOnearinside= True;
- break;
- case '9':
- qh_option ("Q9-pick-furthest", NULL, NULL);
- qh PICKfurthest= True;
- break;
- case 'G':
- i= qh_strtol (s, &t);
- if (qh GOODpoint)
- fprintf (qh ferr, "qhull warning: good point already defined for option 'QGn'. Ignored\n");
- else if (s == t)
- fprintf (qh ferr, "qhull warning: missing good point id for option 'QGn'. Ignored\n");
- else if (i < 0 || *s == '-') {
- qh GOODpoint= i-1;
- qh_option ("QGood-if-dont-see-point", &i, NULL);
- }else {
- qh GOODpoint= i+1;
- qh_option ("QGood-if-see-point", &i, NULL);
- }
- s= t;
- break;
- case 'J':
- if (!isdigit(*s) && *s != '-')
- qh JOGGLEmax= 0.0;
- else {
- qh JOGGLEmax= (realT) qh_strtod (s, &s);
- qh_option ("QJoggle", NULL, &qh JOGGLEmax);
- }
- break;
- case 'R':
- if (!isdigit(*s) && *s != '-')
- fprintf (qh ferr, "qhull warning: missing random seed for option 'QRn'. Ignored\n");
- else {
- qh ROTATErandom= i= qh_strtol(s, &s);
- if (i > 0)
- qh_option ("QRotate-id", &i, NULL );
- else if (i < -1)
- qh_option ("QRandom-seed", &i, NULL );
- }
- break;
- case 'V':
- i= qh_strtol (s, &t);
- if (qh GOODvertex)
- fprintf (qh ferr, "qhull warning: good vertex already defined for option 'QVn'. Ignored\n");
- else if (s == t)
- fprintf (qh ferr, "qhull warning: no good point id given for option 'QVn'. Ignored\n");
- else if (i < 0) {
- qh GOODvertex= i - 1;
- qh_option ("QV-good-facets-not-point", &i, NULL);
- }else {
- qh_option ("QV-good-facets-point", &i, NULL);
- qh GOODvertex= i + 1;
- }
- s= t;
- break;
- default:
- s--;
- fprintf (qh ferr, "qhull warning: unknown 'Q' qhull option %c, rest ignored\n", (int)s[0]);
- while (*++s && !isspace(*s));
- break;
- }
- }
- break;
- case 'T':
- while (*s && !isspace(*s)) {
- if (isdigit(*s) || *s == '-')
- qh IStracing= qh_strtol(s, &s);
- else switch(*s++) {
- case 'c':
- qh_option ("Tcheck-frequently", NULL, NULL);
- qh CHECKfrequently= True;
- break;
- case 's':
- qh_option ("Tstatistics", NULL, NULL);
- qh PRINTstatistics= True;
- break;
- case 'v':
- qh_option ("Tverify", NULL, NULL);
- qh VERIFYoutput= True;
- break;
- case 'z':
- if (!qh fout)
- fprintf (qh ferr, "qhull warning: output file undefined (stdout). Option 'Tz' ignored.\n");
- else {
- qh_option ("Tz-stdout", NULL, NULL);
- qh ferr= qh fout;
- qhmem.ferr= qh fout;
- }
- break;
- case 'C':
- if (!isdigit(*s))
- fprintf (qh ferr, "qhull warning: missing point id for cone for trace option 'TCn'. Ignored\n");
- else {
- i= qh_strtol (s, &s);
- qh_option ("TCone-stop", &i, NULL);
- qh STOPcone= i + 1;
- }
- break;
- case 'F':
- if (!isdigit(*s))
- fprintf (qh ferr, "qhull warning: missing frequency count for trace option 'TFn'. Ignored\n");
- else {
- qh REPORTfreq= qh_strtol (s, &s);
- qh_option ("TFacet-log", &qh REPORTfreq, NULL);
- qh REPORTfreq2= qh REPORTfreq/2; /* for tracemerging() */
- }
- break;
- case 'I':
- if (s[0] != ' ' || s[1] == '\"' || s[1] == '\'' ||isspace (s[1])) {
- s++;
- fprintf (qh ferr, "qhull warning: option 'TI' mistyped.\nUse 'TI', one space, file name, and space or end-of-line.\nDo not use quotes. Option 'FI' ignored.\n");
- }else { /* not a procedure because of qh_option (filename, NULL, NULL); */
- char filename[500], *t= filename;
- s++;
- while (*s) {
- if (t - filename >= sizeof (filename)-2) {
- fprintf (qh ferr, "qhull error: filename for 'TI' too long.\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- if (isspace (*s))
- break;
- *(t++)= *s++;
- }
- *t= '\0';
- if (!freopen (filename, "r", stdin)) {
- fprintf (qh ferr, "qhull error: could not open file \"%s\".", filename);
- qh_errexit (qh_ERRinput, NULL, NULL);
- }else {
- qh_option ("TInput-file", NULL, NULL);
- qh_option (filename, NULL, NULL);
- }
- }
- break;
- case 'O':
- if (s[0] != ' ' || s[1] == '\"' || isspace (s[1])) {
- s++;
- fprintf (qh ferr, "qhull warning: option 'TO' mistyped.\nUse 'TO', one space, file name, and space or end-of-line.\nThe file name may be enclosed in single quotes.\nDo not use double quotes. Option 'FO' ignored.\n");
- }else { /* not a procedure because of qh_option (filename, NULL, NULL); */
- char filename[500], *t= filename;
- boolT isquote= False;
- s++;
- if (*s == '\'') {
- isquote= True;
- s++;
- }
- while (*s) {
- if (t - filename >= sizeof (filename)-2) {
- fprintf (qh ferr, "qhull error: filename for 'TO' too long.\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- if (isquote) {
- if (*s == '\'') {
- s++;
- isquote= False;
- break;
- }
- }else if (isspace (*s))
- break;
- *(t++)= *s++;
- }
- *t= '\0';
- if (isquote)
- fprintf (qh ferr, "qhull error: missing end quote for option 'TO'. Rest of line ignored.\n");
- else if (!freopen (filename, "w", stdout)) {
- fprintf (qh ferr, "qhull error: could not open file \"%s\".", filename);
- qh_errexit (qh_ERRinput, NULL, NULL);
- }else {
- qh_option ("TOutput-file", NULL, NULL);
- qh_option (filename, NULL, NULL);
- }
- }
- break;
- case 'P':
- if (!isdigit(*s))
- fprintf (qh ferr, "qhull warning: missing point id for trace option 'TPn'. Ignored\n");
- else {
- qh TRACEpoint= qh_strtol (s, &s);
- qh_option ("Trace-point", &qh TRACEpoint, NULL);
- }
- break;
- case 'M':
- if (!isdigit(*s))
- fprintf (qh ferr, "qhull warning: missing merge id for trace option 'TMn'. Ignored\n");
- else {
- qh TRACEmerge= qh_strtol (s, &s);
- qh_option ("Trace-merge", &qh TRACEmerge, NULL);
- }
- break;
- case 'R':
- if (!isdigit(*s))
- fprintf (qh ferr, "qhull warning: missing rerun count for trace option 'TRn'. Ignored\n");
- else {
- qh RERUN= qh_strtol (s, &s);
- qh_option ("TRerun", &qh RERUN, NULL);
- }
- break;
- case 'V':
- i= qh_strtol (s, &t);
- if (s == t)
- fprintf (qh ferr, "qhull warning: missing furthest point id for trace option 'TVn'. Ignored\n");
- else if (i < 0) {
- qh STOPpoint= i - 1;
- qh_option ("TV-stop-before-point", &i, NULL);
- }else {
- qh STOPpoint= i + 1;
- qh_option ("TV-stop-after-point", &i, NULL);
- }
- s= t;
- break;
- case 'W':
- if (!isdigit(*s))
- fprintf (qh ferr, "qhull warning: missing max width for trace option 'TWn'. Ignored\n");
- else {
- qh TRACEdist= (realT) qh_strtod (s, &s);
- qh_option ("TWide-trace", NULL, &qh TRACEdist);
- }
- break;
- default:
- s--;
- fprintf (qh ferr, "qhull warning: unknown 'T' trace option %c, rest ignored\n", (int)s[0]);
- while (*++s && !isspace(*s));
- break;
- }
- }
- break;
- default:
- fprintf (qh ferr, "qhull warning: unknown flag %c (%x)\n", (int)s[-1],
- (int)s[-1]);
- break;
- }
- if (s-1 == prev_s && *s && !isspace(*s)) {
- fprintf (qh ferr, "qhull warning: missing space after flag %c (%x); reserved for menu. Skipped.\n",
- (int)*prev_s, (int)*prev_s);
- while (*s && !isspace(*s))
- s++;
- }
- }
- if (isgeom && !qh FORCEoutput && qh PRINTout[1])
- fprintf (qh ferr, "qhull warning: additional output formats are not compatible with Geomview\n");
- /* set derived values in qh_initqhull_globals */
- } /* initflags */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="initqhull_buffers">-</a>
- qh_initqhull_buffers()
- initialize global memory buffers
- notes:
- must match qh_freebuffers()
- */
- void qh_initqhull_buffers (void) {
- int k;
- qh TEMPsize= (qhmem.LASTsize - sizeof (setT))/SETelemsize;
- if (qh TEMPsize <= 0 || qh TEMPsize > qhmem.LASTsize)
- qh TEMPsize= 8; /* e.g., if qh_NOmem */
- qh other_points= qh_setnew (qh TEMPsize);
- qh del_vertices= qh_setnew (qh TEMPsize);
- qh coplanarset= qh_setnew (qh TEMPsize);
- qh NEARzero= (realT *)qh_memalloc(qh hull_dim * sizeof(realT));
- qh lower_threshold= (realT *)qh_memalloc((qh input_dim+1) * sizeof(realT));
- qh upper_threshold= (realT *)qh_memalloc((qh input_dim+1) * sizeof(realT));
- qh lower_bound= (realT *)qh_memalloc((qh input_dim+1) * sizeof(realT));
- qh upper_bound= (realT *)qh_memalloc((qh input_dim+1) * sizeof(realT));
- for(k= qh input_dim+1; k--; ) {
- qh lower_threshold[k]= -REALmax;
- qh upper_threshold[k]= REALmax;
- qh lower_bound[k]= -REALmax;
- qh upper_bound[k]= REALmax;
- }
- qh gm_matrix= (coordT *)qh_memalloc((qh hull_dim+1) * qh hull_dim * sizeof(coordT));
- qh gm_row= (coordT **)qh_memalloc((qh hull_dim+1) * sizeof(coordT *));
- } /* initqhull_buffers */
- /*-<a href="qh-globa.htm#TOC"
- >-------------------------------</a><a name="initqhull_globals">-</a>
- qh_initqhull_globals( points, numpoints, dim, ismalloc )
- initialize globals
- if ismalloc
- points were malloc'd and qhull should free at end
- returns:
- sets qh.first_point, num_points, input_dim, hull_dim and others
- seeds random number generator (seed=1 if tracing)
- modifies qh.hull_dim if ((qh.DELAUNAY and qh.PROJECTdelaunay) or qh.PROJECTinput)
- adjust user flags as needed
- also checks DIM3 dependencies and constants
- notes:
- do not use qh_point() since an input transformation may move them elsewhere
- see:
- qh_initqhull_start() sets default values for non-zero globals
- design:
- initialize points array from input arguments
- test for qh.ZEROcentrum
- (i.e., use opposite vertex instead of cetrum for convexity testing)
- test for qh.PRINTgood (i.e., only print 'good' facets)
- initialize qh.CENTERtype, qh.normal_size,
- qh.center_size, qh.TRACEpoint/level,
- initialize and test random numbers
- check for conflicting print output options
- */
- void qh_initqhull_globals (coordT *points, int numpoints, int dim, boolT ismalloc) {
- int seed, pointsneeded, extra= 0, i, randi, k;
- boolT printgeom= False, printmath= False, printcoplanar= False;
- realT randr;
- realT factorial;
- time_t timedata;
- trace0((qh ferr, "qh_initqhull_globals: for %s | %s\n", qh rbox_command,
- qh qhull_command));
- qh POINTSmalloc= ismalloc;
- qh first_point= points;
- qh num_points= numpoints;
- qh hull_dim= qh input_dim= dim;
- if (!qh NOpremerge && !qh MERGEexact && !qh PREmerge && qh JOGGLEmax > REALmax/2) {
- qh MERGING= True;
- if (qh hull_dim <= 4) {
- qh PREmerge= True;
- qh_option ("_pre-merge", NULL, NULL);
- }else {
- qh MERGEexact= True;
- qh_option ("Qxact_merge", NULL, NULL);
- }
- }else if (qh MERGEexact)
- qh MERGING= True;
- if (!qh NOpremerge && qh JOGGLEmax > REALmax/2) {
- #ifdef qh_NOmerge
- qh JOGGLEmax= 0.0;
- #endif
- }
- if (qh TRIangulate && qh JOGGLEmax < REALmax/2 && qh PRINTprecision)
- fprintf(qh ferr, "qhull warning: joggle ('QJ') always produces simplicial output. Triangulated output ('Qt') does nothing.\n");
- if (qh JOGGLEmax < REALmax/2 && qh DELAUNAY && !qh SCALEinput && !qh SCALElast) {
- qh SCALElast= True;
- qh_option ("Qbbound-last-qj", NULL, NULL);
- }
- if (qh MERGING && !qh POSTmerge && qh premerge_cos > REALmax/2
- && qh premerge_centrum == 0) {
- qh ZEROcentrum= True;
- qh ZEROall_ok= True;
- qh_option ("_zero-centrum", NULL, NULL);
- }
- if (qh JOGGLEmax < REALmax/2 && REALepsilon > 2e-8 && qh PRINTprecision)
- fprintf(qh ferr, "qhull warning: real epsilon, %2.2g, is probably too large for joggle ('QJn')\nRecompile with double precision reals (see user.h).\n",
- REALepsilon);
- #ifdef qh_NOmerge
- if (qh MERGING) {
- fprintf (qh ferr, "qhull input error: merging not installed (qh_NOmerge + 'Qx', 'Cn' or 'An')\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- #endif
- if (!(qh PRINTgood || qh PRINTneighbors)) {
- if (qh KEEParea || qh KEEPminArea < REALmax/2 || qh KEEPmerge || qh DELAUNAY
- || (!qh ONLYgood && (qh GOODvertex || qh GOODpoint))) {
- qh PRINTgood= True;
- qh_option ("Pgood", NULL, NULL);
- }
- }
- if (qh DELAUNAY && qh KEEPcoplanar && !qh KEEPinside) {
- qh KEEPinside= True;
- qh_option ("Qinterior-keep", NULL, NULL);
- }
- if (qh DELAUNAY && qh HALFspace) {
- fprintf (qh ferr, "qhull input error: can not use Delaunay ('d') or Voronoi ('v') with halfspace intersection ('H')\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- if (!qh DELAUNAY && (qh UPPERdelaunay || qh ATinfinity)) {
- fprintf (qh ferr, "qhull input error: use upper-Delaunay ('Qu') or infinity-point ('Qz') with Delaunay ('d') or Voronoi ('v')\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- if (qh UPPERdelaunay && qh ATinfinity) {
- fprintf (qh ferr, "qhull input error: can not use infinity-point ('Qz') with upper-Delaunay ('Qu')\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- if (qh SCALElast && !qh DELAUNAY && qh PRINTprecision)
- fprintf (qh ferr, "qhull input warning: option 'Qbb' (scale-last-coordinate) is normally used with 'd' or 'v'\n");
- qh DOcheckmax= (!qh SKIPcheckmax && qh MERGING );
- qh KEEPnearinside= (qh DOcheckmax && !(qh KEEPinside && qh KEEPcoplanar)
- && !qh NOnearinside);
- if (qh MERGING)
- qh CENTERtype= qh_AScentrum;
- else if (qh VORONOI)
- qh CENTERtype= qh_ASvoronoi;
- if (qh TESTvneighbors && !qh MERGING) {
- fprintf(qh ferr, "qhull input error: test vertex neighbors ('Qv') needs a merge option\n");
- qh_errexit (qh_ERRinput, NULL ,NULL);
- }
- if (qh PROJECTinput || (qh DELAUNAY && qh PROJECTdelaunay)) {
- qh hull_dim -= qh PROJECTinput;
- if (qh DELAUNAY) {
- qh hull_dim++;
- extra= 1;
- }
- }
- if (qh hull_dim <= 1) {
- fprintf(qh ferr, "qhull error: dimension %d must be > 1\n", qh hull_dim);
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- for (k= 2, factorial=1.0; k < qh hull_dim; k++)
- factorial *= k;
- qh AREAfactor= 1.0 / factorial;
- trace2((qh ferr, "qh_initqhull_globals: initialize globals. dim %d numpoints %d malloc? %d projected %d to hull_dim %d\n",
- dim, numpoints, ismalloc, qh PROJECTinput, qh hull_dim));
- qh normal_size= qh hull_dim * sizeof(coordT);
- qh center_size= qh normal_size - sizeof(coordT);
- pointsneeded= qh hull_dim+1;
- if (qh hull_dim > qh_DIMmergeVertex) {
- qh MERGEvertices= False;
- qh_option ("Q3-no-merge-vertices-dim-high", NULL, NULL);
- }
- if (qh GOODpoint)
- pointsneeded++;
- #ifdef qh_NOtrace
- if (qh IStracing) {
- fprintf (qh ferr, "qhull input error: tracing is not installed (qh_NOtrace in user.h)");
- qh_errexit (qh_ERRqhull, NULL, NULL);
- }
- #endif
- if (qh RERUN > 1) {
- qh TRACElastrun= qh IStracing; /* qh_build_withrestart duplicates next conditional */
- if (qh IStracing != -1)
- qh IStracing= 0;
- }else if (qh TRACEpoint != -1 || qh TRACEdist < REALmax/2 || qh TRACEmerge) {
- qh TRACElevel= (qh IStracing? qh IStracing : 3);
- qh IStracing= 0;
- }
- if (qh ROTATErandom == 0 || qh ROTATErandom == -1) {
- seed= time (&timedata);
- if (qh ROTATErandom == -1) {
- seed= -seed;
- qh_option ("QRandom-seed", &seed, NULL );
- }else
- qh_option ("QRotate-random", &seed, NULL);
- qh ROTATErandom= seed;
- }
- seed= qh ROTATErandom;
- if (seed == INT_MIN) /* default value */
- seed= 1;
- else if (seed < 0)
- seed= -seed;
- qh_RANDOMseed_(seed);
- randr= 0.0;
- for (i= 1000; i--; ) {
- randi= qh_RANDOMint;
- randr += randi;
- if (randi > qh_RANDOMmax) {
- fprintf (qh ferr, "\
- qhull configuration error (qh_RANDOMmax in user.h):\n\
- random integer %d > qh_RANDOMmax (%.8g)\n",
- randi, qh_RANDOMmax);
- qh_errexit (qh_ERRinput, NULL, NULL);
- }
- }
- qh_RANDOMseed_(seed);
- randr = randr/1000;
- if (randr < qh_RANDOMmax/10
- || randr > qh_RANDOMmax * 5)
- fprintf (qh ferr, "\
- qhull configuration warning (qh_RANDOMmax in user.h):\n\
- average of 1000 random integers (%.2g) is much different than expected (%.2g).\n\
- Is qh_RANDOMmax (%.2g) wrong?\n",
- randr, qh_RANDOMmax/2.0, qh_RANDOMmax);
- qh RANDOMa= 2.0 * qh RANDOMfactor/qh_RANDOMmax;
- qh RANDOMb= 1.0 - qh RANDOMfactor;
- if (qh_HASHfactor < 1.1) {
- fprintf(qh ferr, "qhull internal error (qh_initqhull_globals): qh_HASHfactor %d must be at least 1.1. Qhull uses linear hash probing\n",
- qh_HASHfactor);
- qh_errexit (qh_ERRqhull, NULL, NULL);
- }
- if (numpoints+extra < pointsneeded) {
- fprintf(qh ferr,"qhull input error: not enough points (%d) to construct initial simplex (need %d)\n",
- numpoints, pointsneeded);
- qh_errexit(qh_ERRinput, NULL, NULL);
- }
- if (qh PRINTtransparent) {
- if (qh hull_dim != 4 || !qh DELAUNAY || qh VORONOI || qh DROPdim >= 0) {
- fprintf(qh ferr,"qhull input error: transparent Delaunay ('Gt') needs 3-d Delaunay ('d') w/o 'GDn'\n");
- qh_errexit(qh_ERRinput, NULL, NULL);
- }
- qh DROPdim = 3;
- qh PRINTridges = True;
- }
- for (i= qh_PRINTEND; i--; ) {
- if (qh PRINTout[i] == qh_PRINTgeom)
- printgeom= True;
- else if (qh PRINTout[i] == qh_PRINTmathematica || qh PRINTout[i] == qh_PRINTmaple)
- printmath= True;
- else if (qh PRINTout[i] == qh_PRINTcoplanars)
- printcoplanar= True;
- else if (qh PRINTout[i] == qh_PRINTpointnearest)
- printcoplanar= True;
- else if (qh PRINTout[i] == qh_PRINTpointintersect && !qh HALFspace) {
- fprintf (qh ferr, "qhull input error: option 'Fp' is only used for \nhalfspace intersection ('Hn,n,n').\n");
- qh_errexit (qh_ERRinput, NULL, NULL);
- }else if (qh PRINTout[i] == qh_PRINTtriangles && (qh HALFspace || qh VORONOI)) {
- fprintf (qh ferr, "qhull input error: option 'Ft' is not available