A Multigrid Method for Nonconforming FE-Discretisations with Application to Non-Matching Grids

Nonconforming finite element discretisations require special care in the construction of the prolongation and restriction in the multigrid process. In this paper, a general scheme is proposed, which guarantees the approximation property. As an example, the technique is applied to the discretisation by non-matching grids (mortar elements). Received: October 15, 1998     

Problem Dependent Generalized Prewavelets

In this paper, we present a new approach to construct robust multilevel algorithms for elliptic differential equations. The multilevel algorithms consist of multiplicative subspace corrections in spaces spanned by problem dependent generalized prewavelets. These generalized prewavelets are constructed by a local orthogonalization of hierarchical basis functions with respect to a so-called local coarse-grid space. Numerical results show that the local orthogonalization leads to a smaller constant in strengthened Cauchy-Schwarz inequality than the original hierarchical basis functions. This holds also for several equations with discontinuous coefficients. Thus, the corresponding multilevel algorithm is a fast and robust iterative solver. Received November 13, 2001; revised October 21, 2002 Published online: December 12, 2002     

Application of an Adaptive Sparse-Grid Technique to a Model Singular Perturbation Problem

In this paper we show how, under minimal conditions, a combination extrapolation can be introduced for an adaptive sparse grid. We apply this technique for the solution of a two-dimensional model singular perturbation problem, defined on the domain exterior of a circle. Received October 18, 1999     

The Streamline–Diffusion Method for Conforming and Nonconforming Finite Elements of Lowest Order Applied to Convection–Diffusion Problems

We consider the streamline-diffusion finite element method with finite elements of lowest order for solving convection-diffusion problems. Our investigations cover both conforming and nonconforming finite element approximations on triangular and quadrilateral meshes. Although the considered finite elements are of the same interpolation order their stability and approximation properties are quite different. We give a detailed overview on the stability and the convergence properties in the L ²- and in the streamline–diffusion norm. Numerical experiments show that often the theoretical predictions on the convergence properties are sharp. Received December 7, 1999; revised October 5, 2000     

The Inf-Sup Condition for the Mapped Q k ?P k?1 disc Element in Arbitrary Space Dimensions

 One of the most popular pairs of finite elements for solving mixed formulations of the Stokes and Navier–Stokes problem is the Q _k −P _k−1 ^disc element. Two possible versions of the discontinuous pressure space can be considered: one can either use an unmapped version of the P _k−1 ^disc space consisting of piecewise polynomial functions of degree at most k−1 on each cell or define a mapped version where the pressure space is defined as the image of a polynomial space on a reference cell. Since the reference transformation is in general not affine but multilinear, the two variants are not equal on arbitrary meshes. It is well-known, that the inf-sup condition is satisfied for the first variant. In the present paper we show that the latter approach satisfies the inf-sup condition as well for k≥2 in any space dimension. Received January 31, 2001; revised May 2, 2002 Published online: July 26, 2002     

Nitsche Type Mortaring for some Elliptic Problem with Corner Singularities

The paper deals with Nitsche type mortaring as a finite element method (FEM) for treating non-matching meshes of triangles at the interface of some domain decomposition. The approach is applied to the Poisson equation with Dirichlet boundary conditions (as a model problem) under the aspect that the interface passes re-entrant corners of the domain. For such problems and non-matching meshes with and without local refinement near the re-entrant corner, some properties of the finite element scheme and error estimates are proved. They show that appropriate mesh grading yields convergence rates as known for the classical FEM in presence of regular solutions. Finally, a numerical example illustrates the approach and the theoretical results. Received July 5, 2001; revised February 5, 2002 Published online April 25, 2002     

Robust Multigrid Methods for Nearly Incompressible Elasticity

We consider multigrid methods for problems in linear elasticity which are robust with respect to the Poisson ratio. Therefore, we consider mixed approximations involving the displacement vector and the pressure, where the pressure is approximated by discontinuous functions. Then, the pressure can be eliminated by static condensation. The method is based on a saddle point smoother which was introduced for the Stokes problem and which is transferred to the elasticity system. The performance and the robustness of the multigrid method are demonstrated on several examples with different discretizations in 2D and 3D. Furthermore, we compare the multigrid method for the saddle point formulation and for the condensed positive definite system. Received February 5, 1999; revised October 5, 1999     

Condition Numbers of Approximate Schur Complements in Two- and Three-Dimensional Discretizations on Hierarchically Ordered Grids

