[Bug-gsl] [bug #24812] gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails

[Brian Gough]

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/