203 |
|
* Jittering final test reduces image artifacts. |
204 |
|
*/ |
205 |
|
wt = sqrt(e1) + sqrt(e2); |
206 |
< |
if (wt > ambacc*(0.9 + 0.2*frandom())) |
206 |
> |
wt *= .9 + .2*frandom(); |
207 |
> |
if (wt > ambacc) |
208 |
|
continue; |
209 |
|
if (wt <= 1e-3) |
210 |
|
wt = 1e3; |
382 |
|
dnew.n = div[0].n + 1; |
383 |
|
dnew.t = div[0].t; dnew.p = div[0].p; |
384 |
|
b2 = bright(dnew.v)/dnew.n - bright(ar.rcol); |
385 |
< |
b2 = b2*b2 + div[0].k*div[0].n; |
386 |
< |
dnew.k = b2/dnew.n; |
385 |
> |
b2 = b2*b2 + div[0].k*(div[0].n*div[0].n); |
386 |
> |
dnew.k = b2/(dnew.n*dnew.n); |
387 |
|
/* reinsert */ |
388 |
|
for (k = 0; k < ne-1 && dnew.k < div[k+1].k; k++) |
389 |
|
bcopy(&div[k+1], &div[k], sizeof(AMBSAMP)); |