From: Rhys Ulerich <rhys.ulerich@gmail.com>
To: gsl-discuss@sourceware.org
Subject: Re: [PATCH] Adding fixed-order Gauss-Legendre integration routines
Date: Mon, 22 Jun 2009 14:00:00 -0000 [thread overview]
Message-ID: <4a00655d0906220700o4422e01s46edd50c8c4a1f95@mail.gmail.com> (raw)
In-Reply-To: <4a00655d0906211331y127e5693uda4805cf347561b3@mail.gmail.com>
> The routines include precomputed full double precision points and
> weights for n = 2, ..., 512, and 1024 Gauss points. If the user selects
> any other number, the necessary points and weights are computed
> on the fly and remain quite good.
A quick comment from Pavel about the numerical precision in the
patch's coefficients...
> ... fixed tables included in [the patch] are way beyond full double precision.
> Full double precision is about 1e-16, fixed tables are of 1e-25.
> This doesn't hurt the double calculations and allows future extensions
> e.g. for 'long double' and multi-precision libraries (like gmp or mpfr).
>
> Tables generated on the fly are of 1e-10, which is pretty much the
> maximum precision we can get using double to calculate weights.
> So, on-the-fly tables are of full double precision, fixed have much
> higher precision.
- Rhys
next prev parent reply other threads:[~2009-06-22 14:00 UTC|newest]
Thread overview: 6+ messages / expand[flat|nested] mbox.gz Atom feed top
[not found] <4a00655d0906211327x742b7224of110fa0621fbbad9@mail.gmail.com>
2009-06-21 20:31 ` Rhys Ulerich
2009-06-22 14:00 ` Rhys Ulerich [this message]
2009-06-23 16:52 ` Brian Gough
2009-06-23 16:58 ` Rhys Ulerich
2009-06-24 16:49 ` Rhys Ulerich
2009-06-28 20:41 ` Brian Gough
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=4a00655d0906220700o4422e01s46edd50c8c4a1f95@mail.gmail.com \
--to=rhys.ulerich@gmail.com \
--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).