| 179 |
|
} |
| 180 |
|
addcolor(acol, dp->v); |
| 181 |
|
} |
| 182 |
< |
if (pg != NULL) |
| 182 |
> |
b = bright(acol); |
| 183 |
> |
if (b > FTINY) |
| 184 |
> |
b = ndivs/b; |
| 185 |
> |
else |
| 186 |
> |
b = 0.0; |
| 187 |
> |
if (pg != NULL) { |
| 188 |
|
posgradient(pg, div, &hemi); |
| 189 |
< |
if (dg != NULL) |
| 189 |
> |
for (i = 0; i < 3; i++) |
| 190 |
> |
pg[i] *= b; |
| 191 |
> |
} |
| 192 |
> |
if (dg != NULL) { |
| 193 |
|
dirgradient(dg, div, &hemi); |
| 194 |
+ |
for (i = 0; i < 3; i++) |
| 195 |
+ |
dg[i] *= b; |
| 196 |
+ |
} |
| 197 |
|
free((char *)div); |
| 198 |
|
} |
| 199 |
|
b = 1.0/ndivs; |
| 332 |
|
yd += mag0*sinp + mag1*cosp; |
| 333 |
|
} |
| 334 |
|
for (i = 0; i < 3; i++) |
| 335 |
< |
gv[i] = xd*hp->ux[i] + yd*hp->uy[i]; |
| 335 |
> |
gv[i] = (xd*hp->ux[i] + yd*hp->uy[i])/PI; |
| 336 |
|
} |
| 337 |
|
|
| 338 |
|
|
| 364 |
|
yd += mag * sin(phi); |
| 365 |
|
} |
| 366 |
|
for (i = 0; i < 3; i++) |
| 367 |
< |
gv[i] = (xd*hp->ux[i] + yd*hp->uy[i])/(hp->nt*hp->np); |
| 367 |
> |
gv[i] = (xd*hp->ux[i] + yd*hp->uy[i])*PI/(hp->nt*hp->np); |
| 368 |
|
} |