Follow-up Comment #2, bug #24812 (project gsl):
From: Andrew Benson <[email protected]>
To: [email protected]
Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1
Date: Tue, 24 Aug 2010 21:00:48 -0700
[1 <text/plain; us-ascii (7bit)>]
Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2)
the gsl_sf_hyperg_2F1 function crashes with the following arguments:
gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725)
The expected answer (computed using Maple for example) is 0.3104735522
I've attached an example case:
$ g++ test.cpp -o test.x -lgsl -lgslcblas
$ test.x
gsl: hyperg_2F1.c:755: ERROR: error
Default GSL error handler invoked.
Abort
The crash seems to occur only when the first argument is less than or equal
to
-10. For example,
gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725)
crashes, but
gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725)
evaluates correctly to 0.321419346301974773.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?24812>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
