Fast and stable randomized low-rank matrix approximation
Abstract
Randomized SVD has become an extremely successful approach for efficiently computing a low-rank approximation of matrices. In particular the paper by Halko, Martinsson, and Tropp (SIREV 2011) contains extensive analysis, and has made it a very popular method. The typical complexity for a rank-r approximation of m× n matrices is O(mn n+(m+n)r2) for dense matrices. The classical Nystr\"om method is much faster, but applicable only to positive semidefinite matrices. This work studies a generalization of Nystr\"om method applicable to general matrices, and shows that (i) it has near-optimal approximation quality comparable to competing methods, (ii) the computational cost is the near-optimal O(mn n+r3) for dense matrices, with small hidden constants, and (iii) crucially, it can be implemented in a numerically stable fashion despite the presence of an ill-conditioned pseudoinverse. Numerical experiments illustrate that generalized Nystr\"om can significantly outperform state-of-the-art methods, especially when r 1, achieving up to a 10-fold speedup. The method is also well suited to updating and downdating the matrix.
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.