This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Bessel function bug or limitation?
- From: Rodolphe Conan <rconan at eso dot org>
- Cc: gsl list <gsl-discuss at sources dot redhat dot com>
- Date: Mon, 08 Apr 2002 10:50:37 +0100
- Subject: Bessel function bug or limitation?
- References: <3CB16447.203@fkf.mpg.de> <3CB166CA.6000202@eso.org>
>
> Hi,
>
> I need to compute bessel function for large arguments (for example:
> gsl_sf_bessel_Jn(0,20000))
> and GSL crashes with the following error message:
> gsl: bessel.c:531: ERROR: error
> Abort (core dumped)
>
> I check the gsl_sf_bessel_Jn with 2 version of GSL 0.9.4 and 1.1.1:
> on GSL 0.9.4, GSL gave me the good value but with the wrong sign and
> now with the GSL 1.1.1 it crashes.
>
> Is it a bug or a known limitation of GSL?
>
>
> R. Conan
>
--
Rodolphe Conan
European Southern Observatory
Instrumentation Division
Adaptive Optics Group
Karl-Schwarzschild-Str. 2
D-85748 Garching bei Muenchen
Germany
Tel: +49-89-32006541 Fax: +49-89-3202362
e-mail : rconan@eso.org