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); |