Abstract:We present a new sublinear time algorithm for approximating the spectral density (eigenvalue distribution) of an $n\times n$ normalized graph adjacency or Laplacian matrix. The algorithm recovers the spectrum up to $\epsilon$ accuracy in the Wasserstein-1 distance in $O(n\cdot \text{poly}(1/\epsilon))$ time given sample access to the graph. This result compliments recent work by David Cohen-Steiner, Weihao Kong, Christian Sohler, and Gregory Valiant (2018), which obtains a solution with runtime independent of $n$, but exponential in $1/\epsilon$. We conjecture that the trade-off between dimension dependence and accuracy is inherent. Our method is simple and works well experimentally. It is based on a Chebyshev polynomial moment matching method that employees randomized estimators for the matrix trace. We prove that, for any Hermitian $A$, this moment matching method returns an $\epsilon$ approximation to the spectral density using just $O({1}/{\epsilon})$ matrix-vector products with $A$. By leveraging stability properties of the Chebyshev polynomial three-term recurrence, we then prove that the method is amenable to the use of coarse approximate matrix-vector products. Our sublinear time algorithm follows from combining this result with a novel sampling algorithm for approximating matrix-vector products with a normalized graph adjacency matrix. Of independent interest, we show a similar result for the widely used \emph{kernel polynomial method} (KPM), proving that this practical algorithm nearly matches the theoretical guarantees of our moment matching method. Our analysis uses tools from Jackson's seminal work on approximation with positive polynomial kernels.

Estimation of Laplacian spectra of direct and strong product graphs

A graph product and its Application

A Local Prime Factor Decomposition Algorithm for Strong Product Graphs

Product Graph Learning from Multi-domain Data with Sparsity and Rank Constraints

Spectral Methods for Matrix Product Factorization

Learning Cartesian Product Graphs with Laplacian Constraints

Efficient Estimation for Random Dot Product Graphs via a One-step Procedure

Analyzing the normalized Laplacian spectrum and spanning tree of the cross of the derivative of linear networks

Eigenvalues of graph Laplacians via rank-one perturbations

Estimates of bands for Laplacians on periodic equilateral metric graphs

Laplacian Spectra of a Class of Small-World Networks and Their Applications

Properties and applications of Laplacian spectra for Koch networks

Laplacian Mixture Modeling for Network Analysis and Unsupervised Learning on Graphs

Sublinear Time Spectral Density Estimation

Comparing Graph Spectra of Adjacency and Laplacian Matrices

A kernel-based analysis of Laplacian Eigenmaps

Network Analysis via Partial Spectral Factorization and Gauss Quadrature

Modeling Higher-Order Interactions in Complex Networks by Edge Product of Graphs

The normalized Laplacian spectrum of subdivisions of a graph

New bounds on the Laplacian spectral ratio of connected graphs

On sum of powers of Laplacian eigenvalues and Laplacian Estrada index of graphs