112 |
|
scalecolor(ctmp, dtmp); |
113 |
|
addcolor(cval, ctmp); |
114 |
|
} |
115 |
< |
|
115 |
> |
|
116 |
|
if (ldot < -FTINY && ltdiff > FTINY) { |
117 |
|
/* |
118 |
|
* Compute diffuse transmission. |
122 |
|
scalecolor(ctmp, dtmp); |
123 |
|
addcolor(cval, ctmp); |
124 |
|
} |
125 |
< |
|
126 |
< |
/* PMAP: skip direct specular via ambient bounce if already accounted for |
127 |
< |
* in photon map */ |
125 |
> |
|
126 |
|
if (ambRayInPmap(np->rp)) |
127 |
< |
return; |
127 |
> |
return; /* specular already in photon map */ |
128 |
|
|
129 |
|
if (ldot > FTINY && (np->specfl&(SP_REFL|SP_PURE)) == SP_REFL) { |
130 |
|
/* |
182 |
|
{ |
183 |
|
NORMDAT nd; |
184 |
|
double fest; |
187 |
– |
double transtest, transdist; |
188 |
– |
double mirtest, mirdist; |
185 |
|
int hastexture; |
186 |
|
double d; |
187 |
|
COLOR ctmp; |
188 |
|
int i; |
189 |
|
|
190 |
|
/* PMAP: skip transmitted shadow ray if accounted for in photon map */ |
191 |
< |
if (shadowRayInPmap(r)) |
192 |
< |
return(1); |
191 |
> |
/* No longer needed? |
192 |
> |
if (shadowRayInPmap(r) || ambRayInPmap(r)) |
193 |
> |
return(1); */ |
194 |
> |
|
195 |
|
/* easy shadow test */ |
196 |
|
if (r->crtype & SHADOW && m->otype != MAT_TRANS) |
197 |
|
return(1); |
231 |
|
if (nd.pdot < .001) |
232 |
|
nd.pdot = .001; /* non-zero for dirnorm() */ |
233 |
|
multcolor(nd.mcolor, r->pcol); /* modify material color */ |
236 |
– |
mirtest = transtest = 0; |
237 |
– |
mirdist = transdist = r->rot; |
234 |
|
nd.rspec = m->oargs.farg[3]; |
235 |
|
/* compute Fresnel approx. */ |
236 |
|
if (nd.specfl & SP_PURE && nd.rspec >= FRESTHRESH) { |
251 |
|
nd.specfl |= SP_TBLT; |
252 |
|
if (!hastexture || r->crtype & (SHADOW|AMBIENT)) { |
253 |
|
VCOPY(nd.prdir, r->rdir); |
258 |
– |
transtest = 2; |
254 |
|
} else { |
255 |
< |
for (i = 0; i < 3; i++) /* perturb */ |
256 |
< |
nd.prdir[i] = r->rdir[i] - r->pert[i]; |
255 |
> |
/* perturb */ |
256 |
> |
VSUB(nd.prdir, r->rdir, r->pert); |
257 |
|
if (DOT(nd.prdir, r->ron) < -FTINY) |
258 |
|
normalize(nd.prdir); /* OK */ |
259 |
|
else |
262 |
|
} |
263 |
|
} else |
264 |
|
nd.tdiff = nd.tspec = nd.trans = 0.0; |
265 |
+ |
/* diffuse reflection */ |
266 |
+ |
nd.rdiff = 1.0 - nd.trans - nd.rspec; |
267 |
|
/* transmitted ray */ |
268 |
< |
|
272 |
< |
/* PMAP: skip indirect specular trans via ambient bounce if already |
273 |
< |
* accounted for in photon map */ |
274 |
< |
if (!ambRayInPmap(r) && |
275 |
< |
(nd.specfl&(SP_TRAN|SP_PURE|SP_TBLT)) == (SP_TRAN|SP_PURE)) { |
268 |
> |
if ((nd.specfl&(SP_TRAN|SP_PURE|SP_TBLT)) == (SP_TRAN|SP_PURE)) { |
269 |
|
RAY lr; |
270 |
|
copycolor(lr.rcoef, nd.mcolor); /* modified by color */ |
271 |
|
scalecolor(lr.rcoef, nd.tspec); |
274 |
|
rayvalue(&lr); |
275 |
|
multcolor(lr.rcol, lr.rcoef); |
276 |
|
addcolor(r->rcol, lr.rcol); |
277 |
< |
transtest *= bright(lr.rcol); |
278 |
< |
transdist = r->rot + lr.rt; |
277 |
> |
if (nd.tspec >= 1.0-FTINY) { |
278 |
> |
/* completely transparent */ |
279 |
> |
multcolor(lr.mcol, lr.rcoef); |
280 |
> |
copycolor(r->mcol, lr.mcol); |
281 |
> |
r->rmt = r->rot + lr.rmt; |
282 |
> |
r->rxt = r->rot + lr.rxt; |
283 |
> |
} else if (nd.tspec > nd.tdiff + nd.rdiff) |
284 |
> |
r->rxt = r->rot + raydistance(&lr); |
285 |
|
} |
286 |
< |
} else |
288 |
< |
transtest = 0; |
286 |
> |
} |
287 |
|
|
288 |
< |
if (r->crtype & SHADOW) { /* the rest is shadow */ |
291 |
< |
r->rt = transdist; |
288 |
> |
if (r->crtype & SHADOW) /* the rest is shadow */ |
289 |
|
return(1); |
293 |
– |
} |
290 |
|
/* get specular reflection */ |
291 |
|
if (nd.rspec > FTINY) { |
292 |
|
nd.specfl |= SP_REFL; |
313 |
|
checknorm(nd.vrefl); |
314 |
|
} |
315 |
|
/* reflected ray */ |
316 |
< |
/* PMAP: skip indirect specular refl via ambient ray if already accounted |
321 |
< |
* for in photon map */ |
322 |
< |
if (!ambRayInPmap(r) && |
323 |
< |
(nd.specfl&(SP_REFL|SP_PURE|SP_RBLT)) == (SP_REFL|SP_PURE)) { |
316 |
> |
if ((nd.specfl&(SP_REFL|SP_PURE|SP_RBLT)) == (SP_REFL|SP_PURE)) { |
317 |
|
RAY lr; |
318 |
|
if (rayorigin(&lr, REFLECTED, r, nd.scolor) == 0) { |
319 |
|
VCOPY(lr.rdir, nd.vrefl); |
320 |
|
rayvalue(&lr); |
321 |
|
multcolor(lr.rcol, lr.rcoef); |
322 |
+ |
copycolor(r->mcol, lr.rcol); |
323 |
|
addcolor(r->rcol, lr.rcol); |
324 |
+ |
r->rmt = r->rot; |
325 |
|
if (nd.specfl & SP_FLAT && |
326 |
< |
!hastexture | (r->crtype & AMBIENT)) { |
327 |
< |
mirtest = 2.*bright(lr.rcol); |
333 |
< |
mirdist = r->rot + lr.rt; |
334 |
< |
} |
326 |
> |
!hastexture | (r->crtype & AMBIENT)) |
327 |
> |
r->rmt += raydistance(&lr); |
328 |
|
} |
329 |
|
} |
337 |
– |
/* diffuse reflection */ |
338 |
– |
nd.rdiff = 1.0 - nd.trans - nd.rspec; |
330 |
|
|
331 |
|
if (nd.specfl & SP_PURE && nd.rdiff <= FTINY && nd.tdiff <= FTINY) |
332 |
|
return(1); /* 100% pure specular */ |
333 |
|
|
334 |
< |
/* PMAP: skip indirect gaussian via ambient bounce if already accounted |
335 |
< |
* for in photon map */ |
345 |
< |
if (!ambRayInPmap(r)) |
346 |
< |
if (!(nd.specfl & SP_PURE)) |
347 |
< |
gaussamp(&nd); /* checks *BLT flags */ |
334 |
> |
if (!(nd.specfl & SP_PURE)) |
335 |
> |
gaussamp(&nd); /* checks *BLT flags */ |
336 |
|
|
337 |
|
if (nd.rdiff > FTINY) { /* ambient from this side */ |
338 |
|
copycolor(ctmp, nd.mcolor); /* modified by material color */ |
362 |
|
} |
363 |
|
/* add direct component */ |
364 |
|
direct(r, dirnorm, &nd); |
377 |
– |
/* check distance */ |
378 |
– |
d = bright(r->rcol); |
379 |
– |
if (transtest > d) |
380 |
– |
r->rt = transdist; |
381 |
– |
else if (mirtest > d) |
382 |
– |
r->rt = mirdist; |
365 |
|
|
366 |
|
return(1); |
367 |
|
} |