Low Rank Approximation and Regression in Input Sparsity Time

TL;DR

This paper introduces a new sparse subspace embedding technique that significantly speeds up low-rank approximation and regression tasks, achieving near input sparsity time with improved accuracy and practical efficiency.

Contribution

The authors develop a novel sparse embedding matrix enabling faster algorithms for low-rank approximation and regression, reducing computational complexity to near input sparsity time.

Findings

Achieves subspace embedding in $ nz(A) + ext{poly}(d, rac{1}{eps})$ time.

Provides the fastest known algorithms for $(1+eps)$-approximate regression and low-rank approximation.

Experimental results indicate practical competitiveness of the proposed methods.

Abstract

We design a new distribution over $\poly(r \eps^{-1}) \times n$ matrices $S$ so that for any fixed $n \times d$ matrix $A$ of rank $r$, with probability at least 9/10, $\norm{SAx}_2 = (1 \pm \eps)\norm{Ax}_2$ simultaneously for all $x \in \mathbb{R}^d$. Such a matrix $S$ is called a \emph{subspace embedding}. Furthermore, $SA$ can be computed in $\nnz(A) + \poly(d \eps^{-1})$ time, where $\nnz(A)$ is the number of non-zero entries of $A$. This improves over all previous subspace embeddings, which required at least $\Omega(nd \log d)$ time to achieve this property. We call our matrices $S$ \emph{sparse embedding matrices}. Using our sparse embedding matrices, we obtain the fastest known algorithms for $(1+\eps)$-approximation for overconstrained least-squares regression, low-rank approximation, approximating all leverage scores, and $\ell_p$-regression. The leading order term in the time complexity of our algorithms is $O(\nnz(A))$ or $O(\nnz(A)\log n)$. We optimize the low-order $\poly(d/\eps)$ terms in our running times (or for rank-$k$ approximation, the $n*\poly(k/eps)$ term), and show various tradeoffs. For instance, we also use our methods to design new preconditioners that improve the dependence on $\eps$ in least squares regression to $\log 1/\eps$. Finally, we provide preliminary experimental results which suggest that our algorithms are competitive in practice.