RBF_INTERP_ND
Multidimensional Interpolation with Radial Basis Functions

RBF_INTERP_ND is a FORTRAN90 library which defines and evaluates radial basis function (RBF) interpolants to multidimensional data.

A radial basis interpolant is a useful, but expensive, technique for definining a smooth function which interpolates a set of function values specified at an arbitrary set of data points.

Given nd multidimensional points xd with function values fd, and a basis function phi(r), the form of the interpolant is

       f(x) = sum ( 1 <= i <= nd ) w(i) * phi(||x-xd(i)||)
      

        sum ( 1 <= i <= nd ) w(i) * phi(||xd(j)-xd(i)||) = fd(j)
      

Four families of radial basis functions are provided.

• phi1(r) = sqrt ( r^2 + r0^2 ) (multiquadric)
• phi2(r) = 1 / sqrt ( r^2 + r0^2 ) (inverse multiquadric)
• phi3(r) = r^2 * log ( r / r0 ) (thin plate spline)
• phi4(r) = exp ( -0.5 r^2 / r0^2 ) (gaussian)

RBF_INTERP_ND needs access to the R8LIB library.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Languages:

RBF_INTERP_ND is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

LAGRANGE_INTERP_ND, a FORTRAN90 library which defines and evaluates the Lagrange polynomial p(x) which interpolates a set of data depending on a multidimensional argument x that was evaluated on a product grid, so that p(x(i)) = z(i).

R8LIB, a FORTRAN90 library which contains many utility routines using double precision real (R8) arithmetic.

RBF_INTERP_1D, a FORTRAN90 library which defines and evaluates radial basis function (RBF) interpolants to 1D data.

RBF_INTERP_2D, a FORTRAN90 library which defines and evaluates radial basis function (RBF) interpolants to 2D data.

SHEPARD_INTERP_ND, a FORTRAN90 library which defines and evaluates Shepard interpolants to multidimensional data, based on inverse distance weighting.

SPARSE_INTERP_ND a FORTRAN90 library which can be used to define a sparse interpolant to a function f(x) of a multidimensional argument.

SPINTERP, a MATLAB library which carries out piecewise multilinear hierarchical sparse grid interpolation; an earlier version of this software is ACM TOMS Algorithm 847, by Andreas Klimke;

TEST_INTERP_ND, a FORTRAN90 library which defines test problems for interpolation of data z(x), depending on an M-dimensional argument.

Reference:

1. Richard Franke,
   Scattered Data Interpolation: Tests of Some Methods,
   Mathematics of Computation,
   Volume 38, Number 157, January 1982, pages 181-200.
2. William Press, Brian Flannery, Saul Teukolsky, William Vetterling,
   Numerical Recipes in FORTRAN: The Art of Scientific Computing,
   Third Edition,
   Cambridge University Press, 2007,
   ISBN13: 978-0-521-88068-8,
   LC: QA297.N866.

Source Code:

• rbf_interp_nd.f90, the source code.

Examples and Tests:

• rbf_interp_nd_test.f90, a sample calling program.
• rbf_interp_nd_test.txt, the output file.

List of Routines:

• DAXPY computes constant times a vector plus a vector.
• DDOT forms the dot product of two vectors.
• DNRM2 returns the euclidean norm of a vector.
• DROT applies a plane rotation.
• DROTG constructs a Givens plane rotation.
• DSCAL scales a vector by a constant.
• DSVDC computes the singular value decomposition of a real rectangular matrix.
• DSWAP interchanges two vectors.
• PHI1 evaluates the multiquadric radial basis function.
• PHI2 evaluates the inverse multiquadric radial basis function.
• PHI3 evaluates the thin-plate spline radial basis function.
• PHI4 evaluates the gaussian radial basis function.
• R8MAT_SOLVE_SVD solves a linear system A*x=b using the SVD.
• RBF_INTERP_ND evaluates a radial basis function interpolant.
• RBF_WEIGHT computes weights for radial basis function interpolation.

You can go up one level to the FORTRAN90 source codes.