Premium
AMLI preconditioner for the p ‐version of the FEM
Author(s) -
Beuchler Sven
Publication year - 2003
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.329
Subject(s) - preconditioner , mathematics , finite element method , discretization , conjugate gradient method , stiffness matrix , condition number , domain decomposition methods , dirichlet distribution , dirichlet boundary condition , boundary value problem , matrix (chemical analysis) , mathematical analysis , linear system , eigenvalues and eigenvectors , mathematical optimization , physics , materials science , quantum mechanics , composite material , thermodynamics
Abstract From the literature it is known that the conjugate gradient method with domain decomposition preconditioners is one of the most efficient methods for solving systems of linear algebraic equations resulting from p ‐version finite element discretizations of elliptic boundary value problems. One ingredient of such a preconditioner is a preconditioner related to the Dirichlet problems. In the case of Poisson's equation, we present a preconditioner for the Dirichlet problems which can be interpreted as the stiffness matrix K h,k resulting from the h ‐version finite element discretization of a special degenerated problem. We construct an AMLI preconditioner C h,k for the matrix K h,k and show that the condition number of C −1 h,kK h,k is independent of the discretization parameter. This proof is based on the strengthened Cauchy inequality. The theoretical result is confirmed by numerical examples. Copyright © 2003 John Wiley & Sons, Ltd.