[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [bug #36152] Incorrect asymptotics of spherical Bessel functio
From: |
Patrick Alken |
Subject: |
[Bug-gsl] [bug #36152] Incorrect asymptotics of spherical Bessel functions |
Date: |
Tue, 23 Jul 2013 15:57:43 +0000 |
User-agent: |
Mozilla/5.0 (X11; Linux i686 on x86_64; rv:22.0) Gecko/20100101 Firefox/22.0 |
Follow-up Comment #2, bug #36152 (project gsl):
The problem with gsl_sf_bessel_j0 is related to gsl_sf_sin (see the following
program). Program output is:
---
y_gsl = 9.155440021830e+41 y = -6.452512852658e-01
---
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_sf_trig.h>
int
main()
{
double y, y_gsl;
const double x = 1.0e20;
y_gsl = gsl_sf_sin(x);
y = sin(x);
fprintf(stderr, "y_gsl = %.12e y = %.12e\n", y_gsl, y);
return 0;
}
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?36152>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Bug-gsl] [bug #36152] Incorrect asymptotics of spherical Bessel functions,
Patrick Alken <=