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.

A splitting-based KPIK method for eddy current optimal control problems in an all-at-once approach

A Gradient-based Kernel Optimization Approach for Parabolic Distributed Parameter Control Systems

Sequential Linear Quadratic Control of Bilinear Parabolic PDEs Based on POD Model Reduction

KIOPS: A fast adaptive Krylov subspace solver for exponential integrators

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Fast algorithm based on reduced precision solution criterion for solving optimal control problem

A new stopping criterion for Krylov solvers applied in Interior Point Methods

An Inversion-Free Iterative Algorithm with a Scalar Tuning Parameter for Coupled Riccati Matrix Equations Arising in LQ Optimal Control of Markov Jump Systems

An Efficient Parallel Krylov-Schur Method for Eigen-Analysis of Large-Scale Power Systems

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

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

A Hybrid Direct-Iterative Method for Solving KKT Linear Systems

Transform the Non-linear Programming Problem to the Initial-value Problem to Solve

A robust all-at-once multigrid method for the Stokes control problem

A new approach to constructing efficient stiffly accurate exponential propagation iterative methods of Runge-Kutta type (EPIRK)

NIPM-MPC: An Efficient Null-Space Method Based Interior-Point Method for Model Predictive Control

Rotated block diagonal preconditioners for Navier-Stokes control problems

APIK: Active Physics-Informed Kriging Model with Partial Differential Equations

Block [math]-Circulant Preconditioners for Parabolic Optimal Control Problems

Two Preconditioners for Time-Harmonic Eddy-Current Optimal Control Problems

Direct Iterative Algorithm to Solve the Extended Levine-Athans Equations