Tricubic interpolation


In the mathematical subfield numerical analysis, tricubic interpolation is a method for obtaining values at arbitrary points in 3D space of a function defined on a regular grid. The approach involves approximating the function locally by an expression of the form
This form has 64 coefficients ; requiring the function to have a given value or given directional derivative at a point places one linear constraint on the 64 coefficients.
The term tricubic interpolation is used in more than one context; some experiments measure both the value of a function and its spatial derivatives, and it is desirable to interpolate preserving the values and the measured derivatives at the grid points. Those provide 32 constraints on the coefficients, and another 32 constraints can be provided by requiring smoothness of higher derivatives.
In other contexts, we can obtain the 64 coefficients by considering a 3×3×3 grid of small cubes surrounding the cube inside which we evaluate the function, and fitting the function at the 64 points on the corners of this grid.
The cubic interpolation article indicates that the method is equivalent to a sequential application of one-dimensional cubic interpolators. Let be the value of a monovariable cubic polynomial evaluated at. In many useful cases, these cubic polynomials have the form for some vector which is a function of alone. The tricubic interpolator is equivalent to:
where and.
At first glance, it might seem more convenient to use the 21 calls to described above instead of the matrix described in Lekien and Marsden. However, a proper implementation using a sparse format for the matrix makes the latter more efficient. This aspect is even much more pronounced when interpolation is needed at several locations inside the same cube. In this case, the matrix is used once to compute the interpolation coefficients for the entire cube. The coefficients are then stored and used for interpolation at any location inside the cube. In comparison, sequential use of one-dimensional integrators performs extremely poorly for repeated interpolations because each computational step must be repeated for each new location.