qhstat id[qhstat next++]=name; \
   qhstat doc[name]= string; qhstat count[name]= cnt; qhstat type[name]= stype

qhstat id[qhstat next++]=name; \
   qhstat doc[name]= string; qhstat count[name]= cnt; qhstat type[name]= stype

00047                 {     /* alphabetical after Z/W */
00048     Zacoplanar,
00049     Wacoplanarmax,
00050     Wacoplanartot,
00051     Zangle,
00052     Wangle,
00053     Wanglemax,
00054     Wanglemin,
00055     Zangletests,
00056     Wareatot,
00057     Wareamax,
00058     Wareamin,
00059     Zavoidold,
00060     Wavoidoldmax,
00061     Wavoidoldtot,
00062     Zback0,
00063     Zbestcentrum,
00064     Zbestdist,
00065     Zcentrumtests,
00066     Zcheckpart,
00067     Zcomputefurthest,
00068     Zconcave,
00069     Wconcavemax,
00070     Wconcavetot,
00071     Zconcaveridges,
00072     Zconcaveridge,
00073     Zcoplanar,
00074     Wcoplanarmax,
00075     Wcoplanartot,
00076     Zcoplanarangle,
00077     Zcoplanarcentrum,
00078     Zcoplanarhorizon,
00079     Zcoplanarinside,
00080     Zcoplanarpart,
00081     Zcoplanarridges,
00082     Wcpu,
00083     Zcyclefacetmax,
00084     Zcyclefacettot,
00085     Zcyclehorizon,
00086     Zcyclevertex,
00087     Zdegen,
00088     Wdegenmax,
00089     Wdegentot,
00090     Zdegenvertex,
00091     Zdelfacetdup, 
00092     Zdelridge,
00093     Zdelvertextot,
00094     Zdelvertexmax,
00095     Zdetsimplex,
00096     Zdistcheck,
00097     Zdistconvex,
00098     Zdistgood,
00099     Zdistio,
00100     Zdistplane,
00101     Zdiststat,
00102     Zdistvertex,
00103     Zdistzero,
00104     Zdoc1,
00105     Zdoc2,
00106     Zdoc3,
00107     Zdoc4,
00108     Zdoc5,
00109     Zdoc6,
00110     Zdoc7,
00111     Zdoc8,
00112     Zdoc9,
00113     Zdoc10,
00114     Zdoc11,
00115     Zdropdegen,
00116     Zdropneighbor,
00117     Zdupflip,
00118     Zduplicate,
00119     Wduplicatemax,
00120     Wduplicatetot,
00121     Zdupridge,
00122     Zdupsame,
00123     Zfindfail,
00124     Zflipped, 
00125     Wflippedmax, 
00126     Wflippedtot, 
00127     Zflippedfacets,
00128     Zgauss0,
00129     Zgoodfacet,
00130     Zhashlookup,
00131     Zhashridge,
00132     Zhashridgetest,
00133     Zhashtests,
00134     Zinsidevisible,
00135     Zintersect,
00136     Zintersectfail,
00137     Zintersectmax,
00138     Zintersectnum,
00139     Zintersecttot,
00140     Zmaxneighbors,
00141     Wmaxout,
00142     Wmaxoutside,
00143     Zmaxridges,
00144     Zmaxvertex,
00145     Zmaxvertices,
00146     Zmaxvneighbors,
00147     Zmemfacets,
00148     Zmempoints,
00149     Zmemridges,
00150     Zmemvertices,
00151     Zmergeflipdup,
00152     Zmergehorizon,
00153     Zmergeinittot,
00154     Zmergeinitmax,
00155     Zmergeinittot2,
00156     Zmergeintohorizon,
00157     Zmergenew,
00158     Zmergesettot,
00159     Zmergesetmax,
00160     Zmergesettot2,
00161     Zmergesimplex,
00162     Zmergevertex,
00163     Wmindenom,
00164     Wminvertex,
00165     Zminnorm,
00166     Zmultiridge,
00167     Znearlysingular,
00168     Zneighbor,
00169     Wnewbalance,
00170     Wnewbalance2,
00171     Znewfacettot,
00172     Znewfacetmax,
00173     Znewvertex,
00174     Wnewvertex,
00175     Wnewvertexmax,
00176     Znoarea,
00177     Znonsimplicial,
00178     Znowsimplicial,
00179     Znotgood,
00180     Znotgoodnew,
00181     Znotmax,
00182     Znumfacets,
00183     Znummergemax,
00184     Znummergetot,
00185     Znumneighbors,
00186     Znumridges,
00187     Znumvertices,
00188     Znumvisibility,
00189     Znumvneighbors,
00190     Zonehorizon,
00191     Zpartcoplanar,
00192     Zpartinside,
00193     Zpartition, 
00194     Zpartitionall,
00195     Zpartnear,
00196     Zpbalance,
00197     Wpbalance,
00198     Wpbalance2, 
00199     Zpostfacets, 
00200     Zpremergetot,
00201     Zprocessed,
00202     Zremvertex,
00203     Zremvertexdel,
00204     Zrenameall,
00205     Zrenamepinch,
00206     Zrenameshare,
00207     Zretry,
00208     Wretrymax,
00209     Zridge,
00210     Wridge,
00211     Wridgemax,
00212     Zridge0,
00213     Wridge0,
00214     Wridge0max,
00215     Zridgemid,
00216     Wridgemid,
00217     Wridgemidmax,
00218     Zridgeok,
00219     Wridgeok,
00220     Wridgeokmax,
00221     Zsearchpoints,
00222     Zsetplane,
00223     Ztestvneighbor,
00224     Ztotcheck,
00225     Ztothorizon,
00226     Ztotmerge,
00227     Ztotpartcoplanar,
00228     Ztotpartition,
00229     Ztotridges,
00230     Ztotvertices,
00231     Ztotvisible,
00232     Wvertexmax,
00233     Wvertexmin,
00234     Zvertexridge,
00235     Zvertexridgetot,
00236     Zvertexridgemax,
00237     Zvertices,
00238     Zvisfacettot,
00239     Zvisfacetmax,
00240     Zvisvertextot,
00241     Zvisvertexmax,
00242     Zwidefacet,
00243     Zwidevertices,
00244     ZEND};

