1 |
|
#ifndef lint |
2 |
< |
static const char RCSid[] = "$Id$"; |
2 |
> |
static const char RCSid[] = "$Id$"; |
3 |
|
#endif |
4 |
|
/* |
5 |
|
* raytrace.c - routines for tracing and shading rays. |
20 |
|
unsigned long raynum = 0; /* next unique ray number */ |
21 |
|
unsigned long nrays = 0; /* number of calls to localhit */ |
22 |
|
|
23 |
< |
static FLOAT Lambfa[5] = {PI, PI, PI, 0.0, 0.0}; |
23 |
> |
static RREAL Lambfa[5] = {PI, PI, PI, 0.0, 0.0}; |
24 |
|
OBJREC Lamb = { |
25 |
|
OVOID, MAT_PLASTIC, "Lambertian", |
26 |
|
{0, 5, NULL, Lambfa}, NULL, |
31 |
|
static int raymove(), checkhit(); |
32 |
|
static void checkset(); |
33 |
|
|
34 |
– |
#ifndef MAXLOOP |
35 |
– |
#define MAXLOOP 0 /* modifier loop detection */ |
36 |
– |
#endif |
37 |
– |
|
34 |
|
#define RAYHIT (-1) /* return value for intercepted ray */ |
35 |
|
|
36 |
|
|
91 |
|
{ |
92 |
|
r->rno = raynum++; |
93 |
|
r->newcset = r->clipset; |
94 |
+ |
r->hitf = rayhit; |
95 |
|
r->robj = OVOID; |
96 |
|
r->ro = NULL; |
97 |
|
r->rox = NULL; |
98 |
|
r->rt = r->rot = FHUGE; |
99 |
|
r->pert[0] = r->pert[1] = r->pert[2] = 0.0; |
100 |
+ |
r->uv[0] = r->uv[1] = 0.0; |
101 |
|
setcolor(r->pcol, 1.0, 1.0, 1.0); |
102 |
|
setcolor(r->rcol, 0.0, 0.0, 0.0); |
103 |
|
} |
154 |
|
{ |
155 |
|
int gotmat; |
156 |
|
register OBJREC *m; |
159 |
– |
#if MAXLOOP |
160 |
– |
static int depth = 0; |
161 |
– |
/* check for infinite loop */ |
162 |
– |
if (depth++ >= MAXLOOP) |
163 |
– |
objerror(r->ro, USER, "possible modifier loop"); |
164 |
– |
#endif |
157 |
|
r->rt = r->rot; /* set effective ray length */ |
158 |
|
for (gotmat = 0; !gotmat && mod != OVOID; mod = m->omod) { |
159 |
|
m = objptr(mod); |
164 |
|
} |
165 |
|
******/ |
166 |
|
/* hack for irradiance calculation */ |
167 |
< |
if (do_irrad && !(r->crtype & ~(PRIMARY|TRANS))) { |
167 |
> |
if (do_irrad && !(r->crtype & ~(PRIMARY|TRANS)) && |
168 |
> |
m->otype != MAT_CLIP && |
169 |
> |
(ofun[m->otype].flags & (T_M|T_X))) { |
170 |
|
if (irr_ignore(m->otype)) { |
177 |
– |
#if MAXLOOP |
178 |
– |
depth--; |
179 |
– |
#endif |
171 |
|
raytrans(r); |
172 |
|
return(1); |
173 |
|
} |
177 |
|
/* materials call raytexture */ |
178 |
|
gotmat = (*ofun[m->otype].funp)(m, r); |
179 |
|
} |
189 |
– |
#if MAXLOOP |
190 |
– |
depth--; |
191 |
– |
#endif |
180 |
|
return(gotmat); |
181 |
|
} |
182 |
|
|
213 |
|
} |
214 |
|
|
215 |
|
|
216 |
+ |
void |
217 |
|
raytexture(r, mod) /* get material modifiers */ |
218 |
|
RAY *r; |
219 |
< |
int mod; |
219 |
> |
OBJECT mod; |
220 |
|
{ |
221 |
|
register OBJREC *m; |
233 |
– |
#if MAXLOOP |
234 |
– |
static int depth = 0; |
235 |
– |
/* check for infinite loop */ |
236 |
– |
if (depth++ >= MAXLOOP) |
237 |
– |
objerror(r->ro, USER, "modifier loop"); |
238 |
– |
#endif |
222 |
|
/* execute textures and patterns */ |
223 |
|
for ( ; mod != OVOID; mod = m->omod) { |
224 |
|
m = objptr(mod); |
234 |
|
objerror(r->ro, USER, errmsg); |
235 |
|
} |
236 |
|
} |
254 |
– |
#if MAXLOOP |
255 |
– |
depth--; /* end here */ |
256 |
– |
#endif |
237 |
|
} |
238 |
|
|
239 |
|
|
254 |
|
/* compute foreground and background */ |
255 |
|
foremat = backmat = 0; |
256 |
|
/* foreground */ |
257 |
< |
copystruct(&fr, r); |
257 |
> |
fr = *r; |
258 |
|
if (coef > FTINY) |
259 |
|
foremat = rayshade(&fr, fore); |
260 |
|
/* background */ |
261 |
< |
copystruct(&br, r); |
261 |
> |
br = *r; |
262 |
|
if (coef < 1.0-FTINY) |
263 |
|
backmat = rayshade(&br, back); |
264 |
|
/* check for transparency */ |
265 |
< |
if (backmat ^ foremat) |
265 |
> |
if (backmat ^ foremat) { |
266 |
|
if (backmat && coef > FTINY) |
267 |
|
raytrans(&fr); |
268 |
|
else if (foremat && coef < 1.0-FTINY) |
269 |
|
raytrans(&br); |
270 |
+ |
} |
271 |
|
/* mix perturbations */ |
272 |
|
for (i = 0; i < 3; i++) |
273 |
|
r->pert[i] = coef*fr.pert[i] + (1.0-coef)*br.pert[i]; |
390 |
|
} |
391 |
|
|
392 |
|
|
393 |
+ |
void |
394 |
+ |
rayhit(oset, r) /* standard ray hit test */ |
395 |
+ |
OBJECT *oset; |
396 |
+ |
RAY *r; |
397 |
+ |
{ |
398 |
+ |
OBJREC *o; |
399 |
+ |
int i; |
400 |
+ |
|
401 |
+ |
for (i = oset[0]; i > 0; i--) { |
402 |
+ |
o = objptr(oset[i]); |
403 |
+ |
if ((*ofun[o->otype].funp)(o, r)) |
404 |
+ |
r->robj = oset[i]; |
405 |
+ |
} |
406 |
+ |
} |
407 |
+ |
|
408 |
+ |
|
409 |
|
int |
410 |
|
localhit(r, scene) /* check for hit in the octree */ |
411 |
|
register RAY *r; |
464 |
|
} |
465 |
|
cxset[0] = 0; |
466 |
|
raymove(curpos, cxset, sflags, r, scene); |
467 |
< |
return(r->ro != NULL & r->ro != &Aftplane); |
467 |
> |
return((r->ro != NULL) & (r->ro != &Aftplane)); |
468 |
|
} |
469 |
|
|
470 |
|
|
561 |
|
OBJECT *cxs; |
562 |
|
{ |
563 |
|
OBJECT oset[MAXSET+1]; |
567 |
– |
register OBJREC *o; |
568 |
– |
register int i; |
564 |
|
|
565 |
|
objset(oset, cu->cutree); |
566 |
< |
checkset(oset, cxs); /* eliminate double-checking */ |
567 |
< |
for (i = oset[0]; i > 0; i--) { |
568 |
< |
o = objptr(oset[i]); |
569 |
< |
if ((*ofun[o->otype].funp)(o, r)) |
570 |
< |
r->robj = oset[i]; |
576 |
< |
} |
577 |
< |
if (r->ro == NULL) |
566 |
> |
checkset(oset, cxs); /* avoid double-checking */ |
567 |
> |
|
568 |
> |
(*r->hitf)(oset, r); /* test for hit in set */ |
569 |
> |
|
570 |
> |
if (r->robj == OVOID) |
571 |
|
return(0); /* no scores yet */ |
572 |
|
|
573 |
|
return(incube(cu, r->rop)); /* hit OK if in current cube */ |