Shape-preserving uni-variate and bi-variate spline interpolation
