[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] gsl_sf_bessel_jl bug

From: Brian Gough
Subject: Re: [Bug-gsl] gsl_sf_bessel_jl bug
Date: Wed, 17 Oct 2007 14:52:00 +0100
User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI)

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.


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 -

reply via email to

[Prev in Thread] Current Thread [Next in Thread]