Abstract:We consider the solution of large stiff systems of ordinary differential equations with explicit exponential Runge--Kutta integrators. These problems arise from semi-discretized semi-linear parabolic partial differential equations on continuous domains or on inherently discrete graph domains. A series of results reduces the requirement of computing linear combinations of $\varphi$-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 increase 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.

A low-rank algorithm for solving Lyapunov operator $φ$-functions within the matrix-valued exponential integrators

Low-rank plus diagonal approximations for Riccati-like matrix differential equations

Projected exponential methods for stiff dynamical low-rank approximation problems

A rank-adaptive robust integrator for dynamical low-rank approximation

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

An increasing rank Riemannian method for generalized Lyapunov equations

Low-Rank Methods for Solving Discrete-Time Projected Lyapunov Equations

Bernstein Operational Matrix Method for Solving the Numerical Solution of High Order Integro-Differential Equation with Weakly Singular

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

A Framework of the Harmonic Arnoldi Method for Evaluating Φ-Functions with Applications to Exponential Integrators

Randomized low-rank Runge-Kutta methods

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

An unconventional robust integrator for dynamical low-rank approximation

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Low rank approximation method for perturbed linear systems with applications to elliptic type stochastic PDEs

Full- and low-rank exponential Euler integrators for the Lindblad equation

A New Method for Computing $\varphi$-functions and Their Condition Numbers of Large Sparse Matrices

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

Positive Operator Based Iterative Algorithms for Solving Lyapunov Equations for Itô Stochastic Systems with Markovian Jumps

Low-synchronization Arnoldi Methods for the Matrix Exponential with Application to Exponential Integrators

A low-rank Lie-Trotter splitting approach for nonlinear fractional complex Ginzburg-Landau equations