bugGNU Scientific Library - Bugs: bug #42042, nan bug in bessel_Jnu

 
 

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

bug #42042: nan bug in bessel_Jnu

Submitted by:  Patrick Alken <psa>
Submitted on:  Thu 03 Apr 2014 09:36:21 PM UTC  
 
Category: Runtime errorSeverity: 3 - Normal
Operating System: Status: None
Assigned to: NoneOpen/Closed: Open
Release: 

Thu 03 Apr 2014 09:36:21 PM UTC, original submission:

Hello everyone,

I have this particular problem where I need

gsl_sf_bessel_Jnu(0.5, 3/2 * pi)

and I tested it with the following code:

#include <stdio.h>
#include <math.h>
#include <gsl/gsl_sf_bessel.h>

int
main (void)
{
double y1 = gsl_sf_bessel_Jnu(0.5, 3.0/2.0 * M_PI+1e-15);
double y2 = gsl_sf_bessel_Jnu(0.5, 3.0/2.0 * M_PI);

printf ("%e vs. %e\n", y1, y2);
return 0;
}

Output:

-3.675526e-01 vs. nan

MATLAB e.g. does return the correct value for 3/2 * pi. Is this a bug in
the implementation of the function? It seems so and I have not found
documentation about irregularities at 3/2 PI.

Patrick Alken <psa>
Project Administrator

 

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

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

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

     

     

    No Changes Have Been Made to This Item

    Back to the top


    Powered by Savane 3.1-cleanup