I like this idea too, though perhaps the function name should be
gsl_interp_extrap_eval since its still using the gsl_interp objects
On 05/02/2013 08:54 PM, Rhys Ulerich wrote:
I also miss the good old days when extrapolation was automatically computed
using the chosen interpolation scheme. Can this behavior be added as a flag or
as a new function?
I could see a gsl_extrap_eval() and friends that were identical in
functionality but that did not perform this x < xmin || x > xmax
check. gsl_interp_eval() could then be implemented as this check
followed by calling gsl_extrap_eval(). People wanting the old
behavior could call the new function which would then self-document
that extrapolation is a possibility.
Thoughts?
- Rhys
