Approaching optimality for solving SDD systems

Abstract

We present an algorithm that on input of an n-vertex m-edge weighted graph G and a value k, produces an incremental sparsifier G with n-1 + m/k edges, such that the condition number of G with G is bounded above by O(k2 n), with probability 1-p. The algorithm runs in time O((m n + n2n)(1/p)). As a result, we obtain an algorithm that on input of an n× n symmetric diagonally dominant matrix A with m non-zero entries and a vector b, computes a vector x satisfying ||x-A+b||A<ε ||A+b||A , in expected time O(m2n(1/ε)). The solver is based on repeated applications of the incremental sparsifier that produces a chain of graphs which is then used as input to a recursive preconditioned Chebyshev iteration.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…