Fast algorithms for interpolation with clamped $L$-splines of order four
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , , |
|---|---|
| Format: | Preprint |
| Publié: |
2026
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866911698832588800 |
|---|---|
| author | Kounchev, O. Render, H. Simeonov, G. Tsachev, Ts. |
| author_facet | Kounchev, O. Render, H. Simeonov, G. Tsachev, Ts. |
| contents | Interpolation and smoothing using cubic and generalized splines are fundamental tools in data analysis and statistical modeling. Recently, fast computational algorithms were developed for natural $L$-splines of order four, which arise as piecewise solutions to the differential operator $L_ξ^2 = (\frac{d^2}{dt^2} - ξ^2)^2$. In this paper, we extend this mathematical framework to the important case of clamped (or complete) boundary conditions, where the first derivatives at the interval endpoints are prescribed. We explicitly construct the governing linear system for the interpolation problem and mathematically prove that the resulting tridiagonal matrix is strictly row diagonally dominant, thereby guaranteeing its invertibility and the numerical stability of the fast algorithm. The proposed method is implemented in MATLAB. Furthermore, the developed clamped $L$-splines provide a foundation for constructing multivariate clamped polysplines, which serve as a promising alternative to Physics-Informed Neural Networks (PINNs) for solving partial differential equations in Mathematical Physics. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2605_20283 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Fast algorithms for interpolation with clamped $L$-splines of order four Kounchev, O. Render, H. Simeonov, G. Tsachev, Ts. Numerical Analysis 65D07, 41A15, 65D05, 65F05 Interpolation and smoothing using cubic and generalized splines are fundamental tools in data analysis and statistical modeling. Recently, fast computational algorithms were developed for natural $L$-splines of order four, which arise as piecewise solutions to the differential operator $L_ξ^2 = (\frac{d^2}{dt^2} - ξ^2)^2$. In this paper, we extend this mathematical framework to the important case of clamped (or complete) boundary conditions, where the first derivatives at the interval endpoints are prescribed. We explicitly construct the governing linear system for the interpolation problem and mathematically prove that the resulting tridiagonal matrix is strictly row diagonally dominant, thereby guaranteeing its invertibility and the numerical stability of the fast algorithm. The proposed method is implemented in MATLAB. Furthermore, the developed clamped $L$-splines provide a foundation for constructing multivariate clamped polysplines, which serve as a promising alternative to Physics-Informed Neural Networks (PINNs) for solving partial differential equations in Mathematical Physics. |
| title | Fast algorithms for interpolation with clamped $L$-splines of order four |
| topic | Numerical Analysis 65D07, 41A15, 65D05, 65F05 |
| url | https://arxiv.org/abs/2605.20283 |