Interpolation, Method -> "Delaunay"

Interpolation[funcvaluelist,Method→"Delaunay",(options)] represents a way to perform interpolation utilizing Delaunay triangulation when funcvaluelist is a list of 2D scattered control points together with the function values in these points.
  • To use Interpolation, Method→"Delaunay", you first need to load the Obtuse Package using Needs["Obtuse`"].
In[1]:=
Click for copyable input
Some 2D points:
In[2]:=
Click for copyable input
Add function values in the 2D point table:
In[3]:=
Click for copyable input
Test interpolation in one point:
In[4]:=
Click for copyable input
Out[4]=
Check with the original function in the same point:
In[5]:=
Click for copyable input
Out[5]=
If the same interpolation should be used in many different points, the function can be precompiled without specifying in which point it will be used. The interpolated function is continuous. According to the ordinary definition of the Delaunay interpolation, the method is not specified outside the convex hull (drawn in white) of the interpolation points. To avoid error messages during the generation of the plots, like that below, I have defined extrapolated values. By drawing the outer bisectors from the convex hull points, and interpolating linearly along lines parallel to the edges of the convex hull, continuous extrapolated values can be found.
In[6]:=
Click for copyable input
Out[6]=