[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] gsl_sf_bessel_jl bug
|
From: |
Koichi Takahashi |
|
Subject: |
[Bug-gsl] gsl_sf_bessel_jl bug |
|
Date: |
Tue, 16 Oct 2007 01:45:59 -0700 |
|
User-agent: |
Thunderbird 2.0.0.6 (X11/20071008) |
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
- [Bug-gsl] gsl_sf_bessel_jl bug,
Koichi Takahashi <=