40 |
|
|
41 |
|
int smooth = 0; /* apply smoothing? */ |
42 |
|
int objout = 0; /* output .OBJ format? */ |
43 |
+ |
int rev = 0; /* invert normal directions? */ |
44 |
|
|
45 |
|
char *modname, *surfname; |
46 |
|
|
59 |
|
|
60 |
|
typedef struct { |
61 |
|
int valid; /* point is valid (vertex number) */ |
62 |
< |
int nvalid; /* normal is valid */ |
62 |
> |
int nvalid; /* normal is valid (normal number) */ |
63 |
|
FVECT p; /* vertex position */ |
64 |
|
FVECT n; /* average normal */ |
65 |
|
RREAL uv[2]; /* (u,v) position */ |
79 |
|
|
80 |
|
|
81 |
|
int |
82 |
< |
main(argc, argv) |
82 |
< |
int argc; |
83 |
< |
char *argv[]; |
82 |
> |
main(int argc, char *argv[]) |
83 |
|
{ |
84 |
|
POINT *row0, *row1, *row2, *rp; |
85 |
|
int i, j, m, n; |
110 |
|
smooth++; |
111 |
|
else if (!strcmp(argv[i], "-o")) |
112 |
|
objout++; |
113 |
+ |
else if (!strcmp(argv[i], "-i")) |
114 |
+ |
rev = 1; |
115 |
|
else |
116 |
|
goto userror; |
117 |
|
|
153 |
|
/* print header */ |
154 |
|
fputs("# ", stdout); |
155 |
|
printargs(argc, argv, stdout); |
156 |
< |
eclock = 0; |
156 |
> |
doptimize(1); |
157 |
> |
eclock++; |
158 |
|
/* initialize */ |
159 |
|
comprow(-1.0/m, row0, n); |
160 |
|
comprow(0.0, row1, n); |
162 |
|
compnorms(row0, row1, row2, n); |
163 |
|
if (objout) { |
164 |
|
printf("\nusemtl %s\n\n", modname); |
165 |
+ |
printf("o %s\n\n", surfname); |
166 |
|
putobjrow(row1, n); |
167 |
|
} |
168 |
|
/* for each row */ |
312 |
|
int n |
313 |
|
) |
314 |
|
{ |
315 |
+ |
static FVECT prevNorm; |
316 |
+ |
|
317 |
|
for ( ; n-- >= 0; rp++) { |
318 |
|
if (!rp->valid) |
319 |
|
continue; |
320 |
|
fputs("v ", stdout); |
321 |
|
pvect(rp->p); |
322 |
< |
if (smooth && !ZEROVECT(rp->n)) { |
322 |
> |
rp->valid = ++nverts; |
323 |
> |
printf("\tvt %.9g %.9g\n", rp->uv[0], rp->uv[1]); |
324 |
> |
if (!smooth || ZEROVECT(rp->n)) |
325 |
> |
rp->nvalid = 0; |
326 |
> |
else if (VABSEQ(rp->n, prevNorm)) |
327 |
> |
rp->nvalid = nnorms; |
328 |
> |
else { |
329 |
|
printf("\tvn %.9g %.9g %.9g\n", |
330 |
|
rp->n[0], rp->n[1], rp->n[2]); |
331 |
|
rp->nvalid = ++nnorms; |
332 |
< |
} else |
333 |
< |
rp->nvalid = 0; |
323 |
< |
printf("\tvt %.9g %.9g\n", rp->uv[0], rp->uv[1]); |
324 |
< |
rp->valid = ++nverts; |
332 |
> |
VCOPY(prevNorm, rp->n); |
333 |
> |
} |
334 |
|
} |
335 |
|
} |
336 |
|
|
360 |
|
int axis; |
361 |
|
FVECT v1, v2, vc1, vc2; |
362 |
|
int ok1, ok2; |
363 |
+ |
|
364 |
+ |
if (rev) { /* reverse normals? */ |
365 |
+ |
POINT *pt = p1; p1 = p2; p2 = pt; |
366 |
+ |
} |
367 |
|
/* compute exact normals */ |
368 |
|
ok1 = (p0->valid && p1->valid && p2->valid); |
369 |
|
if (ok1) { |
544 |
|
fvsum(v2, r1[0].p, r1[-1].p, -1.0); |
545 |
|
else |
546 |
|
fvsum(v2, r1[1].p, r1[-1].p, -1.0); |
547 |
< |
fcross(r1[0].n, v1, v2); |
547 |
> |
if (rev) |
548 |
> |
fcross(r1[0].n, v2, v1); |
549 |
> |
else |
550 |
> |
fcross(r1[0].n, v1, v2); |
551 |
|
normalize(r1[0].n); |
552 |
|
skip: |
553 |
|
r0++; r1++; r2++; |