BUGS


1
2
3
4
5
6
7
8
9
10
11
12
13

- No checks are made for intermediate over- or underflows, which may occur
  in the middle of the algorithms although the final result may be
  representable.
  For instance, if in the computation of Im (cos(x+i*y)) = -sin(x)*sinh(y)
  an overflow occurs in sinh(y), the value sign(Im (cos(x+i*y))) * inf
  is returned, even if sin(x) is sufficiently close to 0 such that
  sin(x)*sinh(y) is representable. If furthermore an underflow occurred
  in sin(x) (which has not been observed in practice), then the return
  value would be NaN*(+-inf)=NaN.
  As another example, tan is computed as sin/cos; if there is an overflow
  in both sin and cos, then NaN+i*NaN is returned even if the result
  may be representable.