A hierarchical low rank schur complement preconditioner for indefinite linear systems

Geoffrey Dillon, Vassilis Kalantzis, Yuanzhe Xi, Yousef Saad

Research output: Contribution to journalArticlepeer-review

7 Scopus citations

Abstract

Nonsymmetric and highly indefinite linear systems can be quite difficult to solve by iterative methods. This paper combines ideas from the multilevel Schur low rank preconditioner developed by Y. Xi, R. Li, and Y. Saad [SIAM J. Matrix Anal., 37 (2016), pp. 235–259] with classic block preconditioning strategies in order to handle this case. The method to be described generates a tree structure T that represents a hierarchical decomposition of the original matrix. This decomposition gives rise to a block structured matrix at each level of T . An approximate inverse of the original matrix based on its block LU factorization is computed at each level via a low rank property that characterizes the difference between the inverses of the Schur complement and another block of the reordered matrix. The low rank correction matrix is computed by several steps of the Arnoldi process. Numerical results illustrate the robustness of the proposed preconditioner with respect to indefiniteness for a few discretized partial differential equations and publicly available test problems.

Original languageEnglish (US)
Pages (from-to)A2234-A2252
JournalSIAM Journal on Scientific Computing
Volume40
Issue number4
DOIs
StatePublished - 2018

Bibliographical note

Publisher Copyright:
© 2018 Society for Industrial and Applied Mathematics.

Keywords

  • Block preconditioner
  • Domain decomposition
  • Krylov subspace methods
  • Low rank approximation
  • Multilevel
  • Nested dissection ordering
  • Schur complements

Fingerprint

Dive into the research topics of 'A hierarchical low rank schur complement preconditioner for indefinite linear systems'. Together they form a unique fingerprint.

Cite this