[bug-gsl] gsl_sf_bessel_jl bug (2007-10-16)
|
From: |
Koichi Takahashi <ktakahashi@..........> |
|
Subject: |
gsl_sf_bessel_jl bug |
|
Date: |
Tue, 16 Oct 2007 01:45:59 -0700 |
|
To: |
bug-gsl@gnu.org |
Hello,
gsl_sf_bessel_jl() crashes with at least one specific set of argument
values.
This;
double a = gsl_sf_bessel_jl( 30, 3875.6138424501978 );
fails in gsl_sf_bessel_J_CF1 () and invokes gsl_error().
Slightly different values, like
gsl_sf_bessel_jl( 30, 3875.610000000000 );
or
gsl_sf_bessel_jl( 29, 3875.6138424501978 );
seem to return correct results.
thanks,
Koichi
