On a class of splines free of Gibbs phenomenon

When interpolating data with certain regularity, spline functions are useful. They are defined as piecewise polynomials that satisfy certain regularity conditions at the joints. In the literature about splines it is possible to find several references that study the apparition of Gibbs phenomenon close to jump discontinuities in the results obtained by spline interpolation. This work is devoted to the construction and analysis of a new nonlinear technique that allows to improve the accuracy of splines near jump discontinuities eliminating the Gibbs phenomenon. The adaption is easily attained through a nonlinear modification of the right hand side of the system of equations of the spline, that contains divided differences. The modification is based on the use of a new limiter specifically designed to attain adaption close to jumps in the function. The new limiter can be seen as a nonlinear weighted mean that has better adaption properties than the linear weighted mean. We will prove that the nonlinear modification introduced in the spline keeps the maximum theoretical accuracy in all the domain except at the intervals that contain a jump discontinuity, where Gibbs oscillations are eliminated. Diffusion is introduced, but this is fine if the discontinuity appears due to a discretization of a high gradient with not enough accuracy. The new technique is introduced for cubic splines, but the theory presented allows to generalize the results very easily to splines of any order. The experiments presented satisfy the theoretical aspects analyzed in the paper.

[1]  Riaz A. Usmani,et al.  Inversion of Jacobi's tridiagonal matrix , 1994 .

[2]  Sergio Amat,et al.  New WENO Smoothness Indicators Computationally Efficient in the Presence of Corner Discontinuities , 2017, J. Sci. Comput..

[3]  Sergio Amat,et al.  Data Compression with ENO Schemes: A Case Study☆☆☆ , 2001 .

[4]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[5]  A. Harten.  Multiresolution representation of data: a general framework , 1996 .

[6]  Ming-Jun Lai,et al.  Polygonal Spline Spaces and the Numerical Solution of the Poisson Equation , 2016, SIAM J. Numer. Anal..

[7]  Jacques Liandrat,et al.  On the stability of the PPH nonlinear multiresolution , 2005 .

[8]  Jacques Liandrat,et al.  On a family of nonlinear cell-average multiresolution schemes for image processing: An experimental study , 2015, Math. Comput. Simul..

[9]  Antonio Marquina,et al.  Power ENO methods: a fifth-order accurate weighted power ENO method , 2004 .

[10]  Sergio Amat,et al.  Adaptive interpolation of images using a new nonlinear cell-average scheme , 2012, Math. Comput. Simul..

[11]  Chi-Wang Shu,et al.  High Order ENO and WENO Schemes for Computational Fluid Dynamics , 1999 .

[12]  Sergio Amat,et al.  On multiresolution schemes using a stencil selection procedure: applications to ENO schemes , 2007, Numerical Algorithms.

[13]  Larry L. Schumaker,et al.  Spline functions - computational methods , 2015 .

[14]  Jacques Liandrat,et al.  On a Power WENO Scheme with Improved Accuracy Near Discontinuities , 2017, SIAM J. Sci. Comput..

[15]  Jacques Liandrat,et al.  On a compact non-extrapolating scheme for adaptive image interpolation , 2012, J. Frankl. Inst..

[16]  Sergio Amat,et al.  Improving the compression rate versus L1 error ratio in cell-average error control algorithms , 2013, Numerical Algorithms.

[17]  Jacques Liandrat,et al.  On a class of L1-stable nonlinear cell-average multiresolution schemes , 2010, J. Comput. Appl. Math..

[18]  Francesc Aràndiga,et al.  Nonlinear multiscale decompositions: The approach of A. Harten , 2000, Numerical Algorithms.

[19]  J. Foster,et al.  Gibbs-Wilbraham splines , 1995 .

[20]  Zhimin Zhang,et al.  Convergence and Gibbs' phenomenon in cubic spline interpolation of discontinuous functions , 1997 .

[21]  ShuChi-Wang,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes, II , 1989 .

[22]  P. Alam.  ‘G’ , 2021, Composites Engineering: An A–Z Guide.

[23]  Sergio Amat,et al.  On an algorithm to adapt spline approximations to the presence of discontinuities , 2018, Numerical Algorithms.

[24]  Helmuth Späth,et al.  Exponential spline interpolation , 1969, Computing.

[25]  Sergio Amat,et al.  Reciprocal polynomial extrapolation vs Richardson extrapolation for singular perturbed boundary problems , 2012, Numerical Algorithms.

[26]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[27]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[28]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[29]  Bernard Bialecki,et al.  Modified Nodal Cubic Spline Collocation For Poisson's Equation , 2007, SIAM J. Numer. Anal..

[30]  Larry L. Schumaker,et al.  Stable Approximation and Interpolation with C1 Quartic Bivariate Splines , 2001, SIAM J. Numer. Anal..

[31]  F. B. Richards,et al.  A Gibbs phenomenon for spline functions , 1991 .

[32]  Nira Dyn,et al.  Approximation of piecewise smooth functions and images by edge-adapted (ENO-EA) nonlinear multiresolution techniques , 2008 .

[33]  Jacques Liandrat,et al.  High order nonlinear interpolatory reconstruction operators and associated multiresolution schemes , 2013, J. Comput. Appl. Math..

[34]  Pep Mulet,et al.  Point-Value WENO Multiresolution Applications to Stable Image Compression , 2010, J. Sci. Comput..

[35]  R. DeVore,et al.  Nonlinear approximation , 1998, Acta Numerica.

[36]  D. Schweikert.  An Interpolation Curve Using a Spline in Tension , 1966 .

[37]  Chi-Wang Shu,et al.  On the Gibbs Phenomenon and Its Resolution , 1997, SIAM Rev..

[38]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[39]  A. K. B. Chand,et al.  Generalized Cubic Spline Fractal Interpolation Functions , 2006, SIAM J. Numer. Anal..