00310 {zdoc,zinc,zadd,zmax,zmin,ZTYPEreal,wadd,wmax,wmin,ZTYPEend};

00034                         {
00035   
00036    /* zdef_(type,name,doc,average) */
00037   zzdef_(zdoc, Zdoc2, "precision statistics", -1);
00038   zdef_(zinc, Znewvertex, NULL, -1);
00039   zdef_(wadd, Wnewvertex, "ave. distance of a new vertex to a facet (not 0s)", Znewvertex);
00040   zzdef_(wmax, Wnewvertexmax, "max. distance of a new vertex to a facet", -1);
00041   zdef_(wmax, Wvertexmax, "max. distance of an output vertex to a facet", -1);
00042   zdef_(wmin, Wvertexmin, "min. distance of an output vertex to a facet", -1);
00043   zdef_(wmin, Wmindenom, "min. denominator in hyperplane computation", -1);
00044 
00045   qhstat precision= qhstat next;  /* call qh_precision for each of these */
00046   zzdef_(zdoc, Zdoc3, "precision problems", -1);
00047   zzdef_(zinc, Zcoplanarridges, "coplanar half ridges in output", -1);
00048   zzdef_(zinc, Zconcaveridges, "concave half ridges in output", -1);
00049   zzdef_(zinc, Zflippedfacets, "flipped facets", -1);
00050   zzdef_(zinc, Zcoplanarhorizon, "coplanar horizon facets for new vertices", -1);
00051   zzdef_(zinc, Zcoplanarpart, "coplanar points during partitioning", -1);
00052   zzdef_(zinc, Zminnorm, "degenerate hyperplanes recomputed with gaussian elimination", -1);
00053   zzdef_(zinc, Znearlysingular, "nearly singular or axis-parallel hyperplanes", -1);
00054   zzdef_(zinc, Zback0, "zero divisors during back substitute", -1);
00055   zzdef_(zinc, Zgauss0, "zero divisors during gaussian elimination", -1);
00056   zzdef_(zinc, Zmultiridge, "ridges with multiple neighbors", -1);
00057 }

00058                         {
00059   zzdef_(zdoc, Zdoc1, "summary information", -1);
00060   zdef_(zinc, Zvertices, "number of vertices in output", -1);
00061   zdef_(zinc, Znumfacets, "number of facets in output", -1);
00062   zdef_(zinc, Znonsimplicial, "number of non-simplicial facets in output", -1);
00063   zdef_(zinc, Znowsimplicial, "number of simplicial facets that were merged", -1);
00064   zdef_(zinc, Znumridges, "number of ridges in output", -1);
00065   zdef_(zadd, Znumridges, "average number of ridges per facet", Znumfacets);
00066   zdef_(zmax, Zmaxridges, "maximum number of ridges", -1);
00067   zdef_(zadd, Znumneighbors, "average number of neighbors per facet", Znumfacets);
00068   zdef_(zmax, Zmaxneighbors, "maximum number of neighbors", -1);
00069   zdef_(zadd, Znumvertices, "average number of vertices per facet", Znumfacets);
00070   zdef_(zmax, Zmaxvertices, "maximum number of vertices", -1);
00071   zdef_(zadd, Znumvneighbors, "average number of neighbors per vertex", Zvertices);
00072   zdef_(zmax, Zmaxvneighbors, "maximum number of neighbors", -1);
00073   zdef_(wadd, Wcpu, "cpu seconds for qhull after input", -1);
00074   zdef_(zinc, Ztotvertices, "vertices created altogether", -1);
00075   zzdef_(zinc, Zsetplane, "facets created altogether", -1);
00076   zdef_(zinc, Ztotridges, "ridges created altogether", -1);
00077   zdef_(zinc, Zpostfacets, "facets before post merge", -1);
00078   zdef_(zadd, Znummergetot, "average merges per facet (at most 511)", Znumfacets);
00079   zdef_(zmax, Znummergemax, "  maximum merges for a facet (at most 511)", -1);
00080   zdef_(zinc, Zangle, NULL, -1);
00081   zdef_(wadd, Wangle, "average angle (cosine) of facet normals for all ridges", Zangle);
00082   zdef_(wmax, Wanglemax, "  maximum angle (cosine) of facet normals across a ridge", -1);
00083   zdef_(wmin, Wanglemin, "  minimum angle (cosine) of facet normals across a ridge", -1);
00084   zdef_(wadd, Wareatot, "total area of facets", -1);
00085   zdef_(wmax, Wareamax, "  maximum facet area", -1);
00086   zdef_(wmin, Wareamin, "  minimum facet area", -1);
00087 }  

