101 |
|
if (ldot > FTINY && (np->specfl&(SP_REFL|SP_BADU)) == SP_REFL) { |
102 |
|
/* |
103 |
|
* Compute specular reflection coefficient using |
104 |
< |
* anisotropic gaussian distribution model. |
104 |
> |
* anisotropic Gaussian distribution model. |
105 |
|
*/ |
106 |
|
/* add source width if flat */ |
107 |
|
if (np->specfl & SP_FLAT) |
119 |
|
dtmp1 *= dtmp1 / au2; |
120 |
|
dtmp2 = DOT(np->v, h); |
121 |
|
dtmp2 *= dtmp2 / av2; |
122 |
< |
/* gaussian */ |
122 |
> |
/* new W-G-M-D model */ |
123 |
|
dtmp = DOT(np->pnorm, h); |
124 |
< |
dtmp = (dtmp1 + dtmp2) / (dtmp*dtmp); |
125 |
< |
dtmp = exp(-dtmp) / (4.0*PI * np->pdot * sqrt(au2*av2)); |
124 |
> |
dtmp *= dtmp; |
125 |
> |
dtmp1 = (dtmp1 + dtmp2) / dtmp; |
126 |
> |
dtmp = exp(-dtmp1) * DOT(h,h) / |
127 |
> |
(PI * dtmp*dtmp * sqrt(au2*av2)); |
128 |
|
/* worth using? */ |
129 |
|
if (dtmp > FTINY) { |
130 |
|
copycolor(ctmp, np->scolor); |
131 |
< |
dtmp *= omega; |
131 |
> |
dtmp *= ldot * omega; |
132 |
|
scalecolor(ctmp, dtmp); |
133 |
|
addcolor(cval, ctmp); |
134 |
|
} |
168 |
|
} |
169 |
|
} else |
170 |
|
dtmp = 0.0; |
171 |
< |
/* gaussian */ |
171 |
> |
/* Gaussian */ |
172 |
|
dtmp = exp(-dtmp) * (1.0/PI) * sqrt(-ldot/(np->pdot*au2*av2)); |
173 |
|
/* worth using? */ |
174 |
|
if (dtmp > FTINY) { |
341 |
|
|
342 |
|
|
343 |
|
static void |
344 |
< |
agaussamp( /* sample anisotropic gaussian specular */ |
344 |
> |
agaussamp( /* sample anisotropic Gaussian specular */ |
345 |
|
RAY *r, |
346 |
|
register ANISODAT *np |
347 |
|
) |