[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [bug #42042] nan bug in bessel_Jnu
From: |
Patrick Alken |
Subject: |
[Bug-gsl] [bug #42042] nan bug in bessel_Jnu |
Date: |
Thu, 03 Apr 2014 21:36:21 +0000 |
User-agent: |
Mozilla/5.0 (X11; Linux i686 on x86_64; rv:28.0) Gecko/20100101 Firefox/28.0 |
URL:
<http://savannah.gnu.org/bugs/?42042>
Summary: nan bug in bessel_Jnu
Project: GNU Scientific Library
Submitted by: psa
Submitted on: Thu 03 Apr 2014 09:36:21 PM GMT
Category: Runtime error
Severity: 3 - Normal
Operating System:
Status: None
Assigned to: None
Open/Closed: Open
Release:
Discussion Lock: Any
_______________________________________________________
Details:
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.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?42042>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] [bug #42042] nan bug in bessel_Jnu,
Patrick Alken <=