Improved Bounds on Sample Size for Implicit Matrix Trace Estimators

This article is concerned with Monte Carlo methods for the estimation of the trace of an implicitly given matrix whose information is only available through matrix-vector products. Such a method approximates the trace by an average of expressions of the form , with random vectors drawn from an appropriate distribution. We prove, discuss and experiment with bounds on the number of realizations required to guarantee a probabilistic bound on the relative error of the trace estimation upon employing Rademacher (Hutchinson), Gaussian and uniform unit vector (with and without replacement) probability distributions. In total, one necessary bound and six sufficient bounds are proved, improving upon and extending similar estimates obtained in the seminal work of Avron and Toledo (JACM 58(2). Article 8, 2011) in several dimensions. We first improve their bound on for the Hutchinson method, dropping a term that relates to and making the bound comparable with that for the Gaussian estimator. We further prove new sufficient bounds for the Hutchinson, Gaussian and unit vector estimators, as well as a necessary bound for the Gaussian estimator, which depend more specifically on properties of matrix . As such, they may suggest the type of matrix for which one distribution or another provides a particularly effective or relatively ineffective stochastic estimation method.