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.

Krylov subspace methods for large multidimensional eigenvalue computation

Tensor Krylov subspace methods via the T-product for large Sylvester tensor equations

Deflated Block Krylov Subspace Methods for Large Scale Eigenvalue Problems.

A rational Krylov methods for large scale linear multidimensional dynamical systems

Spectral computation with third-order tensors using the t-product

A model reduction method for large-scale linear multidimensional dynamical systems

Krylov methods for low-rank commuting generalized Sylvester equations

Tensorized block rational Krylov methods for tensor Sylvester equations

A Refined Arnoldi Type Method for Large Scale Eigenvalue Problems

Aitken Extrapolation Method for Computing the Largest Eigenvalue of Nonnegative Tensors.

Subspace embedding with random Khatri-Rao products and its application to eigensolvers

Tensor Golub Kahan based on Einstein product

Randomized sketched TT-GMRES for linear systems with tensor structure

Estimating the numerical range with a Krylov subspace

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations

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

Computing the k-th Eigenvalue of Symmetric $H^2$-Matrices

Block subspace expansions for eigenvalues and eigenvectors approximation

Third-Order Tensors as Operators on Matrices: A Theoretical and Computational Framework with Applications in Imaging