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.
