From: Brian Gough <bjg@gnu.org>
To: Benjamin Redelings <benjamin_redelings@ncsu.edu>
Cc: gsl-discuss@sourceware.org
Subject: Re: Problems with gsl_cdf_gamma_Pinv(0.01,a,b)
Date: Tue, 23 Oct 2007 22:11:00 -0000 [thread overview]
Message-ID: <m3640xh50k.wl%bjg@network-theory.co.uk> (raw)
In-Reply-To: <471D0516.5080103@ncsu.edu>
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/
prev parent reply other threads:[~2007-10-23 22:11 UTC|newest]
Thread overview: 2+ messages / expand[flat|nested] mbox.gz Atom feed top
2007-10-23 5:50 Benjamin Redelings
2007-10-23 22:11 ` Brian Gough [this message]
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=m3640xh50k.wl%bjg@network-theory.co.uk \
--to=bjg@gnu.org \
--cc=benjamin_redelings@ncsu.edu \
--cc=gsl-discuss@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).