# Deflation in preconditioned conjugate gradient methods for Finite Element Problems

Deflation in preconditioned conjugate gradient methods for Finite Element Problems

Author Faculty Date2002

AbstractWe investigate the influence of the value of deflation vectors at interfaces on the rate of convergence of preconditioned conjugate gradient methods applied to a Finite Element discretization for an elliptic equation. Our set-up is a Poisson problem in two dimensions with continuous or discontinuous coefficients that vary in several orders of magnitude. In the continuous case we are interested in the convergence acceleration of deflation on block preconditioners. For the discontinuous case the Finite Element discretization gives a symmetric matrix with very large condition number and hence many iterations are needed to obtain a solution using conjugate gradients. We use an incomplete Choleski decomposition for the symmetric discretization matrix. Subsequently, deflation is applied to eliminate the disadvantageous effects to convergence caused by the remaining small eigenvalues. Here we investigate the effects of several variants of deflation and we propose an optimal choice for the deflation vectors. Finally, we apply deflation to a parallel preconditioned conjugate method to accellerate the calculation.

Subjectporous media

Laplace equation

finite differences

conjugated gradients

preconditioning

deflation

http://resolver.tudelft.nl/uuid:f48941d8-0c37-444f-8257-b56118dff9ee

PublisherDelft University of Technology, Faculty of Electrical Engineering, Mathematics and Computer Science, Delft Institute of Applied Mathematics

SourceReports of the Department of Applied Mathematical Analysis, 02-10

ISSN1389-6520

Part of collectionInstitutional Repository

Document typereport

Rights(c) 2002 Department of Applied Mathematical Analysis