39 |
|
OBJREC *o; /* mesh object */ |
40 |
|
MESHINST *mi; /* current mesh instance */ |
41 |
|
struct EdgeSide { |
42 |
< |
int4 v1i, v2i; /* vertex indices (lowest first) */ |
42 |
> |
int32 v1i, v2i; /* vertex indices (lowest first) */ |
43 |
|
short signum; /* signed volume */ |
44 |
|
} cache[EDGE_CACHE_SIZ]; |
45 |
|
} edge_cache; |
59 |
|
static int |
60 |
|
signed_volume(r, v1, v2) /* get signed volume for ray and edge */ |
61 |
|
register RAY *r; |
62 |
< |
int4 v1, v2; |
62 |
> |
int32 v1, v2; |
63 |
|
{ |
64 |
|
int reversed = 0; |
65 |
|
register struct EdgeSide *ecp; |
66 |
|
|
67 |
|
if (v1 > v2) { |
68 |
< |
int4 t = v2; v2 = v1; v1 = t; |
68 |
> |
int32 t = v2; v2 = v1; v1 = t; |
69 |
|
reversed = 1; |
70 |
|
} |
71 |
|
ecp = &edge_cache.cache[((v2<<11 ^ v1) & 0x7fffffff) % EDGE_CACHE_SIZ]; |
72 |
|
if (ecp->v1i != v1 || ecp->v2i != v2) { |
73 |
|
MESHVERT tv1, tv2; /* compute signed volume */ |
74 |
+ |
FVECT v2d; |
75 |
|
double vol; |
76 |
|
if (!getmeshvert(&tv1, edge_cache.mi->msh, v1, MT_V) || |
77 |
|
!getmeshvert(&tv2, edge_cache.mi->msh, v2, MT_V)) |
78 |
|
objerror(edge_cache.o, INTERNAL, |
79 |
|
"missing mesh vertex in signed_volume"); |
80 |
+ |
VSUB(v2d, tv2.v, r->rorg); |
81 |
|
vol = (tv1.v[0] - r->rorg[0]) * |
82 |
< |
( (tv2.v[1] - r->rorg[1])*r->rdir[2] - |
81 |
< |
(tv2.v[2] - r->rorg[2])*r->rdir[1] ); |
82 |
> |
(v2d[1]*r->rdir[2] - v2d[2]*r->rdir[1]); |
83 |
|
vol += (tv1.v[1] - r->rorg[1]) * |
84 |
< |
( (tv2.v[2] - r->rorg[2])*r->rdir[0] - |
84 |
< |
(tv2.v[0] - r->rorg[0])*r->rdir[2] ); |
84 |
> |
(v2d[2]*r->rdir[0] - v2d[0]*r->rdir[2]); |
85 |
|
vol += (tv1.v[2] - r->rorg[2]) * |
86 |
< |
( (tv2.v[0] - r->rorg[0])*r->rdir[1] - |
87 |
< |
(tv2.v[1] - r->rorg[1])*r->rdir[0] ); |
86 |
> |
(v2d[0]*r->rdir[1] - v2d[1]*r->rdir[0]); |
87 |
|
if (vol > .0) |
88 |
|
ecp->signum = 1; |
89 |
|
else if (vol < .0) |
102 |
|
OBJECT *oset; |
103 |
|
RAY *r; |
104 |
|
{ |
105 |
< |
int4 tvi[3]; |
105 |
> |
int32 tvi[3]; |
106 |
|
int sv1, sv2, sv3; |
107 |
|
MESHVERT tv[3]; |
108 |
+ |
OBJECT tmod; |
109 |
|
FVECT va, vb, nrm; |
110 |
|
double d; |
111 |
|
int i; |
112 |
|
/* check each triangle */ |
113 |
|
for (i = oset[0]; i > 0; i--) { |
114 |
< |
if (!getmeshtrivid(tvi, curmsh, oset[i])) |
114 |
> |
if (!getmeshtrivid(tvi, &tmod, curmsh, oset[i])) |
115 |
|
objerror(edge_cache.o, INTERNAL, |
116 |
|
"missing triangle vertices in mesh_hit"); |
117 |
|
sv1 = signed_volume(r, tvi[0], tvi[1]); |
152 |
|
RAY rcont; |
153 |
|
int flags; |
154 |
|
MESHVERT tv[3]; |
155 |
< |
FLOAT wt[3]; |
155 |
> |
OBJECT tmod; |
156 |
> |
RREAL wt[3]; |
157 |
|
int i; |
158 |
|
/* get the mesh instance */ |
159 |
|
prep_edge_cache(o); |
171 |
|
return(0); /* missed */ |
172 |
|
if (rcont.rot * curmi->x.f.sca >= r->rot) |
173 |
|
return(0); /* not close enough */ |
173 |
– |
|
174 |
– |
r->robj = objndx(o); /* record new hit */ |
175 |
– |
r->ro = o; |
174 |
|
/* transform ray back */ |
175 |
|
r->rot = rcont.rot * curmi->x.f.sca; |
176 |
|
multp3(r->rop, rcont.rop, curmi->x.f.xfm); |
177 |
|
multv3(r->ron, rcont.ron, curmi->x.f.xfm); |
178 |
|
normalize(r->ron); |
179 |
|
r->rod = -DOT(r->rdir, r->ron); |
180 |
< |
/* compute barycentric weights */ |
181 |
< |
flags = getmeshtri(tv, curmsh, rcont.robj, MT_ALL); |
180 |
> |
/* get triangle */ |
181 |
> |
flags = getmeshtri(tv, &tmod, curmsh, rcont.robj, MT_ALL); |
182 |
|
if (!(flags & MT_V)) |
183 |
|
objerror(o, INTERNAL, "missing mesh vertices in o_mesh"); |
184 |
+ |
r->robj = objndx(o); /* set object and material */ |
185 |
+ |
if (o->omod == OVOID && tmod != OVOID) { |
186 |
+ |
r->ro = getmeshpseudo(curmsh, tmod); |
187 |
+ |
r->rox = &curmi->x; |
188 |
+ |
} else |
189 |
+ |
r->ro = o; |
190 |
+ |
/* compute barycentric weights */ |
191 |
|
if (flags & (MT_N|MT_UV)) |
192 |
|
if (get_baryc(wt, rcont.rop, tv[0].v, tv[1].v, tv[2].v) < 0) { |
193 |
|
objerror(o, WARNING, "bad triangle in o_mesh"); |