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?
At Tue, 15 Dec 2009 11:31:10 -0600,
Rhys Ulerich wrote:
> You're right that each component could be extrapolated independently.
> I wrote it using vectors because doing so allows using BLAS calls for
> the linear algebra and speeds up many component use cases. I'd prefer
> to keep the code vector-capable under the covers under the theory that
> people extrapolating only a scalar at a time aren't all that worried
> about speed.
Ideally I'd like to start with a scalar version using normal C arrays,
similar to the gsl_sum functions, for implementation simplicity.
Also, are the "exact/normtable" arguments essential or just for
convenience? If they can be computed easily by the user it would be
good to have only a minimal number of arguments in the base function.
--
Brian Gough