Abstract:The {\em discrepancy} of a matrix $M \in \mathbb{R}^{d \times n}$ is given by $\mathrm{DISC}(M) := \min_{\boldsymbol{x} \in \{-1,1\}^n} \|M\boldsymbol{x}\|_\infty$. An outstanding conjecture, attributed to Komlós, stipulates that $\mathrm{DISC}(M) = O(1)$, whenever $M$ is a Komlós matrix, that is, whenever every column of $M$ lies within the unit sphere. Our main result asserts that $\mathrm{DISC}(M + R/\sqrt{d}) = O(d^{-1/2})$ holds asymptotically almost surely, whenever $M \in \mathbb{R}^{d \times n}$ is Komlós, $R \in \mathbb{R}^{d \times n}$ is a Rademacher random matrix, $d = \omega(1)$, and $n = \omega(d \log d)$. The factor $d^{-1/2}$ normalising $R$ is essentially best possible and the dependency between $n$ and $d$ is asymptotically best possible. Our main source of inspiration is a result by Bansal, Jiang, Meka, Singla, and Sinha (ICALP 2022). They obtained an assertion similar to the one above in the case that the smoothing matrix is Gaussian. They asked whether their result can be attained with the optimal dependency $n = \omega(d \log d)$ in the case of Bernoulli random noise or any other types of discretely distributed noise; the latter types being more conducive for Smoothed Analysis in other discrepancy settings such as the Beck-Fiala problem. For Bernoulli noise, their method works if $n = \omega(d^2)$. In the case of Rademacher noise, we answer the question posed by Bansal, Jiang, Meka, Singla, and Sinha. Our proof builds upon their approach in a strong way and provides a discrete version of the latter. Breaking the $n = \omega(d^2)$ barrier and reaching the optimal dependency $n = \omega(d \log d)$ for Rademacher noise requires additional ideas expressed through a rather meticulous counting argument, incurred by the need to maintain a high level of precision all throughout the discretisation process.

Smoothed analysis of the condition number under low-rank perturbations

Condition numbers for the Moore-Penrose inverse and the least squares problem involving rank-structured matrices

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

A general error analysis for randomized low-rank approximation with application to data assimilation

Spectral Condition-Number Estimation of Large Sparse Matrices

Smoothed analysis of some condition numbers

Randomized low-rank approximations beyond Gaussian random matrices

PERTURBATION ANALYSIS OF A NONLINEAR MATRIX EQUATION

Estimation of high-dimensional low-rank matrices

On the analysis of optimization with fixed-rank matrices: a quotient geometric view

Outliers and bounded rank perturbation for non-Hermitian random band matrices

Low-Rank Updates of Matrix Square Roots

Condition Numbers and Perturbation Analysis for the Tikhonov Regularization of Discrete Ill-Posed Problems.

On the numerical rank of radial basis function kernels in high dimension

New Tools for Smoothed Analysis: Least Singular Value Bounds for Random Matrices with Dependent Entries

Structured level-2 condition numbers of matrix functions

Matrix Perturbation: Davis-Kahan in the Infinity Norm

Smoothed Analysis of the Komlós Conjecture: Rademacher Noise

Low rank approximation method for perturbed linear systems with applications to elliptic type stochastic PDEs

Finite rank perturbation of non-Hermitian random matrices: heavy tail and sparse regimes

Analysis of singular subspaces under random perturbations