ISSN:
1070-5325
Keywords:
iterative methods for linear systems
;
acceleration of convergence
;
preconditioning
;
Engineering
;
Numerical Methods and Modeling
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Mathematics
Notes:
Considering matrices obtained by the application of a five-point stencil on a 2D rectangular grid, we analyse a preconditioning method introduced by Axelsson and Eijkhout, and by Brand and Heinemann. In this method, one performs a (modified) incomplete factorization with respect to a so-called ‘repeated’ or ‘recursive’ red-black ordering of the unknowns while fill-in is accepted provided that the red unknowns in a same level remain uncoupled.Considering discrete second order elliptic PDEs with isotropic coefficients, we show that the condition number is bounded by O(n½ log2(√(5) -1)) where n is the total number of unknowns (½ log2(√(5) - 1) = 0.153), and thus, that the total arithmetic work for the solution is bounded by O(n1.077). Our condition number estimate, which turns out to be better than standard O(log2 n) estimates for any realistic problem size, is purely algebraic and holds in the presence of Neumann boundary conditions and/or discontinuities in the PDE coefficients.Numerical tests are reported, displaying the efficiency of the method and the relevance of our analysis. © 1997 John Wiley & Sons, Ltd.
Additional Material:
2 Ill.
Type of Medium:
Electronic Resource
Permalink