We investigate multilevel incomplete factorizations of M-matrices arising from finite difference discretizations. The nonzero patterns are based on special orderings of the grid points. Hence, the Schur complements that result from block elimination of unknowns refer to a sequence of hierarchical grids. Having reached the coarsest grid, Gaussian elimination yields a complete decomposition of the last Schur complement. The main focus of this paper is a generalization of the recursive five-point/nine-point factorization method (which can be applied in two-dimensional problems) to matrices that stem from discretizations on three-dimensional cartesian grids. Moreover, we present a local analysis that considers fundamental grid cells. Our analysis allows to derive sharp bounds for the condition number associated with one factorization level (two-grid estimates). A comparison in case of the Laplace operator with Dirichlet boundary conditions shows: Estimating the relative condition number of the multilevel preconditioner by multiplying corresponding two-grid values gives the asymptotic bound O(h ^−0.347) for the two- respectively O(h ^−4/5) for the three-dimensional model problem. Received October 19, 1998; revised September 27, 1999     

Data Mining with Sparse Grids

O (h _n ⁻¹ n ^{d −1}) instead of O(h _n ^−d ) grid points and unknowns are involved. Here d denotes the dimension of the feature space and h _n = 2 ⁻ⁿ gives the mesh size. To be precise, we suggest to use the sparse grid combination technique [42] where the classification problem is discretized and solved on a certain sequence of conventional grids with uniform mesh sizes in each coordinate direction. The sparse grid solution is then obtained from the solutions on these different grids by linear combination. In contrast to other sparse grid techniques, the combination method is simpler to use and can be parallelized in a natural and straightforward way. We describe the sparse grid combination technique for the classification problem in terms of the regularization network approach. We then give implementational details and discuss the complexity of the algorithm. It turns out that the method scales only linearly with the number of instances, i.e. the amount of data to be classified. Finally we report on the quality of the classifier built by our new method. Here we consider standard test problems from the UCI repository and problems with huge synthetical data sets in up to 9 dimensions. It turns out that our new method achieves correctness rates which are competitive to that of the best existing methods. Received April 25, 2001     

Fundamental Systems of Numerical Schemes for Linear Convection–Diffusion Equations and their Relationship to Accuracy

A new approach towards the assessment and derivation of numerical methods for convection dominated problems is presented, based on the comparison of the fundamental systems of the continuous and discrete operators. In two or more space dimensions, the dimension of the fundamental system is infinite, and may be identified with a ball. This set is referred to as the true fundamental locus. The fundamental system for a numerical scheme also forms a locus. As a first application, it is shown that a necessary condition for the uniform convergence of a numerical scheme is that the discrete locus should contain the true locus, and it is then shown it is impossible to satisfy this condition with a finite stencil. This shows that results of Shishkin concerning non-uniform convergence at parabolic boundaries are also generic for outflow boundaries. It is shown that the distance between the loci is related to the accuracy of the schemes provided that the loci are sufficiently close. However, if the loci depart markedly, then the situation is rather more complicated. Under suitable conditions, we develop an explicit numerical lower bound on the attainable relative error in terms of the coefficients in the stencil characterising the scheme and the loci. Received December 10, 1999; revised August 14, 2000     

Multilevel Gauging for Edge Elements

The vector potential of a solenoidal vector field, if it exists, is not unique in general. Any procedure that aims to determine such a vector potential typically involves a decision on how to fix it. This is referred to by the term gauging. Gauging is an important issue in computational electromagnetism, whenever discrete vector potentials have to be computed. In this paper a new gauging algorithm for discrete vector potentials is introduced that relies on a hierarchical multilevel decomposition. With minimum computational effort it yields vector potentials whose L ²-norm does not severely blow up. Thus the new approach compares favorably to the widely used co-tree gauging. Received May 27, 1999; revised October 22, 1999     

A Finite Element-Boundary Element Algorithm for Inhomogeneous Boundary Value Problems

Received January 25, 2001; revised July 17, 2001     

Variable Additive Preconditioning Procedures

Iterative methods with variable preconditioners of additive type are proposed. The scaling factors of each summand in the additive preconditioners are optimized within each iteration step. It is proved that the presented methods converge at least as fast as the Richardson's iterative method with the corresponding additive preconditioner with optimal scaling factors. In the presented numerical experiments the suggested methods need nearly the same number of iterations as the usual preconditioned conjugate gradient method with the corresponding additive preconditioner with numerically determined fixed optimal scaling factors. Received: June 10, 1998; revised October 16, 1998     

