A rational function preconditioner for indefinite sparse linear systems

Yuanzhe XI, Yousef Saad

Research output: Contribution to journalArticlepeer-review

5 Scopus citations

Abstract

This paper introduces a rational function preconditioner for linear systems with indefinite sparse matrices A. By resorting to rational functions of A, the algorithm decomposes the spectrum of A into two disjoint regions and approximates the restriction of A1 on these regions separately. We show a systematic way to construct these rational functions so that they can be applied stably and inexpensively. An attractive feature of the proposed approach is that the construction and application of the preconditioner can exploit two levels of parallelism. Moreover, the proposed preconditioner can be modified at a negligible cost into a preconditioner for a nearby matrix of the form AI, which can be useful in some applications. The efficiency and robustness of the proposed preconditioner are demonstrated on a few tests with challenging model problems, including problems arising from the Helmholtz equation in three dimensions.

Original languageEnglish (US)
Pages (from-to)A1145-A1167
JournalSIAM Journal on Scientific Computing
Volume39
Issue number3
DOIs
StatePublished - 2017

Bibliographical note

Funding Information:
This work was supported by NSF grants DMS-1216366 and DMS-1521573 and by the Minnesota Supercomputing Institute.

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

Keywords

  • Approximate inverse
  • Cauchy integral
  • De ation
  • Helmholtz equation
  • Incomplete LU
  • Rational function

Fingerprint Dive into the research topics of 'A rational function preconditioner for indefinite sparse linear systems'. Together they form a unique fingerprint.

Cite this