| 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)); |
