English

Sparse graph based sketching for fast numerical linear algebra

Numerical Analysis 2021-02-12 v1 Numerical Analysis

Abstract

In recent years, a variety of randomized constructions of sketching matrices have been devised, that have been used in fast algorithms for numerical linear algebra problems, such as least squares regression, low-rank approximation, and the approximation of leverage scores. A key property of sketching matrices is that of subspace embedding. In this paper, we study sketching matrices that are obtained from bipartite graphs that are sparse, i.e., have left degree~s that is small. In particular, we explore two popular classes of sparse graphs, namely, expander graphs and magical graphs. For a given subspace URn\mathcal{U} \subseteq \mathbb{R}^n of dimension kk, we show that the magical graph with left degree s=2s=2 yields a (1±ϵ)(1\pm \epsilon) 2{\ell}_2-subspace embedding for U\mathcal{U}, if the number of right vertices (the sketch size) m=O(k2/ϵ2)m = \mathcal{O}({k^2}/{\epsilon^2}). The expander graph with s=O(logk/ϵ)s = \mathcal{O}({\log k}/{\epsilon}) yields a subspace embedding for m=O(klogk/ϵ2)m = \mathcal{O}({k \log k}/{\epsilon^2}). We also discuss the construction of sparse sketching matrices with reduced randomness using expanders based on error-correcting codes. Empirical results on various synthetic and real datasets show that these sparse graph sketching matrices work very well in practice.

Keywords

Cite

@article{arxiv.2102.05758,
  title  = {Sparse graph based sketching for fast numerical linear algebra},
  author = {Dong Hu and Shashanka Ubaru and Alex Gittens and Kenneth L. Clarkson and Lior Horesh and Vassilis Kalantzis},
  journal= {arXiv preprint arXiv:2102.05758},
  year   = {2021}
}
R2 v1 2026-06-23T23:03:14.952Z