Finite Element Conjugate Gradient Method for the Solution of the Pure Neumann Boundary Value Problem

摘要

An important problem in numerical analysis is the solution of the pure Neumann boundary value problem, for which the solution is unique to within a constant term. In the report the finite element and conjugate gradient method are used to solve the boundary value problem. The resulting singular system of linear equations is preconditioned by an incomplete Cholesky decomposition. The existence of the decomposition is proved for singular irreducible M-matrices. At the end of the report an application to groundwater flow is described. (Copyright (c) 1987 by the Faculty of Mathematics and Informatics, Delft, the Netherlands.)