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.

On the use of Simpson's rule in renal models

Simpson's Rule Revisited

Estimating the error term in the Trapezium Rule using a Runge-Kutta method

A version of Simpson's rule for multiple integrals

A constructive proof of Simpson's Rule

A short proof of the error term in Simpson's rule

A Mathematical Model of Renal Blood Distribution Coupling TGF,MR and Tubular System

An Invitation to Hermite's Integration and Summation: A Comparison between Hermite's and Simpson's Rules

Efficient shifted fractional trapezoidal rule for subdiffusion problems with nonsmooth solutions on uniform meshes

One-dimensional Numerical Simulation of the Dialysis Process in an Artifical Kidney

Explicit and implicit error inhibiting schemes with post-processing

The exponential trapezoidal method for semilinear integro-differential equations

DEVELOPMENT AND TESTING OF THE MODIFIED SIMPSON’S RULE FOR DISCRETE ORDINATES TRANSPORT APPLICATIONS

Improving Newton's method performance by parametrization: the case of Richards equation

Modified Singly-Runge-Kutta-TASE methods for the numerical solution of stiff differential equations

Truncation Error Calculation Based on Richardson Extrapolation for Variable-Step Collaborative Simulation

Accuracy Analysis for Explicit-Implicit Finite Volume Schemes on Cut Cell Meshes

New error bounds for Boole's rule

On the role of the epithelium in a model of sodium exchange in renal tubules

Simple derivation of basic quadrature formulas

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators