| 272 |
|
d = DOT(disp,v->hvec); |
| 273 |
|
d2 = DOT(disp,v->vdir); |
| 274 |
|
ip[0] = 180.0/PI * atan2(d,d2) / v->horiz + 0.5 - v->hoff; |
| 275 |
< |
d = d*d + d2*d2; |
| 276 |
< |
if (d <= FTINY*FTINY) |
| 275 |
> |
d2 = d*d + d2*d2; |
| 276 |
> |
if (d2 <= FTINY*FTINY) |
| 277 |
|
return(VL_BAD); /* at pole */ |
| 278 |
< |
if ((v->vaft > FTINY) & (d >= v->vaft*v->vaft)) |
| 278 |
> |
if ((v->vaft > FTINY) & (d2 >= v->vaft*v->vaft)) |
| 279 |
|
rflags |= VL_BEYOND; |
| 280 |
< |
d = 1.0/sqrt(d); |
| 280 |
> |
d = 1.0/sqrt(d2); |
| 281 |
|
ip[1] = DOT(disp,v->vvec)*d/v->vn2 + 0.5 - v->voff; |
| 282 |
|
ip[2] = VLEN(disp); |
| 283 |
|
ip[2] *= (1.0 - v->vfore*d); |
| 630 |
|
if (rp != NULL && !fgetsresolu(rp, fp)) |
| 631 |
|
mvs.ok = 0; |
| 632 |
|
|
| 633 |
< |
fclose(fp); |
| 633 |
> |
if (fp != stdin) |
| 634 |
> |
fclose(fp); |
| 635 |
|
|
| 636 |
|
return(mvs.ok); |
| 637 |
|
} |