Scattered data approximation with a hybrid scheme

A local hybrid radial–polynomial approximation scheme is introduced here to modify the scattered data fitting method presented in [6], generating C1 or C2 approximating spline surfaces. As for the original method, neither triangulation nor derivative estimate is needed, and the computational complexity is linear. The reported numerical experiments relate to two well known test functions and confirm both the accuracy and the shape recovery capability of the proposed hybrid scheme.