Re: [Bug-gsl] Re: Bug in gsl_sf_hyperg_2F1

[Brian Gough]

At Tue, 24 Aug 2010 21:02:26 -0700,
Andrew Benson wrote:
> Some further investigation shows that this combination of parameters is 
> falling through all of the various methods for computing the result in 
> gsl_sf_hyperg_2F1_e and hitting the "/* We give up. */" block at the end. 
> 
> This seems to be related to http://savannah.gnu.org/bugs/?24812 - if I make 
> the changes suggested there (patch attached) then my test code works 
> successfully:
> 
> $ test.x
> = 0.310473552213130111
> 
> I haven't checked that this doesn't break anything else (beyond running "make 
> check" on the newly compiled code). 

Thanks for the bug report.  I've added a test case for it and applied
the earlier patch (with a further correction because one of the a->ap
changes in the function call arguments was also missing).  

It should work now from revision 4676 in the repository
http://bzr.savannah.gnu.org/r/gsl/trunk.

-- 
Brian Gough