00088                         {
00089   zdef_(zdoc, Zdoc9, "build hull statistics", -1);
00090   zzdef_(zinc, Zprocessed, "points processed", -1);
00091   zzdef_(zinc, Zretry, "retries due to precision problems", -1);
00092   zdef_(wmax, Wretrymax, "  max. random joggle", -1);
00093   zdef_(zmax, Zmaxvertex, "max. vertices at any one time", -1);
00094   zdef_(zinc, Ztotvisible, "ave. visible facets per iteration", Zprocessed);
00095   zdef_(zinc, Zinsidevisible, "  ave. visible facets without an horizon neighbor", Zprocessed);
00096   zdef_(zadd, Zvisfacettot,  "  ave. facets deleted per iteration", Zprocessed);
00097   zdef_(zmax, Zvisfacetmax,  "    maximum", -1);
00098   zdef_(zadd, Zvisvertextot, "ave. visible vertices per iteration", Zprocessed);
00099   zdef_(zmax, Zvisvertexmax, "    maximum", -1);
00100   zdef_(zadd, Zdelvertextot, "  ave. vertices deleted per iteration", Zprocessed);
00101   zdef_(zmax, Zdelvertexmax, "    maximum vertices deleted", -1);
00102   zdef_(zinc, Ztothorizon, "ave. horizon facets per iteration", Zprocessed);
00103   zdef_(zadd, Znewfacettot,  "ave. new or merged facets per iteration", Zprocessed);
00104   zdef_(zmax, Znewfacetmax,  "    maximum (includes initial simplex)", -1);
00105   zdef_(wadd, Wnewbalance, "average new facet balance", Zprocessed);
00106   zdef_(wadd, Wnewbalance2, "  standard deviation", -1);
00107   zdef_(wadd, Wpbalance, "average partition balance", Zpbalance);
00108   zdef_(wadd, Wpbalance2, "  standard deviation", -1);
00109   zdef_(zinc, Zpbalance, "  number of trials", -1);
00110   zdef_(zinc, Zsearchpoints, "searches of all points for initial simplex", -1);
00111   zdef_(zinc, Zdetsimplex, "determinants computed (area & initial hull)", -1);
00112   zdef_(zinc, Znoarea, "determinants not computed because vertex too low", -1);
00113   zdef_(zinc, Znotmax, "points ignored (not above max_outside)", -1);
00114   zdef_(zinc, Znotgood, "points ignored (not above a good facet)", -1);
00115   zdef_(zinc, Znotgoodnew, "points ignored (didn't create a good new facet)", -1);
00116   zdef_(zinc, Zgoodfacet, "good facets found", -1);
00117   zzdef_(zinc, Znumvisibility, "distance tests for facet visibility", -1);
00118   zdef_(zinc, Zdistvertex, "distance tests to report minimum vertex", -1);
00119   zdef_(zinc, Ztotcheck, "points checked for facets' outer planes", -1);
00120   zdef_(zinc, Zcheckpart, "  ave. distance tests per check", Ztotcheck);
00121 }

00122                        {
00123   zdef_(zdoc, Zdoc4, "partitioning statistics", -1);
00124   zdef_(zinc, Zpartinside, "inside points", -1);
00125   zdef_(zinc, Zpartnear, "  inside points kept with a facet", -1);
00126   zdef_(zinc, Zcoplanarinside, "  inside points that were coplanar with a facet", -1);
00127   zdef_(wadd, Wmaxout, "difference in max_outside at final check", -1);
00128   
00129   zzdef_(zinc, Zpartitionall, "distance tests for initial partition", -1);
00130   zdef_(zinc, Ztotpartition, "partitions of a point", -1);
00131   zzdef_(zinc, Zpartition, "distance tests for partitioning", -1);
00132   zzdef_(zinc, Zdistcheck, "distance tests for checking flipped facets", -1); 
00133   zzdef_(zinc, Zdistconvex, "distance tests for checking convexity", -1); 
00134   zdef_(zinc, Zdistgood, "distance tests for checking good point", -1); 
00135   zdef_(zinc, Zdistio, "distance tests for output", -1); 
00136   zdef_(zinc, Zdiststat, "distance tests for statistics", -1); 
00137   zdef_(zinc, Zdistplane, "total number of distance tests", -1);
00138   zdef_(zinc, Ztotpartcoplanar, "partitions of coplanar points or deleted vertices", -1);
00139   zzdef_(zinc, Zpartcoplanar, "   distance tests for these partitions", -1);
00140   zdef_(zinc, Zcomputefurthest, "distance tests for computing furthest", -1);
00141 }

