[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [bug #24812] gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails
From: |
Brian Gough |
Subject: |
[Bug-gsl] [bug #24812] gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails |
Date: |
Fri, 27 Aug 2010 08:42:17 +0000 |
User-agent: |
Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 |
Follow-up Comment #2, bug #24812 (project gsl):
From: Andrew Benson <address@hidden>
To: address@hidden
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] [bug #24812] gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails,
Brian Gough <=