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
