DeepAI AI Chat
Log In Sign Up

Kernel-Independent Sum-of-Exponentials with Application to Convolution Quadrature

by   Zixuan Gao, et al.

We propose an accurate algorithm for a novel sum-of-exponentials (SOE) approximation of kernel functions, and develop a fast algorithm for convolution quadrature based on the SOE, which allows an order N calculation for N time steps of approximating a continuous temporal convolution integral. The SOE method is constructed by a combination of the de la Vallée-Poussin sums for a semi-analytical exponential expansion of a general kernel, and a model reduction technique for the minimization of the number of exponentials under given error tolerance. We employ the SOE expansion for the finite part of the splitting convolution kernel such that the convolution integral can be solved as a system of ordinary differential equations due to the exponential kernels. The remaining part is explicitly approximated by employing the generalized Taylor expansion. The significant features of our algorithm are that the SOE method is efficient and accurate, and works for general kernels with controllable upperbound of positive exponents. We provide numerical analysis for the SOE-based convolution quadrature. Numerical results on the convolution integral and integral equations demonstrate attractive performance of both accuracy and efficiency of the proposed method.


page 1

page 2

page 3

page 4


A kernel-independent sum-of-Gaussians method by de la Vallée-Poussin sums

Approximation of interacting kernels by sum of Gaussians (SOG) is freque...

The Fast and Free Memory Method for the efficient computation of convolution kernels

We introduce the Fast Free Memory method (FFM), a new fast method for th...

Composite Quadrature Methods for Weakly Singular Convolution Integrals

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

Novel Convolution Kernels for Computer Vision and Shape Analysis based on Electromagnetism

Computer vision is a growing field with a lot of new applications in aut...

A fast sparse spectral method for nonlinear integro-differential Volterra equations with general kernels

We present a sparse spectral method for nonlinear integro-differential V...

A fast and oblivious matrix compression algorithm for Volterra integral operators

The numerical solution of dynamical systems with memory requires the eff...

A convolution quadrature method for Maxwell's equations in dispersive media

We study the systematic numerical approximation of Maxwell's equations i...