FirstIndexNext

[bug-gsl] [bug #24704] error in gsl_cdf_chisq_Pinv for very large degrees of freedom (2008-10-30)

From: Brian Gough <INVALID.NOREPLY@gnu.org>
Subject: [bug #24704] error in gsl_cdf_chisq_Pinv for very large degrees of freedom
Date: Thu, 30 Oct 2008 19:11:33 +0000
To: johnhg@............, bug-gsl@gnu.org
Cc:


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
johnhg@........





_______________________________________________________

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/