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

[Yoram Burak]

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));
  }
}





--
Yoram Burak
Kavli Institute for Theoretical Physics
UC Santa Barbara, California 93106