Re: [Bug-gsl] gsl_sf_bessel_jl bug

[Brian Gough]

At Tue, 16 Oct 2007 01:45:59 -0700,
Koichi Takahashi wrote:
> gsl_sf_bessel_jl() crashes with at least one specific set of argument
> values.
> 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.

Hello

Thanks for the email.  Please could you send an example program which
shows the error and details of the version of GSL and your compiler
and operating system (the problem does not occur with gsl-1.10 on x86
with gcc).

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/