Re: [Bug-gsl] incorrect values by bessel_steed_array

[Hu Zhan]

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.

Hu


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

jl0.txt
Description: Text document

jl1.txt
Description: Text document

jl2.txt
Description: Text document

jl3.txt
Description: Text document

jl4.txt
Description: Text document

jl5.txt
Description: Text document