Smoothed prolongation multigrid with rapid coarsening and massive smoothing

We prove that within the frame of smoothed prolongations, rapid coarsening between first two levels can be compensated by massive prolongation smoothing and pre- and post-smoothing derived from the prolongator smoother.     

Uniform convergence of multigrid V-cycle on adaptively refined finite element meshes for second order elliptic problems

In this paper we prove the uniform convergence of the standard multigrid V-cycle algorithm with the Gauss-Seidel relaxation performed only on the new nodes and their "immediate" neighbors for discrete elliptic problems on the adaptively refined finite element meshes using the newest vertex bisection algorithm. The proof depends on sharp estimates on the relationship of local mesh sizes and a new stability estimate for the space decomposition based on the Scott-Zhang interpolation operator. Extensive numerical results are reported, which confirm the theoretical analysis.     

A note on convergence of the multigrid V-cycle

Several recent papers have discussed the convergence of the multigrid V-cycle. In particular there are several results for the symmetric case: where the numbers of smoothings before the fine-to-coarse transfer and after the coarse-to-fine transfer are the same. In most instances, the smoother H=I?E^-1A has been limited to the case where E is positive definite and the eigenvalues h of H satisfy 0?h?1. In this note we extend these results to asymmetric V-cycles and the case where -b?h?1 with 0<b<1.     

Notes on convergence of an algebraic multigrid method

The convergence theory for algebraic multigrid (AMG) algorithms proposed in Chang and Huang [Q.S. Chang, Z.H. Huang, Efficient algebraic multigrid algorithms and their convergence, SIAM J. Sci. Comput. 24 (2002) 597–618] is further discussed and a smaller and elegant upper bound is obtained. On the basis of element-free AMGe [V.E. Henson, P.S. Vassilevski, Element-free AMGe: General algorithms for computing interpolation weights in AMG, SIAM J. Sci. Comput. 23(2) (2001) 629–650] we rewrite the interpolation operator for the classical AMG (cAMG), present a uniform expression and then, by introducing a sparse approximate inverse in the Frobenius norm, give a general convergence theorem which is suited for not only cAMG but also AMG for finite elements and element-free AMGe.     

V-cycle convergence of some multigrid methods for ill-posed problems

For ill-posed linear operator equations we consider some V-cycle multigrid approaches, that, in the framework of Bramble, Pasciak, Wang, and Xu (1991), we prove to yield level independent contraction factor estimates. Consequently, we can incorporate these multigrid operators in a full multigrid method, that, together with a discrepancy principle, is shown to act as an iterative regularization method for the underlying infinite-dimensional ill-posed problem. Numerical experiments illustrate the theoretical results.

     

Algebraic multigrid,mixed‐order interpolation,and incompressible fluid flow

This paper presents the results of numerical experiments on the use of equal‐order and mixed‐order interpolations in algebraic multigrid (AMG) solvers for the fully coupled equations of incompressible fluid flow. Several standard test problems are addressed for Reynolds numbers spanning the laminar range. The range of unstructured meshes spans over two orders of problem size (over one order of mesh bandwidth). Deficiencies in performance are identified for AMG based on equal‐order interpolations (both zero‐order and first‐order). They take the form of poor, fragile, mesh‐dependent convergence rates. The evidence suggests that a degraded representation of the inter‐field coupling in the coarse‐grid approximation is the cause. Mixed‐order interpolation (first‐order for the vectors, zero‐order for the scalars) is shown to address these deficiencies. Convergence is then robust, independent of the number of coarse grids and (almost) of the mesh bandwidth. The AMG algorithms used are reviewed. Copyright © 2009 John Wiley & Sons, Ltd.     

An improved convergence analysis of smoothed aggregation algebraic multigrid

We present an improved analysis of the smoothed aggregation algebraic multigrid method extending the original proof in [Numer. Math. 2001; 88 :559–579] and its modification in [Multilevel Block Factorization Preconditioners. Matrix‐based Analysis and Algorithms for Solving Finite Element Equations. Springer: New York, 2008]. The new result imposes fewer restrictions on the aggregates that makes it easier to verify in practice. Also, we extend a result in [Appl. Math. 2011] that allows us to use aggressive coarsening at all levels. This is due to the properties of the special polynomial smoother that we use and analyze. In particular, we obtain bounds in the multilevel convergence estimates that are independent of the coarsening ratio. Numerical illustration is also provided. Copyright © 2011 John Wiley & Sons, Ltd.     

Uniform convergence of the multigrid V-cycle for an anisotropic problem

