Skip to content
Preprint

A new analysis of the randomly pivoted Cholesky algorithm

Aug 2026 · 0 citations · 30 references
Mathematics Computer Science

Abstract

The randomly pivoted Cholesky algorithm is one of the leading methods for computing a low-rank approximation to a large positive-semidefinite matrix. However, while it consistently achieves accuracy comparable to or better than competing methods of its type in experiments, its theoretical analysis lags somewhat behind other methods. This paper closes this gap, proving that randomly pivoted Cholesky produces an approximation with expected error within a $1+\varepsilon$ factor of the optimal rank-$r$ approximation in $\mathcal{O}(r/\varepsilon + r\sqrt{\log r})$ steps. This result nearly matches the optimal complexity $\Theta(r/\varepsilon)$ for any low-rank approximation method based on a partial Cholesky decomposition (also known as a column Nystr\"om approximation). The paper also presents bounds on the randomly pivoted Cholesky trace and spectral-norm errors that hold with high probability. The mathematical argument is largely due to GPT 5.6-Sol (Pro), with some refinements by the author.

View source

Similar papers

Preprint Jul 2026

RCLUPPr: a new randomized CholeskyQR with LU preconditioning

In this work, we present the comprehensive rounding error analysis of RCLUPPr proposed in \cite{RCLUPP}, which is a novel randomized CholeskyQR-type algorithm performing LU decomposition with partial pivoting (LUPP decomposition) directly on the tall-skinny $X\in\mathbb{R}^{m\times n}$ with $m \ge n$ and $\mbox{rank}(X...

Haoran Guan, Zhenyu Zou, Yufeng Wei et al. · 0 citations
Preprint Aug 2026

The Resultant Distribution Method: Universality for $p$-adic Random Matrices and Polynomials

We prove universality of limiting local eigenvalue statistics for random matrices over $\mathbb{Z}_p$. In previous work of the author and Van Peski (arXiv:2601.06283), the limiting eigenvalue correlation functions of additive Haar random matrices were studied in arbitrary finite extensions of $\mathbb{Q}_p$. The same H...

Jiahe Shen · 0 citations
Preprint Aug 2026

Fast Length-Squared Sampling for Positive-Semidefinite Matrices

We describe a simple rejection-sampling-based algorithm to perform length-squared sampling on an $n \times n$ positive-semidefinite (psd) matrix: that is, to sample a column with probability proportional to its squared $\ell_2$-norm. The algorithm runs in just $O(n)$ expected time, which is significantly sublinear in t...

Rajarshi Bhattacharjee, Ethan N. Epperly, Cameron Musco et al. · 0 citations
Preprint Aug 2026

A Tight Analysis of Khatri-Rao Oblivious Subspace Embeddings

It is proved that sketching dimension m = O(k^{3/2}/\epsilon^2) suffices for subspace embedding with a Khatri-Rao sketching matrix with any fixed order $d$.

Lorenzo Beretta, Cameron Musco · 1 citation
Preprint Aug 2026

Beyond the Bethe Approximation of the Permanent

The canonical Bethe approximation gives a deterministic approximation to the permanent of every nonnegative matrix within a factor of $(\sqrt{2})^n$. We improve the base of this exponential factor: for some absolute constant $c<\sqrt{2}$, there is a deterministic polynomial-time $c^n$-approximation for the permanent of...

Nima Anari · 0 citations
Open access Aug 2026

An Analytic Expression of the Key Rank and a Saddle-Point Approximation

The partition function of the knapsack problem, which is equivalent to the key rank, is derived with an analytic expression and used to prove the mathematical equivalence between the knapsack and histogram-based methods.

M. Noes · 0 citations

We use cookies to run the site and, with your consent, for analytics and to show ads. See our Cookie Policy.