cubic spline interpolation

The cubic spline interpolation, see #807, is now implemented in separate cpp and h.