English

ContHutch++: Stochastic trace estimation for implicit integral operators

Numerical Analysis 2024-09-26 v2 Numerical Analysis

Abstract

Hutchinson's estimator is a randomized algorithm that computes an ϵ\epsilon-approximation to the trace of any positive semidefinite matrix using O(1/ϵ2)\mathcal{O}(1/\epsilon^2) matrix-vector products. An improvement of Hutchinson's estimator, known as Hutch++, only requires O(1/ϵ)\mathcal{O}(1/\epsilon) matrix-vector products. In this paper, we propose a generalization of Hutch++, which we call ContHutch++, that uses operator-function products to efficiently estimate the trace of any trace-class integral operator. Our ContHutch++ estimates avoid spectral artifacts introduced by discretization and are accompanied by rigorous high-probability error bounds. We use ContHutch++ to derive a new high-order accurate algorithm for quantum density-of-states and also show how it can estimate electromagnetic fields induced by incoherent sources.

Cite

@article{arxiv.2311.07035,
  title  = {ContHutch++: Stochastic trace estimation for implicit integral operators},
  author = {Jennifer Zvonek and Andrew Horning and Alex Townsend},
  journal= {arXiv preprint arXiv:2311.07035},
  year   = {2024}
}
R2 v1 2026-06-28T13:18:50.230Z