Estimating the inverse trace using random forests on graphs
Computation
2023-04-24 v1
Abstract
Some data analysis problems require the computation of (regularised) inverse traces, i.e. quantities of the form . For large matrices, direct methods are unfeasible and one must resort to approximations, for example using a conjugate gradient solver combined with Girard's trace estimator (also known as Hutchinson's trace estimator). Here we describe an unbiased estimator of the regularized inverse trace, based on Wilson's algorithm, an algorithm that was initially designed to draw uniform spanning trees in graphs. Our method is fast, easy to implement, and scales to very large matrices. Its main drawback is that it is limited to diagonally dominant matrices .
Cite
@article{arxiv.1905.02086,
title = {Estimating the inverse trace using random forests on graphs},
author = {Simon Barthelmé and Nicolas Tremblay and Alexandre Gaudillière and Luca Avena and Pierre-Olivier Amblard},
journal= {arXiv preprint arXiv:1905.02086},
year = {2023}
}
Comments
Submitted to GRETSI conference