A performance comparison of continuous and discontinuous Galerkin methods with fast multigrid solvers

Abstract

This study presents a fair performance comparison of the continuous finite element method, the symmetric interior penalty discontinuous Galerkin method, and the hybridized discontinuous Galerkin method. Modern implementations of high-order methods with state-of-the-art multigrid solvers for the Poisson equation are considered, including fast matrix-free implementations with sum factorization on quadrilateral and hexahedral elements. For the hybridized discontinuous Galerkin method, a multigrid approach that combines a grid transfer from the trace space to the space of linear finite elements with algebraic multigrid on further levels is developed. Despite similar solver complexity of the matrix-based HDG solver and matrix-free geometric multigrid schemes with continuous and discontinuous Galerkin finite elements, the latter offer up to order of magnitude faster time to solution, even after including the superconvergence effects. This difference is because of vastly better performance of matrix-free operator evaluation as compared to sparse matrix-vector products. A roofline performance model confirms the advantage of the matrix-free implementation.

Cited by in corpus (27)

• The deal.II finite element library: design, features, and insights
• Fast matrix-free evaluation of discontinuous Galerkin finite element operators
• preCICE v2: A Sustainable and User-Friendly Coupling Library
• A Flexible, Parallel, Adaptive Geometric Multigrid method for FEM
• Hybrid multigrid methods for high-order discontinuous Galerkin discretizations
• HDGlab: An open-source implementation of the hybridisable discontinuous Galerkin method in MATLAB
• A Multilevel Approach for Trace System in HDG Discretizations
• On the implementation of a robust and efficient finite element-based parallel solver for the compressible Navier-Stokes equations
• Efficient Explicit Time Stepping of High Order Discontinuous Galerkin Schemes for Waves
• High-order arbitrary Lagrangian-Eulerian discontinuous Galerkin methods for the incompressible Navier-Stokes equations
• Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
• Slate: extending Firedrake's domain-specific abstraction to hybridized solvers for geoscience and beyond
• A scalable and robust vertex-star relaxation for high-order FEM
• A multiscale approach to hybrid RANS/LES wall modeling within a high-order discontinuous Galerkin scheme using function enrichment
• Preconditioning for a pressure-robust HDG discretization of the Stokes equations
• A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
• End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations
• Local Fourier analysis of multigrid for hybridized and embedded discontinuous Galerkin methods
• Matrix-free GPU-accelerated saddle-point solvers for high-order problems in $H(\mathrm{div})$
• Matrix-Free Higher-Order Finite Element Methods for Hyperelasticity
• Geometric Uncertainty of Patient-Specific Blood Vessels and its Impact on Aortic Hemodynamics
• A hybridizable discontinuous Galerkin method for electromagnetics with a view on subsurface applications
• Homogeneous multigrid for embedded discontinuous Galerkin methods
• A discontinuous Galerkin discretization of elliptic problems with improved convergence properties using summation by parts operators
• An asynchronous discontinuous Galerkin method for massively parallel PDE solvers
• SynthAorta: A 3D Mesh Dataset of Parametrized Physiological Healthy Aortas
• Higher-Order Discontinuous Galerkin Splitting Schemes for Fluids with Variable Viscosity