Abstract:Krylov subspace methods are a ubiquitous tool for computing near-optimal rank $k$ approximations of large matrices. While "large block" Krylov methods with block size at least $k$ give the best known theoretical guarantees, block size one (a single vector) or a small constant is often preferred in practice. Despite their popularity, we lack theoretical bounds on the performance of such "small block" Krylov methods for low-rank approximation. We address this gap between theory and practice by proving that small block Krylov methods essentially match all known low-rank approximation guarantees for large block methods. Via a black-box reduction we show, for example, that the standard single vector Krylov method run for $t$ iterations obtains the same spectral norm and Frobenius norm error bounds as a Krylov method with block size $\ell \geq k$ run for $O(t/\ell)$ iterations, up to a logarithmic dependence on the smallest gap between sequential singular values. That is, for a given number of matrix-vector products, single vector methods are essentially as effective as any choice of large block size. By combining our result with tail-bounds on eigenvalue gaps in random matrices, we prove that the dependence on the smallest singular value gap can be eliminated if the input matrix is perturbed by a small random matrix. Further, we show that single vector methods match the more complex algorithm of [Bakshi et al. `22], which combines the results of multiple block sizes to achieve an improved algorithm for Schatten $p$-norm low-rank approximation.

Efficient Vectors for Block Perturbed Consistent Matrices

Efficient vectors for block perturbed consistent matrices

Efficiency analysis for the Perron vector of a reciprocal matrix

Efficiency of the convex hull of the columns of certain triple perturbed consistent matrices

The complete set of efficient vectors for a reciprocal matrix

Efficient vectors in priority setting methodology

Analysis of natural cardinal ranking vectors for pairwise comparisons and the universal efficiency of the Perron geometric mean

Cycle products and efficient vectors in reciprocal matrices

Pseudoinverse-free randomized block iterative algorithms for consistent and inconsistent linear systems

Efficient Randomized Algorithms for the Fixed-Precision Low-Rank Matrix Approximation

A Variation on the Block Arnoldi Method for Large Unsymmetric Matrix Eigenproblems

Efficient Block Preconditioned Eigensolvers for Linear Response Time-Dependent Density Functional Theory.

New results related to cutters and to an extrapolated block-iterative method for finding a common fixed point of a collection of them

On maximum residual block Kaczmarz method for solving large consistent linear systems

Adaptive Factored Incomplete Inverse Matrices

Perturbation and Inverse Problems of Stochastic Matrices

Angle-free cluster robust Ritz value bounds for restarted block eigensolvers

Block [math]-Circulant Preconditioners for Parabolic Optimal Control Problems

On the Unreasonable Effectiveness of Single Vector Krylov Methods for Low-Rank Approximation

Accurate Computations with Block Checkerboard Pattern Matrices

Robust Principal Component Analysis: Exact Recovery of Corrupted Low-Rank Matrices Via Convex Optimization.