Preconditioned iterative methods are widely used to solve linear systems such as those arising from the finite element formulation of boundary value problems and approximate factorizations are widely used as preconditioners. The ordering of the unknowns is therefore an important issue because it has a strong influence on the convergence behaviour of the iteration method while it is also a decisive aspect for their parallel implementation. Consistent orderings are attractive for parallel implementations and it has been shown that some subclasses of these orderings also enhance the convergence behaviour of the associated iteration methods. This has in particular been shown for the so‐called S/P consistent orderings. A wider definition of this class of orderings has recently been proposed and we investigate here how approximate factorizations should be implemented when using such more general orderings (still called S/P consistent) in order to keep their expected high convergence properties. A simple practical conclusion is suggested, supported by both theoretical and numerical arguments. Copyright © 1994 John Wiley & Sons, Ltd
[1]
Yvan Notay,et al.
S/P images of upper triangular M-matrices
,
1994,
Numer. Linear Algebra Appl..
[2]
Y. Notay.
Solving positive (semi) definite linear systems by preconditioned iterative methods
,
1991
.
[3]
J. Pasciak,et al.
Computer solution of large sparse positive definite systems
,
1982
.
[4]
H. V. D. Vorst,et al.
The convergence behaviour of preconditioned CG and CG-S in the presence of rounding errors
,
1991
.
[5]
R. Beauwens,et al.
Conditioning analysis of positive definite matrices by approximate factorizations
,
1989
.
[6]
I. Duff,et al.
The effect of ordering on preconditioned conjugate gradients
,
1989
.
[7]
R. Beauwens.
Approximate factorizations with S/P consistently ordered M-factors
,
1989
.