skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Toward efficient polynomial preconditioning for GMRES

Journal Article · · Numerical Linear Algebra with Applications
DOI:https://doi.org/10.1002/nla.2427· OSTI ID:1838187

Here, we present a polynomial preconditioner for solving large systems of linear equations. The polynomial is derived from the minimum residual polynomial (the GMRES polynomial) and is more straightforward to compute and implement than many previous polynomial preconditioners. Our current implementation of this polynomial using its roots is naturally more stable than previous methods of computing the same polynomial. We implement further stability control using added roots, and this allows for high degree polynomials. We discuss the effectiveness and challenges of root-adding and give an additional check for stability. In this article, we study the polynomial preconditioner applied to GMRES; however it could be used with any Krylov solver. This polynomial preconditioning algorithm can dramatically improve convergence for some problems, especially for difficult problems, and can reduce dot products by an even greater margin.

Research Organization:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA); National Science Foundation (NSF)
Grant/Contract Number:
NA0003525; DMS-1418677
OSTI ID:
1838187
Report Number(s):
SAND2021-16131J; 702364
Journal Information:
Numerical Linear Algebra with Applications, Vol. 29, Issue 4; ISSN 1070-5325
Publisher:
WileyCopyright Statement
Country of Publication:
United States
Language:
English

References (27)

A stable Richardson iteration method for complex linear systems journal March 1989
Chebyshev polynomials in the solution of large-scale linear systems conference January 1952
A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems journal January 2016
A Robust GMRES-Based Adaptive Polynomial Preconditioning Algorithm for Nonsymmetric Linear Systems journal March 1994
Polynomial Preconditioned GMRES and GMRES-DR journal January 2015
A Restarted GMRES Method Augmented with Eigenvectors journal October 1995
Harmonic projection methods for large non-symmetric eigenvalue problems journal January 1998
A Hybrid GMRES Algorithm for Nonsymmetric Linear Systems
  • Nachtigal, Noël M.; Reichel, Lothar; Trefethen, Lloyd N.
  • SIAM Journal on Matrix Analysis and Applications, Vol. 13, Issue 3 https://doi.org/10.1137/0613050
journal July 1992
An Iterative Solution Method for Linear Systems of Which the Coefficient Matrix is a Symmetric M-Matrix journal January 1977
Improved seed methods for symmetric positive definite linear equations with multiple right-hand sides journal July 2013
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems journal July 1986
Iterative Methods for Sparse Linear Systems book January 2003
GMRESR: a family of nested GMRES methods journal July 1994
Ritz and harmonic Ritz values and the convergence of FOM and GMRES journal June 1999
Leapfrog variants of iterative methods for linear algebraic equations journal November 1988
A Flexible Inner-Outer Preconditioned GMRES Algorithm journal March 1993
A harmonic restarted Arnoldi algorithm for calculating eigenvalues and determining multiplicity journal May 2006
Least Squares Polynomials in the Complex Plane and Their Use for Solving Nonsymmetric Linear Systems journal February 1987
An optimum iterative method for solving any linear system with a square matrix journal March 1988
Chebyshev acceleration techniques for solving nonsymmetric eigenvalue problems journal May 1984
A polynomial preconditioner for the GMRES algorithm journal April 1995
A Comparison of Adaptive Chebyshev and Least Squares Polynomial Preconditioning for Hermitian Positive Definite Linear Systems
  • Ashby, Steven F.; Manteuffel, Thomas A.; Otto, James S.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 1 https://doi.org/10.1137/0913001
journal January 1992
Approximate solutions and eigenvalue bounds from Krylov subspaces journal March 1995
A Newton basis GMRES implementation journal January 1994
Avoiding communication in sparse matrix computations
  • Demmel, James; Hoemmen, Mark; Mohiyuddin, Marghoob
  • Distributed Processing Symposium (IPDPS), 2008 IEEE International Symposium on Parallel and Distributed Processing https://doi.org/10.1109/IPDPS.2008.4536305
conference April 2008
Computing interior eigenvalues of large matrices journal August 1991
Polynomial Preconditioned Arnoldi with Stability Control journal January 2021

Similar Records

A parallel additive Schwarz preconditioned Jacobi-Davidson algorithm for polynomial eigenvalue problems in quantum dot simulation
Journal Article · Tue Apr 20 00:00:00 EDT 2010 · Journal of Computational Physics · OSTI ID:1838187

Polynomial Preconditioned Arnoldi with Stability Control
Journal Article · Mon Jan 04 00:00:00 EST 2021 · SIAM Journal on Scientific Computing · OSTI ID:1838187

Performance of preconditioned iterative solvers in MFiX–Trilinos for fluidized beds
Journal Article · Fri May 18 00:00:00 EDT 2018 · Journal of Supercomputing · OSTI ID:1838187