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.
Hello 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 - http://www.gnu.org/software/gsl/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
