KIOPS: A fast adaptive Krylov subspace solver for exponential integrators

Abstract

This paper presents a new algorithm KIOPS for computing linear combinations of $\phi$-functions that appear in exponential integrators. This algorithm is suitable for large-scale problems in computational physics where little or no information about the spectrum or norm of the Jacobian matrix is known a priori. We first show that such problems can be solved efficiently by computing a single exponential of a modified matrix. Then our approach is to compute an appropriate basis for the Krylov subspace using the incomplete orthogonalization procedure and project the matrix exponential on this subspace. We also present a novel adaptive procedure that significantly reduces the computational complexity of exponential integrators. Our numerical experiments demonstrate that KIOPS outperforms the current state-of-the-art adaptive Krylov algorithm phipm.

References in corpus (5)

• 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
• A new approach to constructing efficient stiffly accurate exponential propagation iterative methods of Runge-Kutta type (EPIRK)

Cited by in corpus (14)

• Optimized auxiliary oscillators for the simulation of general open quantum systems
• KIOPS: A fast adaptive Krylov subspace solver for exponential integrators
• High-order numerical solutions to the shallow-water equations on the rotated cubed-sphere grid
• A $\mu$-mode integrator for solving evolution equations in Kronecker form
• Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes
• A Class of Exponential Integrators Based on Spectral Deferred Correction
• LeXInt: Package for Exponential Integrators employing Leja interpolation
• Efficient adaptive step size control for exponential integrators
• Avoiding matrix exponentials for large transition rate matrices
• Leveraging Mixed Precision in Exponential Time Integration Methods
• A scalable exponential-DG approach for nonlinear conservation laws: with application to Burger and Euler equations
• Exploring exponential time integration for strongly magnetized charged particle motion
• Error analysis of a class of semi-discrete schemes for solving the Gross-Pitaevskii equation at low regularity
• Simulating plasma wave propagation on a superconducting quantum chip