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.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.