bug-gsl
[Top][All Lists]
Advanced

[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/




reply via email to

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