[Bug-gsl] [bug #37209] gsl_sf_bessel_jl_e returns NaN for large inputs without setting error code

[Rhys Ulerich]

URL:
  <http://savannah.gnu.org/bugs/?37209>

                 Summary: gsl_sf_bessel_jl_e returns NaN for large inputs
without setting error code
                 Project: GNU Scientific Library
            Submitted by: rhysu
            Submitted on: Tue 28 Aug 2012 02:15:22 PM GMT
                Category: Accuracy problem
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
                 Release: 
         Discussion Lock: Any

    _______________________________________________________

Details:

From: address@hidden Jakub Benda 
To:  address@hidden
Date: Tue, 28 Aug 2012 01:33:59 -0500 
Subject: [Bug-gsl] Spherical Bessel function failure 
 
The spherical Bessel function returns NaN for some (large) inputs without
setting the error code.
Example program:

--------------------------------------------------------------------------------
#include <stdio.h>
#include <gsl/gsl_sf.h>

int main(void)
{
gsl_sf_result res;
int err = gsl_sf_bessel_jl_e(364, 36.62, &res);

printf("err = %d\n", err);
printf("res.val = %g\n", res.val);
printf("res.err = %g\n", res.err);

return 0;
}
--------------------------------------------------------------------------------

Output:

--------------------------------------------------------------------------------
err = 0
res.val = -nan
res.err = nan
--------------------------------------------------------------------------------

Jakub




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?37209>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/