We have a function-like system that produces a output $z$ value for each 2-dimensional $(x,y)$ input. Suppose we have sampled the system using a set of non-uniform $x,y$ inputs (inputs are not in rectangular grids), and then the generated $(x,y) => z$ data is connected into a triangle mesh through some simple approach. How to perform smoothed interpolation on this stuff (given any free $(x,y)$, produce a $z$), like the Lagrange interpolation performed on 1-dimensional $x => y$ data?
And more, how to expand it to a system that has 3D input like $(x,y,z) => w$?