Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

A block variant of the GMRES method on massively parallel processors

Conference ·
 [1]
  1. Cray Research, Inc., Eagan, MN (United States)
This paper presents a block variant of the GMRES method for solving general unsymmetric linear systems. This algorithm generates a transformed Hessenberg matrix by solely using block matrix operations and block data communications. It is shown that this algorithm with block size s, denoted by BVGMRES(s,m), is theoretically equivalent to the GMRES(s*m) method. The numerical results show that this algorithm can be more efficient than the standard GMRES method on a cache based single CPU computer with optimized BLAS kernels. Furthermore, the gain in efficiency is more significant on MPPs due to both efficient block operations and efficient block data communications. Our numerical results also show that in comparison to the standard GMRES method, the more PEs that are used on an MPP, the more efficient the BVGMRES(s,m) algorithm is.
Research Organization:
Front Range Scientific Computations, Inc., Lakewood, CO (United States)
OSTI ID:
433335
Report Number(s):
CONF-9604167--Vol.1; ON: DE96015306
Country of Publication:
United States
Language:
English

Similar Records

On improving linear solver performance: a block variant of GMRES
Journal Article · Mon May 10 00:00:00 EDT 2004 · SIAM Journal on Scientific Computing · OSTI ID:881896

Optimal size of the block in block GMRES on GPUs: computational model and experiments
Journal Article · Mon Dec 12 19:00:00 EST 2022 · Numerical Algorithms · OSTI ID:2311786

The Massively Parallel Processor
Book · Mon Dec 31 23:00:00 EST 1984 · OSTI ID:6916716