[Bug-gsl] [bug #31528] gsl_sf_bessel_lnKnu overflows for large nu

[eyal gruss]

URL:
  <http://savannah.gnu.org/bugs/?31528>

                 Summary: gsl_sf_bessel_lnKnu overflows for large nu
                 Project: GNU Scientific Library
            Submitted by: eyaler
            Submitted on: Tue 02 Nov 2010 05:13:13 PM GMT
                Category: Runtime error
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
                 Release: 1.14
         Discussion Lock: Any

    _______________________________________________________

Details:

gsl_sf_bessel_lnKnu(nu, x) computes the natural logarithm of the modified
bessel function of the second kind, with real arguments.

the function implementation overflows for large order (nu). for example:
gsl_sf_bessel_lnKnu(180.0, 2.2) should return ~ 175.2, but it gives INF
gsl_sf_bessel_lnKnu(3500.5, 1500) should return ~ 1731.2, but it gives INF




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?31528>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/