Premium
On the multilevel preconditioning of Crouzeix–Raviart elliptic problems
Author(s) -
Kraus J.,
Margenov S.,
Synka J.
Publication year - 2008
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.543
Subject(s) - mathematics , linear subspace , basis (linear algebra) , finite element method , type (biology) , pure mathematics , element (criminal law) , dimension (graph theory) , space (punctuation) , linear algebra , constraint (computer aided design) , cauchy distribution , mathematical analysis , geometry , ecology , linguistics , philosophy , physics , biology , political science , law , thermodynamics
Abstract We consider robust hierarchical splittings of finite element spaces related to non‐conforming discretizations using Crouzeix–Raviart type elements. As is well known, this is the key to the construction of efficient two‐ and multilevel preconditioners. The main contribution of this paper is a theoretical and an experimental comparison of three such splittings. Our starting point is the standard method based on differences and aggregates (DA) as introduced in Blaheta et al. ( Numer. Linear Algebra Appl. 2004; 11 :309–326). On this basis we propose a more general (GDA) splitting, which can be viewed as the solution of a constraint optimization problem (based on certain symmetry assumptions). We further consider the locally optimal (ODA) splitting, which is shown to be equivalent to the first reduce (FR) method from Blaheta et al. ( Numer. Linear Algebra Appl. 2004; 11 :309–326). This means that both, the ODA and the FR splitting, generate the same subspaces, and thus the local constant in the strengthened Cauchy–Bunyakowski–Schwarz inequality is minimal for the FR (respectively ODA) splitting. Moreover, since the DA splitting corresponds to a particular choice in the parameter space of the GDA splitting, which itself is an element in the set of all splittings for which the ODA (or equivalently FR) splitting yields the optimum, we conclude that the chain of inequalities γ 2 FR ⩽γ 2 GDA ⩽γ 2 DA ⩽3/4 holds independently of mesh and/or coefficient anisotropy. Apart from the theoretical considerations, the presented numerical results provide a basis for a comparison of these three approaches from a practical point of view. Copyright © 2007 John Wiley & Sons, Ltd.