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/