[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Bug in gsl_cdf_chisq_Pinv for very large degrees of freedom
From: |
John Halley Gotway |
Subject: |
[Bug-gsl] Bug in gsl_cdf_chisq_Pinv for very large degrees of freedom |
Date: |
Wed, 29 Oct 2008 11:14:48 -0600 |
User-agent: |
Thunderbird 2.0.0.17 (X11/20080914) |
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
- [Bug-gsl] Bug in gsl_cdf_chisq_Pinv for very large degrees of freedom,
John Halley Gotway <=