English

Fast randomized numerical rank estimation for numerically low-rank matrices

Numerical Analysis 2024-01-08 v2 Numerical Analysis

Abstract

Matrices with low-rank structure are ubiquitous in scientific computing. Choosing an appropriate rank is a key step in many computational algorithms that exploit low-rank structure. However, estimating the rank has been done largely in an ad-hoc fashion in large-scale settings. In this work we develop a randomized algorithm for estimating the numerical rank of a (numerically low-rank) matrix. The algorithm is based on sketching the matrix with random matrices from both left and right; the key fact is that with high probability, the sketches preserve the orders of magnitude of the leading singular values. We prove a result on the accuracy of the sketched singular values and show that gaps in the spectrum are detected. For an m×nm\times n (mn)(m\geq n) matrix of numerical rank rr, the algorithm runs with complexity O(mnlogn+r3)O(mn\log n+r^3), or less for structured matrices. The steps in the algorithm are required as a part of many low-rank algorithms, so the additional work required to estimate the rank can be even smaller in practice. Numerical experiments illustrate the speed and robustness of our rank estimator.

Keywords

Cite

@article{arxiv.2105.07388,
  title  = {Fast randomized numerical rank estimation for numerically low-rank matrices},
  author = {Maike Meier and Yuji Nakatsukasa},
  journal= {arXiv preprint arXiv:2105.07388},
  year   = {2024}
}

Comments

To be published in Linear Algebra and its Applications

R2 v1 2026-06-24T02:09:07.359Z