pith. sign in

arxiv: 2311.07035 · v2 · pith:4H7NUNTInew · submitted 2023-11-13 · 🧮 math.NA · cs.NA

ContHutch++: Stochastic trace estimation for implicit integral operators

classification 🧮 math.NA cs.NA
keywords conthutchepsilonproductstracealgorithmestimateestimatorhutch
0
0 comments X
read the original abstract

Hutchinson's estimator is a randomized algorithm that computes an $\epsilon$-approximation to the trace of any positive semidefinite matrix using $\mathcal{O}(1/\epsilon^2)$ matrix-vector products. An improvement of Hutchinson's estimator, known as Hutch++, only requires $\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.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.