[Bug-gsl] [bug #36152] Incorrect asymptotics of spherical Bessel functions

[Jussi Lehtola]

URL:
  <http://savannah.gnu.org/bugs/?36152>

                 Summary: Incorrect asymptotics of spherical Bessel functions
                 Project: GNU Scientific Library
            Submitted by: jussilehtola
            Submitted on: Wed 11 Apr 2012 11:43:45 AM GMT
                Category: Accuracy problem
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
                 Release: 
         Discussion Lock: Any

    _______________________________________________________

Details:

There's something seriously wrong with the spherical Bessel functions. j_n (x)
should go to zero as 1/x when x \to \infty. However, the values computed by
GSL actually diverge!

Test case attached. As you can see, the values turn bad for j_0 at around
x=1e18.



    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: Wed 11 Apr 2012 11:43:45 AM GMT  Name: testbessel.c  Size: 320B   By:
jussilehtola

<http://savannah.gnu.org/bugs/download.php?file_id=25629>

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?36152>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/