Monte Carlo Methods for Estimating the Diagonal of a Real Symmetric Matrix
Monte Carlo Methods for Estimating the Diagonal of a Real Symmetric Matrix
复制标题
估计实对称矩阵对角线的蒙特卡罗方法
DOI:
10.1137/22m1476277
复制
发表时间:
2023
影响因子:
1.5
通讯作者:
Saibaba, Arvind K.
中科院分区:
文献类型:
--
作者:
Hallman, Eric;Ipsen, Ilse C.;Saibaba, Arvind K.
For real symmetric matrices that are accessible only through matrix vector products, we present Monte Carlo estimators for computing the diagonal elements. Our probabilistic bounds for normwise absolute and relative errors apply to Monte Carlo estimators based on random Rademacher, sparse Rademacher, and normalized and unnormalized Gaussian vectors and to vectors with bounded fourth moments. The novel use of matrix concentration inequalities in our proofs represents a systematic model for future analyses. Our bounds mostly do not depend explicitly on the matrix dimension, target different error measures than existing work, and imply that the accuracy of the estimators increases with the diagonal dominance of the matrix. Applications to derivative-based global sensitivity metrics and node centrality measures in network science corroborate this, as do numerical experiments on synthetic test matrices. We recommend against the use in practice of sparse Rademacher vectors, which are the basis for many randomized sketching and sampling algorithms, because they tend to deliver barely a digit of accuracy even under large sampling amounts.
登录
查看更多内容
DOI:
10.21220/s2cc7x
发表时间:
2016
期刊:
ArXiv
影响因子:
--
作者:
Jesse Laeuchli
通讯作者:
Jesse Laeuchli
DOI:
--
发表时间:
2022
期刊:
arXiv.org
影响因子:
--
作者:
R. Baston;Y. Nakatsukasa
通讯作者:
Y. Nakatsukasa
DOI:
10.1093/imaiai/ias001
发表时间:
2011
期刊:
Information and Inference: A Journal of the IMA
影响因子:
--
作者:
Richard Y. Chen;Alex Gittens;J. Tropp
通讯作者:
J. Tropp
影响因子:
3.7
作者:
A. Soms
通讯作者:
A. Soms
DOI:
10.1007/978-3-319-97136-0_2
发表时间:
2017
期刊:
ArXiv
影响因子:
--
作者:
Shashanka Ubaru;Y. Saad
通讯作者:
Y. Saad