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.

Improved analyses of the randomized power method and block Lanczos method

Improved Analyses of the Randomized Power Method and Block Lanczos Method

On the Error Estimate of the Randomized Double Block Kaczmarz Method

Efficient Bounds and Estimates for Canonical Angles in Randomized Subspace Approximations

Analysis of eigenvalue condition numbers for a class of randomized numerical methods for singular matrix pencils

Randomized matrix-free quadrature: unified and uniform bounds for stochastic Lanczos quadrature and the kernel polynomial method

An Accelerated Block Randomized Kaczmarz Method

Paved with good intentions: Analysis of a randomized block Kaczmarz method

A posteriori error bounds for the block-Lanczos method for matrix function approximation

Fast and Accurate Randomized Algorithms for Linear Systems and Eigenvalue Problems

On pseudoinverse-free randomized methods for linear systems: Unified framework and acceleration

Randomized block Kaczmarz method with projection for solving least squares

An analysis on stochastic Lanczos quadrature with asymmetric quadrature nodes

A linear algebra perspective on the random multi-block ADMM: the QP case

On the alternating randomized block Kaczmarz method

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

On the adaptive deterministic block Kaczmarz method with momentum for solving large-scale consistent linear systems

On global randomized block Kaczmarz algorithm for solving large-scale matrix equations