Efficient adaptive step size control for exponential integrators
arXiv:2102.02524 · doi:10.1016/j.camwa.2022.07.011
Abstract
Traditional step size controllers make the tacit assumption that the cost of a time step is independent of the step size. This is reasonable for explicit and implicit integrators that use direct solvers. In the context of exponential integrators, however, an iterative approach, such as Krylov methods or polynomial interpolation, to compute the action of the required matrix functions is usually employed. In this case, the assumption of constant cost is not valid. This is, in particular, a problem for higher-order exponential integrators, which are able to take relatively large time steps based on accuracy considerations. In this paper, we consider an adaptive step size controller for exponential Rosenbrock methods that determines the step size based on the premise of minimizing computational cost. The largest allowed step size, given by accuracy considerations, merely acts as a constraint. We test this approach on a range of nonlinear partial differential equations. Our results show significant improvements (up to a factor of 4 reduction in the computational cost) over the traditional step size controller for a wide range of tolerances.
arXiv comment23 pages, 14 figures, matches with the published version
References in corpus (10)
- A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators
- KIOPS: A fast adaptive Krylov subspace solver for exponential integrators
- The Leja method revisited: backward error analysis for the matrix exponential
- On the performance of exponential integrators for problems in magnetohydrodynamics
- Exponential methods for solving hyperbolic problems with application to kinetic equations
- An exponential integrator for the drift-kinetic model
- An adaptive step size controller for iterative implicit methods
- Exponential Integrators for Resistive Magnetohydrodynamics: Matrix-free Leja Interpolation and Efficient Adaptive Time Stepping
- Explicit Exponential Rosenbrock Methods and their Application in Visual Computing
- EPIRK-W and EPIRK-K time discretization methods