[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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

From: Brian Gough
Subject: Re: [Bug-gsl] Re: Bug in gsl_sf_hyperg_2F1
Date: Fri, 27 Aug 2010 10:12:21 +0100
User-agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO)

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

Brian Gough

reply via email to

[Prev in Thread] Current Thread [Next in Thread]