static double sc(basis_s *v,simplex *s, int k, int j) {
/* amount by which to scale up vector, for reduce_inner */
double labound;
static int lscale;
static double max_scale,
ldetbound,
Sb;
if (j<10) {
labound = logb(v->sqa)/2;
max_scale = exact_bits - labound - 0.66*(k-2)-1 -DELIFT;
if (max_scale<1) {
warning(-10, overshot exact arithmetic);
max_scale = 1;
}
if (j==0) {
int i;
neighbor *sni;
basis_s *snib;
ldetbound = DELIFT;
Sb = 0;
for (i=k-1,sni=s->neigh+k-1;i>0;i--,sni--) {
snib = sni->basis;
Sb += snib->sqb;
ldetbound += logb(snib->sqb)/2 + 1;
ldetbound -= snib->lscale;
}
}
}
if (ldetbound - v->lscale + logb(v->sqb)/2 + 1 < 0) {
DEBS(-2)
DEBTR(-2) DEBEXP(-2, ldetbound)
print_simplex_f(s, DFILE, &print_neighbor_full);
print_basis(DFILE,v);
EDEBS
return 0;
} else {
int main(int argc, char **argv) {
short pine = 0,
output = 1,
hist = 0,
vol = 0,
ahull = 0,
ofn = 0,
ifn = 0;
int option;
double alpha = 0;
int main_out_form=0, alpha_out_form=0;
simplex *root;
fg *faces_gr;
mult_up = 1;
while ((option = getopt(argc, argv, "i:m:rs:do:X:a:Af:")) != EOF) {
switch (option) {
case 'm' :
sscanf(optarg,"%lf",&mult_up);
DEBEXP(-4,mult_up);
break;
case 'a' :
vd = ahull = 1;
switch(optarg[0]) {
case 'a': sscanf(optarg+1,"%lf",&alpha); break;
case '\0': break;
default: tell_options();
}
break;
default :
tell_options();
exit(1);
}
}
INFILE = stdin;
OUTFILE = stdout;
DFILE = stderr;
read_next_site(-1);
if (dim > MAXDIM) panic("dimension bound MAXDIM exceeded");
point_size = site_size = sizeof(Coord)*dim;
shuf = &noshuffle;
get_site_n = read_next_site;
root = build_convex_hull(get_next_site, site_numm, dim, vd);
out_func* aof = out_funcs[alpha_out_form];
if (alpha==0) alpha=find_alpha(root);
alph_test(0,0,&alpha);
make_output(root, visit_outside_ashape, afacets_print, aof, OUTFILE);
free_hull_storage();
exit(0);
}
