Sketched inverse iteration applied to the sketching error gives a top-k eigensolver whose convergence rate is proportional to the quality of a Nyström preconditioner and depends only on the final spectral gap.
An implementation of a randomized algorithm for principal component analysis
1 Pith paper cite this work. Polarity classification is still indexing.
abstract
Recent years have witnessed intense development of randomized methods for low-rank approximation. These methods target principal component analysis (PCA) and the calculation of truncated singular value decompositions (SVD). The present paper presents an essentially black-box, fool-proof implementation for Mathworks' MATLAB, a popular software platform for numerical computation. As illustrated via several tests, the randomized algorithms for low-rank approximation outperform or at least match the classical techniques (such as Lanczos iterations) in basically all respects: accuracy, computational efficiency (both speed and memory usage), ease-of-use, parallelizability, and reliability. However, the classical procedures remain the methods of choice for estimating spectral norms, and are far superior for calculating the least singular values and corresponding singular vectors (or singular subspaces).
fields
math.NA 1years
2025 1verdicts
CONDITIONAL 1representative citing papers
citing papers explorer
-
What is a Sketch-and-Precondition Derivation for Low-Rank Approximation? Inverse Power Error or Inverse Power Estimation?
Sketched inverse iteration applied to the sketching error gives a top-k eigensolver whose convergence rate is proportional to the quality of a Nyström preconditioner and depends only on the final spectral gap.