The standard implementation of the hybrid GMRES algorithm for solving large nonsymmetric linear systems involves a Gram-Schmidt process which is a potential source of significant numerical error. An alternative implem...The standard implementation of the hybrid GMRES algorithm for solving large nonsymmetric linear systems involves a Gram-Schmidt process which is a potential source of significant numerical error. An alternative implementation is outlined here in which orthogonalization by Householder transformations replaces the Gram-Schmidt process. Numerical experiments show that the new implementation is more stable.展开更多
文摘The standard implementation of the hybrid GMRES algorithm for solving large nonsymmetric linear systems involves a Gram-Schmidt process which is a potential source of significant numerical error. An alternative implementation is outlined here in which orthogonalization by Householder transformations replaces the Gram-Schmidt process. Numerical experiments show that the new implementation is more stable.