Analysis of stochastic Lanczos quadrature for spectrum approximation

Abstract

The cumulative empirical spectral measure (CESM) [A] : R [0,1] of a n× n symmetric matrix A is defined as the fraction of eigenvalues of A less than a given threshold, i.e., [A](x) := Σi=1n 1n x1D7D9[ λi[A]≤ x]. Spectral sums tr(f[A]) can be computed as the Riemann--Stieltjes integral of f against [A], so the task of estimating CESM arises frequently in a number of applications, including machine learning. We present an error analysis for stochastic Lanczos quadrature (SLQ). We show that SLQ obtains an approximation to the CESM within a Wasserstein distance of t \: | λmax[A] - λmin[A] | with probability at least 1-η, by applying the Lanczos algorithm for 12 t-1 + 12 iterations to 4 ( n+2 )-1t-2 (2nη-1) vectors sampled independently and uniformly from the unit sphere. We additionally provide (matrix-dependent) a posteriori error bounds for the Wasserstein and Kolmogorov--Smirnov distances between the output of this algorithm and the true CESM. The quality of our bounds is demonstrated using numerical experiments.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…