Move assert() from 1D constructor to its second_derivative method.
This will assure that there are enough number of points for second order derivative. Also lets the interpolator to run the value() method with less points (order+1)
This will assure that there are enough number of points for second order derivative. Also lets the interpolator to run the value() method with less points (order+1)