[Bug-gsl] Bug in gsl_cdf_chisq_Pinv for very large degrees of freedom

[John Halley Gotway]

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