00142                        {
00143   zdef_(zdoc, Zdoc5, "statistics for matching ridges", -1);
00144   zdef_(zinc, Zhashlookup, "total lookups for matching ridges of new facets", -1);
00145   zdef_(zinc, Zhashtests, "average number of tests to match a ridge", Zhashlookup);
00146   zdef_(zinc, Zhashridge, "total lookups of subridges (duplicates and boundary)", -1);
00147   zdef_(zinc, Zhashridgetest, "average number of tests per subridge", Zhashridge);
00148   zdef_(zinc, Zdupsame, "duplicated ridges in same merge cycle", -1);
00149   zdef_(zinc, Zdupflip, "duplicated ridges with flipped facets", -1);
00150 
00151   zdef_(zdoc, Zdoc6, "statistics for determining merges", -1);
00152   zdef_(zinc, Zangletests, "angles computed for ridge convexity", -1);
00153   zdef_(zinc, Zbestcentrum, "best merges used centrum instead of vertices",-1);
00154   zzdef_(zinc, Zbestdist, "distance tests for best merge", -1);
00155   zzdef_(zinc, Zcentrumtests, "distance tests for centrum convexity", -1);
00156   zzdef_(zinc, Zdistzero, "distance tests for checking simplicial convexity", -1);
00157   zdef_(zinc, Zcoplanarangle, "coplanar angles in getmergeset", -1);
00158   zdef_(zinc, Zcoplanarcentrum, "coplanar centrums in getmergeset", -1);
00159   zdef_(zinc, Zconcaveridge, "concave ridges in getmergeset", -1);
00160 }

00161                        {
00162   zdef_(zdoc, Zdoc7, "statistics for merging", -1);
00163   zdef_(zinc, Zpremergetot, "merge iterations", -1);
00164   zdef_(zadd, Zmergeinittot, "ave. initial non-convex ridges per iteration", Zpremergetot);
00165   zdef_(zadd, Zmergeinitmax, "  maximum", -1);
00166   zdef_(zadd, Zmergesettot, "  ave. additional non-convex ridges per iteration", Zpremergetot);
00167   zdef_(zadd, Zmergesetmax, "  maximum additional in one pass", -1);
00168   zdef_(zadd, Zmergeinittot2, "initial non-convex ridges for post merging", -1);
00169   zdef_(zadd, Zmergesettot2, "  additional non-convex ridges", -1);
00170   zdef_(wmax, Wmaxoutside, "max distance of vertex or coplanar point above facet (w/roundoff)", -1);
00171   zdef_(wmin, Wminvertex, "max distance of merged vertex below facet (or roundoff)", -1);
00172   zdef_(zinc, Zwidefacet, "centrums frozen due to a wide merge", -1);
00173   zdef_(zinc, Zwidevertices, "centrums frozen due to extra vertices", -1);
00174   zzdef_(zinc, Ztotmerge, "total number of facets or cycles of facets merged", -1);
00175   zdef_(zinc, Zmergesimplex, "merged a simplex", -1);
00176   zdef_(zinc, Zonehorizon, "simplices merged into coplanar horizon", -1);
00177   zzdef_(zinc, Zcyclehorizon, "cycles of facets merged into coplanar horizon", -1);
00178   zzdef_(zadd, Zcyclefacettot, "  ave. facets per cycle", Zcyclehorizon);
00179   zdef_(zmax, Zcyclefacetmax, "  max. facets", -1);
00180   zdef_(zinc, Zmergeintohorizon, "new facets merged into horizon", -1);
00181   zdef_(zinc, Zmergenew, "new facets merged", -1);
00182   zdef_(zinc, Zmergehorizon, "horizon facets merged into new facets", -1);
00183   zdef_(zinc, Zmergevertex, "vertices deleted by merging", -1);
00184   zdef_(zinc, Zcyclevertex, "vertices deleted by merging into coplanar horizon", -1);
00185   zdef_(zinc, Zdegenvertex, "vertices deleted by degenerate facet", -1);
00186   zdef_(zinc, Zmergeflipdup, "merges due to flipped facets in duplicated ridge", -1);
00187   zdef_(zinc, Zneighbor, "merges due to redundant neighbors", -1);
00188   zdef_(zadd, Ztestvneighbor, "non-convex vertex neighbors", -1); 
00189 }

