clopen.solutions

Sharp worst-case factors for randomly pivoted Cholesky and LU

Recorded by the registry's curators from GitHub; the authors did not submit it here.

Matthew J. Colbrook

For every positive integer $r$, an explicit one-parameter family of real,
entrywise-positive, positive-definite matrices of order $r+1$ makes the
expected trace error of $r$ randomly pivoted Cholesky steps approach $2^r$
times the optimal rank-$r$ trace error. The same family makes the expected
squared Frobenius error of $r$ randomly pivoted LU steps approach $4^r$ times
the optimal rank-$r$ squared Frobenius error. Thus both previously established
universal factors are sharp as suprema. The same suprema hold on
positive-definite correlation matrices when the dimension is unrestricted. In
particular, a polynomial-in-$r$ factor for same-rank RPCholesky and a $2^r$
factor for same-rank RPLU are impossible without additional hypotheses. The
proof combines scale-separated Cauchy–Binet expansions with an exact
cancellation over pivot histories. A separate $2\times2$ rational example
already disproves the literal $2^r$ RPLU bound at $r=1$.

Claims

LEAN VERIFIED OP-00115 — Polynomial trace-error factor after exactly the target rank of pivots
Disproof · fidelity not yet established
LEAN VERIFIED OP-00116 — Improve the randomized LU squared-error factor to 2^k
Disproof · fidelity not yet established

Files and sources

Paper
OpenProblemsInNLA @ 0e916df, manuscripts/sharp_random_pivoting.pdf ↗
Posted on GitHub
11 September 2026

No verification run is on record for this version.