Piecewise Hermite interpolation in one and two variables with applications to partial differential equations