Sublinear randomized algorithms for skeleton decompositions
Abstract
Let $A$ be a $n$ by $n$ matrix. A skeleton decomposition is any factorization of the form $CUR$ where $C$ comprises columns of $A$, and $R$ comprises rows of $A$. In this paper, we consider uniformly sampling $ł\simeq k \log n$ rows and columns to produce a skeleton decomposition. The algorithm runs in $O(ł^3)$ time, and has the following error guarantee. Let $\norm{\cdot}$ denote the 2-norm. Suppose $A\simeq X B Y^T$ where $X,Y$ each have $k$ orthonormal columns. Assuming that $X,Y$ are incoherent, we show that with high probability, the approximation error $\norm{A-CUR}$ will scale with $(n/ł)\norm{A-X B Y^T}$ or better. A key step in this algorithm involves regularization. This step is crucial for a nonsymmetric $A$ as empirical results suggest. Finally, we use our proof framework to analyze two existing algorithms in an intuitive way.
- Publication:
-
arXiv e-prints
- Pub Date:
- October 2011
- DOI:
- 10.48550/arXiv.1110.4193
- arXiv:
- arXiv:1110.4193
- Bibcode:
- 2011arXiv1110.4193C
- Keywords:
-
- Mathematics - Numerical Analysis;
- Computer Science - Numerical Analysis