📄 unix.c
字号:
/*<html><pre> -<a href="qh-c.htm#qhull"
>-------------------------------</a><a name="TOP">-</a>
unix.c
Unix version of Qhull
see qh-c.htm
copyright (c) 1993-1999, The Geometry Center
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <math.h>
#include "qhull.h"
#include "mem.h"
#include "qset.h"
#if __MWERKS__ && __POWERPC__
#include <SIOUX.h>
#include <Files.h>
#include <console.h>
#include <Desk.h>
#elif __cplusplus
extern "C" {
int isatty (int);
}
#elif _MSC_VER
#include <io.h>
#define isatty _isatty
#else
int isatty (int); /* returns 1 if stdin is a tty
if "Undefined symbol" this can be deleted along with call in main() */
#endif
/*-<a href="qh-c.htm#qhull"
>-------------------------------</a><a name="prompt">-</a>
qh_prompt and qh_version
version and long prompt for qhull
notes:
change date: Changes.txt, Announce.txt, README.txt, qhull.man
qhull-news.html, Eudora signatures
change version: README.txt, qhull.html, file_id.diz, Makefile
change year: Copying.txt
check download size
recompile user_eg.c, rbox.c, qhull.c
see:
concise prompt below
*/
char qh_version[]= "version 2.6 1999/04/19";
char qh_prompta[]= "\n\
qhull- compute convex hulls and related structures.\n\
http://www.geom.umn.edu/locate/qhull %s\n\
input (stdin):\n\
first lines: dimension and number of points (or vice-versa).\n\
other lines: point coordinates, best if one point per line\n\
comments: start with a non-numeric character\n\
halfspaces: use dim plus one and put offset after coefficients.\n\
May be preceeded by a single interior point ('H').\n\
\n\
options:\n\
d - Delaunay triangulation by lifting points to a paraboloid\n\
v - Voronoi diagram from the Delaunay triangulation\n\
Hn,n,... - halfspace intersection about [n,n,0,...]\n\
d Qu - furthest-site Delaunay triangulation (upper convex hull)\n\
v Qu - furthest-site Voronoi diagram\n\
QJ - joggle input instead of merging facets\n\
\n\
Qopts- Qhull control options:\n\
Qbk:n - scale coord k so that low bound is n\n\
QBk:n - scale coord k so that upper bound is n (QBk is %2.2g)\n\
QbB - scale input to unit cube\n\
Qbb - scale last coordinate to [0,m] for Delaunay triangulations\n\
Qbk:0Bk:0 - remove k-th coordinate from input\n\
QJn - randomly joggle input in range [-n,n]\n\
QRn - random rotation (n=seed, n=0 time, n=-1 time/no rotate)\n\
%s%s%s%s"; /* split up qh_prompt for Visual C++ */
char qh_promptb[]= "\
Qc - keep coplanar points with nearest facet\n\
Qf - partition point to furthest outside facet\n\
Qg - only build good facets (needs 'QGn', 'QVn', or 'PdD')\n\
Qm - only process points that would increase max_outside\n\
Qi - keep interior points with nearest facet\n\
Qr - process random outside points instead of furthest ones\n\
Qs - search all points for the initial simplex\n\
Qu - for 'd', compute upper hull without point at-infinity\n\
returns furthest-site Delaunay triangulation\n\
Qv - test vertex neighbors for convexity\n\
Qx - exact pre-merges (skips coplanar and angle coplanar facets)\n\
Qz - add point-at-infinity to Delaunay triangulation\n\
QGn - good facet if visible from point n, -n for not visible\n\
QVn - good facet if it includes point n, -n if not\n\
Q0 - turn off default premerge with 'C-0'/'Qx'\n\
Q1 - sort merges by type instead of angle\n\
Q2 - merge all non-convex at once instead of independent sets\n\
Q3 - do not merge redundant vertices\n\
Q4 - avoid old->new merges\n\
Q5 - do not correct outer planes at end of qhull\n\
Q6 - do not pre-merge concave or coplanar facets\n\
Q7 - depth-first processing instead of breadth-first\n\
Q8 - do not process near-inside points\n\
Q9 - process furthest of furthest points\n\
";
char qh_promptc[]= "\
Topts- Trace options:\n\
T4 - trace at level n, 4=all, 5=mem/gauss, -1= events\n\
Tc - check frequently during execution\n\
Ts - print statistics\n\
Tv - verify result: structure, convexity, and point inclusion\n\
Tz - send all output to stdout\n\
TFn - report summary when n or more facets created\n\
TO file - output results to file, may be enclosed in single quotes\n\
TPn - turn on tracing when point n added to hull\n\
TMn - turn on tracing at merge n\n\
TWn - trace merge facets when width > n\n\
TRn - rerun qhull n times. Use with 'QJn'\n\
TVn - stop qhull after adding point n, -n for before (see TCn)\n\
TCn - stop qhull after building cone for point n (see TVn)\n\
\n\
Precision options (default detects precision problems without correction):\n\
Cn - radius of centrum (roundoff added). Merge facets if non-convex\n\
An - cosine of maximum angle. Merge facets if cosine > n or non-convex\n\
C-0 roundoff, A-0.99/C-0.01 pre-merge, A0.99/C0.01 post-merge\n\
En - max roundoff error for distance computation\n\
Rn - randomly perturb computations by a factor of [1-n,1+n]\n\
Vn - min distance above plane for a visible facet (default 3C-n or En)\n\
Un - max distance below plane for a new, coplanar point (default Vn)\n\
Wn - min facet width for outside point (before roundoff, default 2Vn)\n\
\n\
Output formats (may be combined; if none, produces a summary to stdout):\n\
f - all fields of all facets\n\
i - vertices incident to each facet\n\
m - Mathematica output (2-d and 3-d)\n\
o - OFF file format (dim, points and facets; Voronoi regions)\n\
n - normals with offsets\n\
p - point coordinates (Voronoi vertices)\n\
s - print summary to stderr\n\
";
char qh_promptd[]= "\
Fopts- additional input/output formats:\n\
Fa - print area for each facet\n\
FA - compute total area and volume for option 's'\n\
Fc - print count and coplanar points for each facet\n\
FC - print centrum or Voronoi center for each facet\n\
Fd - use cdd format for input (homogeneous with offset first)\n\
FD - use cdd format for normals (offset first)\n\
FF - print facets w/o ridges\n\
Fi - print inner planes for each facet\n\
Fi - print inner ridge normals (bounded) for Voronoi diagram\n\
FI - print id for each facet\n\
Fm - print merge count for each facet (511 max)\n\
Fn - print count and neighbors for each facet\n\
FN - print count and vertex neighbors for each point\n\
Fo - print outer planes (or max_outside) for each facet\n\
Fo - print outer ridge normals (unbounded)for Voronoi diagram\n\
FO - print all options\n\
Fp - print point for each halfspace intersection\n\
FP - print nearest vertex for each coplanar point (v,p,f,d)\n\
FQ - print Qhull and input command\n\
Fs - print summary: #int (6), dimension, #points,\n\
tot vertices, tot facets, #vertices in output, #facets\n\
#real (2), max outer plane and min vertex\n\
FS - print sizes: #int (0), #real(2) tot area, tot vol\n\
Ft - print Delaunay triangulation with added centrums\n\
Fv - print count and vertices for each facet\n\
Fv - print Voronoi diagram as vertices for each pair of sites\n\
FV - print average vertex (feasible point for 'H')\n\
Fx - print extreme points by point id (2-d, in order)\n\
";
char qh_prompte[]= "\
Gopts- Geomview output (2-d, 3-d and 4-d; 2-d Voronoi)\n\
Ga - all points as dots\n\
Gp - coplanar points and vertices as radii\n\
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -