A power schur complement low-rank correction preconditioner for general sparse linear systems

Qingqing Zheng, Yuanzhe Xi, Yousef Saad

Research output: Contribution to journalArticlepeer-review

7 Scopus citations

Abstract

A parallel preconditioner is proposed for general large sparse linear systems that combines a power series expansion method with low-rank correction techniques. To enhance convergence, a power series expansion is added to a basic Schur complement iterative scheme by exploiting a standard matrix splitting of the Schur complement. One of the goals of the power series approach is to improve the eigenvalue separation of the preconditioner thus allowing an effective application of a low-rank correction technique. Experiments indicate that this combination can be quite robust when solving highly indefinite linear systems. The preconditioner exploits a domain-decomposition approach, and its construction starts with the use of a graph partitioner to reorder the original coefficient matrix. In this framework, unknowns corresponding to interface variables are obtained by solving a linear system whose coefficient matrix is the Schur complement. Unknowns associated with the interior variables are obtained by solving a block diagonal linear system where parallelism can be easily exploited. Numerical examples are provided to illustrate the effectiveness of the proposed preconditioner, with an emphasis on highlighting its robustness properties in the indefinite case.

Original languageEnglish (US)
Pages (from-to)659-682
Number of pages24
JournalSIAM Journal on Matrix Analysis and Applications
Volume42
Issue number4
DOIs
StatePublished - 2021

Bibliographical note

Publisher Copyright:
© 2021 Society for Industrial and Applied Mathematics

Keywords

  • Domain decomposition
  • Krylov subspace method
  • Low-rank correction
  • Parallel preconditioner
  • Power series expansion
  • Schur complement

Fingerprint

Dive into the research topics of 'A power schur complement low-rank correction preconditioner for general sparse linear systems'. Together they form a unique fingerprint.

Cite this