skip to main content

Title: Hierarchical probing for estimating the trace of the matrix inverse on toroidal lattices

The standard approach for computing the trace of the inverse of a very large, sparse matrix $A$ is to view the trace as the mean value of matrix quadratures, and use the Monte Carlo algorithm to estimate it. This approach is heavily used in our motivating application of Lattice QCD. Often, the elements of $$A^{-1}$$ display certain decay properties away from the non zero structure of $A$, but random vectors cannot exploit this induced structure of $$A^{-1}$$. Probing is a technique that, given a sparsity pattern of $A$, discovers elements of $A$ through matrix-vector multiplications with specially designed vectors. In the case of $$A^{-1}$$, the pattern is obtained by distance-$k$ coloring of the graph of $A$. For sufficiently large $k$, the method produces accurate trace estimates but the cost of producing the colorings becomes prohibitively expensive. More importantly, it is difficult to search for an optimal $k$ value, since none of the work for prior choices of $k$ can be reused.
 [1] ;  [1] ;  [2]
  1. College of William and Mary, Williamsburg, VA
  2. College of William and Mary, Williamsburg, VA; Jefferson Lab
Publication Date:
OSTI Identifier:
Report Number(s):
JLAB-THY-13-1725; DOE/OR/23177-2756; arXiv:1302.4018
Journal ID: 1095-7197; NSF grant No. CCF-1218349
DOE Contract Number:
AC05-06OR23177; FC02-12ER4189
Resource Type:
Resource Relation:
Journal Name: SIAM J. Sci. Comput.; Journal Volume: 35; Journal Issue: 5; Conference: Twelfth Copper Mountain Conference on Iterative Methods, Copper Mountain, Colorado, USA, March 25–30, 2012
Research Org:
Thomas Jefferson National Accelerator Facility, Newport News, VA (United States)
Sponsoring Org:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21)
Country of Publication:
United States