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: Rhys Ulerich <rhys dot ulerich at gmail dot com>
- To: Brian Gough <bjg at gnu dot org>
- Cc: gsl-discuss at sourceware dot org
- Date: Sat, 6 Feb 2010 09:25:47 -0600
- 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>
Sorry for the long hiatus on this Richardson extrapolation thread:
> I think the best interface would be the analog of the existing levin
> sum routines:
>
> ...
>
> int gsl_extrap_richardson (const double *array,
> ? ? ? ? ? ? ? ? ? ? ? ? ? const size_t n,
> ? ? ? ? ? ? ? ? ? ? ? ? ? const double t,
> ? ? ? ? ? ? ? ? ? ? ? ? ? const size_t k
> ? ? ? ? ? ? ? ? ? ? ? ? ? gsl_extrap_richardson_workspace * w,
> ? ? ? ? ? ? ? ? ? ? ? ? ? double *result, double *abserr);
>
> It would be useful to have both a result and abserr, to take into
> account cancellation error.
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.
Brian, if this is an easier thing to talk through on the phone and
you're willing to do so, please contact me directly.
- Rhys