Abstract:
The 3-term Lanczos process leads, for a symmetric matrix, to bases for Krylov subspaces of increasing dimension. The Lanczos basis, together with the recurrence coefficients, can be used for the solution of symmetric indefinite linear systems, by solving the reduced system in one way or another. This leads to well-known methods: MINRES, GMRES, and SYMMLQ. We will discuss in what way and to what extent these approaches differ in their sensitivity to rounding errors. In our analysis we will assume that the Lanczos basis is generated in exactly the same way for the different methods, and we will not consider the errors in the Lanczos process itself. We will show that the method of solution may lead, under certain circumstances, to large additional errors, that are not corrected by continuing the iteration process. Our findings are supported and illustrated by numerical examples. 1 Introduction We will consider iterative methods for the construction of approximate solutions, starting with...
Citations
|
898
|
GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems
– Saad, Schultz
- 1986
|
|
558
|
The Symmetric Eigenvalue Problem
– Parlett
|
|
506
|
Accuracy and Stability of Numerical Algorithms
– Higham
- 2002
|
|
417
|
Methods of conjugate gradients for solving linear systems
– Hestenes, Stiefel
- 1952
|
|
388
|
der Vorst. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods
– Barrett, Berry, et al.
- 1994
|
|
203
|
Iterative Methods for Solving Linear Systems
– Greenbaum
- 1997
|
|
176
|
Solution of sparse indefinite systems of linear equations
– PAIGE, SAUNDERS
|
|
104
|
Matrix Computations, Third Edition
– Golub
- 1996
|
|
62
|
Condition numbers and equilibration of matrices
– SLUIS
- 1969
|
|
51
|
der Vorst, Approximate solutions and eigenvalue bounds from Krylov subspaces
– Paige, Parlett, et al.
- 1995
|
|
48
|
A theoretical comparison of the Arnoldi and GMRES algorithms
– Brown
- 1991
|
|
39
|
Polynomial Based Iteration Methods for Symmetric Linear Systems. Ad-vances in Numerical Mathematics
– Fischer
- 1998
|
|
32
|
Error analysis of the Lanczos algorithm for tridiagonalizing a symmetric matrix
– Paige
- 1976
|
|
29
|
A survey of preconditioned iterative methods
– Bruaset
- 1995
|
|
29
|
The superlinear convergence behaviour of GMRES
– Vorst, Vuik
- 1993
|
|
28
|
Behavior of slightly perturbed Lanczos and conjugate-gradient
– Greenbaum
- 1989
|
|
24
|
Accuracy and effectiveness of the Lanczos algorithm for the symmetric eigenproblem
– Paige
- 1980
|
|
21
|
s, Predicting the behavior of finite precision Lanczos and conjugate gradient computations
– Greenbaum, Strako
- 1992
|
|
17
|
E#cient high accuracy solutions with GMRES(m
– Turner, Walker
- 1992
|
|
17
|
der Vorst, Reliable updated residuals in hybrid Bi-CG methods
– Sleijpen, Van
- 1994
|
|
16
|
Relaxationsmethoden bester Strategie zur Losung linearer Gleichungssysteme
– Stiefel
- 1955
|
|
12
|
The superlinear convergence behaviour of
– Vorst, Vuik
- 1993
|
|
4
|
A Survey of Preconditioned Iterative
– Bruaset
- 1995
|
|
3
|
Rozlo zn' ik, and Z. Strako s, Numerical stability of GMRES
– a, Greenbaum, et al.
- 1995
|
|
3
|
Accuracy and e#ectiveness of the Lanczos algorithm for the symmetric eigenproblem
– Paige
- 1980
|
|
2
|
Condition numbers and equilibration of matrices
– Sluis
- 1969
|
|
1
|
Rozlo zn k, and Z. Strako s, Numerical stability of GMRES
– a, Greenbaum, et al.
- 1995
|