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 */ |
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 |
> |
eclock = 1; |
157 |
|
/* initialize */ |
158 |
|
comprow(-1.0/m, row0, n); |
159 |
|
comprow(0.0, row1, n); |
311 |
|
int n |
312 |
|
) |
313 |
|
{ |
314 |
+ |
static FVECT prevNorm; |
315 |
+ |
|
316 |
|
for ( ; n-- >= 0; rp++) { |
317 |
|
if (!rp->valid) |
318 |
|
continue; |
319 |
|
fputs("v ", stdout); |
320 |
|
pvect(rp->p); |
321 |
< |
if (smooth && !ZEROVECT(rp->n)) { |
321 |
> |
rp->valid = ++nverts; |
322 |
> |
printf("\tvt %.9g %.9g\n", rp->uv[0], rp->uv[1]); |
323 |
> |
if (!smooth || ZEROVECT(rp->n)) |
324 |
> |
rp->nvalid = 0; |
325 |
> |
else if (VABSEQ(rp->n, prevNorm)) |
326 |
> |
rp->nvalid = nnorms; |
327 |
> |
else { |
328 |
|
printf("\tvn %.9g %.9g %.9g\n", |
329 |
|
rp->n[0], rp->n[1], rp->n[2]); |
330 |
|
rp->nvalid = ++nnorms; |
331 |
< |
} else |
332 |
< |
rp->nvalid = 0; |
322 |
< |
printf("\tvt %.9g %.9g\n", rp->uv[0], rp->uv[1]); |
323 |
< |
rp->valid = ++nverts; |
331 |
> |
VCOPY(prevNorm, rp->n); |
332 |
> |
} |
333 |
|
} |
334 |
|
} |
335 |
|
|
359 |
|
int axis; |
360 |
|
FVECT v1, v2, vc1, vc2; |
361 |
|
int ok1, ok2; |
362 |
+ |
|
363 |
+ |
if (rev) { /* reverse normals? */ |
364 |
+ |
POINT *pt = p1; p1 = p2; p2 = pt; |
365 |
+ |
} |
366 |
|
/* compute exact normals */ |
367 |
|
ok1 = (p0->valid && p1->valid && p2->valid); |
368 |
|
if (ok1) { |
543 |
|
fvsum(v2, r1[0].p, r1[-1].p, -1.0); |
544 |
|
else |
545 |
|
fvsum(v2, r1[1].p, r1[-1].p, -1.0); |
546 |
< |
fcross(r1[0].n, v1, v2); |
546 |
> |
if (rev) |
547 |
> |
fcross(r1[0].n, v2, v1); |
548 |
> |
else |
549 |
> |
fcross(r1[0].n, v1, v2); |
550 |
|
normalize(r1[0].n); |
551 |
|
skip: |
552 |
|
r0++; r1++; r2++; |