[Top][All Lists]

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

[Bug-gsl] [bug #28718] gsl_sf_hyperg_1F1(-1.5, 1.5, -100) catastrophic c

From: Brian Gough
Subject: [Bug-gsl] [bug #28718] gsl_sf_hyperg_1F1(-1.5, 1.5, -100) catastrophic cancellation in series
Date: Sat, 23 Jan 2010 14:09:26 +0000
User-agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv: Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11


                 Summary: gsl_sf_hyperg_1F1(-1.5,1.5,-100) catastrophic
cancellation in series
                 Project: GNU Scientific Library
            Submitted by: bjg
            Submitted on: Sat 23 Jan 2010 02:09:25 PM GMT
                Category: Accuracy problem
                Severity: 3 - Normal
        Operating System: 
                  Status: Confirmed
             Assigned to: bjg
             Open/Closed: Open
                 Release: 1.13
         Discussion Lock: Any



From: Lorenzo Moneta <address@hidden>
To: <address@hidden>
CC: Weibin Li <address@hidden>, Brian Gough
Subject: Re: [Bug-gsl] gsl_sf_hyperg_1F1
Date: Thu, 21 Jan 2010 16:12:45 +0100


 I have found a similar problem, probably related to the inaccuracy in the
series calculation of 
the hypergeometric function for 
gsl_sf_hyperg_1F1( a < 0, b > 0, x << 0). 

For example : 

  gsl_sf_hyperg_1F1(-1.5, 1.5, - 99.9)  returns the correct value of 

while gsl_sf_hyperg_1F1(-1.5, 1.5, -100.) returns -1.13760e+20 when the
correct value is 456.44  (from Mathematica)

I think the Kummer formula (    1F1(a,b,x) = exp(-x) * 1F1(b-a,b,-x) could be
used for larger value of |x| ( at least works for |x| = 700), while 
the current cut-off is for |x| < 100. 

 Best Regards



Reply to this item at:


  Message sent via/by Savannah

reply via email to

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