Re: [Bug-gsl] Error in evaluating scaled regular modified bessel functions of large arguments

[Brian Gough]

Yoram Burak writes:
 > Attempt to evaluate gsl_sf_bessel_In_scaled(n,x) results in an error
 > for n>=2 and x larger than about 10^7.
 > 
 > #include <stdio.h>
 > #include <gsl/gsl_sf_bessel.h>
 > 
 > int main()
 > {
 >   double logx, x;
 > 
 >   for(logx = 0; logx < 20; logx++) {
 >     x = exp(logx);
 >     printf("%13le  %13le\n", x, gsl_sf_bessel_In_scaled(2,x));
 >   }
 > }
 > 

Hello, 

Thanks for the bug report -- can you send the output of the program on
your system for comparison as well.

-- 
Brian Gough