A rational function preconditioner for indefinite sparse linear systems

Yuanzhe XI, Yousef Saad

Research output: Contribution to journalArticlepeer-review

10 Scopus citations


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
Issue number3
StatePublished - 2017

Bibliographical note

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


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


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

Cite this