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
- scientific article; zbMATH DE number 2212009 (Why is no real title available?)
- 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
- 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)- Accelerating the simulation of kinetic shear Alfvén waves with a dynamical low-rank approximation
- Fast and accurate computation of divided differences for analytic functions, with an application to the exponential function
- BAMPHI: matrix-free and transpose-free action of linear combinations of -functions from exponential integrators
- 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
- Approximation of the matrix exponential for matrices with a skinny field of values
- Further development of efficient and accurate time integration schemes for meteorological models
- A posteriori error estimates for exponential midpoint integrator finite element method for parabolic equations
- 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
- Accelerating Exponential Integrators to Efficiently Solve Semilinear Advection-Diffusion-Reaction Equations
- A study of defect-based error estimates for the Krylov approximation of \(\varphi\)-functions
- Exploring exponential time integration for strongly magnetized charged particle motion
- Convergence of a low-rank Lie-Trotter splitting for stiff matrix differential equations
- Should exponential integrators be used for advection-dominated problems?
- Numerical low-rank approximation of matrix differential equations
- Efficient third-order tensor-oriented directional splitting for exponential integrators
- A low-rank Lie-Trotter splitting approach for nonlinear fractional complex Ginzburg-Landau equations
- Robust linear stability analysis and a new method for computing the action of the matrix exponential
- Computing the action of trigonometric and hyperbolic matrix functions
- Interpolating discrete advection--diffusion propagators at Leja sequences
- On the algorithm by Al-Mohy and Higham for computing the action of the matrix exponential: a posteriori roundoff error estimation
- An adaptive step size controller for iterative implicit methods
- A scaling and recovering algorithm for the matrix -functions
- On the performance of exponential integrators for problems in magnetohydrodynamics
- Limited‐memory polynomial methods for large‐scale matrix functions
- A \(\mu\)-mode approach for exponential integrators: actions of \(\varphi\)-functions of Kronecker sums
- On the stability of exponential integrators for non-diffusive equations
- Efficient adaptive step size control for exponential integrators
- Exponential polynomial block methods
- Low-synchronization Arnoldi algorithms with application to exponential integrators
- On a nonlinear Schrödinger equation for nucleons in one space dimension
- Numerical study of the Serre-Green-Naghdi equations and a fully dispersive counterpart
- A nonlocal osmosis model for enhanced multi-image fusion
- A residual based error estimate for Leja interpolation of matrix functions
- GPU acceleration of splitting schemes applied to differential matrix equations
- Efficient exponential Runge-Kutta methods of high order: construction and implementation
- A low-rank algorithm for strongly damped wave equations with visco-elastic damping and mass terms
- Comparison of software for computing the action of the matrix exponential
- Error analysis of the truncated Taylor series expansion method for computing matrix exponential
- Efficient exponential time integration for simulating nonlinear coupled oscillators
- Comparing Leja and Krylov Approximations of Large Scale Matrix Exponentials
- Exponential Rosenbrock methods and their application in visual computing
- A Monte Carlo method for computing the action of a matrix exponential on a vector
- Exponential methods for solving hyperbolic problems with application to collisionless kinetic equations
- An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
- KIOPS: a fast adaptive Krylov subspace solver for exponential integrators
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)
