On an New Algorithm for Function Approximation with Full Accuracy in the Presence of Discontinuities Based on the Immersed Interface Method

This paper is devoted to the construction and analysis of an adapted and nonlinear multiresolution algorithm designed for interpolation or approximation of discontinuous univariate functions. The adaption attained allows to avoid numerical artifacts that appear when using linear algorithms and, at the same time, to obtain a high order of accuracy close to the singularities. It is known that linear algorithms are stable and convergent for smooth functions, but diffusion and Gibbs effect appear if the functions are piecewise continuous. Our aim is to develop an algorithm for function approximation with full accuracy that is capable to adapt to corners (kinks) and jump discontinuities, that uses a centered stencil and that does not use extrapolation. In order to reach this goal, we will need some information about the jumps in the function that we want to approximate and its derivatives. If this information is available, the algorithm is the most compact possible in the sense that the stencil is fixed and we do not need a stencil selection procedure as other algorithms do, such as ENO subcell resolution (ENO-SR). If the information about the jumps is not available, we will show a technique to approximate it. The algorithm is based on linear interpolation plus correction terms that provide the desired accuracy close to corners or jump discontinuities.

[1]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[2]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[3]  Jacques Liandrat,et al.  Analysis of a class of nonlinear subdivision schemes and associated multiresolution transforms , 2011, Adv. Comput. Math..

[4]  Rafael C. González,et al.  Local Determination of a Moving Contrast Edge , 1985, IEEE Transactions on Pattern Analysis and Machine Intelligence.

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

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

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

[8]  J. C. Trillo,et al.  On a Nonlinear Cell-Average Multiresolution Scheme for Image Compression , 2012 .

[9]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

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

[11]  Jacques Liandrat,et al.  Analysis of a New Nonlinear Subdivision Scheme. Applications in Image Processing , 2006, Found. Comput. Math..

[12]  Jacques Liandrat,et al.  On a nonlinear subdivision scheme avoiding Gibbs oscillations and converging towards Cs functions with s>1 , 2011, Math. Comput..

[13]  A. Cohen,et al.  Quasilinear subdivision schemes with applications to ENO interpolation , 2003 .

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

[15]  Zhilin Li,et al.  The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Frontiers in Applied Mathematics) , 2006 .

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

[17]  A. Harten.  ENO schemes with subcell resolution , 1989 .

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

[19]  Nira Dyn,et al.  Interpolation and Approximation of Piecewise Smooth Functions , 2005, SIAM J. Numer. Anal..

[20]  Pep Mulet,et al.  Adaptive interpolation of images , 2003, Signal Process..

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

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

[23]  Zhilin Li,et al.  The immersed interface method for the Navier-Stokes equations with singular forces , 2001 .

[24]  Albert Cohen,et al.  Tensor product multiresolution analysis with error control for compact image representation , 2002, Signal Process..

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

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

[27]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

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

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