Re: Where a generalized Richardson extrapolation routine would fit in GSL?

[Brian Gough <bjg at gnu dot org>]

At Wed, 16 Dec 2009 18:43:44 -0600,
Rhys Ulerich wrote:
> > Ideally I'd like to start with a scalar version using normal C arrays,
> > similar to the gsl_sum functions, for implementation simplicity.
> 
> For API simplicity?  Or simplicity of the underlying implementation?

Where possible it is good to have an API that does not depend on
gsl_vector and gsl_matrix as it makes it simpler to use (or reuse) the
code from C programs with a minimal dependence on GSL.  Implementing
each routine to be minimal (without extra features) also makes life
simpler.

> There's just enough to the argument processing and general normtable
> handling that it's easier to get good test coverage if there's one
> implementation that does the whole enchilada and then everything else
> sits as convenience wrappers.  That prevents changes to one
> convenience function from accidentally changing the interface's
> semantics (especially for the 'k' handling).
> 
> All that said, I'm happy to make the interface conform.  Would you be
> willing to take draft public API in gsl_extrap.h (inline below) and
> change the functions and their signatures to resemble what you'd like?

I will send that once I've looked at it a bit more.