This is the mail archive of the
mailing list for the GSL project.
Floating-point exceptions in the GSL
- From: Thanh Vo <thanhvv at gmail dot com>
- To: gsl-discuss at sourceware dot org, bug-gsl at gnu dot org
- Cc: "Earl T. Barr" <etbarr at ucdavis dot edu>, Zhendong Su <su at cs dot ucdavis dot edu>, Brian Gough <bjg at gnu dot org>
- Date: Thu, 3 Mar 2011 10:12:06 -0800
- Subject: Floating-point exceptions in the GSL
My name's Thanh Vo. My research group is developing a tool to detect
floating point exceptions in numerical code. On Aug 11th 2010, we contacted
the list to report exception triggering inputs we had found in the GSL.
Brian Gough confirmed the exceptions were valid warnings. We have continued
our work since then. In this email, we report more, potentially problematic
For each exception, we describe the name of function, the type of exception
thrown, the inputs that trigger the exception, and the source line that
throws the exception.
Invalid: nu = 0.000000e+00 and x = 0.000000e+00
Location: line 359 at bessel.c
Divide ByZero: x = 0.000000e+00
Location: line 228 at sinint.c
Overflow: x = -5.996154e+153
Location: line 95 at legendre_poly.c
Underflow: x = 8.988466e+307
Location: line 1342 at gamma.c
You can see the full list of potential bugs at
If anyone could confirm whether any of these exceptions are useful warnings
or actual bugs, we would appreciate it.
We are finishing the tool and will release it as open source after
publishing the results.
Looking forward to hearing from you,