160 |
|
static float * |
161 |
|
getambdiffs(AMBHEMI *hp) |
162 |
|
{ |
163 |
+ |
const double normf = 1./bright(hp->acoef); |
164 |
|
float *earr = (float *)calloc(hp->ns*hp->ns, sizeof(float)); |
165 |
|
float *ep; |
166 |
|
AMBSAMP *ap; |
174 |
|
for (j = 0; j < hp->ns; j++, ap++, ep++) { |
175 |
|
b = bright(ap[0].v); |
176 |
|
if (i) { /* from above */ |
177 |
< |
d2 = b - bright(ap[-hp->ns].v); |
177 |
> |
d2 = normf*(b - bright(ap[-hp->ns].v)); |
178 |
|
d2 *= d2; |
179 |
|
ep[0] += d2; |
180 |
|
ep[-hp->ns] += d2; |
181 |
|
} |
182 |
|
if (!j) continue; |
183 |
|
/* from behind */ |
184 |
< |
d2 = b - bright(ap[-1].v); |
184 |
> |
d2 = normf*(b - bright(ap[-1].v)); |
185 |
|
d2 *= d2; |
186 |
|
ep[0] += d2; |
187 |
|
ep[-1] += d2; |
188 |
|
if (!i) continue; |
189 |
|
/* diagonal */ |
190 |
< |
d2 = b - bright(ap[-hp->ns-1].v); |
190 |
> |
d2 = normf*(b - bright(ap[-hp->ns-1].v)); |
191 |
|
d2 *= d2; |
192 |
|
ep[0] += d2; |
193 |
|
ep[-hp->ns-1] += d2; |
231 |
|
goto done; /* nothing left to do */ |
232 |
|
nss = *ep/e2rem*cnt + frandom(); |
233 |
|
for (n = 1; n <= nss && ambsample(hp,i,j,n); n++) |
234 |
< |
--cnt; |
234 |
> |
if (!--cnt) goto done; |
235 |
|
e2rem -= *ep++; /* update remainder */ |
236 |
|
} |
237 |
|
done: |
249 |
|
AMBHEMI *hp; |
250 |
|
double d; |
251 |
|
int n, i, j; |
252 |
+ |
/* insignificance check */ |
253 |
+ |
if (bright(rcol) <= FTINY) |
254 |
+ |
return(NULL); |
255 |
|
/* set number of divisions */ |
256 |
|
if (ambacc <= FTINY && |
257 |
|
wt > (d = 0.8*intens(rcol)*r->rweight/(ambdiv*minweight))) |
674 |
|
for (a1 = ang-ang_res; a1 <= ang+ang_res; a1 += ang_step) |
675 |
|
flgs |= 1L<<(int)(16/PI*(a1 + 2.*PI*(a1 < 0))); |
676 |
|
} |
673 |
– |
/* add low-angle incident (< 20deg) */ |
674 |
– |
if (fabs(hp->rp->rod) <= 0.342) { |
675 |
– |
u = -DOT(hp->rp->rdir, uv[0]); |
676 |
– |
v = -DOT(hp->rp->rdir, uv[1]); |
677 |
– |
if ((r0*r0*u*u + r1*r1*v*v) > hp->rp->rot*hp->rp->rot) { |
678 |
– |
ang = atan2a(v, u); |
679 |
– |
ang += 2.*PI*(ang < 0); |
680 |
– |
ang *= 16/PI; |
681 |
– |
if ((ang < .5) | (ang >= 31.5)) |
682 |
– |
flgs |= 0x80000001; |
683 |
– |
else |
684 |
– |
flgs |= 3L<<(int)(ang-.5); |
685 |
– |
} |
686 |
– |
} |
677 |
|
return(flgs); |
678 |
|
} |
679 |
|
|