Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations
From MaRDI portal
Publication:2222294
DOI10.1016/J.JCP.2019.06.001zbMATH Open1452.65322arXiv1805.11930OpenAlexW2953067464WikidataQ127756098 ScholiaQ127756098MaRDI QIDQ2222294FDOQ2222294
Authors: P. Bastian, Eike H. Müller, Steffen Müthing, Marian Piatkowski 
Publication date: 26 January 2021
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: Efficient and suitably preconditioned iterative solvers for elliptic partial differential equations (PDEs) of the convection-diffusion type are used in all fields of science and engineering. To achieve optimal performance, solvers have to exhibit high arithmetic intensity and need to exploit every form of parallelism available in modern manycore CPUs. The computationally most expensive components of the solver are the repeated applications of the linear operator and the preconditioner. For discretisations based on higher-order Discontinuous Galerkin methods, sum-factorisation results in a dramatic reduction of the computational complexity of the operator application while, at the same time, the matrix-free implementation can run at a significant fraction of the theoretical peak floating point performance. Multigrid methods for high order methods often rely on block-smoothers to reduce high-frequency error components within one grid cell. Traditionally, this requires the assembly and expensive dense matrix solve in each grid cell, which counteracts any improvements achieved in the fast matrix-free operator application. To overcome this issue, we present a new matrix-free implementation of block-smoothers. Inverting the block matrices iteratively avoids storage and factorisation of the matrix and makes it is possible to harness the full power of the CPU. We implemented a hybrid multigrid algorithm with matrix-free block-smoothers in the high order DG space combined with a low order coarse grid correction using algebraic multigrid where only low order components are explicitly assembled. The effectiveness of this approach is demonstrated by solving a set of representative elliptic PDEs of increasing complexity, including a convection dominated problem and the stationary SPE10 benchmark.
Full work available at URL: https://arxiv.org/abs/1805.11930
Recommendations
- Efficient low-order refined preconditioners for high-order matrix-free continuous and discontinuous Galerkin methods
- Smoothed aggregation multigrid solvers for high-order discontinuous Galerkin methods for elliptic problems
- Multigrid algorithms for high order discontinuous Galerkin methods
- Hybrid multigrid methods for high-order discontinuous Galerkin discretizations
- Fast tensor product Schwarz smoothers for high-order discontinuous Galerkin methods
Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55) Boundary value problems for second-order elliptic systems (35J57) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)
Cites Work
- The \texttt{deal.II} library, version 8.4
- A generic grid interface for parallel and adaptive scientific computing. I: Abstract framework
- A generic grid interface for parallel and adaptive scientific computing. II: Implementation and tests in DUNE
- GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
- Methods of conjugate gradients for solving linear systems
- Parallel iterative methods for sparse linear systems
- Numerical recipes. The art of scientific computing.
- Estimation of penalty parameters for symmetric interior penalty Galerkin methods
- Sur l'approximation de la solution des équations de Navier-Stokes par la méthode des pas fractionnaires. II
- An Interior Penalty Finite Element Method with Discontinuous Elements
- An Elliptic Collocation-Finite Element Method with Interior Penalties
- Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
- Numerical Solution of the Navier-Stokes Equations
- Spectral/hp Element Methods for Computational Fluid Dynamics
- A discontinuous \(hp\) finite element method for convection-diffusion problems
- Approximation of the velocity by coupling discontinuous Galerkin and mixed finite element methods for flow problems
- A discontinuous Galerkin method with weighted averages for advection-diffusion equations with locally small and anisotropic diffusivity
- A review of algebraic multigrid
- An overlapping Schwarz method for spectral element solution of the incompressible Navier-Stokes equations
- Title not available (Why is that?)
- A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure
- A discontinuous \(hp\) finite element method for diffusion problems
- A Flexible Inner-Outer Preconditioned GMRES Algorithm
- From \(h\) to \(p\) efficiently: implementing finite and spectral/hp element methods to achieve optimal performance for low- and high-order discretisations
- Improved energy estimates for interior penalty, constrained and discontinuous Galerkin methods for elliptic problems. I
- A generic interface for parallel cell-based finite element operator application
- A conservative discontinuous Galerkin semi-implicit formulation for the Navier-Stokes equations in nonhydrostatic mesoscale modeling
- An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier-Stokes equations
- Spectral methods for problems in complex geometries
- Title not available (Why is that?)
- Algebraic multigrid for discontinuous Galerkin discretizations of heterogeneous elliptic problems
- Error estimation in preconditioned conjugate gradients
- Fast simplicial finite element algorithms using Bernstein polynomials
- Fast simplicial quadrature-based finite element operators using Bernstein polynomials
- Semi-automatic sparse preconditioners for high-order finite element methods on non-uniform meshes
- A parallel algebraic multigrid method for elliptic problems with highly discontinuous coefficients
- Efficient nonlinear solvers for nodal high-order finite elements in 3D
- Algebraic multigrid techniques for discontinuous Galerkin methods with varying polynomial order
- A stable and high-order accurate discontinuous Galerkin based splitting method for the incompressible Navier-Stokes equations
Cited In (18)
- Nektar++: design and implementation of an implicit, spectral/\(hp\) element, compressible flow solver using a Jacobian-free Newton Krylov approach
- Efficient low-order refined preconditioners for high-order matrix-free continuous and discontinuous Galerkin methods
- High-performance implementation of discontinuous Galerkin methods with application in fluid flow
- Fast tensor product Schwarz smoothers for high-order discontinuous Galerkin methods
- The \textsc{Dune} framework: basic concepts and recent developments
- High level implementation of geometric multigrid solvers for finite element problems: applications in atmospheric modelling
- A matrix-free high-order solver for the numerical solution of cardiac electrophysiology
- Hybrid multigrid methods for high-order discontinuous Galerkin discretizations
- Low-Order Preconditioning for the High-Order Finite Element de Rham Complex
- Multigrid preconditioners for the hybridised discontinuous Galerkin discretisation of the shallow water equations
- Advantages of static condensation in implicit compressible Navier-Stokes DGSEM solvers
- Scaling to the stars -- a linearly scaling elliptic solver for \(p\)-multigrid
- High-order matrix-free incompressible flow solvers with GPU acceleration and low-order refined preconditioners
- A multigrid preconditioner for Jacobian-free Newton-Krylov methods
- A Robust Multigrid Solver for Isogeometric Analysis Based on Multiplicative Schwarz Smoothers
- Textbook efficiency: massively parallel matrix-free multigrid for the Stokes system
- Manycore parallel computing for a hybridizable discontinuous Galerkin nested multigrid method
- Efficient multigrid solution of elliptic interface problems using viscosity-upwinded local discontinuous Galerkin methods
Uses Software
This page was built for publication: Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222294)
