toms886, a FORTRAN90 code 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