299 |
|
/* modify threshold */ |
300 |
|
ourthresh = shadthresh / r->rweight; |
301 |
|
/* test for shadows */ |
302 |
< |
nhits = 0; |
303 |
< |
for (sn = 0; sn < ncnts; sn++) { |
302 |
> |
for (nhits = 0, hwt = 0.0, sn = 0; sn < ncnts; |
303 |
> |
hwt += (double)source[scp->sno].nhits / |
304 |
> |
(double)source[scp->sno].ntests, |
305 |
> |
sn++) { |
306 |
|
/* check threshold */ |
307 |
|
if ((sn+nshadcheck>=ncnts ? cntord[sn].brt : |
308 |
|
cntord[sn].brt-cntord[sn+nshadcheck].brt) |
331 |
|
nhits++; |
332 |
|
source[scp->sno].nhits++; |
333 |
|
} |
334 |
< |
/* surface hit rate */ |
335 |
< |
if (sn > 0) |
336 |
< |
hwt = (double)nhits / (double)sn; |
334 |
> |
/* source hit rate */ |
335 |
> |
if (hwt > FTINY) |
336 |
> |
hwt = (double)nhits / hwt; |
337 |
|
else |
338 |
|
hwt = 0.5; |
339 |
|
#ifdef DEBUG |
340 |
< |
sprintf(errmsg, "%d tested, %d untested, %f hit rate\n", |
340 |
> |
sprintf(errmsg, "%d tested, %d untested, %f conditional hit rate\n", |
341 |
|
sn, ncnts-sn, hwt); |
342 |
|
eputs(errmsg); |
343 |
|
#endif |
346 |
|
scp = srccnt + cntord[sn].sndx; |
347 |
|
prob = hwt * (double)source[scp->sno].nhits / |
348 |
|
(double)source[scp->sno].ntests; |
349 |
+ |
if (prob > 1.0) |
350 |
+ |
prob = 1.0; |
351 |
|
scalecolor(scp->val, prob); |
352 |
|
addcolor(r->rcol, scp->val); |
353 |
|
} |