English

Well-Conditioned Oblivious Perturbations in Linear Space

Data Structures and Algorithms 2026-04-28 v1 Machine Learning Numerical Analysis Numerical Analysis Probability Machine Learning

Abstract

Perturbing a deterministic nn-dimensional matrix with small Gaussian noise is a cornerstone of smoothed analysis of algorithms [Spielman and Teng, JACM 2004], as it reduces the condition number of the input to O(n)O(n), and with it the complexity of many matrix algorithms. However, when deployed algorithmically, these perturbations are expensive due to the cost of generating and storing n2n^2 Gaussian random variables. We propose a perturbation that requires generating and storing O(n)O(n) random numbers in O(logn)O(\log n) bits of precision, and reduces the condition number of any deterministic matrix to O(n)O(n), matching Gaussian perturbations. Our result in particular implies a better complexity for the perturbed conjugate gradient algorithm, showing that we can solve an n×nn\times n linear system in linear space to within an arbitrarily small constant backward error using O(n)O(n) matrix-vector products. In our construction, we introduce the concept of a pattern matrix, which is a dense deterministic matrix that maps all sparse vectors into dense vectors, and we combine it with a sparse perturbation whose entries are dependent and located in a non-uniform fashion. In order to analyze this construction, we develop new techniques for lower bounding the smallest singular value of a random matrix with dependent entries.

Keywords

Cite

@article{arxiv.2604.23193,
  title  = {Well-Conditioned Oblivious Perturbations in Linear Space},
  author = {Shabarish Chenakkod and Michał Dereziński and Xiaoyu Dong and Mark Rudelson},
  journal= {arXiv preprint arXiv:2604.23193},
  year   = {2026}
}
R2 v1 2026-07-01T12:34:55.138Z