| 50 |
|
rw = 1.0; |
| 51 |
|
setscolor(r->rcoef, 1., 1., 1.); |
| 52 |
|
} else { |
| 53 |
< |
rw = sintens((COLORV *)rc); |
| 53 |
> |
rw = sintens(rc); |
| 54 |
|
if (rw > 1.0) |
| 55 |
|
rw = 1.0; /* avoid calculation growth */ |
| 56 |
|
if (rc != r->rcoef) |
| 261 |
|
RAY *r |
| 262 |
|
) |
| 263 |
|
{ |
| 264 |
< |
SCOLOR ce, ca; |
| 264 |
> |
COLOR ce; |
| 265 |
|
double re, ge, be; |
| 266 |
|
|
| 267 |
|
if (intens(r->cext) <= 1./FHUGE) |
| 274 |
|
ge *= 1. - colval(r->albedo,GRN); |
| 275 |
|
be *= 1. - colval(r->albedo,BLU); |
| 276 |
|
} |
| 277 |
< |
setscolor(ce, re<=FTINY ? 1. : re>92. ? 0. : exp(-re), |
| 277 |
> |
setcolor(ce, re<=FTINY ? 1. : re>92. ? 0. : exp(-re), |
| 278 |
|
ge<=FTINY ? 1. : ge>92. ? 0. : exp(-ge), |
| 279 |
|
be<=FTINY ? 1. : be>92. ? 0. : exp(-be)); |
| 280 |
< |
smultscolor(r->rcol, ce); /* path extinction */ |
| 280 |
> |
smultcolor(r->rcol, ce); /* path extinction */ |
| 281 |
|
if (r->crtype & SHADOW || intens(r->albedo) <= FTINY) |
| 282 |
|
return; /* no scattering */ |
| 283 |
|
|
| 284 |
|
/* PMAP: indirect inscattering accounted for by volume photons? */ |
| 285 |
|
if (!volumePhotonMapping) { |
| 286 |
+ |
SCOLOR ca; |
| 287 |
|
setscolor(ca, |
| 288 |
< |
colval(r->albedo,RED)*colval(ambval,RED)*(1.-scolval(ce,RED)), |
| 289 |
< |
colval(r->albedo,GRN)*colval(ambval,GRN)*(1.-scolval(ce,GRN)), |
| 290 |
< |
colval(r->albedo,BLU)*colval(ambval,BLU)*(1.-scolval(ce,BLU))); |
| 291 |
< |
saddscolor(r->rcol, ca); /* ambient in scattering */ |
| 288 |
> |
colval(r->albedo,RED)*colval(ambval,RED)*(1.-colval(ce,RED)), |
| 289 |
> |
colval(r->albedo,GRN)*colval(ambval,GRN)*(1.-colval(ce,GRN)), |
| 290 |
> |
colval(r->albedo,BLU)*colval(ambval,BLU)*(1.-colval(ce,BLU))); |
| 291 |
> |
saddscolor(r->rcol, ca); /* ambient in scattering */ |
| 292 |
|
} |
| 293 |
|
|
| 294 |
|
srcscatter(r); /* source in scattering */ |