bugGNU Scientific Library - Bugs: bug #36152, Incorrect asymptotics of spherical...

 
 
Show feedback again

You are not allowed to post comments on this tracker with your current authentification level.

bug #36152: Incorrect asymptotics of spherical Bessel functions

Submitted by:  Jussi Lehtola <jussilehtola>
Submitted on:  Wed 11 Apr 2012 11:43:45 AM UTC  
 
Category: Accuracy problemSeverity: 3 - Normal
Operating System: Status: None
Assigned to: NoneOpen/Closed: Open
Release: 

Wed 11 Apr 2012 11:43:45 AM UTC, original submission:

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.

Jussi Lehtola <jussilehtola>

 

Attached Files
file #25629:  testbessel.c added by jussilehtola (320B - text/x-csrc)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by jussilehtola (Submitted the item)
  •  

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follows 1 latest change.

    Date Changed By Updated Field Previous Value => Replaced By
    Wed 11 Apr 2012 11:43:45 AM UTCjussilehtolaAttached File-=>Added testbessel.c, #25629
    Show feedback again

    Back to the top


    Powered by Savane 3.1-cleanup