At Sat, 27 Oct 2007 19:47:12 -0400,
Chris Mihelich wrote:
> I wish to report and correct the following bugs in sys/ldfrexp.c:
>
> (1) gsl_ldexp computes a power of 2 (p2) that may overflow or
> underflow even though the correct answer is representable;
> consequently the return value x * p2 may be spuriously zero or
> infinity.
>
> (2) gsl_frexp can hang (this is what bit me in the rear at work
> yesterday) if its call to gsl_ldexp spuriously returns zero or
> infinity as in (1). The first while loop divides infinity by 2.0 over
> and over and wonders why it doesn't get any smaller; the second while
> loop has an analogous problem with zero.
Thanks for the bug report and code. I will fix those problems for
the next release.
--
best regards,
Brian Gough
(GSL Maintainer)
Network Theory Ltd,
Commercial support for GSL ---
http://www.network-theory.com/gsl/
