Problems with gsl_cdf_gamma_Pinv(0.01,a,b)

[Brian Gough]

At Mon, 22 Oct 2007 16:16:22 -0400,
Benjamin Redelings wrote:
> Q1: I'm very much interested in advice on how to handle this kind of 
> situation.  Should I expect this to work in GSL, or is there some reason 
> why GSL does not want this to work (e.g. speed issues)?  Should I write 
> some code myself to revert to a gaussian when a > 10?

Hello,

Thanks for the bug report, I have added a test case for that value.
You should expect this sort of thing to work in GSL -- definitely.

I added the hard failures in 1.10 as the problem was reported when I
was making the release and I didn't have time to fix it fully.  Some
of the other inverse functions were also affected so I thought it best
to report an error.

I would welcome any patches or other suggestions on the best way to
calculate this.  High precision test values to add to test.c are also
useful.

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/