[Top][All Lists]

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

Re: [Bug-gsl] incorrect values by bessel_steed_array

From: Hu Zhan
Subject: Re: [Bug-gsl] incorrect values by bessel_steed_array
Date: Wed, 26 Jan 2005 09:23:38 -0800

I'm using gcc & gsl 1.6. The problem existed in gsl 1.5 as well. The
files (jl[0-5].txt) are in the attachment. The error gets progressively
smaller from l = 0 to 5 for x < 0.7.


On Wed, 2005-01-26 at 06:02, Brian Gough wrote:
> Hu Zhan writes:
>  > When I use gsl_sf_bessel_jl_steed_array(LMAX, x, wrk) with a large LMAX,
>  > e.g. LMAX = 80, the function returns wrong values at small l and x. My
>  > system is Debian Linux with Opteron CPU.
>  > 
>  > Attached is a program demonstrating this problem. It writes five files
>  > for l = 0 to 5 in the format of 
>  > 
>  > x  value_by_gsl_sf_bessel_jl_steed_array  value_by_gsl_sf_bessel_jl.
>  > 
>  > One can easily see the difference between the two functions. 
> Thanks for the email.
> Can you send the output of the program on your system for reference,
> so we can be sure we're dealing with the same problem. Thanks

Attachment: jl0.txt
Description: Text document

Attachment: jl1.txt
Description: Text document

Attachment: jl2.txt
Description: Text document

Attachment: jl3.txt
Description: Text document

Attachment: jl4.txt
Description: Text document

Attachment: jl5.txt
Description: Text document

reply via email to

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