| 179 |
|
} |
| 180 |
|
addcolor(acol, dp->v); |
| 181 |
|
} |
| 182 |
< |
if (pg != NULL) |
| 183 |
< |
posgradient(pg, div, &hemi); |
| 184 |
< |
if (dg != NULL) |
| 185 |
< |
dirgradient(dg, div, &hemi); |
| 182 |
> |
b = bright(acol); |
| 183 |
> |
if (b > FTINY) { |
| 184 |
> |
b = ndivs/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 |
> |
} |
| 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; |
| 343 |
|
yd += mag0*sinp + mag1*cosp; |
| 344 |
|
} |
| 345 |
|
for (i = 0; i < 3; i++) |
| 346 |
< |
gv[i] = xd*hp->ux[i] + yd*hp->uy[i]; |
| 346 |
> |
gv[i] = (xd*hp->ux[i] + yd*hp->uy[i])/PI; |
| 347 |
|
} |
| 348 |
|
|
| 349 |
|
|
| 375 |
|
yd += mag * sin(phi); |
| 376 |
|
} |
| 377 |
|
for (i = 0; i < 3; i++) |
| 378 |
< |
gv[i] = (xd*hp->ux[i] + yd*hp->uy[i])/(hp->nt*hp->np); |
| 378 |
> |
gv[i] = (xd*hp->ux[i] + yd*hp->uy[i])*PI/(hp->nt*hp->np); |
| 379 |
|
} |
