The backward error analysis is a great tool which allows selecting in an effective way the scaling parameter s and the polynomial degree of approximation m when the action of the matrix exponential exp(A)v has to be approximated by (p(m)(s(-1)A)(s)v=exp(A+Delta A)v. We propose here a rigorous bound for the relative backward error Delta A(2)/A(2), which is of particular interest for matrices whose field of values is skinny, such as the discretization of the advection-diffusion or the Schrodinger operators. The numerical results confirm the superiority of the new approach with respect to methods based on the classical power series expansion of the backward error for the matrices of our interest, both in terms of computational cost and achieved accuracy.

Approximation of the matrix exponential for matrices with a skinny field of values

Caliari, M
;
Cassini, F;Zivcovich, F
2020-01-01

Abstract

The backward error analysis is a great tool which allows selecting in an effective way the scaling parameter s and the polynomial degree of approximation m when the action of the matrix exponential exp(A)v has to be approximated by (p(m)(s(-1)A)(s)v=exp(A+Delta A)v. We propose here a rigorous bound for the relative backward error Delta A(2)/A(2), which is of particular interest for matrices whose field of values is skinny, such as the discretization of the advection-diffusion or the Schrodinger operators. The numerical results confirm the superiority of the new approach with respect to methods based on the classical power series expansion of the backward error for the matrices of our interest, both in terms of computational cost and achieved accuracy.
2020
Backward error analysis
Action of matrix exponential
Leja-Hermite interpolation
Skinny field of values
File in questo prodotto:
File Dimensione Formato  
CCZ20.pdf

non disponibili

Tipologia: Versione dell'editore
Licenza: Accesso ristretto
Dimensione 467.04 kB
Formato Adobe PDF
467.04 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11562/1030250
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 5
social impact