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/
