Abstract:In this paper, we examine a number of additive and multiplicative multilevel iterative methods and preconditioners in the setting of two-dimensional local mesh refinement. While standard multilevel methods are effective for uniform refinement-based discretizations of elliptic equations, they tend to be less effective for algebraic systems which arise from discretizations on locally refined meshes, losing their optimal behavior in both storage and computational complexity. Our primary focus here is on BPX-style additive and multiplicative multilevel preconditioners, and on various stabilizations of the additive and multiplicative hierarchical basis method (HB), and their use in the local mesh refinement setting. In this article, we describe in detail the implementation of these types of algorithms, including detailed discussions of the datastructures and traversal algorithms we employ for obtaining optimal storage and computational complexity in our implementations. We show how each of the algorithms can be implemented using standard datatypes available in languages such as C and FORTRAN, so that the resulting algorithms have optimal (linear) storage requirements, and so that the resulting multilevel method or preconditioner can be applied with optimal (linear) computational costs. Our implementations are performed in both C and MATLAB using the Finite Element ToolKit (FETK), an open source finite element software package. We finish the paper with a sequence of numerical experiments illustrating the effectiveness of a number of BPX and stabilized HB variants for several examples requiring local refinement.

Substructuring Preconditioners for an h-p Nitsche-type method

Optimal preconditioners for Nitsche-XFEM discretizations of interface problems

Substructuring preconditioners with novel interface solvers for general elliptic-type equations in three dimensions

A substructuring preconditioner with vertex-related interface solvers for elliptic-type equations in three dimensions

Boundary Element Methods for the Laplace Hypersingular Integral Equation on Multiscreens: a two-level Substructuring Preconditioner

Substructuring the Hiptmair-Xu preconditioner for positive Maxwell problems

Preconditioners for Elliptic Problems Via Non-Uniform Meshes.

Efficient additive Schwarz preconditioning for hypersingular integral equations on locally refined triangulations

Subspace and auxiliary space preconditioners for high-order interior penalty discretizations in $H(\mathrm{div})$

Nonoverlapping Domain Decomposition Methods with A Simple Coarse Space for Elliptic Problems

Substructuring Preconditioners with a Simple Coarse Space for 2-D 3-T Radiation Diffusion Equations

A PRECONDITIONER BASED ON NONOVERLAPPING DOMAIN DECOMPOSITION FOR SOLVING H(curl) ELLIPTIC PROBLEM

Double Saddle-Point Preconditioning for Krylov Methods in the Inexact Sequential Homotopy Method

Local Refinement and Multilevel Preconditioning: Implementation and Numerical Experiments

Efficient Parallel Preconditioners for High-Order Finite Element Discretizations of H(grad) and H(curl) Problems

An AMG saddle point preconditioner with application to mixed Poisson problems on adaptive quad/cube meshes

A One-level Additive Schwarz Preconditioner for a Discontinuous Petrov-Galerkin Method

Cell-by-cell approximate Schur complement technique in preconditioning of meshfree discretized piezoelectric equations

A deflated Schur complement method for the iterative solution of a high-order discontinuous element discretization of the Poisson equation

Nodal auxiliary space preconditioners for mixed virtual element methods

Nitsche stabilized Virtual element approximations for a Brinkman problem with mixed boundary conditions