| 64 |
|
RREAL uv[2]; /* (u,v) position */ |
| 65 |
|
} POINT; |
| 66 |
|
|
| 67 |
+ |
int nverts = 0; /* vertex output count */ |
| 68 |
+ |
int nnorms = 0; /* normal output count */ |
| 69 |
|
|
| 70 |
|
void loaddata(char *file, int m, int n, int pointsize); |
| 71 |
|
double l_dataval(char *nam); |
| 72 |
|
void putobjrow(POINT *rp, int n); |
| 73 |
+ |
void putobjvert(POINT *p); |
| 74 |
|
void putsquare(POINT *p0, POINT *p1, POINT *p2, POINT *p3); |
| 75 |
|
void comprow(double s, POINT *row, int siz); |
| 76 |
|
void compnorms(POINT *r0, POINT *r1, POINT *r2, int siz); |
| 86 |
|
int i, j, m, n; |
| 87 |
|
char stmp[256]; |
| 88 |
|
|
| 89 |
+ |
esupport |= E_VARIABLE|E_FUNCTION|E_RCONST; |
| 90 |
+ |
esupport &= ~(E_OUTCHAN|E_INCHAN); |
| 91 |
|
varset("PI", ':', PI); |
| 92 |
|
funset("hermite", 5, ':', l_hermite); |
| 93 |
|
funset("bezier", 5, ':', l_bezier); |
| 99 |
|
for (i = 8; i < argc; i++) |
| 100 |
|
if (!strcmp(argv[i], "-e")) |
| 101 |
|
scompile(argv[++i], NULL, 0); |
| 102 |
< |
else if (!strcmp(argv[i], "-f")) |
| 103 |
< |
fcompile(argv[++i]); |
| 104 |
< |
else if (!strcmp(argv[i], "-s")) |
| 102 |
> |
else if (!strcmp(argv[i], "-f")) { |
| 103 |
> |
char *fpath = getpath(argv[++i], getrlibpath(), 0); |
| 104 |
> |
if (fpath == NULL) { |
| 105 |
> |
fprintf(stderr, "%s: cannot find file '%s'\n", |
| 106 |
> |
argv[0], argv[i]); |
| 107 |
> |
quit(1); |
| 108 |
> |
} |
| 109 |
> |
fcompile(fpath); |
| 110 |
> |
} else if (!strcmp(argv[i], "-s")) |
| 111 |
|
smooth++; |
| 112 |
|
else if (!strcmp(argv[i], "-o")) |
| 113 |
|
objout++; |
| 116 |
|
|
| 117 |
|
modname = argv[1]; |
| 118 |
|
surfname = argv[2]; |
| 119 |
< |
m = atoi(argv[6]); |
| 120 |
< |
n = atoi(argv[7]); |
| 119 |
> |
m = eval(argv[6]) + .5; |
| 120 |
> |
n = eval(argv[7]) + .5; |
| 121 |
|
if (m <= 0 || n <= 0) |
| 122 |
|
goto userror; |
| 123 |
|
if (!strcmp(argv[5], "-") || access(argv[5], 4) == 0) { /* file? */ |
| 209 |
|
{ |
| 210 |
|
FILE *fp; |
| 211 |
|
char word[64]; |
| 212 |
< |
register int size; |
| 213 |
< |
register RREAL *dp; |
| 212 |
> |
int size; |
| 213 |
> |
RREAL *dp; |
| 214 |
|
|
| 215 |
|
datarec.flags = HASBORDER; /* assume border values */ |
| 216 |
|
datarec.m = m+1; |
| 266 |
|
) |
| 267 |
|
{ |
| 268 |
|
double u, v; |
| 269 |
< |
register int i, j; |
| 270 |
< |
register RREAL *dp; |
| 269 |
> |
int i, j; |
| 270 |
> |
RREAL *dp; |
| 271 |
|
double d00, d01, d10, d11; |
| 272 |
|
/* compute coordinates */ |
| 273 |
|
u = argument(1); v = argument(2); |
| 305 |
|
|
| 306 |
|
void |
| 307 |
|
putobjrow( /* output vertex row to .OBJ */ |
| 308 |
< |
register POINT *rp, |
| 308 |
> |
POINT *rp, |
| 309 |
|
int n |
| 310 |
|
) |
| 311 |
|
{ |
| 301 |
– |
static int nverts = 0; |
| 302 |
– |
static int nnorms = 0; |
| 303 |
– |
|
| 312 |
|
for ( ; n-- >= 0; rp++) { |
| 313 |
|
if (!rp->valid) |
| 314 |
|
continue; |
| 327 |
|
|
| 328 |
|
|
| 329 |
|
void |
| 330 |
+ |
putobjvert( /* put out OBJ vertex index triplet */ |
| 331 |
+ |
POINT *p |
| 332 |
+ |
) |
| 333 |
+ |
{ |
| 334 |
+ |
int pti = p->valid ? p->valid-nverts-1 : 0; |
| 335 |
+ |
int ni = p->nvalid ? p->nvalid-nnorms-1 : 0; |
| 336 |
+ |
|
| 337 |
+ |
printf(" %d/%d/%d", pti, pti, ni); |
| 338 |
+ |
} |
| 339 |
+ |
|
| 340 |
+ |
|
| 341 |
+ |
void |
| 342 |
|
putsquare( /* put out a square */ |
| 343 |
|
POINT *p0, |
| 344 |
|
POINT *p1, |
| 370 |
|
return; |
| 371 |
|
if (objout) { /* output .OBJ faces */ |
| 372 |
|
if (ok1 & ok2 && fdot(vc1,vc2) >= 1.0-FTINY*FTINY) { |
| 373 |
< |
printf("f %d/%d/%d %d/%d/%d %d/%d/%d %d/%d/%d\n", |
| 374 |
< |
p0->valid, p0->valid, p0->nvalid, |
| 375 |
< |
p1->valid, p1->valid, p1->nvalid, |
| 376 |
< |
p3->valid, p3->valid, p3->nvalid, |
| 357 |
< |
p2->valid, p2->valid, p2->nvalid); |
| 373 |
> |
putc('f', stdout); |
| 374 |
> |
putobjvert(p0); putobjvert(p1); |
| 375 |
> |
putobjvert(p3); putobjvert(p2); |
| 376 |
> |
putc('\n', stdout); |
| 377 |
|
return; |
| 378 |
|
} |
| 379 |
< |
if (ok1) |
| 380 |
< |
printf("f %d/%d/%d %d/%d/%d %d/%d/%d\n", |
| 381 |
< |
p0->valid, p0->valid, p0->nvalid, |
| 382 |
< |
p1->valid, p1->valid, p1->nvalid, |
| 383 |
< |
p2->valid, p2->valid, p2->nvalid); |
| 384 |
< |
if (ok2) |
| 385 |
< |
printf("f %d/%d/%d %d/%d/%d %d/%d/%d\n", |
| 386 |
< |
p2->valid, p2->valid, p2->nvalid, |
| 387 |
< |
p1->valid, p1->valid, p1->nvalid, |
| 388 |
< |
p3->valid, p3->valid, p3->nvalid); |
| 379 |
> |
if (ok1) { |
| 380 |
> |
putc('f', stdout); |
| 381 |
> |
putobjvert(p0); putobjvert(p1); putobjvert(p2); |
| 382 |
> |
putc('\n', stdout); |
| 383 |
> |
} |
| 384 |
> |
if (ok2) { |
| 385 |
> |
putc('f', stdout); |
| 386 |
> |
putobjvert(p2); putobjvert(p1); putobjvert(p3); |
| 387 |
> |
putc('\n', stdout); |
| 388 |
> |
} |
| 389 |
|
return; |
| 390 |
|
} |
| 391 |
|
/* compute normal interpolation */ |
| 456 |
|
void |
| 457 |
|
comprow( /* compute row of values */ |
| 458 |
|
double s, |
| 459 |
< |
register POINT *row, |
| 459 |
> |
POINT *row, |
| 460 |
|
int siz |
| 461 |
|
) |
| 462 |
|
{ |
| 463 |
|
double st[2]; |
| 464 |
|
int end; |
| 465 |
|
int checkvalid; |
| 466 |
< |
register int i; |
| 466 |
> |
int i; |
| 467 |
|
|
| 468 |
|
if (smooth) { |
| 469 |
|
i = -1; /* compute one past each end */ |
| 497 |
|
|
| 498 |
|
void |
| 499 |
|
compnorms( /* compute row of averaged normals */ |
| 500 |
< |
register POINT *r0, |
| 501 |
< |
register POINT *r1, |
| 502 |
< |
register POINT *r2, |
| 500 |
> |
POINT *r0, |
| 501 |
> |
POINT *r1, |
| 502 |
> |
POINT *r2, |
| 503 |
|
int siz |
| 504 |
|
) |
| 505 |
|
{ |
| 540 |
|
|
| 541 |
|
int |
| 542 |
|
norminterp( /* compute normal interpolation */ |
| 543 |
< |
register FVECT resmat[4], |
| 543 |
> |
FVECT resmat[4], |
| 544 |
|
POINT *p0, |
| 545 |
|
POINT *p1, |
| 546 |
|
POINT *p2, |
| 550 |
|
#define u ((ax+1)%3) |
| 551 |
|
#define v ((ax+2)%3) |
| 552 |
|
|
| 553 |
< |
register int ax; |
| 553 |
> |
int ax; |
| 554 |
|
MAT4 eqnmat; |
| 555 |
|
FVECT v1; |
| 556 |
< |
register int i, j; |
| 556 |
> |
int i, j; |
| 557 |
|
|
| 558 |
|
if (!smooth) /* no interpolation if no smoothing */ |
| 559 |
|
return(-1); |
