Nearly Linear Time Algorithms For Preconditioning And Solving Symmetric, Diagonally Dominant Linear Systems
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS(2014)
摘要
We present a randomized algorithm that on input a symmetric, weakly diagonally dominant n-by-n matrix A with m nonzero entries and an n-vector b produces an (x) over tilde such that parallel to(x) over tilde - A(dagger) b parallel to(A) <= epsilon parallel to A(dagger) b parallel to(A) in expected time O(m log(c) n log(1/epsilon)) for some constant c. By applying this algorithm inside the inverse power method, we compute approximate Fiedler vectors in a similar amount of time. The algorithm applies subgraph preconditioners in a recursive fashion. These preconditioners improve upon the subgraph preconditioners first introduced by Vaidya in 1990. For any symmetric, weakly diagonally dominant matrix A with nonpositive off-diagonal entries and k >= 1, we construct in time O(m log(c) n) a preconditioner B of A with at most 2(n - 1) + O((m/k) log(39) n) nonzero off-diagonal entries such that the finite generalized condition number kappa(f) (A, B) is at most k, for some other constant c. In the special case when the nonzero structure of the matrix is planar the corresponding linear system solver runs in expected time O(n log(2) n+n log n log log n log(1/epsilon)). We hope that our introduction of algorithms of low asymptotic complexity will lead to the development of algorithms that are also fast in practice.
更多查看译文
关键词
support theory,preconditioning,linear equation solvers,symmetric,diagonally dominant matrices,Laplacians
AI 理解论文
溯源树
样例
生成溯源树,研究论文发展脉络