We obtain an expression for the error in the approximation of f (A)b and bT f (A)b with rational Krylov methods, where A is a symmetric matrix, b is a vector and the function f admits an integral representation. The error expression is obtained by linking the matrix function error with the error in the approximate solution of shifted linear systems using the same rational Krylov subspace, and it can be exploited to derive both a priori and a posteriori error bounds. The error bounds are a generalization of the ones given in Chen et al. for the Lanczos method for matrix functions. A technique that we employ in the rational Krylov context can also be applied to refine the bounds for the Lanczos case.
Error bounds for the approximation of matrix functions with rational Krylov methods
Simunec, Igor
2024
Abstract
We obtain an expression for the error in the approximation of f (A)b and bT f (A)b with rational Krylov methods, where A is a symmetric matrix, b is a vector and the function f admits an integral representation. The error expression is obtained by linking the matrix function error with the error in the approximate solution of shifted linear systems using the same rational Krylov subspace, and it can be exploited to derive both a priori and a posteriori error bounds. The error bounds are a generalization of the ones given in Chen et al. for the Lanczos method for matrix functions. A technique that we employ in the rational Krylov context can also be applied to refine the bounds for the Lanczos case.File | Dimensione | Formato | |
---|---|---|---|
Numerical Linear Algebra App - 2024 - Simunec - Error bounds for the approximation of matrix functions with rational Krylov.pdf
accesso aperto
Tipologia:
Published version
Licenza:
Creative Commons
Dimensione
2.48 MB
Formato
Adobe PDF
|
2.48 MB | Adobe PDF |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.