A Solution Method for a Certain Nonlinear Interface Problem in Unbounded Domains

In this paper, we consider a kind of nonlinear interface problem in unbounded domains. To solve this problem, we discuss a new coupling of finite element and boundary element by adding an auxiliary circle. We first derive the optimal error estimate of finite element approximation to the coupled FEM-BEM problem. Then we introduce a preconditioning steepest descent method for solving the discrete system by constructing a cheap domain decomposition preconditioner. Moreover, we give a complete analysis to the convergence speed of this iterative method. Received March 30, 2000; revised November 29, 2000     

A Class of Smoothers for Saddle Point Problems

In this paper smoothing properties are shown for a class of iterative methods for saddle point problems with smoothing rates of the order 1/m, where m is the number of smoothing steps. This generalizes recent results by Braess and Sarazin, who could prove this rates for methods where, in the context of the Stokes problem, the pressure correction equation is solved exactly, which is not needed here. Received December 4, 1998; revised April 14, 2000     

A Cascadic Multigrid Algorithm for Semilinear Indefinite Elliptic Problems

We propose a cascadic multigrid algorithm for a semilinear indefinite elliptic problem. We use a standard finite element discretization with piecewise linear finite elements. The arising nonlinear equations are solved by a cascadic organization of Newton's method with frozen derivative on a sequence of nested grids. This gives a simple version of a multigrid method without projections on coarser grids. The cascadic multigrid algorithm starts on a comparatively coarse grid where the number of unknowns is small enough to obtain an approximate solution within sufficiently high precision without substantial computational effort. On each finer grid we perform exactly one Newton step taking the approximate solution from the coarsest grid as initial guess. The linear Newton systems are solved iteratively by a Jacobi-type iteration with special parameters using the approximate solution from the previous grid as initial guess. We prove that for a sufficiently fine initial grid and for a sufficiently good start approximation the algorithm yields an approximate solution within the discretization error on the finest grid and that the method has multigrid complexity with logarithmic multiplier. Received February 1999, revised July 13, 1999     

Incomplete Cross Approximation in the Mosaic-Skeleton Method

The mosaic-skeleton method was bred in a simple observation that rather large blocks in very large matrices coming from integral formulations can be approximated accurately by a sum of just few rank-one matrices (skeletons). These blocks might correspond to a region where the kernel is smooth enough, and anyway it can be a region where the kernel is approximated by a short sum of separable functions (functional skeletons). Since the effect of approximations is like that of having small-rank matrices, we find it pertinent to say about mosaic ranks of a matrix which turn out to be pretty small for many nonsingular matrices. On the first stage, the method builds up an appropriate mosaic partitioning using the concept of a tree of clusters and some extra information rather than the matrix entries (related to the mesh). On the second stage, it approximates every allowed block by skeletons using the entries of some rather small cross which is chosen by an adaptive procedure. We focus chiefly on some aspects of practical implementation and numerical examples on which the approximation time was found to grow almost linearly in the matrix size. Received February 13, 1999; revised October 26, 1999     

A Sparse Matrix Arithmetic Based on Cal H-Matrices. Part I: Introduction to {Cal H}-Matrices

A class of matrices (-matrices) is introduced which have the following properties. (i) They are sparse in the sense that only few data are needed for their representation. (ii) The matrix-vector multiplication is of almost linear complexity. (iii) In general, sums and products of these matrices are no longer in the same set, but their truncations to the -matrix format are again of almost linear complexity. (iv) The same statement holds for the inverse of an -matrix. This paper is the first of a series and is devoted to the first introduction of the -matrix concept. Two concret formats are described. The first one is the simplest possible. Nevertheless, it allows the exact inversion of tridiagonal matrices. The second one is able to approximate discrete integral operators. Received: July 30, 1998; revised December 28, 1998     

Solution of Large Scale Algebraic Matrix Riccati Equations by Use of Hierarchical Matrices

In previous papers, a class of hierarchical matrices (ℋ-matrices) is introduced which are data-sparse and allow an approximate matrix arithmetic of almost optimal complexity. Here, we investigate a new approach to exploit the ℋ-matrix structure for the solution of large scale Lyapunov and Riccati equations as they typically arise for optimal control problems where the constraint is a partial differential equation of elliptic type. This approach leads to an algorithm of linear-logarithmic complexity in the size of the matrices. Received July 30, 2002; revised December 16, 2002 Published online: April 22, 2003