We describe how to perform the backward error analysis for the approximation of exp(A)v by p(s −1 A) s v, for any given polynomial p(x). The result of this analysis is an optimal choice of the scaling parameter s which assures a bound on the backward error, i.e. the equivalence of the approximation with the exponential of a slightly perturbed matrix. Thanks to the SageMath package expbea we have developed, one can optimize the performance of the given polynomial approximation. On the other hand, we employ the package for the analysis of polynomials interpolating the exponential function at so called Leja–Hermite points. The resulting method for the action of the matrix exponential can be considered an extension of both Taylor series approximation and Leja point interpolation. We illustrate the behavior of the new approximation with several numerical examples.

### Backward error analysis of polynomial approximations for computing the action of the matrix exponential

#### Abstract

We describe how to perform the backward error analysis for the approximation of exp(A)v by p(s −1 A) s v, for any given polynomial p(x). The result of this analysis is an optimal choice of the scaling parameter s which assures a bound on the backward error, i.e. the equivalence of the approximation with the exponential of a slightly perturbed matrix. Thanks to the SageMath package expbea we have developed, one can optimize the performance of the given polynomial approximation. On the other hand, we employ the package for the analysis of polynomials interpolating the exponential function at so called Leja–Hermite points. The resulting method for the action of the matrix exponential can be considered an extension of both Taylor series approximation and Leja point interpolation. We illustrate the behavior of the new approximation with several numerical examples.
##### Scheda breve Scheda completa Scheda completa (DC)
2018
Bacward error analysis, action of matrix exponential, Leja-Hermite interpolation, Taylor series
File in questo prodotto:
File
CKZ17.pdf

accesso aperto

Tipologia: Documento in Pre-print
Licenza: Dominio pubblico
Dimensione 307.97 kB
Utilizza questo identificativo per citare o creare un link a questo documento: `https://hdl.handle.net/11562/987259`