1 |
#ifndef lint |
2 |
static const char RCSid[] = "$Id: objutil.c,v 2.16 2021/04/07 16:40:28 greg Exp $"; |
3 |
#endif |
4 |
/* |
5 |
* Basic .OBJ scene handling routines. |
6 |
* |
7 |
* Created by Greg Ward on Wed Feb 11 2004. |
8 |
*/ |
9 |
|
10 |
#include <stdlib.h> |
11 |
#include <ctype.h> |
12 |
#include "rtio.h" |
13 |
#include "rtmath.h" |
14 |
#include "rterror.h" |
15 |
#include "lookup.h" |
16 |
#include "objutil.h" |
17 |
|
18 |
/* Find an existing name in a list of names */ |
19 |
int |
20 |
findName(const char *nm, const char **nmlist, int n) |
21 |
{ |
22 |
int i; |
23 |
|
24 |
for (i = n; i-- > 0; ) |
25 |
if (!strcmp(nmlist[i], nm)) |
26 |
break; |
27 |
return(i); |
28 |
} |
29 |
|
30 |
/* Clear face selection */ |
31 |
void |
32 |
clearSelection(Scene *sc, int set) |
33 |
{ |
34 |
Face *f; |
35 |
|
36 |
for (f = sc->flist; f != NULL; f = f->next) |
37 |
if (set) |
38 |
f->flags |= FACE_SELECTED; |
39 |
else |
40 |
f->flags &= ~FACE_SELECTED; |
41 |
} |
42 |
|
43 |
/* Select faces by object name (modifies current) */ |
44 |
void |
45 |
selectGroup(Scene *sc, const char *gname, int invert) |
46 |
{ |
47 |
int gid = getGroupID(sc, gname); |
48 |
Face *f; |
49 |
|
50 |
if (gid < 0) |
51 |
return; |
52 |
for (f = sc->flist; f != NULL; f = f->next) |
53 |
if (f->grp == gid) { |
54 |
if (invert) |
55 |
f->flags &= ~FACE_SELECTED; |
56 |
else |
57 |
f->flags |= FACE_SELECTED; |
58 |
} |
59 |
} |
60 |
|
61 |
/* Select faces by material name (modifies current) */ |
62 |
void |
63 |
selectMaterial(Scene *sc, const char *mname, int invert) |
64 |
{ |
65 |
int mid = getMaterialID(sc, mname); |
66 |
Face *f; |
67 |
|
68 |
if (mid < 0) |
69 |
return; |
70 |
for (f = sc->flist; f != NULL; f = f->next) |
71 |
if (f->mat == mid) { |
72 |
if (invert) |
73 |
f->flags &= ~FACE_SELECTED; |
74 |
else |
75 |
f->flags |= FACE_SELECTED; |
76 |
} |
77 |
} |
78 |
|
79 |
/* Invert face selection */ |
80 |
void |
81 |
invertSelection(Scene *sc) |
82 |
{ |
83 |
Face *f; |
84 |
|
85 |
for (f = sc->flist; f != NULL; f = f->next) |
86 |
f->flags ^= FACE_SELECTED; |
87 |
} |
88 |
|
89 |
/* Count selected faces */ |
90 |
int |
91 |
numberSelected(Scene *sc) |
92 |
{ |
93 |
int nsel = 0; |
94 |
Face *f; |
95 |
|
96 |
for (f = sc->flist; f != NULL; f = f->next) |
97 |
nsel += ((f->flags & FACE_SELECTED) != 0); |
98 |
return(nsel); |
99 |
} |
100 |
|
101 |
/* Execute callback on indicated faces */ |
102 |
int |
103 |
foreachFace(Scene *sc, int (*cb)(Scene *, Face *, void *), |
104 |
int flreq, int flexc, void *c_data) |
105 |
{ |
106 |
const int fltest = flreq | flexc; |
107 |
int sum = 0; |
108 |
Face *f; |
109 |
|
110 |
if ((sc == NULL) | (cb == NULL)) |
111 |
return(0); |
112 |
for (f = sc->flist; f != NULL; f = f->next) |
113 |
if ((f->flags & fltest) == flreq) { |
114 |
int res = (*cb)(sc, f, c_data); |
115 |
if (res < 0) |
116 |
return(res); |
117 |
sum += res; |
118 |
} |
119 |
return(sum); |
120 |
} |
121 |
|
122 |
/* Callback for removeTexture() */ |
123 |
static int |
124 |
remFaceTexture(Scene *sc, Face *f, void *dummy) |
125 |
{ |
126 |
int hadTexture = 0; |
127 |
int i; |
128 |
|
129 |
for (i = f->nv; i-- > 0; ) { |
130 |
if (f->v[i].tid < 0) |
131 |
continue; |
132 |
f->v[i].tid = -1; |
133 |
hadTexture = 1; |
134 |
} |
135 |
return(hadTexture); |
136 |
} |
137 |
|
138 |
/* Remove texture coordinates from the indicated faces */ |
139 |
int |
140 |
removeTexture(Scene *sc, int flreq, int flexc) |
141 |
{ |
142 |
return(foreachFace(sc, remFaceTexture, flreq, flexc, NULL)); |
143 |
} |
144 |
|
145 |
/* Callback for removeNormals() */ |
146 |
static int |
147 |
remFaceNormal(Scene *sc, Face *f, void *dummy) |
148 |
{ |
149 |
int hadNormal = 0; |
150 |
int i; |
151 |
|
152 |
for (i = f->nv; i-- > 0; ) { |
153 |
if (f->v[i].nid < 0) |
154 |
continue; |
155 |
f->v[i].nid = -1; |
156 |
hadNormal = 1; |
157 |
} |
158 |
return(hadNormal); |
159 |
} |
160 |
|
161 |
/* Remove surface normals from the indicated faces */ |
162 |
int |
163 |
removeNormals(Scene *sc, int flreq, int flexc) |
164 |
{ |
165 |
return(foreachFace(sc, remFaceNormal, flreq, flexc, NULL)); |
166 |
} |
167 |
|
168 |
/* Callback for changeGroup() */ |
169 |
static int |
170 |
chngFaceGroup(Scene *sc, Face *f, void *ptr) |
171 |
{ |
172 |
int grp = *(int *)ptr; |
173 |
|
174 |
if (f->grp == grp) |
175 |
return(0); |
176 |
f->grp = grp; |
177 |
return(1); |
178 |
} |
179 |
|
180 |
/* Change group for the indicated faces */ |
181 |
int |
182 |
changeGroup(Scene *sc, const char *gname, int flreq, int flexc) |
183 |
{ |
184 |
int grp = getGroupID(sc, gname); |
185 |
|
186 |
if (grp < 0) { |
187 |
sc->grpname = chunk_alloc(char *, sc->grpname, sc->ngrps); |
188 |
sc->grpname[grp=sc->ngrps++] = savqstr((char *)gname); |
189 |
} |
190 |
return(foreachFace(sc, chngFaceGroup, flreq, flexc, (void *)&grp)); |
191 |
} |
192 |
|
193 |
/* Callback for changeMaterial() */ |
194 |
static int |
195 |
chngFaceMaterial(Scene *sc, Face *f, void *ptr) |
196 |
{ |
197 |
int mat = *(int *)ptr; |
198 |
|
199 |
if (f->mat == mat) |
200 |
return(0); |
201 |
f->mat = mat; |
202 |
return(1); |
203 |
} |
204 |
|
205 |
/* Change material for the indicated faces */ |
206 |
int |
207 |
changeMaterial(Scene *sc, const char *mname, int flreq, int flexc) |
208 |
{ |
209 |
int mat = getMaterialID(sc, mname); |
210 |
|
211 |
if (mat < 0) { |
212 |
sc->matname = chunk_alloc(char *, sc->matname, sc->nmats); |
213 |
sc->matname[mat=sc->nmats++] = savqstr((char *)mname); |
214 |
} |
215 |
return(foreachFace(sc, chngFaceMaterial, flreq, flexc, (void *)&mat)); |
216 |
} |
217 |
|
218 |
/* Compare floating point values for (near) equality */ |
219 |
static int |
220 |
fdiffer(double f1, double f2, double eps) |
221 |
{ |
222 |
if (f2 != .0) |
223 |
f1 = f1/f2 - 1.; |
224 |
return((f1 > eps) | (f1 < -eps)); |
225 |
} |
226 |
|
227 |
/* Compare two texture coordinates for (near) equality */ |
228 |
static int |
229 |
tex_diff(const TexCoord *t1, const TexCoord *t2, double eps) |
230 |
{ |
231 |
if (fdiffer(t1->u, t2->u, eps)) |
232 |
return(1); |
233 |
if (fdiffer(t1->v, t2->v, eps)) |
234 |
return(1); |
235 |
return(0); |
236 |
} |
237 |
|
238 |
/* Compare two surface normals for (near) equality */ |
239 |
static int |
240 |
norm_diff(const Normal n1, const Normal n2, double eps) |
241 |
{ |
242 |
if (fabs(n1[0]-n2[0]) > eps) |
243 |
return(1); |
244 |
if (fabs(n1[1]-n2[1]) > eps) |
245 |
return(1); |
246 |
if (fabs(n1[2]-n2[2]) > eps) |
247 |
return(1); |
248 |
return(0); |
249 |
} |
250 |
|
251 |
/* Replace the given vertex with an equivalent one */ |
252 |
static int |
253 |
replace_vertex(Scene *sc, int prev, int repl, double eps) |
254 |
{ |
255 |
int repl_tex[10]; |
256 |
int repl_norm[10]; |
257 |
Face *f, *flast=NULL; |
258 |
int i, j=0; |
259 |
/* check to see if it's even used */ |
260 |
if (sc->vert[prev].vflist == NULL) |
261 |
return(0); |
262 |
/* get replacement textures */ |
263 |
repl_tex[0] = -1; |
264 |
for (f = sc->vert[repl].vflist; f != NULL; f = f->v[j].fnext) { |
265 |
/* make sure prev isn't in there */ |
266 |
for (j = f->nv; j-- > 0; ) |
267 |
if (f->v[j].vid == prev) |
268 |
return(0); |
269 |
for (j = f->nv; j-- > 0; ) |
270 |
if (f->v[j].vid == repl) |
271 |
break; |
272 |
if (j < 0) |
273 |
goto linkerr; |
274 |
if (f->v[j].tid < 0) |
275 |
continue; |
276 |
/* see if it's new */ |
277 |
for (i = 0; repl_tex[i] >= 0; i++) { |
278 |
if (repl_tex[i] == f->v[j].tid) |
279 |
break; |
280 |
if (!tex_diff(&sc->tex[repl_tex[i]], |
281 |
&sc->tex[f->v[j].tid], eps)) { |
282 |
f->v[j].tid = repl_tex[i]; |
283 |
break; /* consolidate */ |
284 |
} |
285 |
} |
286 |
if (repl_tex[i] >= 0) |
287 |
continue; /* have this one already */ |
288 |
/* else add it */ |
289 |
repl_tex[i++] = f->v[j].tid; |
290 |
repl_tex[i] = -1; |
291 |
if (i >= 9) |
292 |
break; /* that's all we have room for */ |
293 |
} |
294 |
/* get replacement normals */ |
295 |
repl_norm[0] = -1; |
296 |
for (f = sc->vert[repl].vflist; f != NULL; f = f->v[j].fnext) { |
297 |
for (j = f->nv; j-- > 0; ) |
298 |
if (f->v[j].vid == repl) |
299 |
break; |
300 |
if (f->v[j].nid < 0) |
301 |
continue; |
302 |
/* see if it's new */ |
303 |
for (i = 0; repl_norm[i] >= 0; i++) { |
304 |
if (repl_norm[i] == f->v[j].nid) |
305 |
break; |
306 |
if (!norm_diff(sc->norm[repl_norm[i]], |
307 |
sc->norm[f->v[j].nid], eps)) { |
308 |
f->v[j].nid = repl_norm[i]; |
309 |
break; /* consolidate */ |
310 |
} |
311 |
} |
312 |
if (repl_norm[i] >= 0) |
313 |
continue; /* have this one already */ |
314 |
/* else add it */ |
315 |
repl_norm[i++] = f->v[j].nid; |
316 |
repl_norm[i] = -1; |
317 |
if (i >= 9) |
318 |
break; /* that's all we have room for */ |
319 |
} |
320 |
/* replace occurrences of vertex */ |
321 |
for (f = sc->vert[prev].vflist; f != NULL; f = f->v[j].fnext) { |
322 |
for (j = f->nv; j-- > 0; ) |
323 |
if (f->v[j].vid == prev) |
324 |
break; |
325 |
if (j < 0) |
326 |
goto linkerr; |
327 |
/* XXX doesn't allow for multiple references to prev in face */ |
328 |
f->v[j].vid = repl; /* replace vertex itself */ |
329 |
if (faceArea(sc, f, NULL) <= FTINY*FTINY) |
330 |
f->flags |= FACE_DEGENERATE; |
331 |
if (f->v[j].tid >= 0) /* replace texture if appropriate */ |
332 |
for (i = 0; repl_tex[i] >= 0; i++) { |
333 |
if (repl_tex[i] == f->v[j].tid) |
334 |
break; |
335 |
if (!tex_diff(&sc->tex[repl_tex[i]], |
336 |
&sc->tex[f->v[j].tid], eps)) { |
337 |
f->v[j].tid = repl_tex[i]; |
338 |
break; |
339 |
} |
340 |
} |
341 |
if (f->v[j].nid >= 0) /* replace normal if appropriate */ |
342 |
for (i = 0; repl_norm[i] >= 0; i++) { |
343 |
if (repl_norm[i] == f->v[j].nid) |
344 |
break; |
345 |
if (!norm_diff(sc->norm[repl_norm[i]], |
346 |
sc->norm[f->v[j].nid], eps)) { |
347 |
f->v[j].nid = repl_norm[i]; |
348 |
break; |
349 |
} |
350 |
} |
351 |
flast = f; |
352 |
} |
353 |
/* transfer face list */ |
354 |
flast->v[j].fnext = sc->vert[repl].vflist; |
355 |
sc->vert[repl].vflist = sc->vert[prev].vflist; |
356 |
sc->vert[prev].vflist = NULL; |
357 |
return(1); |
358 |
linkerr: |
359 |
error(CONSISTENCY, "Link error in replace_vertex()"); |
360 |
return(0); /* shouldn't return */ |
361 |
} |
362 |
|
363 |
/* Eliminate duplicate vertices, return # joined */ |
364 |
int |
365 |
coalesceVertices(Scene *sc, double eps) |
366 |
{ |
367 |
int nelim = 0; |
368 |
LUTAB myLookup; |
369 |
LUENT *le; |
370 |
char vertfmt[32], vertbuf[64]; |
371 |
double d; |
372 |
int i; |
373 |
|
374 |
if (eps >= 1.) |
375 |
return(0); |
376 |
if (sc->nverts <= 3) |
377 |
return(0); |
378 |
/* create hash table */ |
379 |
myLookup.hashf = lu_shash; |
380 |
myLookup.keycmp = (lut_keycmpf_t *)strcmp; |
381 |
myLookup.freek = (lut_free_t *)freeqstr; |
382 |
myLookup.freed = NULL; |
383 |
if (!lu_init(&myLookup, sc->nverts)) |
384 |
return(0); |
385 |
if (eps <= 5e-15) |
386 |
strcpy(vertfmt, "%.15e %.15e %.15e"); |
387 |
else { |
388 |
int nsigdig = 0; |
389 |
for (d = eps; d < 0.5; d *= 10.) |
390 |
++nsigdig; |
391 |
sprintf(vertfmt, "%%.%de %%.%de %%.%de", |
392 |
nsigdig, nsigdig, nsigdig); |
393 |
} |
394 |
/* find coicident vertices */ |
395 |
for (i = 0; i < sc->nverts; i++) { |
396 |
if (verbose && !((i+1) & 0x3fff)) |
397 |
fprintf(stderr, "%3d%% complete\r", 100*i/sc->nverts); |
398 |
/* check for match */ |
399 |
sprintf(vertbuf, vertfmt, sc->vert[i].p[0], |
400 |
sc->vert[i].p[1], sc->vert[i].p[2]); |
401 |
le = lu_find(&myLookup, vertbuf); |
402 |
if (le->data != NULL) { /* coincident vertex */ |
403 |
nelim += replace_vertex(sc, i, |
404 |
(Vertex *)le->data - sc->vert, eps); |
405 |
continue; |
406 |
} |
407 |
if (le->key == NULL) /* else create new entry */ |
408 |
le->key = savqstr(vertbuf); |
409 |
le->data = (char *)&sc->vert[i]; |
410 |
} |
411 |
lu_done(&myLookup); /* clean up */ |
412 |
return(nelim); |
413 |
} |
414 |
|
415 |
/* Identify duplicate faces */ |
416 |
int |
417 |
findDuplicateFaces(Scene *sc) |
418 |
{ |
419 |
int nchecked = 0; |
420 |
int nfound = 0; |
421 |
Face *f, *f1; |
422 |
int vid; |
423 |
int i, j; |
424 |
/* start fresh */ |
425 |
for (f = sc->flist; f != NULL; f = f->next) |
426 |
f->flags &= ~FACE_DUPLICATE; |
427 |
/* check each face */ |
428 |
for (f = sc->flist; f != NULL; f = f->next) { |
429 |
nchecked++; |
430 |
if (verbose && !(nchecked & 0x3fff)) |
431 |
fprintf(stderr, "%3d%% complete\r", |
432 |
100*nchecked/sc->nfaces); |
433 |
if (f->flags & FACE_DUPLICATE) |
434 |
continue; /* already identified */ |
435 |
vid = f->v[0].vid; |
436 |
/* look for duplicates */ |
437 |
for (f1 = sc->vert[vid].vflist; f1 != NULL; |
438 |
f1 = f1->v[j].fnext) { |
439 |
for (j = f1->nv; j-- > 0; ) |
440 |
if (f1->v[j].vid == vid) |
441 |
break; |
442 |
if (j < 0) |
443 |
break; /* missing link! */ |
444 |
if (f1 == f) |
445 |
continue; /* shouldn't happen */ |
446 |
if (f1->flags & FACE_DUPLICATE) |
447 |
continue; /* already marked */ |
448 |
if (f1->nv != f->nv) |
449 |
continue; /* couldn't be dup. */ |
450 |
for (i = f->nv; --i > 0; ) |
451 |
if (f->v[i].vid != f1->v[(j+i)%f1->nv].vid) |
452 |
break; /* vertex mismatch */ |
453 |
if (i) { |
454 |
#if DUP_CHECK_REVERSE /* check reverse direction */ |
455 |
for (i = f->nv; --i > 0; ) |
456 |
if (f1->v[(j+f1->nv-i)%f1->nv].vid |
457 |
!= f->v[i].vid) |
458 |
break; |
459 |
if (i) /* no match */ |
460 |
#endif |
461 |
continue; |
462 |
} |
463 |
f1->flags |= FACE_DUPLICATE; |
464 |
++nfound; |
465 |
} |
466 |
} |
467 |
if (verbose) |
468 |
fprintf(stderr, "Found %d duplicate faces\n", nfound); |
469 |
return(nfound); |
470 |
} |
471 |
|
472 |
/* Delete indicated faces */ |
473 |
int |
474 |
deleteFaces(Scene *sc, int flreq, int flexc) |
475 |
{ |
476 |
const int fltest = flreq | flexc; |
477 |
int orig_nfaces = sc->nfaces; |
478 |
Face fhead; |
479 |
Face *f, *ftst; |
480 |
|
481 |
fhead.next = sc->flist; |
482 |
f = &fhead; |
483 |
while ((ftst = f->next) != NULL) |
484 |
if ((ftst->flags & fltest) == flreq) { |
485 |
Face *vf; /* remove from vertex lists */ |
486 |
int vid, i, j; |
487 |
for (i = 0; i < ftst->nv; i++) { |
488 |
vid = ftst->v[i].vid; |
489 |
vf = sc->vert[vid].vflist; |
490 |
if (vf == ftst) { |
491 |
sc->vert[vid].vflist = ftst->v[i].fnext; |
492 |
continue; |
493 |
} |
494 |
while (vf != NULL) { |
495 |
for (j = vf->nv; j-- > 0; ) |
496 |
if (vf->v[j].vid == vid) |
497 |
break; |
498 |
if (j < 0) |
499 |
break; /* error */ |
500 |
if (vf->v[j].fnext == ftst) { |
501 |
vf->v[j].fnext = |
502 |
ftst->v[i].fnext; |
503 |
break; |
504 |
} |
505 |
vf = vf->v[j].fnext; |
506 |
} |
507 |
} |
508 |
f->next = ftst->next; /* remove from scene list */ |
509 |
efree((char *)ftst); |
510 |
sc->nfaces--; |
511 |
} else |
512 |
f = f->next; |
513 |
sc->flist = fhead.next; |
514 |
return(orig_nfaces - sc->nfaces); |
515 |
} |
516 |
|
517 |
/* Compute face area (and normal) */ |
518 |
double |
519 |
faceArea(const Scene *sc, const Face *f, Normal nrm) |
520 |
{ |
521 |
FVECT fnrm; |
522 |
double area; |
523 |
FVECT v1, v2, v3; |
524 |
double *p0; |
525 |
int i; |
526 |
|
527 |
if (f->flags & FACE_DEGENERATE) |
528 |
return(.0); /* should we check this? */ |
529 |
fnrm[0] = fnrm[1] = fnrm[2] = .0; |
530 |
p0 = sc->vert[f->v[0].vid].p; |
531 |
VSUB(v1, sc->vert[f->v[1].vid].p, p0); |
532 |
for (i = 2; i < f->nv; i++) { |
533 |
VSUB(v2, sc->vert[f->v[i].vid].p, p0); |
534 |
fcross(v3, v1, v2); |
535 |
VADD(fnrm, fnrm, v3); |
536 |
VCOPY(v1, v2); |
537 |
} |
538 |
area = 0.5*normalize(fnrm); |
539 |
if (nrm != NULL) { |
540 |
if (area != 0.) |
541 |
VCOPY(nrm, fnrm); |
542 |
else |
543 |
nrm[0] = nrm[1] = nrm[2] = 0.; |
544 |
} |
545 |
return(area); |
546 |
} |
547 |
|
548 |
/* Add a descriptive comment */ |
549 |
void |
550 |
addComment(Scene *sc, const char *comment) |
551 |
{ |
552 |
sc->descr = chunk_alloc(char *, sc->descr, sc->ndescr); |
553 |
sc->descr[sc->ndescr++] = savqstr((char *)comment); |
554 |
} |
555 |
|
556 |
/* Find index for comment containing the given string (starting from n) */ |
557 |
int |
558 |
findComment(Scene *sc, const char *match, int n) |
559 |
{ |
560 |
if (n >= sc->ndescr) |
561 |
return(-1); |
562 |
n *= (n > 0); |
563 |
while (n < sc->ndescr) |
564 |
if (strstr(sc->descr[n], match) != NULL) |
565 |
return(n); |
566 |
return(-1); |
567 |
} |
568 |
|
569 |
/* Clear comments */ |
570 |
void |
571 |
clearComments(Scene *sc) |
572 |
{ |
573 |
while (sc->ndescr > 0) |
574 |
freeqstr(sc->descr[--sc->ndescr]); |
575 |
efree((char *)sc->descr); |
576 |
sc->ndescr = 0; |
577 |
} |
578 |
|
579 |
/* Add a face to a vertex face list */ |
580 |
static void |
581 |
add2facelist(Scene *sc, Face *f, int i) |
582 |
{ |
583 |
int vid = f->v[i].vid; |
584 |
Face *fp = sc->vert[vid].vflist; |
585 |
int j; |
586 |
|
587 |
f->v[i].fnext = NULL; /* will put at end */ |
588 |
if (fp == NULL) { /* new list */ |
589 |
sc->vert[vid].vflist = f; |
590 |
return; |
591 |
} |
592 |
for ( ; ; ) { /* else find position */ |
593 |
if (fp == f) |
594 |
return; /* already in list */ |
595 |
for (j = fp->nv; j-- > 0; ) |
596 |
if (fp->v[j].vid == vid) |
597 |
break; |
598 |
if (j < 0) |
599 |
error(CONSISTENCY, "Link error in add2facelist()"); |
600 |
if (fp->v[j].fnext == NULL) |
601 |
break; /* reached the end */ |
602 |
fp = fp->v[j].fnext; |
603 |
} |
604 |
fp->v[j].fnext = f; /* append new face */ |
605 |
} |
606 |
|
607 |
/* Add a vertex to a scene */ |
608 |
int |
609 |
addVertex(Scene *sc, double x, double y, double z) |
610 |
{ |
611 |
sc->vert = chunk_alloc(Vertex, sc->vert, sc->nverts); |
612 |
sc->vert[sc->nverts].p[0] = x; |
613 |
sc->vert[sc->nverts].p[1] = y; |
614 |
sc->vert[sc->nverts].p[2] = z; |
615 |
sc->vert[sc->nverts].vflist = NULL; |
616 |
return(sc->nverts++); |
617 |
} |
618 |
|
619 |
/* Add a texture coordinate to a scene */ |
620 |
int |
621 |
addTexture(Scene *sc, double u, double v) |
622 |
{ |
623 |
sc->tex = chunk_alloc(TexCoord, sc->tex, sc->ntex); |
624 |
sc->tex[sc->ntex].u = u; |
625 |
sc->tex[sc->ntex].v = v; |
626 |
return(sc->ntex++); |
627 |
} |
628 |
|
629 |
/* Add a surface normal to a scene */ |
630 |
int |
631 |
addNormal(Scene *sc, double xn, double yn, double zn) |
632 |
{ |
633 |
FVECT nrm; |
634 |
|
635 |
nrm[0] = xn; nrm[1] = yn; nrm[2] = zn; |
636 |
if (normalize(nrm) == .0) |
637 |
return(-1); |
638 |
sc->norm = chunk_alloc(Normal, sc->norm, sc->nnorms); |
639 |
VCOPY(sc->norm[sc->nnorms], nrm); |
640 |
return(sc->nnorms++); |
641 |
} |
642 |
|
643 |
/* Set current (last) group */ |
644 |
void |
645 |
setGroup(Scene *sc, const char *nm) |
646 |
{ |
647 |
sc->lastgrp = getGroupID(sc, nm); |
648 |
if (sc->lastgrp >= 0) |
649 |
return; |
650 |
sc->grpname = chunk_alloc(char *, sc->grpname, sc->ngrps); |
651 |
sc->grpname[sc->lastgrp=sc->ngrps++] = savqstr((char *)nm); |
652 |
} |
653 |
|
654 |
/* Set current (last) material */ |
655 |
void |
656 |
setMaterial(Scene *sc, const char *nm) |
657 |
{ |
658 |
sc->lastmat = getMaterialID(sc, nm); |
659 |
if (sc->lastmat >= 0) |
660 |
return; |
661 |
sc->matname = chunk_alloc(char *, sc->matname, sc->nmats); |
662 |
sc->matname[sc->lastmat=sc->nmats++] = savqstr((char *)nm); |
663 |
} |
664 |
|
665 |
/* Add new face to a scene */ |
666 |
Face * |
667 |
addFace(Scene *sc, VNDX vid[], int nv) |
668 |
{ |
669 |
Face *f; |
670 |
int i; |
671 |
|
672 |
if (nv < 3) |
673 |
return(NULL); |
674 |
f = (Face *)emalloc(sizeof(Face)+sizeof(VertEnt)*(nv-3)); |
675 |
f->flags = 0; |
676 |
f->nv = nv; |
677 |
f->grp = sc->lastgrp; |
678 |
f->mat = sc->lastmat; |
679 |
for (i = 0; i < nv; i++) { /* add each vertex */ |
680 |
int j; |
681 |
f->v[i].vid = vid[i][0]; |
682 |
f->v[i].tid = vid[i][1]; |
683 |
f->v[i].nid = vid[i][2]; |
684 |
f->v[i].fnext = NULL; |
685 |
for (j = i; j-- > 0; ) |
686 |
if (f->v[j].vid == vid[i][0]) |
687 |
break; |
688 |
if (j < 0) { /* first occurrence? */ |
689 |
f->v[i].fnext = sc->vert[vid[i][0]].vflist; |
690 |
sc->vert[vid[i][0]].vflist = f; |
691 |
} else if (nv == 3) /* degenerate triangle? */ |
692 |
f->flags |= FACE_DEGENERATE; |
693 |
} |
694 |
f->next = sc->flist; /* push onto face list */ |
695 |
sc->flist = f; |
696 |
sc->nfaces++; |
697 |
/* check face area */ |
698 |
if (!(f->flags & FACE_DEGENERATE) && faceArea(sc, f, NULL) <= FTINY*FTINY) |
699 |
f->flags |= FACE_DEGENERATE; |
700 |
return(f); |
701 |
} |
702 |
|
703 |
/* Allocate an empty scene */ |
704 |
Scene * |
705 |
newScene(void) |
706 |
{ |
707 |
Scene *sc = (Scene *)ecalloc(1, sizeof(Scene)); |
708 |
/* default group & material */ |
709 |
sc->grpname = chunk_alloc(char *, sc->grpname, sc->ngrps); |
710 |
sc->grpname[sc->ngrps++] = savqstr("DEFAULT_GROUP"); |
711 |
sc->matname = chunk_alloc(char *, sc->matname, sc->nmats); |
712 |
sc->matname[sc->nmats++] = savqstr("DEFAULT_MATERIAL"); |
713 |
|
714 |
return(sc); |
715 |
} |
716 |
|
717 |
/* Duplicate a scene, optionally selecting faces */ |
718 |
Scene * |
719 |
dupScene(const Scene *osc, int flreq, int flexc) |
720 |
{ |
721 |
int fltest = flreq | flexc; |
722 |
Scene *sc; |
723 |
const Face *fo; |
724 |
Face *f; |
725 |
int i; |
726 |
|
727 |
if (osc == NULL) |
728 |
return(NULL); |
729 |
sc = newScene(); |
730 |
for (i = 0; i < osc->ndescr; i++) |
731 |
addComment(sc, osc->descr[i]); |
732 |
if (osc->ngrps > 1) { |
733 |
sc->grpname = (char **)erealloc((char *)sc->grpname, |
734 |
sizeof(char *) * (osc->ngrps+(CHUNKSIZ-1))); |
735 |
for (i = 1; i < osc->ngrps; i++) |
736 |
sc->grpname[i] = savqstr(osc->grpname[i]); |
737 |
sc->ngrps = osc->ngrps; |
738 |
} |
739 |
if (osc->nmats > 1) { |
740 |
sc->matname = (char **)erealloc((char *)sc->matname, |
741 |
sizeof(char *) * (osc->nmats+(CHUNKSIZ-1))); |
742 |
for (i = 1; i < osc->nmats; i++) |
743 |
sc->matname[i] = savqstr(osc->matname[i]); |
744 |
sc->nmats = osc->nmats; |
745 |
} |
746 |
if (osc->nverts) { |
747 |
sc->vert = (Vertex *)emalloc(sizeof(Vertex) * |
748 |
(osc->nverts+(CHUNKSIZ-1))); |
749 |
memcpy(sc->vert, osc->vert, sizeof(Vertex)*osc->nverts); |
750 |
sc->nverts = osc->nverts; |
751 |
for (i = 0; i < sc->nverts; i++) |
752 |
sc->vert[i].vflist = NULL; |
753 |
} |
754 |
if (osc->ntex) { |
755 |
sc->tex = (TexCoord *)emalloc(sizeof(TexCoord) * |
756 |
(osc->ntex+(CHUNKSIZ-1))); |
757 |
memcpy(sc->tex, osc->tex, sizeof(TexCoord)*osc->ntex); |
758 |
sc->ntex = osc->ntex; |
759 |
} |
760 |
if (osc->nnorms) { |
761 |
sc->norm = (Normal *)emalloc(sizeof(Normal) * |
762 |
(osc->nnorms+CHUNKSIZ)); |
763 |
memcpy(sc->norm, osc->norm, sizeof(Normal)*osc->nnorms); |
764 |
sc->nnorms = osc->nnorms; |
765 |
} |
766 |
for (fo = osc->flist; fo != NULL; fo = fo->next) { |
767 |
if ((fo->flags & fltest) != flreq) |
768 |
continue; |
769 |
f = (Face *)emalloc(sizeof(Face) + sizeof(VertEnt)*(fo->nv-3)); |
770 |
memcpy(f, fo, sizeof(Face) + sizeof(VertEnt)*(fo->nv-3)); |
771 |
for (i = 0; i < f->nv; i++) |
772 |
add2facelist(sc, f, i); |
773 |
f->next = sc->flist; |
774 |
sc->flist = f; |
775 |
sc->nfaces++; |
776 |
} |
777 |
deleteUnreferenced(sc); /* jetsam */ |
778 |
return(sc); |
779 |
} |
780 |
|
781 |
/* Add one scene to another, not checking for redundancies */ |
782 |
int |
783 |
addScene(Scene *scdst, const Scene *scsrc) |
784 |
{ |
785 |
VNDX my_vlist[4]; |
786 |
int *vert_map = NULL; |
787 |
int tex_off = 0; |
788 |
int norm_off = 0; |
789 |
VNDX *vlist = my_vlist; |
790 |
int vllen = sizeof(my_vlist)/sizeof(VNDX); |
791 |
int cur_mat = 0; |
792 |
int cur_grp = 0; |
793 |
int fcnt = 0; |
794 |
const Face *f; |
795 |
int i; |
796 |
|
797 |
if ((scdst == NULL) | (scsrc == NULL)) |
798 |
return(-1); |
799 |
if (scsrc->nfaces <= 0) |
800 |
return(0); |
801 |
/* map vertices */ |
802 |
vert_map = (int *)emalloc(sizeof(int)*scsrc->nverts); |
803 |
for (i = 0; i < scsrc->nverts; i++) { |
804 |
const Vertex *v = scsrc->vert + i; |
805 |
if (v->vflist == NULL) { |
806 |
vert_map[i] = -1; |
807 |
continue; |
808 |
} |
809 |
vert_map[i] = addVertex(scdst, v->p[0], v->p[1], v->p[2]); |
810 |
} |
811 |
tex_off = scdst->ntex; /* append texture coords */ |
812 |
if (scsrc->ntex > 0) { |
813 |
scdst->tex = (TexCoord *)erealloc((char *)scdst->tex, |
814 |
sizeof(TexCoord)*(tex_off+scsrc->ntex+(CHUNKSIZ-1))); |
815 |
memcpy(scdst->tex+tex_off, scsrc->tex, |
816 |
sizeof(TexCoord)*scsrc->ntex); |
817 |
} |
818 |
norm_off = scdst->nnorms; /* append normals */ |
819 |
if (scsrc->nnorms > 0) { |
820 |
scdst->norm = (Normal *)erealloc((char *)scdst->norm, |
821 |
sizeof(Normal)*(norm_off+scsrc->nnorms+(CHUNKSIZ-1))); |
822 |
memcpy(scdst->norm+norm_off, scsrc->norm, |
823 |
sizeof(Normal)*scsrc->nnorms); |
824 |
} |
825 |
/* add faces */ |
826 |
scdst->lastgrp = scdst->lastmat = 0; |
827 |
for (f = scsrc->flist; f != NULL; f = f->next) { |
828 |
if (f->grp != cur_grp) |
829 |
setGroup(scdst, scsrc->grpname[cur_grp = f->grp]); |
830 |
if (f->mat != cur_mat) |
831 |
setMaterial(scdst, scsrc->matname[cur_mat = f->mat]); |
832 |
if (f->nv > vllen) { |
833 |
vlist = (VNDX *)( vlist == my_vlist ? |
834 |
emalloc(sizeof(VNDX)*f->nv) : |
835 |
erealloc((char *)vlist, sizeof(VNDX)*f->nv) ); |
836 |
vllen = f->nv; |
837 |
} |
838 |
memset(vlist, 0xff, sizeof(VNDX)*f->nv); |
839 |
for (i = f->nv; i-- > 0; ) { |
840 |
if (f->v[i].vid >= 0) |
841 |
vlist[i][0] = vert_map[f->v[i].vid]; |
842 |
if (f->v[i].tid >= 0) |
843 |
vlist[i][1] = f->v[i].tid + tex_off; |
844 |
if (f->v[i].nid >= 0) |
845 |
vlist[i][2] = f->v[i].nid + norm_off; |
846 |
} |
847 |
fcnt += (addFace(scdst, vlist, f->nv) != NULL); |
848 |
} |
849 |
/* clean up */ |
850 |
if (vlist != my_vlist) efree((char *)vlist); |
851 |
efree((char *)vert_map); |
852 |
return(fcnt); |
853 |
} |
854 |
|
855 |
#define MAXAC 100 |
856 |
|
857 |
/* Transform entire scene */ |
858 |
int |
859 |
xfScene(Scene *sc, int xac, char *xav[]) |
860 |
{ |
861 |
char comm[24+MAXAC*8]; |
862 |
char *cp; |
863 |
XF myxf; |
864 |
FVECT vec; |
865 |
int i; |
866 |
|
867 |
if ((sc == NULL) | (xac <= 0) | (xav == NULL)) |
868 |
return(0); |
869 |
/* compute matrix */ |
870 |
if (xf(&myxf, xac, xav) < xac) |
871 |
return(0); |
872 |
/* transform vertices */ |
873 |
for (i = 0; i < sc->nverts; i++) { |
874 |
VCOPY(vec, sc->vert[i].p); |
875 |
multp3(vec, vec, myxf.xfm); |
876 |
VCOPY(sc->vert[i].p, vec); |
877 |
} |
878 |
/* transform normals */ |
879 |
for (i = 0; i < sc->nnorms; i++) { |
880 |
VCOPY(vec, sc->norm[i]); |
881 |
multv3(vec, vec, myxf.xfm); |
882 |
vec[0] /= myxf.sca; vec[1] /= myxf.sca; vec[2] /= myxf.sca; |
883 |
VCOPY(sc->norm[i], vec); |
884 |
} |
885 |
/* add comment */ |
886 |
cp = strcpy(comm, "Transformed by:"); |
887 |
for (i = 0; i < xac; i++) { |
888 |
while (*cp) cp++; |
889 |
*cp++ = ' '; |
890 |
strcpy(cp, xav[i]); |
891 |
} |
892 |
addComment(sc, comm); |
893 |
return(xac); /* all done */ |
894 |
} |
895 |
|
896 |
/* Ditto, using transform string rather than pre-parsed words */ |
897 |
int |
898 |
xfmScene(Scene *sc, const char *xfm) |
899 |
{ |
900 |
char *xav[MAXAC+1]; |
901 |
int xac, i; |
902 |
|
903 |
if ((sc == NULL) | (xfm == NULL)) |
904 |
return(0); |
905 |
/* skip spaces at beginning */ |
906 |
while (isspace(*xfm)) |
907 |
xfm++; |
908 |
if (!*xfm) |
909 |
return(0); |
910 |
/* parse string into words */ |
911 |
xav[0] = strcpy((char *)emalloc(strlen(xfm)+1), xfm); |
912 |
xac = 1; i = 0; |
913 |
for ( ; ; ) { |
914 |
while (!isspace(xfm[++i])) |
915 |
if (!xfm[i]) |
916 |
break; |
917 |
while (isspace(xfm[i])) |
918 |
xav[0][i++] = '\0'; |
919 |
if (!xfm[i]) |
920 |
break; |
921 |
if (xac >= MAXAC-1) { |
922 |
free(xav[0]); |
923 |
return(0); |
924 |
} |
925 |
xav[xac++] = xav[0] + i; |
926 |
} |
927 |
xav[xac] = NULL; |
928 |
i = xfScene(sc, xac, xav); |
929 |
efree((char *)xav[0]); |
930 |
return(i); |
931 |
} |
932 |
#undef MAXAC |
933 |
|
934 |
/* Delete unreferenced vertices, normals, texture coords */ |
935 |
void |
936 |
deleteUnreferenced(Scene *sc) |
937 |
{ |
938 |
int *vmap; |
939 |
Face *f; |
940 |
int nused, i; |
941 |
/* allocate index map */ |
942 |
if (!sc->nverts) |
943 |
return; |
944 |
i = sc->nverts; |
945 |
if (sc->ntex > i) |
946 |
i = sc->ntex; |
947 |
if (sc->nnorms > i) |
948 |
i = sc->nnorms; |
949 |
vmap = (int *)emalloc(sizeof(int)*i); |
950 |
/* remap positions */ |
951 |
for (i = nused = 0; i < sc->nverts; i++) { |
952 |
if (sc->vert[i].vflist == NULL) { |
953 |
vmap[i] = -1; |
954 |
continue; |
955 |
} |
956 |
if (nused != i) |
957 |
sc->vert[nused] = sc->vert[i]; |
958 |
vmap[i] = nused++; |
959 |
} |
960 |
if (nused == sc->nverts) |
961 |
goto skip_pos; |
962 |
sc->vert = (Vertex *)erealloc((char *)sc->vert, |
963 |
sizeof(Vertex)*(nused+(CHUNKSIZ-1))); |
964 |
sc->nverts = nused; |
965 |
for (f = sc->flist; f != NULL; f = f->next) |
966 |
for (i = f->nv; i--; ) |
967 |
if ((f->v[i].vid = vmap[f->v[i].vid]) < 0) |
968 |
error(CONSISTENCY, |
969 |
"Link error in del_unref_verts()"); |
970 |
skip_pos: |
971 |
/* remap texture coord's */ |
972 |
if (!sc->ntex) |
973 |
goto skip_tex; |
974 |
memset((void *)vmap, 0, sizeof(int)*sc->ntex); |
975 |
for (f = sc->flist; f != NULL; f = f->next) |
976 |
for (i = f->nv; i--; ) |
977 |
if (f->v[i].tid >= 0) |
978 |
vmap[f->v[i].tid] = 1; |
979 |
for (i = nused = 0; i < sc->ntex; i++) { |
980 |
if (!vmap[i]) |
981 |
continue; |
982 |
if (nused != i) |
983 |
sc->tex[nused] = sc->tex[i]; |
984 |
vmap[i] = nused++; |
985 |
} |
986 |
if (nused == sc->ntex) |
987 |
goto skip_tex; |
988 |
sc->tex = (TexCoord *)erealloc((char *)sc->tex, |
989 |
sizeof(TexCoord)*(nused+(CHUNKSIZ-1))); |
990 |
sc->ntex = nused; |
991 |
for (f = sc->flist; f != NULL; f = f->next) |
992 |
for (i = f->nv; i--; ) |
993 |
if (f->v[i].tid >= 0) |
994 |
f->v[i].tid = vmap[f->v[i].tid]; |
995 |
skip_tex: |
996 |
/* remap normals */ |
997 |
if (!sc->nnorms) |
998 |
goto skip_norms; |
999 |
memset((void *)vmap, 0, sizeof(int)*sc->nnorms); |
1000 |
for (f = sc->flist; f != NULL; f = f->next) |
1001 |
for (i = f->nv; i--; ) |
1002 |
if (f->v[i].nid >= 0) |
1003 |
vmap[f->v[i].nid] = 1; |
1004 |
for (i = nused = 0; i < sc->nnorms; i++) { |
1005 |
if (!vmap[i]) |
1006 |
continue; |
1007 |
if (nused != i) |
1008 |
memcpy(sc->norm[nused], sc->norm[i], sizeof(Normal)); |
1009 |
vmap[i] = nused++; |
1010 |
} |
1011 |
if (nused == sc->nnorms) |
1012 |
goto skip_norms; |
1013 |
sc->norm = (Normal *)erealloc((char *)sc->norm, |
1014 |
sizeof(Normal)*(nused+(CHUNKSIZ-1))); |
1015 |
sc->nnorms = nused; |
1016 |
for (f = sc->flist; f != NULL; f = f->next) |
1017 |
for (i = f->nv; i--; ) |
1018 |
if (f->v[i].nid >= 0) |
1019 |
f->v[i].nid = vmap[f->v[i].nid]; |
1020 |
skip_norms: |
1021 |
/* clean up */ |
1022 |
efree((char *)vmap); |
1023 |
} |
1024 |
|
1025 |
/* Free a scene */ |
1026 |
void |
1027 |
freeScene(Scene *sc) |
1028 |
{ |
1029 |
int i; |
1030 |
Face *f; |
1031 |
|
1032 |
if (sc == NULL) |
1033 |
return; |
1034 |
clearComments(sc); |
1035 |
for (i = sc->ngrps; i-- > 0; ) |
1036 |
freeqstr(sc->grpname[i]); |
1037 |
efree((char *)sc->grpname); |
1038 |
for (i = sc->nmats; i-- > 0; ) |
1039 |
freeqstr(sc->matname[i]); |
1040 |
efree((char *)sc->matname); |
1041 |
efree((char *)sc->vert); |
1042 |
efree((char *)sc->tex); |
1043 |
efree((char *)sc->norm); |
1044 |
while ((f = sc->flist) != NULL) { |
1045 |
sc->flist = f->next; |
1046 |
efree((char *)f); |
1047 |
} |
1048 |
efree((char *)sc); |
1049 |
} |