View metadata, citation and similar papers at core.ac.uk brought to you by CORE
provided by Elsevier - Publisher Connector
Linear Algebra and its Applications 434 (2011) 2308–2324
Contents lists available at ScienceDirect
Linear Algebra and its Applications
journal homepage: www.elsevier.com/locate/laa
On element-by-element Schur complement approximations Maya Neytcheva
Uppsala University, Department of Information Technology, Box 337, SE-751 05 Uppsala, Sweden
ARTICLE INFO ABSTRACT
Article history: We discuss a methodology to construct sparse approximations of Received 19 September 2009 Schur complements of two-by-two block matrices arising in finite Accepted 22 February 2010 element discretizations of partial differential equations. Earlier re- sults from [2] are extended to more general symmetric positive Submitted by V. Mehrmann definite matrices of two-by-two block form. The applicability of the method for general symmetric and nonsymmetric matrices is Keywords: analysed. The paper demonstrates the applicability of the presented Iterative methods method providing extensive numerical experiments. Two-by-two block structured matrices © 2010 Elsevier Inc. All rights reserved. Symmetric Nonsymmetric Schur complements Sparse matrix approximations based on finite element discretizations
1. Introduction
This paper deals with general real symmetric or nonsymmetric sparse matrices which admit block- partitioning of the form A A }N A = 11 12 1 (1) A21 A22 }N2 = − −1 and the associated Schur complement matrix SA A22 A21A11 A12 , under the assumption that the block A11 is invertible. There is a plethora of applications where the arising matrices admit the form (1), such as discretized compressible and incompressible flow problems, constraint optimization problems, linear elasticity, to name a few. The task to construct a good approximation of SA is an important issue when constructing preconditioners for matrices of the above type. It is known that preconditioners which utilize the block structure of A are among the best performing ones, however, for their construction we must
E-mail address: [email protected]
0024-3795/$ - see front matter © 2010 Elsevier Inc. All rights reserved. doi:10.1016/j.laa.2010.03.031 M. Neytcheva / Linear Algebra and its Applications 434 (2011) 2308–2324 2309 approximate the arising Schur complement. The literature on how to solve and how to precondition such matrices, in particular saddle point problems, is nearly inexhaustible. We refer to [15]asarich and up-today source. One can look upon two-by-two block structured matrices from a different point of view. Namely, if we impose such as structure on a matrix, then its exact block-factorization − A 0 I A 1A A = 11 1 11 12 A21 SA 0 I2 provides possibilities to construct highly efficient preconditioners of block-factorized form by approx- imating some of the matrix blocks. Thus, the need to approximate Schur complements appears again. Furthermore, the block-factorized form is easily reframed as a multilevel procedure, where we would need to approximate Schur complements of Schur complements in a recursive manner, provided that the Schur complement itself can be split into two-by-two block form. One established class of methods of this type is the so-called Algebraic Multilevel Iteration (AMLI) (cf., e.g. [6] and more recently [28,20], to follow their development). Theoretical results as well as numerical tests have shown that good preconditioners require high quality approximations of SA. At the same time, Schur complements are in general dense matrices and could have large dimensions. For some special cases, good preconditioners of SA are known, e.g. when A11 is diagonal or could be approximated by a diagonal matrix, the pressure mass matrix for the Schur −1 complement of the Stokes matrix, etc. The search for accurate approximations of SA or directly of SA has lead to very efficient but highly computationally demanding preconditioners, for example based on the properties of the underlying partial differential operators (see some examples in [27,14] and the references therein). Still, the known recipes are problem- or discretization-dependent and, in general, the question how to approximate Schur complement matrices remains open. We consider here a discretization-based approach how to construct a high quality sparse approxi- mation of SA, namely, applicable for matrices A, arising in the framework of partial differential equations (PDE) discretized by a standard (conforming) finite element method (FEM). There exists already some research on this topic in [19,24,23,2,25]. This paper is based on [2,25]. It can be viewed as an upgrade of [2]. There, only finite element matrices corresponding to second order selfadjoint elliptic problems are considered. Those are dis- cretized on a fine mesh, obtained by a regular refinement of some given coarse mesh and the matrices are partitioned in the sets ‘added fine mesh nodepoints’ and ‘coarse nodepoints’. Here we consider a more general setting and discuss, test and develop further some ideas, described briefly in [25], also for more general nonsymmetric matrices. The paper has the following structure. Section 2 introduces some notations and the context in which we construct approximate Schur complement matrices. In Section 3 we derive some estimates to quantify the quality of the approximated Schur complements for some classes of matrices. Section 4 contains numerical illustrations of the applicability of the approach. Some conclusions and questions to be addressed in a future research are given in Section 5.
2. Notation and preliminaries
As already stated above, we consider matrices partitioned in a two-by-two block form. These matrices are assumed to arise from FEM discretizations of scalar or vector partial differential equations, or systems of equations. After discretization, a linear system of equations arises, denoted by Au = f, × where A ∈ RN N is supposedly a large and sparse matrix and therefore the system is to be solved using a preconditioned iterative method. Our final goal is to construct a robust and computationally efficient preconditioner for A based on some block-factorized approximation of it. For our purposes, we view the system matrix A as an assembly of local finite element matrices, which themselves admit a two-by-two block structure, namely, 2310 M. Neytcheva / Linear Algebra and its Applications 434 (2011) 2308–2324
Fig. 1. Two macroelements in a mesh.