108 |
|
/* worth using? */ |
109 |
|
if (dtmp > FTINY) { |
110 |
|
copycolor(ctmp, np->scolor); |
111 |
< |
dtmp *= omega / np->pdot; |
111 |
> |
dtmp *= omega * sqrt(ldot/np->pdot); |
112 |
|
scalecolor(ctmp, dtmp); |
113 |
|
addcolor(cval, ctmp); |
114 |
|
} |
134 |
|
/* worth using? */ |
135 |
|
if (dtmp > FTINY) { |
136 |
|
copycolor(ctmp, np->mcolor); |
137 |
< |
dtmp *= np->tspec * omega / np->pdot; |
137 |
> |
dtmp *= np->tspec * omega * sqrt(ldot/np->pdot); |
138 |
|
scalecolor(ctmp, dtmp); |
139 |
|
addcolor(cval, ctmp); |
140 |
|
} |
192 |
|
colval(nd.scolor,i) += (1.0-colval(nd.scolor,i))*dtmp; |
193 |
|
nd.rspec += (1.0-nd.rspec)*dtmp; |
194 |
|
/* check threshold */ |
195 |
< |
if (specthresh > FTINY && |
196 |
< |
((specthresh >= 1.-FTINY || |
197 |
< |
specthresh + (.05 - .1*frandom()) > nd.rspec))) |
195 |
> |
if (!(nd.specfl & SP_PURE) && |
196 |
> |
specthresh > FTINY && |
197 |
> |
(specthresh >= 1.-FTINY || |
198 |
> |
specthresh > nd.rspec)) |
199 |
|
nd.specfl |= SP_RBLT; |
200 |
|
/* compute reflected ray */ |
201 |
|
for (i = 0; i < 3; i++) |
222 |
|
if (nd.tspec > FTINY) { |
223 |
|
nd.specfl |= SP_TRAN; |
224 |
|
/* check threshold */ |
225 |
< |
if (specthresh > FTINY && |
226 |
< |
((specthresh >= 1.-FTINY || |
227 |
< |
specthresh + |
227 |
< |
(.05 - .1*frandom()) > nd.tspec))) |
225 |
> |
if (!(nd.specfl & SP_PURE) && specthresh > FTINY && |
226 |
> |
(specthresh >= 1.-FTINY || |
227 |
> |
specthresh > nd.tspec)) |
228 |
|
nd.specfl |= SP_TBLT; |
229 |
|
if (r->crtype & SHADOW || |
230 |
|
DOT(r->pert,r->pert) <= FTINY*FTINY) { |
310 |
|
double rv[2]; |
311 |
|
double d, sinp, cosp; |
312 |
|
register int i; |
313 |
+ |
/* quick test */ |
314 |
+ |
if ((np->specfl & (SP_REFL|SP_RBLT)) != SP_REFL && |
315 |
+ |
(np->specfl & (SP_TRAN|SP_TBLT)) != SP_TRAN) |
316 |
+ |
return; |
317 |
|
/* set up sample coordinates */ |
318 |
|
v[0] = v[1] = v[2] = 0.0; |
319 |
|
for (i = 0; i < 3; i++) |