| 160 |
|
if (dtmp > FTINY*FTINY) { |
| 161 |
|
dtmp1 = DOT(h,np->pnorm); |
| 162 |
|
dtmp = 1.0 - dtmp1*dtmp1/dtmp; |
| 163 |
< |
if (dtmp > FTINY*FTINY) { |
| 164 |
< |
dtmp1 = DOT(h,np->u); |
| 165 |
< |
dtmp1 *= dtmp1 / au2; |
| 166 |
< |
dtmp2 = DOT(h,np->v); |
| 167 |
< |
dtmp2 *= dtmp2 / av2; |
| 168 |
< |
dtmp = (dtmp1 + dtmp2) / dtmp; |
| 169 |
< |
} |
| 163 |
> |
} |
| 164 |
> |
if (dtmp > FTINY*FTINY) { |
| 165 |
> |
dtmp1 = DOT(h,np->u); |
| 166 |
> |
dtmp1 *= dtmp1 / au2; |
| 167 |
> |
dtmp2 = DOT(h,np->v); |
| 168 |
> |
dtmp2 *= dtmp2 / av2; |
| 169 |
> |
dtmp = (dtmp1 + dtmp2) / dtmp; |
| 170 |
> |
dtmp = exp(-dtmp); |
| 171 |
|
} else |
| 172 |
< |
dtmp = 0.0; |
| 172 |
> |
dtmp = 1.0; |
| 173 |
|
/* Gaussian */ |
| 174 |
< |
dtmp = exp(-dtmp) * (1.0/PI) * sqrt(-ldot/(np->pdot*au2*av2)); |
| 174 |
> |
dtmp *= (1.0/PI) * sqrt(-ldot/(np->pdot*au2*av2)); |
| 175 |
|
/* worth using? */ |
| 176 |
|
if (dtmp > FTINY) { |
| 177 |
|
copyscolor(sctmp, np->mcolor); |
