Authors

Abstract

We implement a linear-time algorithm for interpolation on a regular multidimensional gridin the Julia language. The algorithm is an approximate Laplace interpolation (Press, 1992)when no parameters are given; and when parameters m Z and > 0 are set, the interpolantapproximates a Matrn kernel, of which radial basis functions and polyharmonic splines are aspecial case. We implement, in addition, Neumann, Dirichlet (trivial), and average boundaryconditions with potentially different aspect ratios in the different dimensions. The interpolantfunctions in arbitrary dimensions.