Philipp Hennig

Stochastic Lanczos Quadrature for Computational Uncertainty in Linear Algebra

It is shown that a single Lanczos run yields both a matrix-free approximation of the inverse and a Gaussian distribution over the exact result, whose covariance measures what the truncation discarded and concentrates as the number of steps grows.