Efficient low-order refined preconditioners for high-order matrix-free continuous and discontinuous Galerkin methods

Abstract

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on FEM-SEM equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order ($p=1$) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in $h$ and $p$. The method is extended to interior penalty and BR2 discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

arXiv comment26 pages, 9 figures

References in corpus (2)

• Nodal Discontinuous Galerkin Methods on Graphics Processors
• Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods

Cited by in corpus (9)

• Efficient Exascale Discretizations: High-Order Finite Element Methods
• Low-order preconditioning for the high-order finite element de Rham complex
• High-performance finite elements with MFEM
• A Family of Independent Variable Eddington Factor Methods with Efficient Preconditioned Iterative Solvers
• A scalable and robust vertex-star relaxation for high-order FEM
• End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations
• Matrix-free GPU-accelerated saddle-point solvers for high-order problems in $H(\mathrm{div})$
• Conservative and accurate solution transfer between high-order and low-order refined finite element spaces
• An implementation of tensor product patch smoothers on GPU