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 |
|
} |