FlexTrace: Exchangeable Randomized Trace Estimation for Matrix Functions
Abstract
We consider the task of estimating the trace of a matrix function, , of a large symmetric positive semi-definite matrix . 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 , which can be very expensive. In this article, we introduce a novel trace estimator, FlexTrace, an exchangeable, single-pass method that estimates solely using matvecs with . We consider the case where is an operator monotone matrix function with , which includes functions such as and , 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 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