| 88 |
|
if (source[si->sn].sflags & SFLAT) { |
| 89 |
|
si->dom *= sflatform(si->sn, dir); |
| 90 |
|
if (si->dom <= FTINY) { /* behind source */ |
| 91 |
< |
si->sp = si->np; |
| 91 |
> |
si->np = 0; |
| 92 |
|
goto tryagain; |
| 93 |
|
} |
| 94 |
|
si->dom *= (double)(size[SU]*size[SV])/(MAXSPART*MAXSPART); |
| 158 |
|
/* first check point location */ |
| 159 |
|
clrpart(si->spt); |
| 160 |
|
sp = &source[si->sn]; |
| 161 |
< |
rad2 = 1.273 * DOT(sp->ss[SV],sp->ss[SV]); |
| 161 |
> |
rad2 = 1.365 * DOT(sp->ss[SV],sp->ss[SV]); |
| 162 |
|
v[0] = ro[0] - sp->sloc[0]; |
| 163 |
|
v[1] = ro[1] - sp->sloc[1]; |
| 164 |
|
v[2] = ro[2] - sp->sloc[2]; |