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/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
