This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: problem in gsl_sf_bessel_Jn
- To: Arvid Requate <arvid at Mathematik dot Uni-Bielefeld dot DE>
- Subject: Re: problem in gsl_sf_bessel_Jn
- From: Brian Gough <bjg at network-theory dot co dot uk>
- Date: Wed, 11 Jul 2001 09:58:03 +0100 (BST)
- Cc: gsl-discuss at sources dot redhat dot com
- References: <20010709212817.A21572@Mathematik.Uni-Bielefeld.DE>
Arvid Requate writes:
> there seem to be problems within gsl_sf_bessel_Jn and gsl_sf_bessel_Jn_array:
> - for n<168 they do not return exactly the same results.
> - for n=168 only gsl_sf_bessel_Jn returns successfull
> - for n>168 both routines return with errors (e.g. in airy.c)
>
> some example code (c++) is attached.
Thanks for the bug report. We'll look into it