00190                        {
00191   zdef_(zinc, Zacoplanar, "merges due to angle coplanar facets", -1);
00192   zdef_(wadd, Wacoplanartot, "  average merge distance", Zacoplanar);
00193   zdef_(wmax, Wacoplanarmax, "  maximum merge distance", -1);
00194   zdef_(zinc, Zcoplanar, "merges due to coplanar facets", -1);
00195   zdef_(wadd, Wcoplanartot, "  average merge distance", Zcoplanar);
00196   zdef_(wmax, Wcoplanarmax, "  maximum merge distance", -1);
00197   zdef_(zinc, Zconcave, "merges due to concave facets", -1);
00198   zdef_(wadd, Wconcavetot, "  average merge distance", Zconcave);
00199   zdef_(wmax, Wconcavemax, "  maximum merge distance", -1);
00200   zdef_(zinc, Zavoidold, "coplanar/concave merges due to avoiding old merge", -1);
00201   zdef_(wadd, Wavoidoldtot, "  average merge distance", Zavoidold);
00202   zdef_(wmax, Wavoidoldmax, "  maximum merge distance", -1);
00203   zdef_(zinc, Zdegen, "merges due to degenerate facets", -1);
00204   zdef_(wadd, Wdegentot, "  average merge distance", Zdegen);
00205   zdef_(wmax, Wdegenmax, "  maximum merge distance", -1);
00206   zdef_(zinc, Zflipped, "merges due to removing flipped facets", -1);
00207   zdef_(wadd, Wflippedtot, "  average merge distance", Zflipped);
00208   zdef_(wmax, Wflippedmax, "  maximum merge distance", -1);
00209   zdef_(zinc, Zduplicate, "merges due to duplicated ridges", -1);
00210   zdef_(wadd, Wduplicatetot, "  average merge distance", Zduplicate);
00211   zdef_(wmax, Wduplicatemax, "  maximum merge distance", -1);
00212 }

00213                        {
00214   zdef_(zdoc, Zdoc8, "renamed vertex statistics", -1);
00215   zdef_(zinc, Zrenameshare, "renamed vertices shared by two facets", -1);
00216   zdef_(zinc, Zrenamepinch, "renamed vertices in a pinched facet", -1);
00217   zdef_(zinc, Zrenameall, "renamed vertices shared by multiple facets", -1);
00218   zdef_(zinc, Zfindfail, "rename failures due to duplicated ridges", -1);
00219   zdef_(zinc, Zdupridge, "  duplicate ridges detected", -1);
00220   zdef_(zinc, Zdelridge, "deleted ridges due to renamed vertices", -1);
00221   zdef_(zinc, Zdropneighbor, "dropped neighbors due to renamed vertices", -1);
00222   zdef_(zinc, Zdropdegen, "degenerate facets due to dropped neighbors", -1);
00223   zdef_(zinc, Zdelfacetdup, "  facets deleted because of no neighbors", -1);
00224   zdef_(zinc, Zremvertex, "vertices removed from facets due to no ridges", -1);
00225   zdef_(zinc, Zremvertexdel, "  deleted", -1);
00226   zdef_(zinc, Zintersectnum, "vertex intersections for locating redundant vertices", -1);
00227   zdef_(zinc, Zintersectfail, "intersections failed to find a redundant vertex", -1);
00228   zdef_(zinc, Zintersect, "intersections found redundant vertices", -1);
00229   zdef_(zadd, Zintersecttot, "   ave. number found per vertex", Zintersect);
00230   zdef_(zmax, Zintersectmax, "   max. found for a vertex", -1);
00231   zdef_(zinc, Zvertexridge, NULL, -1);
00232   zdef_(zadd, Zvertexridgetot, "  ave. number of ridges per tested vertex", Zvertexridge);
00233   zdef_(zmax, Zvertexridgemax, "  max. number of ridges per tested vertex", -1);
00234 
00235   zdef_(zdoc, Zdoc10, "memory usage statistics (in bytes)", -1);
00236   zdef_(zadd, Zmemfacets, "for facets and their normals, neighbor and vertex sets", -1);
00237   zdef_(zadd, Zmemvertices, "for vertices and their neighbor sets", -1);
00238   zdef_(zadd, Zmempoints, "for input points and outside and coplanar sets",-1);
00239   zdef_(zadd, Zmemridges, "for ridges and their vertex sets", -1);
00240 } /* allstat */

00242                        {
00243   qhstat vridges= qhstat next;
00244   zzdef_(zdoc, Zdoc11, "Voronoi ridge statistics", -1);
00245   zzdef_(zinc, Zridge, "non-simplicial Voronoi vertices for all ridges", -1);
00246   zzdef_(wadd, Wridge, "  ave. distance to ridge", Zridge);
00247   zzdef_(wmax, Wridgemax, "  max. distance to ridge", -1);
00248   zzdef_(zinc, Zridgemid, "bounded ridges", -1);
00249   zzdef_(wadd, Wridgemid, "  ave. distance of midpoint to ridge", Zridgemid);
00250   zzdef_(wmax, Wridgemidmax, "  max. distance of midpoint to ridge", -1);
00251   zzdef_(zinc, Zridgeok, "bounded ridges with ok normal", -1);
00252   zzdef_(wadd, Wridgeok, "  ave. angle to ridge", Zridgeok);
00253   zzdef_(wmax, Wridgeokmax, "  max. angle to ridge", -1);
00254   zzdef_(zinc, Zridge0, "bounded ridges with near-zero normal", -1);
00255   zzdef_(wadd, Wridge0, "  ave. angle to ridge", Zridge0);
00256   zzdef_(wmax, Wridge0max, "  max. angle to ridge", -1);
00257 } /* allstat */

