This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: Problems with gsl_cdf_gamma_Pinv(0.01,a,b)
- From: Brian Gough <bjg at gnu dot org>
- To: Benjamin Redelings <benjamin_redelings at ncsu dot edu>
- Cc: gsl-discuss at sourceware dot org
- Date: Tue, 23 Oct 2007 22:53:31 +0100
- Subject: Re: Problems with gsl_cdf_gamma_Pinv(0.01,a,b)
- References: <471D0516.5080103@ncsu.edu>
- Reply-to: gsl-discuss at sourceware dot org
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/