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

Comparing ray/src/rt/o_mesh.c (file contents):
Revision 2.10 by schorsch, Tue Mar 30 16:13:01 2004 UTC vs.
Revision 2.11 by greg, Sat Feb 25 19:49:16 2006 UTC

# Line 58 | Line 58 | OBJREC *o;
58  
59  
60   static int
61 < signed_volume(r, v1, v2)        /* get signed volume for ray and edge */
61 > volume_sign(r, v1, v2)  /* get signed volume for ray and edge */
62   register RAY    *r;
63   int32           v1, v2;
64   {
# Line 70 | Line 70 | int32          v1, v2;
70                  reversed = 1;
71          }
72          ecp = &edge_cache.cache[((v2<<11 ^ v1) & 0x7fffffff) % EDGE_CACHE_SIZ];
73 <        if (ecp->v1i != v1 || ecp->v2i != v2) {
73 >        if ((ecp->v1i != v1) | (ecp->v2i != v2)) {
74                  MESHVERT        tv1, tv2;       /* compute signed volume */
75                  FVECT           v2d;
76                  double          vol;
77 <                if (!getmeshvert(&tv1, edge_cache.mi->msh, v1, MT_V) ||
77 >                if (!getmeshvert(&tv1, edge_cache.mi->msh, v1, MT_V) |
78                              !getmeshvert(&tv2, edge_cache.mi->msh, v2, MT_V))
79                          objerror(edge_cache.o, INTERNAL,
80 <                                "missing mesh vertex in signed_volume");
80 >                                "missing mesh vertex in volume_sign");
81                  VSUB(v2d, tv2.v, r->rorg);
82                  vol = (tv1.v[0] - r->rorg[0]) *
83                                  (v2d[1]*r->rdir[2] - v2d[2]*r->rdir[1]);
# Line 85 | Line 85 | int32          v1, v2;
85                                  (v2d[2]*r->rdir[0] - v2d[0]*r->rdir[2]);
86                  vol += (tv1.v[2] - r->rorg[2]) *
87                                  (v2d[0]*r->rdir[1] - v2d[1]*r->rdir[0]);
88 <                if (vol > .0)
89 <                        ecp->signum = 1;
90 <                else if (vol < .0)
91 <                        ecp->signum = -1;
92 <                else
93 <                        ecp->signum = 0;
88 >                                                /* don't generate 0 */
89 >                ecp->signum = vol > .0 ? 1 : -1;
90                  ecp->v1i = v1;
91                  ecp->v2i = v2;
92          }
# Line 115 | Line 111 | RAY    *r;
111                  if (!getmeshtrivid(tvi, &tmod, curmsh, oset[i]))
112                          objerror(edge_cache.o, INTERNAL,
113                                  "missing triangle vertices in mesh_hit");
114 <                sv1 = signed_volume(r, tvi[0], tvi[1]);
115 <                sv2 = signed_volume(r, tvi[1], tvi[2]);
116 <                sv3 = signed_volume(r, tvi[2], tvi[0]);
117 <                if (sv1 != sv2 || sv2 != sv3)   /* compare volume signs */
118 <                        if (sv1 && sv2 && sv3)
119 <                                continue;
114 >                sv1 = volume_sign(r, tvi[0], tvi[1]);
115 >                sv2 = volume_sign(r, tvi[1], tvi[2]);
116 >                sv3 = volume_sign(r, tvi[2], tvi[0]);
117 >                                                /* compare volume signs */
118 >                if ((sv1 != sv2) | (sv2 != sv3))
119 >                        continue;
120                                                  /* compute intersection */
121                  getmeshvert(&tv[0], curmsh, tvi[0], MT_V);
122                  getmeshvert(&tv[1], curmsh, tvi[1], MT_V);

Diff Legend

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