English

Ultrasparse Ultrasparsifiers and Faster Laplacian System Solvers

Data Structures and Algorithms 2023-04-04 v2 Optimization and Control

Abstract

In this paper we provide an O(m(loglogn)O(1)log(1/ϵ))O(m (\log \log n)^{O(1)} \log(1/\epsilon))-expected time algorithm for solving Laplacian systems on nn-node mm-edge graphs, improving improving upon the previous best expected runtime of O(mlogn(loglogn)O(1)log(1/ϵ))O(m \sqrt{\log n} (\log \log n)^{O(1)} \log(1/\epsilon)) achieved by (Cohen, Kyng, Miller, Pachocki, Peng, Rao, Xu 2014). To obtain this result we provide efficient constructions of p\ell_p-stretch graph approximations with improved stretch and sparsity bounds. Additionally, as motivation for this work, we show that for every set of vectors in Rd\mathbb{R}^d (not just those induced by graphs) and all k>1k > 1 there exist ultrasparsifiers with d1+O(d/k)d-1 + O(d/\sqrt{k}) re-weighted vectors of relative condition number at most kk. For small kk, this improves upon the previous best known relative condition number of O~(klogd)\tilde{O}(\sqrt{k \log d}), which is only known for the graph case.

Keywords

Cite

@article{arxiv.2011.08806,
  title  = {Ultrasparse Ultrasparsifiers and Faster Laplacian System Solvers},
  author = {Arun Jambulapati and Aaron Sidford},
  journal= {arXiv preprint arXiv:2011.08806},
  year   = {2023}
}

Comments

56 pages. Updated version includes slightly improved running time and new sparsity bounds for graph ultrasparsifiers