English

Analysis of probing techniques for sparse approximation and trace estimation of decaying matrix functions

Numerical Analysis 2021-02-09 v2 Numerical Analysis

Abstract

The computation of matrix functions f(A)f(A), or related quantities like their trace, is an important but challenging task, in particular for large and sparse matrices AA. In recent years, probing methods have become an often considered tool in this context, as they allow to replace the computation of f(A)f(A) or tr(f(A))\text{tr}(f(A)) by the evaluation of (a small number of) quantities of the form f(A)vf(A)v or vTf(A)vv^Tf(A)v, respectively. These tasks can then efficiently be solved by standard techniques like, e.g., Krylov subspace methods. It is well-known that probing methods are particularly efficient when f(A)f(A) is approximately sparse, e.g., when the entries of f(A)f(A) show a strong off-diagonal decay, but a rigorous error analysis is lacking so far. In this paper we develop new theoretical results on the existence of sparse approximations for f(A)f(A) and error bounds for probing methods based on graph colorings. As a by-product, by carefully inspecting the proofs of these error bounds, we also gain new insights into when to stop the Krylov iteration used for approximating f(A)vf(A)v or vTf(A)vv^Tf(A)v, thus allowing for a practically efficient implementation of the probing methods.

Keywords

Cite

@article{arxiv.2009.01589,
  title  = {Analysis of probing techniques for sparse approximation and trace estimation of decaying matrix functions},
  author = {Andreas Frommer and Claudia Schimmel and Marcel Schweitzer},
  journal= {arXiv preprint arXiv:2009.01589},
  year   = {2021}
}
R2 v1 2026-06-23T18:17:27.847Z