Compact Finite Differences and Cubic Splines

Robert M. Corless

DOI: https://doi.org/10.48550/arXiv.1805.07659

2019-11-22

Abstract:In this paper I uncover and explain---using contour integrals and residues---a connection between cubic splines and a popular compact finite difference formula. The connection is that on a uniform mesh the simplest Padé scheme for generating fourth-order accurate compact finite differences gives \textsl{exactly} the derivatives at the interior nodes needed to guarantee twice-continuous differentiability for cubic splines. %I found this connection surprising, because the two problems being solved are different. I also introduce an apparently new spline-like interpolant that I call a compact cubic interpolant; this is similar to one introduced in 1972 by Swartz and Varga, but has higher order accuracy at the edges. I argue that for mildly nonuniform meshes the compact cubic approach offers some potential advantages, and even for uniform meshes offers a simple way to treat the edge conditions, relieving the user of the burden of deciding to use one of the three standard options: free (natural), complete (clamped), or "not-a-knot" conditions. Finally, I establish that the matrices defining the compact cubic splines (equivalently, the fourth-order compact finite difference formulæ) are positive definite, and in fact totally nonnegative, if all mesh widths are the same sign.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to reveal and explain the connection between the fourth - order compact finite - difference formula and the cubic spline on a uniform grid. Specifically, the author uses the methods of complex integration and residues to show that the fourth - order compact finite - difference formula generated by the simplest Padé scheme can accurately give the derivatives of the interior nodes, thus ensuring the second - order continuous differentiability of the cubic spline.

In addition, the author also introduces a new spline interpolation method, called the compact cubic interpolant. This method has higher precision in handling boundary conditions and simplifies the burden on users when choosing boundary conditions. Finally, the author proves that the matrix defining the compact cubic spline (or the fourth - order compact finite - difference formula) is positive definite when all grid widths have the same sign, and is actually completely non - negative.

### Key problem summary:
1. **Reveal the relationship between compact finite - difference and cubic spline**: Through the methods of complex integration and residues, show the equivalence between the two on a uniform grid.
2. **Introduce compact cubic interpolation**: Propose a new interpolation method with higher boundary precision.
3. **Prove matrix properties**: Prove the positive definiteness and complete non - negativity of the relevant matrix.

### Formula representation:
- Compact finite - difference formula:
\[
\frac{1}{6} f'(t - h)+\frac{2}{3} f'(t)+\frac{1}{6} f'(t + h)=\frac{1}{2h}(f(t + h)-f(t - h))+O(h^{4})
\]
- Second - order derivative condition of cubic spline:
\[
p''_{k - 1}(\tau_{k}^{-})=p''_{k}(\tau_{k}^{+}),\quad k = 1,2,\ldots,n - 1
\]
- New method for handling boundary conditions:
\[
a_{0} f'(\tau_{0})+b_{0} f'(\tau_{1})=c_{0} f(\tau_{0})+c_{1} f(\tau_{1})+c_{2} f(\tau_{2})+c_{3} f(\tau_{3})
\]

Through these studies, the author not only reveals the internal connection between the two numerical methods, but also provides more efficient and accurate solutions for practical applications.