English

Solving Sparse, Symmetric, Diagonally-Dominant Linear Systems in Time $O (m^{1.31})$

Data Structures and Algorithms 2025-10-20 v2 Numerical Analysis Numerical Analysis

Abstract

We present a linear-system solver that, given an nn-by-nn symmetric positive semi-definite, diagonally dominant matrix AA with mm non-zero entries and an nn-vector \bb\bb , produces a vector \xxt\xxt within relative distance ϵ\epsilon of the solution to A\xx=\bbA \xx = \bb in time O(m1.31log(nκf(A)/ϵ)O(1))O (m^{1.31} \log (n \kappa_{f} (A)/\epsilon)^{O (1)}), where κf(A)\kappa_{f} (A) is the log of the ratio of the largest to smallest non-zero eigenvalue of AA. In particular, log(κf(A))=O(blogn)\log (\kappa_{f} (A)) = O (b \log n), where bb is the logarithm of the ratio of the largest to smallest non-zero entry of AA. If the graph of AA has genus m2θm^{2\theta} or does not have a KmθK_{m^{\theta}} minor, then the exponent of mm can be improved to the minimum of 1+5θ1 + 5 \theta and (9/8)(1+θ)(9/8) (1+\theta). The key contribution of our work is an extension of Vaidya's techniques for constructing and analyzing combinatorial preconditioners.

Keywords

Cite

@article{arxiv.cs/0310036,
  title  = {Solving Sparse, Symmetric, Diagonally-Dominant Linear Systems in Time $O (m^{1.31})$},
  author = {Daniel A. Spielman and Shang-Hua Teng},
  journal= {arXiv preprint arXiv:cs/0310036},
  year   = {2025}
}

Comments

fixed a typo on page 9

R2 v1 2026-07-22T12:21:30.145Z