Fast accurate approximation of convolutions with weakly singular kernel and its applications

by   Awanish Kumar Tiwari, et al.

In this article, we present an O(N log N) rapidly convergent algorithm for the numerical approximation of the convolution integral with radially symmetric weakly singular kernels and compactly supported densities. To achieve the reduced computational complexity, we utilize the Fast Fourier Transform (FFT) on a uniform grid of size N for approximating the convolution. To facilitate this and maintain the accuracy, we primarily rely on a periodic Fourier extension of the density with a suitably large period depending on the support of the density. The rate of convergence of the method increases with increasing smoothness of the periodic extension and, in fact, approximations exhibit super-algebraic convergence when the extension is infinitely differentiable. Furthermore, when the density has jump discontinuities, we utilize a certain Fourier smoothing technique to accelerate the convergence to achieve the quadratic rate in the overall approximation. Finally, we apply the integration scheme for numerical solution of certain partial differential equations. Moreover, we apply the quadrature to obtain a fast and high-order Nystöm solver for the solution of the Lippmann-Schwinger integral equation. We validate the performance of the proposed scheme in terms of accuracy as well as computational efficiency through a variety of numerical experiments.


page 13

page 16

page 18

page 19

page 20


Fourier smoothed pre-corrected trapezoidal rule for solution of Lippmann-Schwinger integral equation

For the numerical solution of the Lippmann-Schwinger equation, while the...

Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification

This paper deals with the kernel-based approximation of a multivariate p...

High-order BDF convolution quadrature for subdiffusion models with a singular source term

Anomalous diffusion is often modelled in terms of the subdiffusion equat...

Composite Quadrature Methods for Weakly Singular Convolution Integrals

The well-known Caputo fractional derivative and the corresponding Caputo...

Numerical methods for stochastic Volterra integral equations with weakly singular kernels

In this paper, we first establish the existence, uniqueness and Hölder c...

Efficient high-order accurate Fresnel diffraction via areal quadrature and the nonuniform FFT

We present a fast algorithm for computing the diffracted field from arbi...

On optimal zero-padding of kernel truncation method

The kernel truncation method (KTM) is a commonly-used algorithm to compu...

Please sign up or login with your details

Forgot password? Click here to reset