14 |
|
#include "source.h" |
15 |
|
#include "func.h" |
16 |
|
#include "random.h" |
17 |
+ |
#include "pmapmat.h" |
18 |
|
|
19 |
|
#ifndef MAXITER |
20 |
|
#define MAXITER 10 /* maximum # specular ray attempts */ |
98 |
|
scalecolor(ctmp, dtmp); |
99 |
|
addcolor(cval, ctmp); |
100 |
|
} |
101 |
+ |
|
102 |
+ |
if ((ldot < -FTINY) & (np->tdiff > FTINY)) { |
103 |
+ |
/* |
104 |
+ |
* Compute diffuse transmission. |
105 |
+ |
*/ |
106 |
+ |
copycolor(ctmp, np->mcolor); |
107 |
+ |
dtmp = -ldot * omega * np->tdiff * (1.0/PI); |
108 |
+ |
scalecolor(ctmp, dtmp); |
109 |
+ |
addcolor(cval, ctmp); |
110 |
+ |
} |
111 |
+ |
|
112 |
+ |
/* PMAP: skip direct specular refl via ambient bounce if already |
113 |
+ |
* accounted for in photon map */ |
114 |
+ |
if (ambRayInPmap(np->rp)) |
115 |
+ |
return; |
116 |
+ |
|
117 |
|
if (ldot > FTINY && np->specfl&SP_REFL) { |
118 |
|
/* |
119 |
|
* Compute specular reflection coefficient using |
147 |
|
addcolor(cval, ctmp); |
148 |
|
} |
149 |
|
} |
150 |
< |
if ((ldot < -FTINY) & (np->tdiff > FTINY)) { |
134 |
< |
/* |
135 |
< |
* Compute diffuse transmission. |
136 |
< |
*/ |
137 |
< |
copycolor(ctmp, np->mcolor); |
138 |
< |
dtmp = -ldot * omega * np->tdiff * (1.0/PI); |
139 |
< |
scalecolor(ctmp, dtmp); |
140 |
< |
addcolor(cval, ctmp); |
141 |
< |
} |
150 |
> |
|
151 |
|
if (ldot < -FTINY && np->specfl&SP_TRAN) { |
152 |
|
/* |
153 |
|
* Compute specular transmission. Specular transmission |
276 |
|
|
277 |
|
getacoords(&nd); /* set up coordinates */ |
278 |
|
|
279 |
< |
if (nd.specfl & (SP_REFL|SP_TRAN)) |
279 |
> |
/* PMAP: skip indirect specular via ambient bounce if already accounted |
280 |
> |
* for in photon map */ |
281 |
> |
if (!ambRayInPmap(r) && nd.specfl & (SP_REFL|SP_TRAN)) |
282 |
|
agaussamp(&nd); |
283 |
|
|
284 |
|
if (nd.rdiff > FTINY) { /* ambient from this side */ |
289 |
|
multambient(ctmp, r, nd.pnorm); |
290 |
|
addcolor(r->rcol, ctmp); /* add to returned color */ |
291 |
|
} |
292 |
+ |
|
293 |
|
if (nd.tdiff > FTINY) { /* ambient from other side */ |
294 |
|
FVECT bnorm; |
295 |
|
|
312 |
|
return(1); |
313 |
|
} |
314 |
|
|
303 |
– |
|
315 |
|
static void |
316 |
|
getacoords( /* set up coordinate system */ |
317 |
|
ANISODAT *np |
333 |
|
if (normalize(np->v) == 0.0) { |
334 |
|
if (fabs(np->u_alpha - np->v_alpha) > 0.001) |
335 |
|
objerror(np->mp, WARNING, "illegal orientation vector"); |
336 |
< |
getperpendicular(np->u, np->pnorm); /* punting */ |
336 |
> |
getperpendicular(np->u, np->pnorm, 1); /* punting */ |
337 |
|
fcross(np->v, np->pnorm, np->u); |
338 |
|
np->u_alpha = np->v_alpha = sqrt( 0.5 * |
339 |
|
(np->u_alpha*np->u_alpha + np->v_alpha*np->v_alpha) ); |