00265                              {
00266   int i;
00267   
00268   for (i=ZEND; i--; ) 
00269     qhstat printed[i]= False;
00270 } /* allstatistics */

00280                                  {
00281   facetT *facet, *neighbor, **neighborp;
00282   vertexT *vertex, **vertexp;
00283   realT dotproduct, dist;
00284   int sizneighbors, sizridges, sizvertices, i;
00285   
00286   qh old_randomdist= qh RANDOMdist;
00287   qh RANDOMdist= False;
00288   zval_(Zmempoints)= qh num_points * qh normal_size + 
00289                              sizeof (qhT) + sizeof (qhstatT);
00290   zval_(Zmemfacets)= 0;
00291   zval_(Zmemridges)= 0;
00292   zval_(Zmemvertices)= 0;
00293   zval_(Zangle)= 0;
00294   wval_(Wangle)= 0.0;
00295   zval_(Znumridges)= 0;
00296   zval_(Znumfacets)= 0;
00297   zval_(Znumneighbors)= 0;
00298   zval_(Znumvertices)= 0;
00299   zval_(Znumvneighbors)= 0;
00300   zval_(Znummergetot)= 0;
00301   zval_(Znummergemax)= 0;
00302   zval_(Zvertices)= qh num_vertices - qh_setsize (qh del_vertices);
00303   if (qh MERGING || qh APPROXhull || qh JOGGLEmax < REALmax/2)
00304     wmax_(Wmaxoutside, qh max_outside);
00305   if (qh MERGING)
00306     wmin_(Wminvertex, qh min_vertex);
00307   FORALLfacets
00308     facet->seen= False;
00309   if (qh DELAUNAY) {
00310     FORALLfacets {
00311       if (facet->upperdelaunay != qh UPPERdelaunay)
00312         facet->seen= True; /* remove from angle statistics */
00313     }
00314   }
00315   FORALLfacets {
00316     if (facet->visible && qh NEWfacets)
00317       continue;
00318     sizvertices= qh_setsize (facet->vertices);
00319     sizneighbors= qh_setsize (facet->neighbors);
00320     sizridges= qh_setsize (facet->ridges);
00321     zinc_(Znumfacets);
00322     zadd_(Znumvertices, sizvertices);
00323     zmax_(Zmaxvertices, sizvertices);
00324     zadd_(Znumneighbors, sizneighbors);
00325     zmax_(Zmaxneighbors, sizneighbors);
00326     zadd_(Znummergetot, facet->nummerge);
00327     i= facet->nummerge; /* avoid warnings */
00328     zmax_(Znummergemax, i); 
00329     if (!facet->simplicial) {
00330       if (sizvertices == qh hull_dim) {
00331         zinc_(Znowsimplicial);
00332       }else {
00333         zinc_(Znonsimplicial);
00334       }
00335     }
00336     if (sizridges) {
00337       zadd_(Znumridges, sizridges);
00338       zmax_(Zmaxridges, sizridges);
00339     }
00340     zadd_(Zmemfacets, sizeof (facetT) + qh normal_size + 2*sizeof (setT) 
00341        + SETelemsize * (sizneighbors + sizvertices));
00342     if (facet->ridges) {
00343       zadd_(Zmemridges,
00344          sizeof (setT) + SETelemsize * sizridges + sizridges * 
00345          (sizeof (ridgeT) + sizeof (setT) + SETelemsize * (qh hull_dim-1))/2);
00346     }
00347     if (facet->outsideset)
00348       zadd_(Zmempoints, sizeof (setT) + SETelemsize * qh_setsize (facet->outsideset));
00349     if (facet->coplanarset)
00350       zadd_(Zmempoints, sizeof (setT) + SETelemsize * qh_setsize (facet->coplanarset));
00351     if (facet->seen) /* Delaunay upper envelope */
00352       continue;
00353     facet->seen= True;
00354     FOREACHneighbor_(facet) {
00355       if (neighbor == qh_DUPLICATEridge || neighbor == qh_MERGEridge
00356           || neighbor->seen || !facet->normal || !neighbor->normal)
00357         continue;
00358       dotproduct= qh_getangle(facet->normal, neighbor->normal);
00359       zinc_(Zangle);
00360       wadd_(Wangle, dotproduct);
00361       wmax_(Wanglemax, dotproduct)
00362       wmin_(Wanglemin, dotproduct)
00363     }
00364     if (facet->normal) {
00365       FOREACHvertex_(facet->vertices) {
00366         zinc_(Zdiststat);
00367         qh_distplane(vertex->point, facet, &dist);
00368         wmax_(Wvertexmax, dist);
00369         wmin_(Wvertexmin, dist);
00370       }
00371     }
00372   }
00373   FORALLvertices {
00374     if (vertex->deleted)
00375       continue;
00376     zadd_(Zmemvertices, sizeof (vertexT));
00377     if (vertex->neighbors) {
00378       sizneighbors= qh_setsize (vertex->neighbors);
00379       zadd_(Znumvneighbors, sizneighbors);
00380       zmax_(Zmaxvneighbors, sizneighbors);
00381       zadd_(Zmemvertices, sizeof (vertexT) + SETelemsize * sizneighbors);
00382     }
00383   }
00384   qh RANDOMdist= qh old_randomdist;
00385 } /* collectstatistics */

