This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: Where a generalized Richardson extrapolation routine would fit in GSL?
- From: Brian Gough <bjg at gnu dot org>
- To: Rhys Ulerich <rhys dot ulerich at gmail dot com>
- Cc: gsl-discuss at sourceware dot org
- Date: Mon, 08 Feb 2010 16:56:37 +0000
- Subject: Re: Where a generalized Richardson extrapolation routine would fit in GSL?
- References: <4a00655d0908201247g7d7bd9a1t466f4a66f08df4@mail.gmail.com><4a00655d0911291538y9f29830v984d1a796fdd5d1@mail.gmail.com><m34ooa7qvr.wl%bjg@network-theory.co.uk><4a00655d0912131549w19638273nb51d723e9ddd9273@mail.gmail.com><87d42ggnsv.wl%bjg@network-theory.co.uk><4a00655d0912150931q4b34fd24p94594ed08857254f@mail.gmail.com><87zl5ig7i0.wl%bjg@network-theory.co.uk><4a00655d0912161643k63868ac4i2577f618f6e2196@mail.gmail.com><87zl5hnrkh.wl%bjg@network-theory.co.uk><87ljh0nqp4.wl%bjg@network-theory.co.uk><4a00655d1002060725q21794722n12f69160c26e381a@mail.gmail.com>
At Sat, 6 Feb 2010 09:25:47 -0600,
Rhys Ulerich wrote:
> I'm good on moving the interface away from gsl_matrix and gsl_vector.
> No problem there. But, from looking a bit at the levin sum
> documentation, I'm not sure how to apply the abserr concept to
> Richardson extrapolation. Any guidance would be much appreciated.
The error term in the Levin sum is calculated by simple error
propagation assuming the individual terms are accurate to double
precision. Therefore it is just a measure of the cumulative
cancellation error. Since there are a lot of subtractions, it's worth
tracking this.
--
Brian Gough