Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering and Quasi-Interpolation: A Unified View

Filtering plays a crucial role in postprocessing and analyzing data in scientific and engineering applications. Various application-specific filtering schemes have been proposed based on particular design criteria. In this paper, we focus on establishing the theoretical connection between quasi-interpolation and a class of kernels (based on B-splines) that are specifically designed for the postprocessing of the discontinuous Galerkin (DG) method called smoothness-increasing accuracy-conserving (SIAC) filtering. SIAC filtering, as the name suggests, aims to increase the smoothness of the DG approximation while conserving the inherent accuracy of the DG solution (superconvergence). Superconvergence properties of SIAC filtering has been studied in the literature. In this paper, we present the theoretical results that establish the connection between SIAC filtering to long-standing concepts in approximation theory such as quasi-interpolation and polynomial reproduction. This connection bridges the gap between the two related disciplines and provides a decisive advancement in designing new filters and mathematical analysis of their properties. In particular, we derive a closed formulation for convolution of SIAC kernels with polynomials. We also compare and contrast cardinal spline functions as an example of filters designed for image processing applications with SIAC filters of the same order, and study their properties.

[1]  Jennifer K. Ryan,et al.  Investigation of Smoothness-Increasing Accuracy-Conserving Filters for Improving Streamline Integration through Discontinuous Fields , 2008, IEEE Transactions on Visualization and Computer Graphics.

[2]  Jennifer K. Ryan,et al.  Efficient Implementation of Smoothness-Increasing Accuracy-Conserving (SIAC) Filters for Discontinuous Galerkin Solutions , 2012, J. Sci. Comput..

[3]  Michael Unser,et al.  B-spline signal processing. I. Theory , 1993, IEEE Trans. Signal Process..

[4]  Michael Unser,et al.  Fast B-spline Transforms for Continuous Image Representation and Interpolation , 1991, IEEE Trans. Pattern Anal. Mach. Intell..

[5]  Jennifer K. Ryan,et al.  Smoothness-Increasing Accuracy-Conserving Filters for Discontinuous Galerkin Solutions over Unstructured Triangular Meshes , 2013, SIAM J. Sci. Comput..

[6]  Jennifer K. Ryan,et al.  Accuracy-enhancement of discontinuous Galerkin solutions for convection-diffusion equations in multiple-dimensions , 2012, Math. Comput..

[7]  Robert Haimes,et al.  One-Sided Smoothness-Increasing Accuracy-Conserving Filtering for Enhanced Streamline Integration through Discontinuous Fields , 2008, J. Sci. Comput..

[8]  Jennifer K. Ryan,et al.  On a One-Sided Post-Processing Technique for the Discontinuous Galerkin Methods , 2003 .

[9]  Hartmut Prautzsch,et al.  Box Splines , 2002, Handbook of Computer Aided Geometric Design.

[10]  M. Marsden.  An identity for spline functions with applications to variation-diminishing spline approximation☆ , 1970 .

[11]  Bernardo Cockburn.  Discontinuous Galerkin methods , 2003 .

[12]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[13]  J. Daniel,et al.  Splines with Nonnegative B-Spline Coefficients , 1974 .

[14]  Gershon Elber,et al.  Geometric modeling with splines - an introduction , 2001 .

[15]  Jennifer K. Ryan,et al.  Position-Dependent Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering for Improving Discontinuous Galerkin Solutions , 2011, SIAM J. Sci. Comput..

[16]  Akram Aldroubi,et al.  B-SPLINE SIGNAL PROCESSING: PART I-THEORY , 1993 .

[17]  M. Unser.  Sampling-50 years after Shannon , 2000, Proceedings of the IEEE.

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

[19]  Michael Unser,et al.  Splines: a perfect fit for signal and image processing , 1999, IEEE Signal Process. Mag..

[20]  Robert M. Kirby,et al.  Smoothness-increasing accuracy-conserving filters (siac) for discontinuous galerkin solutions , 2013 .

[21]  Alireza Entezari,et al.  Quasi Interpolation With Voronoi Splines , 2011, IEEE Transactions on Visualization and Computer Graphics.

[22]  Endre Süli,et al.  Enhanced accuracy by post-processing for finite element methods for hyperbolic equations , 2003, Math. Comput..

[23]  Jennifer K. Ryan,et al.  Smoothness-Increasing Accuracy-Conserving (SIAC) Postprocessing for Discontinuous Galerkin Solutions over Structured Triangular Meshes , 2011, SIAM J. Numer. Anal..

[24]  Carl de Boor,et al.  Approximation Order without Quasi-Interpolants , 1993 .

[25]  Jennifer K. Ryan,et al.  Superconvergent error estimates for position-dependent smoothness-increasing accuracy-conserving (SIAC) post-processing of discontinuous Galerkin solutions , 2014, Math. Comput..

[26]  I. J. Schoenberg.  Cardinal Spline Interpolation , 1987 .

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

[28]  Jennifer K. Ryan,et al.  Quantification of Errors Introduced in the Numerical Approximation and Implementation of Smoothness-Increasing Accuracy Conserving (SIAC) Filtering of Discontinuous Galerkin (DG) Fields , 2010, J. Sci. Comput..

[29]  K. Höllig.  Finite element methods with B-splines , 1987 .

[30]  G. Strang,et al.  A Fourier Analysis of the Finite Element Variational Method , 2011 .

[31]  Jennifer K. Ryan,et al.  Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering for Discontinuous Galerkin Solutions: Improved Errors Versus Higher-Order Accuracy , 2012, J. Sci. Comput..

[32]  Vidar Thomée,et al.  High order local approximations to derivatives in the finite element method , 1977 .

[33]  Mahsa Mirzargar,et al.  A reconstruction framework for common sampling lattices , 2012 .

[34]  J. Bramble,et al.  Higher order local accuracy by averaging in the finite element method , 1977 .

[35]  Eduard Gröller,et al.  On Visualization and Reconstruction from Non‐Uniform Point Sets using B‐splines , 2009, Comput. Graph. Forum.

[36]  Jennifer K. Ryan,et al.  Smoothness-Increasing Accuracy-Conserving (SIAC) Filters for Discontinuous Galerkin Solutions: Application to Structured Tetrahedral Meshes , 2014, J. Sci. Comput..