00394                               {
00395 
00396 #if qh_QHpointer
00397   free (qh_qhstat);
00398   qh_qhstat= NULL;
00399 #endif
00400 } /* freestatistics */

00411                               {
00412   int i;
00413   realT realx;
00414   int intx;
00415 
00416 #if qh_QHpointer
00417   if (!(qh_qhstat= (qhstatT *)malloc (sizeof(qhstatT)))) {
00418     fprintf (qhmem.ferr, "qhull error (qh_initstatistics): insufficient memory\n");
00419     exit (1);  /* can not use qh_errexit() */
00420   }
00421 #endif
00422   
00423   qhstat next= 0;
00424   qh_allstatA();
00425   qh_allstatB();
00426   qh_allstatC();
00427   qh_allstatD();
00428   qh_allstatE();
00429   qh_allstatF();
00430   qh_allstatG();
00431   qh_allstatH();
00432   qh_allstatI();
00433   if (qhstat next > sizeof(qhstat id)) {
00434     fprintf (qhmem.ferr, "qhull error (qh_initstatistics): increase size of qhstat.id[].\n\
00435       qhstat.next %d should be <= sizeof(qhstat id) %d\n", qhstat next, sizeof(qhstat id));
00436 #if 0 /* for locating error, Znumridges should be duplicated */
00437     for (i=0; i < ZEND; i++) {
00438       int j;
00439       for (j=i+1; j < ZEND; j++) {
00440         if (qhstat id[i] == qhstat id[j]) {
00441           fprintf (qhmem.ferr, "qhull error (qh_initstatistics): duplicated statistic %d at indices %d and %d\n", 
00442               qhstat id[i], i, j);
00443         }
00444       }
00445     }
00446 #endif 
00447     exit (1);  /* can not use qh_errexit() */
00448   }
00449   qhstat init[zinc].i= 0;
00450   qhstat init[zadd].i= 0;
00451   qhstat init[zmin].i= INT_MAX;
00452   qhstat init[zmax].i= INT_MIN;
00453   qhstat init[wadd].r= 0;
00454   qhstat init[wmin].r= REALmax;
00455   qhstat init[wmax].r= -REALmax;
00456   for (i=0; i < ZEND; i++) {
00457     if (qhstat type[i] > ZTYPEreal) {
00458       realx= qhstat init[(unsigned char)(qhstat type[i])].r;
00459       qhstat stats[i].r= realx;
00460     }else if (qhstat type[i] != zdoc) {
00461       intx= qhstat init[(unsigned char)(qhstat type[i])].i;
00462       qhstat stats[i].i= intx;
00463     }
00464   }
00465 } /* initstatistics */

00476                                               {
00477   boolT isnew= False;
00478   int start, i;
00479 
00480   if (qhstat type[qhstat id[index]] == zdoc) 
00481     start= index+1;
00482   else
00483     start= index;
00484   for (i= start; i < qhstat next && qhstat type[qhstat id[i]] != zdoc; i++) {
00485     if (!qh_nostatistic(qhstat id[i]) && !qhstat printed[qhstat id[i]])
00486         isnew= True;
00487   }
00488   *nextindex= i;
00489   return isnew;
00490 } /* newstats */

00498                              {
00499   
00500   if ((qhstat type[i] > ZTYPEreal
00501        &&qhstat stats[i].r == qhstat init[(unsigned char)(qhstat type[i])].r)
00502       || (qhstat type[i] < ZTYPEreal
00503           &&qhstat stats[i].i == qhstat init[(unsigned char)(qhstat type[i])].i))
00504     return True;
00505   return False;
00506 } /* nostatistic */

00515                                                     {
00516 
00517   qh_allstatistics();
00518   qh_collectstatistics();
00519   qh_printstatistics (fp, string);
00520   qh_memstatistics (fp);
00521 }

