Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 1, Page 744-770, March 2024. We consider the solution of large stiff systems of ODEs with explicit exponential Runge–Kutta integrators. These problems arise from semidiscretized semilinear parabolic PDEs on continuous domains or on inherently discrete graph domains. A series of results reduces the requirement of computing linear combinations of [math]-functions in exponential integrators to the approximation of the action of a smaller number of matrix exponentials on certain vectors. State-of-the-art computational methods use polynomial Krylov subspaces of adaptive size for this task. They have the drawback that the required number of Krylov subspace iterations to obtain a desired tolerance increases drastically with the spectral radius of the discrete linear differential operator, e.g., the problem size. We present an approach that leverages rational Krylov subspace methods promising superior approximation qualities. We prove a novel a posteriori error estimate of rational Krylov approximations to the action of the matrix exponential on vectors for single time points, which allows for an adaptive approach similar to existing polynomial Krylov techniques. We discuss pole selection and the efficient solution of the arising sequences of shifted linear systems by direct and preconditioned iterative solvers. Numerical experiments show that our method outperforms the state of the art for sufficiently large spectral radii of the discrete linear differential operators. The key to this are approximately constant numbers of rational Krylov iterations, which enable a near-linear scaling of the runtime with respect to the problem size.

Matrix Pencil Optimal Iterative Algorithms and Restarted Versions for Linear Matrix Equation and Pseudoinverse

Iterative Methods by Space Decomposition and Subspace Correction

Low-Rank Matrix Iteration Using Polynomial-Filtered Subspace Extraction.

Sketched and truncated polynomial Krylov methods: Evaluation of matrix functions

An Indefinite Variant of LOBPCG for Definite Matrix Pencils

A Relaxed Algorithm for Online Matrix Inversion

Latest Inversion-Free Iterative Scheme for Solving a Pair of Nonlinear Matrix Equations

Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations

Generalized Pseudospectral Shattering and Inverse-Free Matrix Pencil Diagonalization

A positivity preserving inexact Noda iteration for computing the smallest eigenpair of a large irreducible M -matrix

Restarted Generalized Krylov Subspace Methods for Solving Large-Scale Polynomial Eigenvalue Problems

A Riemannian optimization method to compute the nearest singular pencil

Randomized Complete Pivoting for Solving Symmetric Indefinite Linear Systems

The shift-and-invert Arnoldi method for singular matrix pencils

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

Solutions to a Quadratic Inverse Eigenvalue Problem

A preconditioned Krylov subspace method for linear inverse problems with general-form Tikhonov regularization

Efficient scaling and squaring method for the matrix exponential

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems