[Top][All Lists]

[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


                 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
         Discussion Lock: Any



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>

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;


-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:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]