toms886_test, a FORTRAN90 code which calls toms886(), which implements an interpolation procedure based on "Padua points", defined in the square [-1,+1]^2, whose interpolating power is especially good. It is possible to map these points to the general rectangle, ellipse or triangle to do interpolation on these regions as well.
The computer code and data files made available on this web page are distributed under the MIT license
toms886, a FORTRAN90 library which implements an interpolation procedure based on "Padua points", defined in the square [-1,+1]^2, whose interpolating power is especially good. It is possible to map these points to the general rectangle, ellipse or triangle to do interpolation on these regions as well.
ELLIPSE applies the procedure to an ellipse.
RECTANGLE applies the procedure to a rectangle.
TRIANGLE applies the procedure to a triangle.