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.

Libkrylov: A modular open‐source software library for extremely large on‐the‐fly matrix computations

Preconditioned Krylov solvers for kernel regression

Pipelined, Flexible Krylov Subspace Methods

Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations

KIOPS: A fast adaptive Krylov subspace solver for exponential integrators

Krylov Solvers for Interior Point Methods with Applications in Radiation Therapy and Support Vector Machines

Efficient Block Preconditioned Eigensolvers for Linear Response Time-Dependent Density Functional Theory.

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Sketched and truncated polynomial Krylov methods: Evaluation of matrix functions

Deflated Block Krylov Subspace Methods for Large Scale Eigenvalue Problems.

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

Comparative study of inner-outer Krylov solvers for linear systems in structured and high-order unstructured CFD problems

LancBiO: dynamic Lanczos-aided bilevel optimization via Krylov subspace

On a vectorized basic linear algebra package for prototyping codes in MATLAB

Krylov subspace methods for large multidimensional eigenvalue computation

A Hybrid Direct-Iterative Method for Solving KKT Linear Systems

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

LAW: A Tool for Improved Productivity with High-Performance Linear Algebra Codes. Design and Applications

XAMG: A library for solving linear systems with multiple right-hand side vectors

Parallel-in-time Multilevel Krylov Methods: A Prototype

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix