ContHutch++: Stochastic trace estimation for implicit integral operators
Abstract
Hutchinson's estimator is a randomized algorithm that computes an -approximation to the trace of any positive semidefinite matrix using matrix-vector products. An improvement of Hutchinson's estimator, known as Hutch++, only requires 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}
}