Preconditioning and two-level multigrid methods of arbitrary degree of approximation

Mathematics of Computation - Tập 40 Số 161 - Trang 219-242
Owe Axelsson, Ivar Gustafsson

Tóm tắt

Let h be a mesh parameter corresponding to a finite element mesh for an elliptic problem. We describe preconditioning methods for two-level meshes which, for most problems solved in practice, behave as methods of optimal order in both storage and computational complexity. Namely, per mesh point, these numbers are bounded above by relatively small constants for all h ⩾ h 0 h \geqslant {h_0} , where h 0 {h_0} is small enough to cover all but excessively fine meshes. We note that, in practice, multigrid methods are actually solved on a finite, often even a fixed number of grid levels, in which case also these methods are not asymptotically optimal as h → 0 h \to 0 . Numerical tests indicate that the new methods are about as fast as the best implementations of multigrid methods applied on general problems (variable coefficients, general domains and boundary conditions) for all but excessively fine meshes. Furthermore, most of the latter methods have been implemented only for difference schemes of second order of accuracy, whereas our methods are applicable to higher order approximations. We claim that our scheme could be added fairly easily to many existing finite element codes.

Từ khóa


Tài liệu tham khảo

Axelsson, O., 1976, A class of iterative methods for finite element equations, Comput. Methods Appl. Mech. Engrg., 9, 123, 10.1016/0045-7825(76)90056-6

O. Axelsson & I. Gustafsson, A Preconditioned Conjugate Gradient Method for Finite Element Equations, Which is Stable for Rounding Errors, Information Processing 80 (S. H. Lavington, ed.), North-Holland, Amsterdam, 1980, pp. 723-728.

R. Bank & T. Dupont, Analysis of a Two-Level Scheme for Solving Finite Element Equations, Report CNA-159, Center for Numerical Analysis, The University of Texas at Austin, 1980.

Brandt, Achi, 1977, Multi-level adaptive solutions to boundary-value problems, Math. Comp., 31, 333, 10.2307/2006422

I. Fried, "Bounds on the extremal eigenvalues of the finite element stiffness and mass matrices and their spectral condition numbers," J. Sound Vibration, v. 22, 1972, pp. 407-418.

I. Gustafsson, Stability and Rate of Convergence of Modified Incomplete Cholesky Factorization Methods, Thesis, Report 79.02R, Department of Computer Sciences, Chalmers University of Technology, Göteborg, Sweden, 1979.

Hemker, P. W., 1980, Introduction to multigrid methods, 59

Kershaw, David S., 1978, The incomplete Cholesky-conjugate gradient method for the iterative solution of systems of linear equations, J. Comput. Phys., 26, 43, 10.1016/0021-9991(78)90098-0

Meijerink, J. A., 1977, An iterative solution method for linear systems of which the coefficient matrix is a symmetric 𝑀-matrix, Math. Comp., 31, 148, 10.2307/2005786

U. Trottenberg, Private communication, 1981.

Thông tin
Thông tin xuất bản

Mathematics of Computation

Tập 40 Số 161

219-242

Thông tin tác giả