Valérie Frayssé, Luc Giraud, Serge Gratton, Julien Langou

In this article we describe our implementations of the GMRES algorithm for both real and complex, single and double precision arithmetics suitable for serial, shared memory and distributed memory computers. For the sake of portability, simplicity, flexibility and efficiency the GMRES solvers have been implemented in Fortran 77 using the reverse communication mechanism for the matrix-vector product, the preconditioning and the dot product computations. For distributed memory computation, several orthogonalization procedures have been implemented to reduce the cost of the dot product calculation, which is a well-known bottleneck of efficiency for the Krylov methods. Either implicit or explicit calculation of the residual at restart are possible depending on the actual cost of the matrix-vector product. Finally the implemented stopping criterion is based on a normwise backward error.

ArticleIn this article we describe our implementations of the GMRES algorithm for both real and complex, single and double precision arithmetics suitable for serial, shared memory and distributed memory computers. For the sake of portability, simplicity, flexibility and efficiency the GMRES solvers have been implemented in Fortran 77 using the reverse communication mechanism for the matrix-vector product, the preconditioning and the dot product computations. For distributed memory computation, several orthogonalization procedures have been implemented to reduce the cost of the dot product calculation, which is a well-known bottleneck of efficiency for the Krylov methods. Either implicit or explicit calculation of the residual at restart are possible depending on the actual cost of the matrix-vector product. Finally the implemented stopping criterion is based on a normwise backward error.

Valérie Frayssé,
Luc Giraud,
Serge Gratton,
Julien Langou,
et al.
" Algorithm 842: A set of GMRES routines for real and complex arithmetics on high performance computers."
*Journal ACM Transactions on Mathematical Software (TOMS)*.
doi:10.1145/1067967.1067970.
Retrieved 06/20/2019 from researchcompendia.org/compendia/2013.4/

**Compendium Type**: Published Papers

**Primary Research Field**: Computer and Information Sciences

**Secondary Research Field**: Mathematics

**Content License**: Public Domain Mark

**Code License**: MIT License

created 12/12/2013

modified 01/16/2014

blog comments powered by Disqus