Fast and oblivious convolution quadrature
Abstract
We give an algorithm to compute N steps of a convolution quadrature approximation to a continuous temporal convolution using only O(N N) multiplications and O( N) active memory. The method does not require evaluations of the convolution kernel, but instead O( N) evaluations of its Laplace transform, which is assumed sectorial. The algorithm can be used for the stable numerical solution with quasi-optimal complexity of linear and nonlinear integral and integro-differential equations of convolution type. In a numerical example we apply it to solve a subdiffusion equation with transparent boundary conditions.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.