English

Error bounds for the approximation of matrix functions with rational Krylov methods

Numerical Analysis 2023-11-07 v1 Numerical Analysis

Abstract

We obtain an expression for the error in the approximation of f(A)bf(A) \boldsymbol{b} and bTf(A)b\boldsymbol{b}^T f(A) \boldsymbol{b} with rational Krylov methods, where AA is a symmetric matrix, b\boldsymbol{b} is a vector and the function ff 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 [T. Chen, A. Greenbaum, C. Musco, C. Musco, SIAM J. Matrix Anal. Appl., 43 (2022), pp. 787--811] (arXiv:2106.09806) 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.

Keywords

Cite

@article{arxiv.2311.02701,
  title  = {Error bounds for the approximation of matrix functions with rational Krylov methods},
  author = {Igor Simunec},
  journal= {arXiv preprint arXiv:2311.02701},
  year   = {2023}
}

Comments

25 pages, 6 figures

R2 v1 2026-06-28T13:12:04.843Z