My second question is no doubt more complex. Suppose I have scattered
data points on some unstructured grid on the plane, and I have a
function value at each point. That is, the only data I have is a bunch
of (x,y,z) triplets that can be given to me in any order. They should
represent samplings from a smooth function, hopefully a solution to a
shallow water equation without shocks, once I get my code working. How
can I form a quasi-smooth surface which interpolates with small
polygons the data points that I give? I realise this may be a bit of a
complicated problem in computational geometry, but it also seems like
it should be a problem common enough for someone to have solved it
very well already.