A method of bivariate interpolation and smooth surface fitting based on local procedures | Synapse