English

FlexTrace: Exchangeable Randomized Trace Estimation for Matrix Functions

Numerical Analysis 2026-03-09 v1 Numerical Analysis

Abstract

We consider the task of estimating the trace of a matrix function, tr(f(A)){\rm tr}(f({\bf A})), of a large symmetric positive semi-definite matrix A{\bf A}. This problem arises in multiple applications, including kernel methods and inverse problems. A key challenge across existing trace estimation methods is the need for matrix-vector products (matvecs) with f(A)f({\bf A}), which can be very expensive. In this article, we introduce a novel trace estimator, FlexTrace, an exchangeable, single-pass method that estimates tr(f(A)){\rm tr}(f({\bf A})) solely using matvecs with A{\bf A}. We consider the case where ff is an operator monotone matrix function with f(0)=0f(0)=0, which includes functions such as log(1+x)\log(1+x) and x1/2x^{1/2}, and derive probabilistic bounds showcasing the theoretical advantages of FlexTrace. Numerical experiments across synthetic examples and application domains demonstrate that FlexTrace provides substantially more accurate estimates of the trace of f(A)f({\bf A}) compared to existing methods.

Keywords

Cite

@article{arxiv.2603.05721,
  title  = {FlexTrace: Exchangeable Randomized Trace Estimation for Matrix Functions},
  author = {Madhusudan Madhavan and Alen Alexanderian and Arvind K. Saibaba},
  journal= {arXiv preprint arXiv:2603.05721},
  year   = {2026}
}

Comments

35 pages, 10 figures

R2 v1 2026-07-01T11:05:49.950Z