In this paper, we consider the linear systems arising from the standard finite element discretizations of certain second order anisotropic problems with variable coefficients on a rectangle. We study the performance of a V-cycle multigrid method applied to the finite element equations. Since the usual ``regularity and approximation' assumption does not hold for the anisotropic finite element problems, the standard multigrid convergence theory cannot be applied directly. In this paper, a modification of the theory of Braess and Hackbusch will be presented. We show that the V-cycle multigrid iteration with a line smoother is a uniform contraction in the energy norm. In the verification of the hypotheses in our theory, we use a weighted -norm estimate for the error in the Galerkin finite element approximation and a smoothing property of the line smoothers which is proved in this paper.

     

First order convergence of matroids

The model theory based notion of the first order convergence unifies the notions of the left-convergence for dense structures and the Benjamini–Schramm convergence for sparse structures. It is known that every first order convergent sequence of graphs with bounded tree-depth can be represented by an analytic limit object called a limit modeling. We establish the matroid counterpart of this result: every first order convergent sequence of matroids with bounded branch-depth representable over a fixed finite field has a limit modeling, i.e., there exists an infinite matroid with the elements forming a probability space that has asymptotically the same first order properties. We show that neither of the bounded branch-depth assumption nor the representability assumption can be removed.     

The convergence of V-cycle multigrid algorithms for axisymmetric Laplace and Maxwell equations

We investigate some simple finite element discretizations for the axisymmetric Laplace equation and the azimuthal component of the axisymmetric Maxwell equations as well as multigrid algorithms for these discretizations. Our analysis is targeted at simple model problems and our main result is that the standard V-cycle with point smoothing converges at a rate independent of the number of unknowns. This is contrary to suggestions in the existing literature that line relaxations and semicoarsening are needed in multigrid algorithms to overcome difficulties caused by the singularities in the axisymmetric Maxwell problems. Our multigrid analysis proceeds by applying the well known regularity based multigrid theory. In order to apply this theory, we prove regularity results for the axisymmetric Laplace and Maxwell equations in certain weighted Sobolev spaces. These, together with some new finite element error estimates in certain weighted Sobolev norms, are the main ingredients of our analysis.

     

Analysis of the convergence of the multigrid method for rigid problems

Translated from Matematicheskie Zametki, Vol. 48, No. 2, pp. 53–63, August, 1990.     

Nearly optimal convergence result for multigrid with aggressive coarsening and polynomial smoothing

We analyze a general multigrid method with aggressive coarsening and polynomial smoothing. We use a special polynomial smoother that originates in the context of the smoothed aggregation method. Assuming the degree of the smoothing polynomial is, on each level k, at least Ch _k+1/h _k , we prove a convergence result independent of h _k+1/h _k . The suggested smoother is cheaper than the overlapping Schwarz method that allows to prove the same result. Moreover, unlike in the case of the overlapping Schwarz method, analysis of our smoother is completely algebraic and independent of geometry of the problem and prolongators (the geometry of coarse spaces).     

An optimal order multigrid method for biharmonic,C 1 finite element equations

Summary In this paper, we study a special multigrid method for solving large linear systems which arise from discretizing biharmonic problems by the Hsieh-Clough-Tocher,C ¹ macro finite elements or several otherC ¹ finite elements. Since the multipleC ¹ finite element spaces considered are not nested, the nodal interpolation operator is used to transfer functions between consecutive levels in the multigrid method. This method converges with the optimal computational order.     

Remark on convergence of algebraic multigrid in the form of matrix decomposition

We introduce the convergence of algebraic multigrid in the form of matrix decomposition. The convergence is proved in block versions of the multi-elimination incomplete LU (BILUM) factorization technique and the approximation of their inverses to preserve sparsity. The convergence theorem can be applied to general interpolation operator. Furthermore, we discuss the error caused by the error matrix.     

Duality for unbounded order convergence and applications

Unbounded order convergence has lately been systematically studied as a generalization of almost everywhere convergence to the abstract setting of vector and Banach lattices. This paper presents a duality theory for unbounded order convergence. We define the unbounded order dual (or uo-dual) \({X_{uo}^\sim }\) of a Banach lattice X and identify it as the order continuous part of the order continuous dual \({X_n^\sim }\). The result allows us to characterize the Banach lattices that have order continuous preduals and to show that an order continuous predual is unique when it exists. Applications to the Fenchel–Moreau duality theory of convex functionals are given. The applications are of interest in the theory of risk measures in Mathematical Finance.     

Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004