Numerical integration rules with improved accuracy close to singularities

by   Sergio Amat, et al.

Sometimes it is necessary to obtain a numerical integration using only discretised data. In some cases, the data contains singularities which position is known but does not coincide with a discretisation point, and the jumps in the function and its derivatives are available at these positions. The motivation of this paper is to use the previous information to obtain numerical quadrature formulas that allow approximating the integral of the discrete data over certain intervals accurately. This work is devoted to the construction and analysis of a new nonlinear technique that allows to obtain accurate numerical integrations of any order using data that contains singularities, and when the integrand is only known at grid points. The novelty of the technique consists in the inclusion of correction terms with a closed expression that depends on the size of the jumps of the function and its derivatives at the singularities, that are supposed to be known. The addition of these terms allows recovering the accuracy of classical numerical integration formulas even close to the singularities, as these correction terms account for the error that the classical integration formulas commit up to their accuracy at smooth zones. Thus, the correction terms can be added during the integration or as post-processing, which is useful if the main calculation of the integral has been already done using classical formulas. The numerical experiments performed allow us to confirm the theoretical conclusions reached in this paper.


page 1

page 2

page 3

page 4


Expansion of generalized Stieltjes constants in terms of derivatives of Hurwitz zeta-functions

Generalized Stieltjes constants γ n (a) are the coecients in the Laurent...

Low-regularity integrators for nonlinear Dirac equations

In this work, we consider the numerical integration of the nonlinear Dir...

Analysis of resolution of tomographic-type reconstruction from discrete data for a class of conormal distributions

Let f(x), x∈R^2, be a piecewise smooth function with a jump discontinuit...

Cell-average WENO with progressive order of accuracy close to discontinuities with applications to signal processing

In this paper we translate to the cell-average setting the algorithm for...

Optimal quadrature formulas for CT image reconstruction in the Sobolev space of non-periodic functions

In the present paper, optimal quadrature formulas in the sense of Sard a...

Theoretical error estimates for computing the matrix logarithm by Padé-type approximants

In this article, we focus on the error that is committed when computing ...

Fast parallel calculation of modified Bessel function of the second kind and its derivatives

There are three main types of numerical computations for the Bessel func...

Please sign up or login with your details

Forgot password? Click here to reset