On Solving Linear Systems in Sublinear Time
Abstract
We study \emph{sublinear} algorithms that solve linear systems locally. In the classical version of this problem the input is a matrix and a vector in the range of , and the goal is to output satisfying . For the case when the matrix is symmetric diagonally dominant (SDD), the breakthrough algorithm of Spielman and Teng [STOC 2004] approximately solves this problem in near-linear time (in the input size which is the number of non-zeros in ), and subsequent papers have further simplified, improved, and generalized the algorithms for this setting. Here we focus on computing one (or a few) coordinates of , which potentially allows for sublinear algorithms. Formally, given an index together with and as above, the goal is to output an approximation for , where is a fixed solution to . Our results show that there is a qualitative gap between SDD matrices and the more general class of positive semidefinite (PSD) matrices. For SDD matrices, we develop an algorithm that approximates a single coordinate in time that is polylogarithmic in , provided that is sparse and has a small condition number (e.g., Laplacian of an expander graph). The approximation guarantee is additive for accuracy parameter . We further prove that the condition-number assumption is necessary and tight. In contrast to the SDD matrices, we prove that for certain PSD matrices , the running time must be at least polynomial in . This holds even when one wants to obtain the same additive approximation, and has bounded sparsity and condition number.
Cite
@article{arxiv.1809.02995,
title = {On Solving Linear Systems in Sublinear Time},
author = {Alexandr Andoni and Robert Krauthgamer and Yosef Pogrow},
journal= {arXiv preprint arXiv:1809.02995},
year = {2026}
}