bug-gsl
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] addition information for reported bug


From: Yoram Burak
Subject: [Bug-gsl] addition information for reported bug
Date: Fri, 23 Dec 2005 16:35:59 -0800

gsl version: 1.7
computer system: PowerMac G4
operating system: Mac OS X 10.3.9
compiler: gcc 3.3

---------- Forwarded message ----------
From: Yoram Burak <address@hidden>
Date: Dec 23, 2005 4:23 PM
Subject: Error in evaluating scaled regular modified bessel functions
of large arguments
To: address@hidden


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


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

Voice: +1 805 893-7310
Fax: +1 805 893-2431
E-mail: address@hidden




reply via email to

[Prev in Thread] Current Thread [Next in Thread]