Computer Science Technical Reports
CS at VT

A New Adaptive GMRES Algorithm for Achieving High Accuracy

Sosonkina, Maria and Watson, Layne T and Kapania, Rakesh K (1996) A New Adaptive GMRES Algorithm for Achieving High Accuracy. Technical Report ncstrl.vatech_cs//TR-96-09, Computer Science, Virginia Polytechnic Institute and State University.

Full text available as:
Postscript - Requires a viewer, such as GhostView (380796)


GMRES(k) is widely used for solving nonsymmetric linear systems. However, it is inadequate either when it converges only for k close to the problem size or when numerical error in the modified Gram-Schmidt process used in the GMRES orthogonalization phase dramatically affects the algorithm performance. An adaptive version of GMRES(k) which tunes the restart value k based on criteria estimating the GMRES conversion rate for the given problem is proposed here. This adaptive GMRES(k) procedure outperforms standard GMRES(k), several other GMRES-like methods, and QMR on actual large scale sparse structural mechanics postbuckling and analog circuit simulation problems. There are some applications, such as homotopy methods for high Reynolds number viscous flows, solid mechanics postbuckling analysis, and analog circuit simulation, where very high accuracy in the linear system solutions is essential. In this context, the modified Gram-Schmidt process in GMRES can fail causing the entire GMRES iteration to fail. It is shown that the adaptive GMRES(k) with the orthogonalization performed by Householder transformations succeeds whenever GMRES(k) with the orthogonolization performed by the modified Gram-Schmidt process fails, and the extra cost of computing Householder transformations is justified for these applications.

Item Type:Departmental Technical Report
Subjects:Computer Science > Historical Collection(Till Dec 2001)
ID Code:448
Deposited By:User autouser
Deposited On:05 December 2001
Alternative Locations: URL:, URL: