Allow various interpolation methods for spline-coefficients of the kernel, and q0.
Right now the interpolation is linear for the kernel, and hard coded cubic splines are used for q0.
Right now the interpolation is linear for the kernel, and hard coded cubic splines are used for q0.