180 |
|
addcolor(acol, dp->v); |
181 |
|
} |
182 |
|
b = bright(acol); |
183 |
< |
if (b > FTINY) |
183 |
> |
if (b > FTINY) { |
184 |
|
b = ndivs/b; |
185 |
< |
else |
186 |
< |
b = 0.0; |
187 |
< |
if (pg != NULL) { |
188 |
< |
posgradient(pg, div, &hemi); |
189 |
< |
for (i = 0; i < 3; i++) |
190 |
< |
pg[i] *= b; |
185 |
> |
if (pg != NULL) { |
186 |
> |
posgradient(pg, div, &hemi); |
187 |
> |
for (i = 0; i < 3; i++) |
188 |
> |
pg[i] *= b; |
189 |
> |
} |
190 |
> |
if (dg != NULL) { |
191 |
> |
dirgradient(dg, div, &hemi); |
192 |
> |
for (i = 0; i < 3; i++) |
193 |
> |
dg[i] *= b; |
194 |
> |
} |
195 |
> |
} else { |
196 |
> |
if (pg != NULL) |
197 |
> |
for (i = 0; i < 3; i++) |
198 |
> |
pg[i] = 0.0; |
199 |
> |
if (dg != NULL) |
200 |
> |
for (i = 0; i < 3; i++) |
201 |
> |
dg[i] = 0.0; |
202 |
|
} |
192 |
– |
if (dg != NULL) { |
193 |
– |
dirgradient(dg, div, &hemi); |
194 |
– |
for (i = 0; i < 3; i++) |
195 |
– |
dg[i] *= b; |
196 |
– |
} |
203 |
|
free((char *)div); |
204 |
|
} |
205 |
|
b = 1.0/ndivs; |
255 |
|
dp = da; |
256 |
|
for (i = 0; i < hp->nt; i++) |
257 |
|
for (j = 0; j < hp->np; j++) { |
258 |
+ |
#ifdef DEBUG |
259 |
+ |
if (dp->t != i || dp->p != j) |
260 |
+ |
error(CONSISTENCY, |
261 |
+ |
"division order in comperrs"); |
262 |
+ |
#endif |
263 |
|
b = bright(dp[0].v); |
264 |
|
if (i > 0) { /* from above */ |
265 |
|
b2 = bright(dp[-hp->np].v) - b; |