From: Peter Vanroose <Peter.V<anroose@es...>  20020208 20:05:57

> cos_angle to return (1+eps). vcl_acos(1+eps) returns NaN. > > Alternatives > 1. rewrite cos_angle so that it is immune to this problem, However, we can't > see how  any suggestions welcome. Maybe reimplement vcl_acos such that it always returns 0 when its argument is > 1 (and Pi when its argument is < 1) ? Or replace the line return vcl_acos( abs_r( cos_angle(a, b) ) ); in vcl_vector.txx by abs_r c = abs_r( cos_angle(a, b) ); return c<1 ? vcl_acos(c) : vcl_acos(1.0); I've tried out the latter and it seems to work fine. Peter. 