The Leja method revisited: backward error analysis for the matrix exponential

Abstract

The Leja method is a polynomial interpolation procedure that can be used to compute matrix functions. In particular, computing the action of the matrix exponential on a given vector is a typical application. This quantity is required, e.g., in exponential integrators. The Leja method essentially depends on three parameters: the scaling parameter, the location of the interpolation points, and the degree of interpolation. We present here a backward error analysis that allows us to determine these three parameters as a function of the prescribed accuracy. Additional aspects that are required for an efficient and reliable implementation are discussed. Numerical examples that illustrate the performance of our Matlab code are included.

References in corpus (1)

• A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators

Cited by in corpus (20)

• KIOPS: A fast adaptive Krylov subspace solver for exponential integrators
• Further development of efficient and accurate time integration schemes for meteorological models
• Anisotropic osmosis filtering for shadow removal in images
• On the performance of exponential integrators for problems in magnetohydrodynamics
• Exponential methods for solving hyperbolic problems with application to kinetic equations
• A low-rank Lie-Trotter splitting approach for nonlinear fractional complex Ginzburg-Landau equations
• Magnus integrators on multicore CPUs and GPUs
• High-order numerical solutions to the shallow-water equations on the rotated cubed-sphere grid
• Limited-memory polynomial methods for large-scale matrix functions
• An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
• A $\mu$-mode integrator for solving evolution equations in Kronecker form
• An adaptive step size controller for iterative implicit methods
• Computing the Action of Trigonometric and Hyperbolic Matrix Functions
• Efficient adaptive step size control for exponential integrators
• Numerical study of the Serre-Green-Naghdi equations and a fully dispersive counterpart
• ParaExp using Leapfrog as Integrator for High-Frequency Electromagnetic Simulations
• Explicit Exponential Rosenbrock Methods and their Application in Visual Computing
• Exploring exponential time integration for strongly magnetized charged particle motion
• A study of defect-based error estimates for the Krylov approximation of $\phi$-functions
• On a nonlinear Schr{ö}dinger equation for nucleons in one space dimension