Journal
JOURNAL OF THE ACM
Volume 58, Issue 2, Pages -Publisher
ASSOC COMPUTING MACHINERY
DOI: 10.1145/1944345.1944349
Keywords
Algorithms; Trace estimation; implicit linear operators
Categories
Funding
- Israel Science Foundation [1045/09]
Ask authors/readers for more resources
We analyze the convergence of randomized trace estimators. Starting at 1989, several algorithms have been proposed for estimating the trace of a matrix by 1/M Sigma(M)(i=1) z(i)(T) A(zi), where the z(i) are random vectors; different estimators use different distributions for the z(i)s, all of which lead to E (1/M Sigma(M)(i=1)z(i)(T)Az(i)) = trace(A). These algorithms are useful in applications in which there is no explicit representation of Abut rather an efficient method compute z(T) Az given z. Existing results only analyze the variance of the different estimators. In contrast, we analyze the number of samples M required to guarantee that with probability at least 1-delta, the relative error in the estimate is at most epsilon. We argue that such bounds are much more useful in applications than the variance. We found that these bounds rank the estimators differently than the variance; this suggests that minimum-variance estimators may not be the best. We also make two additional contributions to this area. The first is a specialized bound for projection matrices, whose trace (rank) needs to be computed in electronic structure calculations. The second is a new estimator that uses less randomness than all the existing estimators.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available