Stable approximation and interpolation with C1 quartic bivariate splines

We show how two recent algorithms for computing C1 quartic interpolating splines can be stabilized to ensure that, for smooth functions, they provide full approximation power with approximation constants depending only on the smallest angle in the triangulation