2025/04/29 by Vladu, Adrian
#Data Structures and Algorithms (cs.DS) #FOS: Computer and information sciences #FOS: Mathematics #Optimization and Control (math.OC)
paper · doi:10.48550/arxiv.2504.20619
We study two fundamental optimization problems: (1) scaling a symmetric positive definite matrix by a positive diagonal matrix so that the resulting matrix has row and column sums equal to 1; and (2) minimizing a quadratic function subject to hard non-negativity constraints. Both problems lend themselves to efficient algorithms based on interior point methods (IPMs). For general instances, standard self-concordance theory places a limit on the iteration complexity of these methods at \widetildeO(n1/2), where n denotes the matrix dimension. We show via an amortized analysis that, when the input matrix is an M-matrix, an IPM with adaptive step sizes solves both problems in only \widetildeO(n1/3) iterations. As a corollary, using fast Laplacian solvers, we obtain an ℓ2 flow diffusion algorithm with depth \widetildeO(n1/3) and work \widetildeO(n1/3\cdotnnz). This result marks a significant instance in which a standard log-barrier IPM permits provably fewer than Θ(n1/2) iterations.