The Leja method revisited: backward error analysis for the matrix exponential
From MaRDI portal
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.
Recommendations
- A residual based error estimate for Leja interpolation of matrix functions
- Comparison of software for computing the action of the matrix exponential
- Comparing Leja and Krylov Approximations of Large Scale Matrix Exponentials
- Computing the action of the matrix exponential, with an application to exponential integrators
- The scaling and squaring method for the matrix exponential revisited
Cites work
- A residual based error estimate for Leja interpolation of matrix functions
- Accurate evaluation of divided differences for polynomial interpolation of exponential propagators
- Algorithm 919: A Krylov subspace algorithm for evaluating the -functions appearing in exponential integrators
- Comparison of software for computing the action of the matrix exponential
- Computing the action of the matrix exponential, with an application to exponential integrators
- Error Analysis of Direct Methods of Matrix Inversion
- Expokit
- Exponential integrators
- Functions of Matrices
- scientific article; zbMATH DE number 2212009 (Why is no real title available?)
- Interpolating discrete advection--diffusion propagators at Leja sequences
- Meshfree exponential integrators
- Newton interpolation at Leja points
- Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later
- The University of Florida sparse matrix collection
Cited in
(47)- On-the-fly backward error estimate for matrix exponential approximation by Taylor algorithm
- Backward error analysis of polynomial approximations for computing the action of the matrix exponential
- Numerical low-rank approximation of matrix differential equations
- An adaptive step size controller for iterative implicit methods
- On the performance of exponential integrators for problems in magnetohydrodynamics
- Interpolating discrete advection--diffusion propagators at Leja sequences
- KIOPS: a fast adaptive Krylov subspace solver for exponential integrators
- Efficient exponential Runge-Kutta methods of high order: construction and implementation
- Efficient adaptive step size control for exponential integrators
- BAMPHI: matrix-free and transpose-free action of linear combinations of -functions from exponential integrators
- On the stability of exponential integrators for non-diffusive equations
- An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
- Exponential methods for solving hyperbolic problems with application to collisionless kinetic equations
- A study of defect-based error estimates for the Krylov approximation of \(\varphi\)-functions
- A low-rank Lie-Trotter splitting approach for nonlinear fractional complex Ginzburg-Landau equations
- A \(\mu\)-mode integrator for solving evolution equations in Kronecker form
- High-order numerical solutions to the shallow-water equations on the rotated cubed-sphere grid
- Numerical study of the Serre-Green-Naghdi equations and a fully dispersive counterpart
- Approximation of the matrix exponential for matrices with a skinny field of values
- A residual based error estimate for Leja interpolation of matrix functions
- A Monte Carlo method for computing the action of a matrix exponential on a vector
- GPU acceleration of splitting schemes applied to differential matrix equations
- Further development of efficient and accurate time integration schemes for meteorological models
- On the algorithm by Al-Mohy and Higham for computing the action of the matrix exponential: a posteriori roundoff error estimation
- Comparison of software for computing the action of the matrix exponential
- Efficient exponential time integration for simulating nonlinear coupled oscillators
- Robust linear stability analysis and a new method for computing the action of the matrix exponential
- Exponential polynomial block methods
- On a nonlinear Schrödinger equation for nucleons in one space dimension
- Exponential Rosenbrock methods and their application in visual computing
- Convergence of a low-rank Lie-Trotter splitting for stiff matrix differential equations
- Comparing Leja and Krylov Approximations of Large Scale Matrix Exponentials
- Computing the action of trigonometric and hyperbolic matrix functions
- Error analysis of the truncated Taylor series expansion method for computing matrix exponential
- Limited‐memory polynomial methods for large‐scale matrix functions
- Accelerating the simulation of kinetic shear Alfvén waves with a dynamical low-rank approximation
- Accelerating Exponential Integrators to Efficiently Solve Semilinear Advection-Diffusion-Reaction Equations
- Fast and accurate computation of divided differences for analytic functions, with an application to the exponential function
- Exploring exponential time integration for strongly magnetized charged particle motion
- A \(\mu\)-mode approach for exponential integrators: actions of \(\varphi\)-functions of Kronecker sums
- Low-synchronization Arnoldi algorithms with application to exponential integrators
- A nonlocal osmosis model for enhanced multi-image fusion
- A low-rank algorithm for strongly damped wave equations with visco-elastic damping and mass terms
- A posteriori error estimates for exponential midpoint integrator finite element method for parabolic equations
- Efficient third-order tensor-oriented directional splitting for exponential integrators
- A scaling and recovering algorithm for the matrix -functions
- Should exponential integrators be used for advection-dominated problems?
Describes a project that uses
Uses Software
This page was built for publication: The Leja method revisited: backward error analysis for the matrix exponential
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2811992)
