ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/src/common/mesh.c
(Generate patch)

Comparing src/common/mesh.c (file contents):
Revision 2.6 by greg, Wed Apr 23 00:52:33 2003 UTC vs.
Revision 2.24 by greg, Tue Nov 6 01:04:23 2012 UTC

# Line 5 | Line 5 | static const char RCSid[] = "$Id$";
5   * Mesh support routines
6   */
7  
8 < #include "standard.h"
8 > #include <string.h>
9 >
10 > #include "rtio.h"
11 > #include "rtmath.h"
12 > #include "rterror.h"
13 > #include "paths.h"
14   #include "octree.h"
15   #include "object.h"
16   #include "otypes.h"
# Line 14 | Line 19 | static const char RCSid[] = "$Id$";
19   /* An encoded mesh vertex */
20   typedef struct {
21          int             fl;
22 <        uint4           xyz[3];
23 <        int4            norm;
24 <        uint4           uv[2];
22 >        uint32          xyz[3];
23 >        int32           norm;
24 >        uint32          uv[2];
25   } MCVERT;
26  
27   #define  MPATCHBLKSIZ   128             /* patch allocation block size */
# Line 27 | Line 32 | static MESH    *mlist = NULL;          /* list of loaded meshes *
32  
33  
34   static unsigned long
35 < cvhash(cvp)                             /* hash an encoded vertex */
31 < MCVERT  *cvp;
35 > cvhash(const char *p)                   /* hash an encoded vertex */
36   {
37 +        const MCVERT    *cvp = (const MCVERT *)p;
38          unsigned long   hval;
39          
40          if (!(cvp->fl & MT_V))
# Line 44 | Line 49 | MCVERT *cvp;
49  
50  
51   static int
52 < cvcmp(v1, v2)                           /* compare encoded vertices */
48 < register MCVERT *v1, *v2;
52 > cvcmp(const char *vv1, const char *vv2)         /* compare encoded vertices */
53   {
54 +        const MCVERT    *v1 = (const MCVERT *)vv1, *v2 = (const MCVERT *)vv2;
55          if (v1->fl != v2->fl)
56                  return(1);
57          if (v1->xyz[0] != v2->xyz[0])
# Line 68 | Line 73 | register MCVERT        *v1, *v2;
73  
74  
75   MESH *
76 < getmesh(mname, flags)                   /* get mesh data */
77 < char    *mname;
78 < int     flags;
76 > getmesh(                                /* get new mesh data reference */
77 >        char    *mname,
78 >        int     flags
79 > )
80   {
81          char  *pathname;
82 <        register MESH  *ms;
82 >        MESH  *ms;
83  
84          flags &= IO_LEGAL;
85          for (ms = mlist; ms != NULL; ms = ms->next)
86                  if (!strcmp(mname, ms->name)) {
87 <                        if ((ms->ldflags & flags) == flags) {
88 <                                ms->nref++;
83 <                                return(ms);             /* loaded */
84 <                        }
85 <                        break;                  /* load the rest */
87 >                        ms->nref++;     /* increase reference count */
88 >                        break;
89                  }
90 <        if (ms == NULL) {
90 >        if (ms == NULL) {               /* load first time */
91                  ms = (MESH *)calloc(1, sizeof(MESH));
92                  if (ms == NULL)
93                          error(SYSTEM, "out of memory in getmesh");
# Line 94 | Line 97 | int    flags;
97                  ms->next = mlist;
98                  mlist = ms;
99          }
100 <        if ((pathname = getpath(mname, getlibpath(), R_OK)) == NULL) {
100 >        if ((pathname = getpath(mname, getrlibpath(), R_OK)) == NULL) {
101                  sprintf(errmsg, "cannot find mesh file \"%s\"", mname);
102                  error(USER, errmsg);
103          }
# Line 106 | Line 109 | int    flags;
109  
110  
111   MESHINST *
112 < getmeshinst(o, flags)                   /* create mesh instance */
113 < OBJREC  *o;
114 < int     flags;
112 > getmeshinst(                            /* create mesh instance */
113 >        OBJREC  *o,
114 >        int     flags
115 > )
116   {
117 <        register MESHINST  *ins;
117 >        MESHINST  *ins;
118  
119          flags &= IO_LEGAL;
120          if ((ins = (MESHINST *)o->os) == NULL) {
# Line 128 | Line 132 | int    flags;
132                  ins->msh = NULL;
133                  o->os = (char *)ins;
134          }
135 <        if (ins->msh == NULL || (ins->msh->ldflags & flags) != flags)
135 >        if (ins->msh == NULL)
136                  ins->msh = getmesh(o->oargs.sarg[0], flags);
137 +        else if ((flags &= ~ins->msh->ldflags))
138 +                readmesh(ins->msh,
139 +                        getpath(o->oargs.sarg[0], getrlibpath(), R_OK),
140 +                                flags);
141          return(ins);
142   }
143  
144  
145   int
146 < getmeshtrivid(tvid, mo, mp, ti)         /* get triangle vertex ID's */
147 < int4    tvid[3];
148 < OBJECT  *mo;
149 < MESH    *mp;
142 < OBJECT  ti;
146 > nextmeshtri(                            /* get next triangle ID */
147 >        OBJECT *tip,
148 >        MESH *mp
149 > )
150   {
151 +        int             pn;
152 +        MESHPATCH       *pp;
153 +
154 +        if (*tip == OVOID) {                    /* check for first index */
155 +                *tip = 0;
156 +                return(mp->npatches > 0);       /* assumes 1 local triangle */
157 +        }
158 +        pn = *tip >> 10;
159 +        if (pn >= mp->npatches)                 /* past end? */
160 +                return(0);
161 +        pp = &mp->patch[pn];
162 +        if (!(*tip & 0x200)) {                  /* local triangle? */
163 +                if ((*tip & 0x1ff) < pp->ntris-1) {
164 +                        ++*tip;
165 +                        return(1);
166 +                }
167 +                *tip &= ~0x1ff;                 /* move on to single-joiners */
168 +                *tip |= 0x200;
169 +                if (pp->nj1tris)                /* is there at least one? */
170 +                        return(1);
171 +        }
172 +        if (!(*tip & 0x100)) {                  /* single joiner? */
173 +                if ((*tip & 0xff) < pp->nj1tris-1) {
174 +                        ++*tip;
175 +                        return(1);
176 +                }
177 +                *tip &= ~0xff;                  /* move on to double-joiners */
178 +                *tip |= 0x100;
179 +                if (pp->nj2tris)                /* is there one? */
180 +                        return(1);
181 +        }
182 +        if ((*tip & 0xff) < pp->nj2tris-1) {    /* double-joiner? */
183 +                ++*tip;
184 +                return(1);
185 +        }
186 +        *tip = ++pn << 10;                      /* first in next patch */
187 +        return(pn < mp->npatches);
188 + }
189 +
190 + int
191 + getmeshtrivid(                          /* get triangle vertex ID's */
192 +        int32   tvid[3],
193 +        OBJECT  *mo,
194 +        MESH    *mp,
195 +        OBJECT  ti
196 + )
197 + {
198          int             pn = ti >> 10;
199          MESHPATCH       *pp;
200  
# Line 196 | Line 250 | OBJECT ti;
250  
251  
252   int
253 < getmeshvert(vp, mp, vid, what)  /* get triangle vertex from ID */
254 < MESHVERT        *vp;
255 < MESH            *mp;
256 < int4            vid;
257 < int             what;
253 > getmeshvert(                            /* get triangle vertex from ID */
254 >        MESHVERT        *vp,
255 >        MESH            *mp,
256 >        int32           vid,
257 >        int             what
258 > )
259   {
260          int             pn = vid >> 8;
261          MESHPATCH       *pp;
262          double          vres;
263 <        register int    i;
263 >        int     i;
264          
265          vp->fl = 0;
266          if (pn >= mp->npatches)
# Line 240 | Line 295 | int            what;
295  
296  
297   OBJREC *
298 < getmeshpseudo(mp, mo)           /* get mesh pseudo object for material */
299 < MESH    *mp;
300 < OBJECT  mo;
298 > getmeshpseudo(                  /* get mesh pseudo object for material */
299 >        MESH    *mp,
300 >        OBJECT  mo
301 > )
302   {
303          if (mo < mp->mat0 || mo >= mp->mat0 + mp->nmats)
304                  error(INTERNAL, "modifier out of range in getmeshpseudo");
305          if (mp->pseudo == NULL) {
306 <                register int    i;
306 >                int     i;
307                  mp->pseudo = (OBJREC *)calloc(mp->nmats, sizeof(OBJREC));
308                  if (mp->pseudo == NULL)
309                          error(SYSTEM, "out of memory in getmeshpseudo");
# Line 262 | Line 318 | OBJECT mo;
318  
319  
320   int
321 < getmeshtri(tv, mo, mp, ti, wha) /* get triangle vertices */
322 < MESHVERT        tv[3];
323 < OBJECT          *mo;
324 < MESH            *mp;
325 < OBJECT          ti;
326 < int             wha;
321 > getmeshtri(                     /* get triangle vertices */
322 >        MESHVERT        tv[3],
323 >        OBJECT          *mo,
324 >        MESH            *mp,
325 >        OBJECT          ti,
326 >        int             wha
327 > )
328   {
329 <        int4    tvid[3];
329 >        int32   tvid[3];
330  
331          if (!getmeshtrivid(tvid, mo, mp, ti))
332                  return(0);
# Line 282 | Line 339 | int            wha;
339   }
340  
341  
342 < int4
343 < addmeshvert(mp, vp)             /* find/add a mesh vertex */
344 < register MESH   *mp;
345 < MESHVERT        *vp;
342 > int32
343 > addmeshvert(                    /* find/add a mesh vertex */
344 >        MESH    *mp,
345 >        MESHVERT        *vp
346 > )
347   {
348          LUENT           *lvp;
349          MCVERT          cv;
350 <        register int    i;
350 >        int     i;
351  
352          if (!(vp->fl & MT_V))
353                  return(-1);
# Line 299 | Line 357 | MESHVERT       *vp;
357                          return(-1);
358                  if (vp->v[i] >= mp->mcube.cuorg[i] + mp->mcube.cusize)
359                          return(-1);
360 <                cv.xyz[i] = (uint4)(4294967296. *
360 >                cv.xyz[i] = (uint32)(4294967296. *
361                                  (vp->v[i] - mp->mcube.cuorg[i]) /
362                                  mp->mcube.cusize);
363          }
364 <        if (vp->fl & MT_N)
364 >        if (vp->fl & MT_N)              /* assumes normalized! */
365                  cv.norm = encodedir(vp->n);
366          if (vp->fl & MT_UV)
367                  for (i = 0; i < 2; i++) {
# Line 311 | Line 369 | MESHVERT       *vp;
369                                  return(-1);
370                          if (vp->uv[i] >= mp->uvlim[1][i])
371                                  return(-1);
372 <                        cv.uv[i] = (uint4)(4294967296. *
372 >                        cv.uv[i] = (uint32)(4294967296. *
373                                          (vp->uv[i] - mp->uvlim[0][i]) /
374                                          (mp->uvlim[1][i] - mp->uvlim[0][i]));
375                  }
# Line 328 | Line 386 | MESHVERT       *vp;
386          if (lvp == NULL)
387                  goto nomem;
388          if (lvp->key == NULL) {
389 <                lvp->key = (char *)malloc(sizeof(MCVERT)+sizeof(int4));
390 <                bcopy((void *)&cv, (void *)lvp->key, sizeof(MCVERT));
389 >                lvp->key = (char *)malloc(sizeof(MCVERT)+sizeof(int32));
390 >                memcpy((void *)lvp->key, (void *)&cv, sizeof(MCVERT));
391          }
392          if (lvp->data == NULL) {        /* new vertex */
393 <                register MESHPATCH      *pp;
393 >                MESHPATCH       *pp;
394                  if (mp->npatches <= 0) {
395                          mp->patch = (MESHPATCH *)calloc(MPATCHBLKSIZ,
396                                          sizeof(MESHPATCH));
# Line 345 | Line 403 | MESHVERT       *vp;
403                                                  (void *)mp->patch,
404                                          (mp->npatches + MPATCHBLKSIZ)*
405                                                  sizeof(MESHPATCH));
406 <                                bzero((void *)(mp->patch + mp->npatches),
406 >                                memset((void *)(mp->patch + mp->npatches), '\0',
407                                          MPATCHBLKSIZ*sizeof(MESHPATCH));
408                          }
409                          if (mp->npatches++ >= 1L<<22)
# Line 353 | Line 411 | MESHVERT       *vp;
411                  }
412                  pp = &mp->patch[mp->npatches-1];
413                  if (pp->xyz == NULL) {
414 <                        pp->xyz = (uint4 (*)[3])calloc(256, 3*sizeof(int4));
414 >                        pp->xyz = (uint32 (*)[3])calloc(256, 3*sizeof(int32));
415                          if (pp->xyz == NULL)
416                                  goto nomem;
417                  }
# Line 361 | Line 419 | MESHVERT       *vp;
419                          pp->xyz[pp->nverts][i] = cv.xyz[i];
420                  if (cv.fl & MT_N) {
421                          if (pp->norm == NULL) {
422 <                                pp->norm = (int4 *)calloc(256, sizeof(int4));
422 >                                pp->norm = (int32 *)calloc(256, sizeof(int32));
423                                  if (pp->norm == NULL)
424                                          goto nomem;
425                          }
# Line 369 | Line 427 | MESHVERT       *vp;
427                  }
428                  if (cv.fl & MT_UV) {
429                          if (pp->uv == NULL) {
430 <                                pp->uv = (uint4 (*)[2])calloc(256,
431 <                                                2*sizeof(uint4));
430 >                                pp->uv = (uint32 (*)[2])calloc(256,
431 >                                                2*sizeof(uint32));
432                                  if (pp->uv == NULL)
433                                          goto nomem;
434                          }
# Line 379 | Line 437 | MESHVERT       *vp;
437                  }
438                  pp->nverts++;
439                  lvp->data = lvp->key + sizeof(MCVERT);
440 <                *(int4 *)lvp->data = (mp->npatches-1) << 8 | (pp->nverts-1);
440 >                *(int32 *)lvp->data = (mp->npatches-1) << 8 | (pp->nverts-1);
441          }
442 <        return(*(int4 *)lvp->data);
442 >        return(*(int32 *)lvp->data);
443   nomem:
444          error(SYSTEM, "out of memory in addmeshvert");
445          return(-1);
# Line 389 | Line 447 | nomem:
447  
448  
449   OBJECT
450 < addmeshtri(mp, tv, mo)          /* add a new mesh triangle */
451 < MESH            *mp;
452 < MESHVERT        tv[3];
453 < OBJECT          mo;
450 > addmeshtri(                     /* add a new mesh triangle */
451 >        MESH            *mp,
452 >        MESHVERT        tv[3],
453 >        OBJECT          mo
454 > )
455   {
456 <        int4                    vid[3], t;
456 >        int32                   vid[3], t;
457          int                     pn[3], i;
458 <        register MESHPATCH      *pp;
458 >        MESHPATCH       *pp;
459  
460          if (!(tv[0].fl & tv[1].fl & tv[2].fl & MT_V))
461                  return(OVOID);
# Line 407 | Line 466 | OBJECT         mo;
466                  pn[i] = vid[i] >> 8;
467          }
468                                  /* normalize material index */
469 <        if (mo != OVOID)
469 >        if (mo != OVOID) {
470                  if ((mo -= mp->mat0) >= mp->nmats)
471                          mp->nmats = mo+1;
472                  else if (mo < 0)
473                          error(INTERNAL, "modifier range error in addmeshtri");
474 +        }
475                                  /* assign triangle */
476          if (pn[0] == pn[1] && pn[1] == pn[2]) { /* local case */
477                  pp = &mp->patch[pn[0]];
# Line 428 | Line 488 | OBJECT         mo;
488                          if (pp->ntris == 0)
489                                  pp->solemat = mo;
490                          else if (pp->trimat == NULL && mo != pp->solemat) {
491 <                                pp->trimat = (int2 *)malloc(
492 <                                                512*sizeof(int2));
491 >                                pp->trimat = (int16 *)malloc(
492 >                                                512*sizeof(int16));
493                                  if (pp->trimat == NULL)
494                                          goto nomem;
495                                  for (i = pp->ntris; i--; )
# Line 485 | Line 545 | nomem:
545  
546  
547   char *
548 < checkmesh(mp)                           /* validate mesh data */
489 < register MESH   *mp;
548 > checkmesh(MESH *mp)                     /* validate mesh data */
549   {
550          static char     embuf[128];
551          int             nouvbounds = 1;
552 <        register int    i;
552 >        int     i;
553                                          /* basic checks */
554          if (mp == NULL)
555                  return("NULL mesh pointer");
# Line 530 | Line 589 | register MESH  *mp;
589                  if (mp->npatches <= 0)
590                          error(WARNING, "no patches in mesh");
591                  for (i = 0; i < mp->npatches; i++) {
592 <                        register MESHPATCH      *pp = &mp->patch[i];
592 >                        MESHPATCH       *pp = &mp->patch[i];
593                          if (pp->nverts <= 0)
594                                  error(WARNING, "no vertices in patch");
595                          else {
# Line 539 | Line 598 | register MESH  *mp;
598                                  if (nouvbounds && pp->uv != NULL)
599                                          return("unreferenced uv coordinates");
600                          }
542 –                        if (pp->ntris + pp->nj1tris + pp->nj2tris <= 0)
543 –                                error(WARNING, "no triangles in patch");
601                          if (pp->ntris > 0 && pp->tri == NULL)
602                                  return("missing patch triangle list");
603                          if (pp->nj1tris > 0 && pp->j1tri == NULL)
# Line 554 | Line 611 | register MESH  *mp;
611  
612  
613   static void
614 < tallyoctree(ot, ecp, lcp, ocp)  /* tally octree size */
615 < OCTREE  ot;
616 < int     *ecp, *lcp, *ocp;
614 > tallyoctree(                    /* tally octree size */
615 >        OCTREE  ot,
616 >        int     *ecp,
617 >        int     *lcp,
618 >        int     *ocp
619 > )
620   {
621          int     i;
622  
# Line 577 | Line 637 | int    *ecp, *lcp, *ocp;
637  
638  
639   void
640 < printmeshstats(ms, fp)          /* print out mesh statistics */
641 < MESH    *ms;
642 < FILE    *fp;
640 > printmeshstats(                 /* print out mesh statistics */
641 >        MESH    *ms,
642 >        FILE    *fp
643 > )
644   {
645          int     lfcnt=0, lecnt=0, locnt=0;
646          int     vcnt=0, ncnt=0, uvcnt=0;
# Line 589 | Line 650 | FILE   *fp;
650          
651          tallyoctree(ms->mcube.cutree, &lecnt, &lfcnt, &locnt);
652          for (i = 0; i < ms->npatches; i++) {
653 <                register MESHPATCH      *pp = &ms->patch[i];
653 >                MESHPATCH       *pp = &ms->patch[i];
654                  vcnt += pp->nverts;
655                  if (pp->norm != NULL) {
656                          for (j = pp->nverts; j--; )
# Line 608 | Line 669 | FILE   *fp;
669                  t2cnt += pp->nj2tris;
670          }
671          fprintf(fp, "Mesh statistics:\n");
672 <        fprintf(fp, "\t%d materials\n", ms->nmats);
672 >        fprintf(fp, "\t%ld materials\n", (long)ms->nmats);
673          fprintf(fp, "\t%d patches (%.2f MBytes)\n", ms->npatches,
674                          (ms->npatches*sizeof(MESHPATCH) +
675 <                        vcnt*3*sizeof(uint4) +
676 <                        nscnt*sizeof(int4) +
677 <                        uvscnt*2*sizeof(uint4) +
675 >                        vcnt*3*sizeof(uint32) +
676 >                        nscnt*sizeof(int32) +
677 >                        uvscnt*2*sizeof(uint32) +
678                          tcnt*sizeof(struct PTri) +
679                          t1cnt*sizeof(struct PJoin1) +
680                          t2cnt*sizeof(struct PJoin2))/(1024.*1024.));
# Line 631 | Line 692 | FILE   *fp;
692  
693  
694   void
695 < freemesh(ms)                    /* free mesh data */
635 < register MESH   *ms;
695 > freemesh(MESH *ms)              /* free mesh data */
696   {
697          MESH    mhead;
698          MESH    *msp;
# Line 660 | Line 720 | register MESH  *ms;
720          octfree(ms->mcube.cutree);
721          lu_done(&ms->lut);
722          if (ms->npatches > 0) {
723 <                register MESHPATCH      *pp = ms->patch + ms->npatches;
723 >                MESHPATCH       *pp = ms->patch + ms->npatches;
724                  while (pp-- > ms->patch) {
725                          if (pp->j2tri != NULL)
726                                  free((void *)pp->j2tri);
# Line 684 | Line 744 | register MESH  *ms;
744  
745  
746   void
747 < freemeshinst(o)                 /* free mesh instance */
688 < OBJREC  *o;
747 > freemeshinst(OBJREC *o)         /* free mesh instance */
748   {
749          if (o->os == NULL)
750                  return;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines