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: [email protected] Jakub Benda
To: [email protected]
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/
