A block variant of the GMRES method on massively parallel processors
- 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
Optimal size of the block in block GMRES on GPUs: computational model and experiments
The Massively Parallel Processor
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