156 |
|
} |
157 |
|
|
158 |
|
|
159 |
< |
/* Estimate variance based on relative ambient division differences */ |
159 |
> |
/* Estimate variance based on ambient division differences */ |
160 |
|
static float * |
161 |
|
getambdiffs(AMBHEMI *hp) |
162 |
|
{ |
175 |
|
b = bright(ap[0].v); |
176 |
|
if (i) { /* from above */ |
177 |
|
b1 = bright(ap[-hp->ns].v); |
178 |
< |
d2 = (b - b1)/(b + b1); |
179 |
< |
d2 *= d2*normf; |
178 |
> |
d2 = b - b1; |
179 |
> |
d2 *= d2*normf/(b + b1); |
180 |
|
ep[0] += d2; |
181 |
|
ep[-hp->ns] += d2; |
182 |
|
} |
183 |
|
if (!j) continue; |
184 |
|
/* from behind */ |
185 |
|
b1 = bright(ap[-1].v); |
186 |
< |
d2 = (b - b1)/(b + b1); |
187 |
< |
d2 *= d2*normf; |
186 |
> |
d2 = b - b1; |
187 |
> |
d2 *= d2*normf/(b + b1); |
188 |
|
ep[0] += d2; |
189 |
|
ep[-1] += d2; |
190 |
|
if (!i) continue; |
191 |
|
/* diagonal */ |
192 |
|
b1 = bright(ap[-hp->ns-1].v); |
193 |
< |
d2 = (b - b1)/(b + b1); |
194 |
< |
d2 *= d2*normf; |
193 |
> |
d2 = b - b1; |
194 |
> |
d2 *= d2*normf/(b + b1); |
195 |
|
ep[0] += d2; |
196 |
|
ep[-hp->ns-1] += d2; |
197 |
|
} |