| 21 |
|
|
| 22 |
|
#include "copyright.h" |
| 23 |
|
|
| 24 |
< |
#include "ray.h" |
| 24 |
> |
#include <string.h> |
| 25 |
|
|
| 26 |
+ |
#include "ray.h" |
| 27 |
|
#include "mesh.h" |
| 27 |
– |
|
| 28 |
|
#include "tmesh.h" |
| 29 |
|
|
| 30 |
|
|
| 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; |
| 52 |
|
/* get mesh instance */ |
| 53 |
|
edge_cache.mi = getmeshinst(edge_cache.o = o, IO_ALL); |
| 54 |
|
/* clear edge cache */ |
| 55 |
< |
bzero((void *)edge_cache.cache, sizeof(edge_cache.cache)); |
| 55 |
> |
memset((void *)edge_cache.cache, '\0', sizeof(edge_cache.cache)); |
| 56 |
|
} |
| 57 |
|
|
| 58 |
|
|
| 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]; |
| 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); |
| 160 |
|
/* copy and transform ray */ |
| 161 |
< |
copystruct(&rcont, r); |
| 161 |
> |
rcont = *r; |
| 162 |
|
multp3(rcont.rorg, r->rorg, curmi->x.b.xfm); |
| 163 |
|
multv3(rcont.rdir, r->rdir, curmi->x.b.xfm); |
| 164 |
|
for (i = 0; i < 3; i++) |
| 171 |
|
return(0); /* missed */ |
| 172 |
|
if (rcont.rot * curmi->x.f.sca >= r->rot) |
| 173 |
|
return(0); /* not close enough */ |
| 172 |
– |
|
| 173 |
– |
r->robj = objndx(o); /* record new hit */ |
| 174 |
– |
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"); |