00534                                                  {
00535   int i, k;
00536   realT ave;
00537   
00538   if (qh num_points != qh num_vertices) {
00539     wval_(Wpbalance)= 0;
00540     wval_(Wpbalance2)= 0;
00541   }else
00542     wval_(Wpbalance2)= qh_stddev (zval_(Zpbalance), wval_(Wpbalance), 
00543                                  wval_(Wpbalance2), &ave);
00544   wval_(Wnewbalance2)= qh_stddev (zval_(Zprocessed), wval_(Wnewbalance), 
00545                                  wval_(Wnewbalance2), &ave);
00546   fprintf (fp, "\n\
00547 %s\n\
00548  qhull invoked by: %s | %s\n%s with options:\n%s\n", string, qh rbox_command, 
00549      qh qhull_command, qh_version, qh qhull_options);
00550   fprintf (fp, "\nprecision constants:\n\
00551  %6.2g max. abs. coordinate in the (transformed) input ('Qbd:n')\n\
00552  %6.2g max. roundoff error for distance computation ('En')\n\
00553  %6.2g max. roundoff error for angle computations\n\
00554  %6.2g min. distance for outside points ('Wn')\n\
00555  %6.2g min. distance for visible facets ('Vn')\n\
00556  %6.2g max. distance for coplanar facets ('Un')\n\
00557  %6.2g max. facet width for recomputing centrum and area\n\
00558 ", 
00559   qh MAXabs_coord, qh DISTround, qh ANGLEround, qh MINoutside, 
00560         qh MINvisible, qh MAXcoplanar, qh WIDEfacet);
00561   if (qh KEEPnearinside)
00562     fprintf(fp, "\
00563  %6.2g max. distance for near-inside points\n", qh NEARinside);
00564   if (qh premerge_cos < REALmax/2) fprintf (fp, "\
00565  %6.2g max. cosine for pre-merge angle\n", qh premerge_cos);
00566   if (qh PREmerge) fprintf (fp, "\
00567  %6.2g radius of pre-merge centrum\n", qh premerge_centrum);
00568   if (qh postmerge_cos < REALmax/2) fprintf (fp, "\
00569  %6.2g max. cosine for post-merge angle\n", qh postmerge_cos);
00570   if (qh POSTmerge) fprintf (fp, "\
00571  %6.2g radius of post-merge centrum\n", qh postmerge_centrum);
00572   fprintf (fp, "\
00573  %6.2g max. distance for merging two simplicial facets\n\
00574  %6.2g max. roundoff error for arithmetic operations\n\
00575  %6.2g min. denominator for divisions\n\
00576   zero diagonal for Gauss: ", qh ONEmerge, REALepsilon, qh MINdenom);
00577   for (k=0; k < qh hull_dim; k++)
00578     fprintf (fp, "%6.2e ", qh NEARzero[k]);
00579   fprintf (fp, "\n\n");
00580   for (i=0 ; i < qhstat next; ) 
00581     qh_printstats (fp, i, &i);
00582 } /* printstatistics */

00594                                                      {
00595 #define NULLfield "       "
00596 
00597   if (id >= ZEND || qhstat printed[id])
00598     return;
00599   if (qhstat type[id] == zdoc) {
00600     fprintf (fp, "%s\n", qhstat doc[id]);
00601     return;
00602   }
00603   start= 0; /* not used */
00604   if (qh_nostatistic(id) || !qhstat doc[id])
00605     return;
00606   qhstat printed[id]= True;
00607   if (qhstat count[id] != -1 
00608       && qhstat stats[(unsigned char)(qhstat count[id])].i == 0)
00609     fprintf (fp, " *0 cnt*");
00610   else if (qhstat type[id] >= ZTYPEreal && qhstat count[id] == -1)
00611     fprintf (fp, "%7.2g", qhstat stats[id].r);
00612   else if (qhstat type[id] >= ZTYPEreal && qhstat count[id] != -1)
00613     fprintf (fp, "%7.2g", qhstat stats[id].r/ qhstat stats[(unsigned char)(qhstat count[id])].i);
00614   else if (qhstat type[id] < ZTYPEreal && qhstat count[id] == -1)
00615     fprintf (fp, "%7d", qhstat stats[id].i);
00616   else if (qhstat type[id] < ZTYPEreal && qhstat count[id] != -1)
00617     fprintf (fp, "%7.3g", (realT) qhstat stats[id].i / qhstat stats[(unsigned char)(qhstat count[id])].i);
00618   fprintf (fp, " %s\n", qhstat doc[id]);
00619 } /* printstatlevel */

00631                                                          {
00632   int j, nexti;
00633 
00634   if (qh_newstats (index, &nexti)) {
00635     fprintf (fp, "\n");
00636     for (j=index; j<nexti; j++)
00637       qh_printstatlevel (fp, qhstat id[j], 0);
00638   }
00639   if (nextindex)
00640     *nextindex= nexti;
00641 } /* printstats */

00659                                                              {
00660   realT stddev;
00661 
00662   *ave= tot/num;
00663   stddev= sqrt (tot2/num - *ave * *ave);
00664   return stddev;
00665 } /* stddev */