Hello, gsl_sf_bessel_jl() crashes with at least one specific set of argument values.
This; 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. thanks, Koichi _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
