[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Error in evaluating scaled regular modified bessel functio
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Error in evaluating scaled regular modified bessel functions of large arguments |
Date: |
Wed, 28 Dec 2005 14:22:57 +0000 |
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