[Bug-gsl] [bug #24704] error in gsl_cdf_chisq_Pinv for very large degrees of freedom

[Brian Gough]

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

                 Summary: error in gsl_cdf_chisq_Pinv for very large degrees
of freedom
                 Project: GNU Scientific Library
            Submitted by: bjg
            Submitted on: Thu 30 Oct 2008 07:11:29 PM GMT
                Category: Runtime error
                Severity: 3 - Normal
        Operating System: 
                  Status: Confirmed
             Assigned to: None
             Open/Closed: Open
                 Release: 1.11
         Discussion Lock: Any

    _______________________________________________________

Details:

Dear GSL,

I've encountered a potential bug in the routine "gsl_cdf_chisq_Pinv" defined
in the file gsl-1.11/cdf/chisqinv.c.

When calling this routine with the degrees of freedom > 1263131, I get the
following error message:

gsl: gamma_inc.c:99: ERROR: error
Default GSL error handler invoked.
Abort

For degrees of freedom <= 1263131, it seems to perform fine.

I'm running GSL version 1.11 and calling the routine as follows:
   v = gsl_cdf_chisq_Pinv(0.05, nu);
where nu = the degrees of freedom.

Thanks,
John Halley Gotway
address@hidden





    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: Thu 30 Oct 2008 07:11:29 PM GMT  Name: chisq3.c  Size: 508B   By: bjg

<http://savannah.gnu.org/bugs/download.php?file_id=16760>

    _______________________________________________________

Reply to this item at:

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

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