Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations
From MaRDI portal
Publication:2222294
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.
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
Cites work
- A conservative discontinuous Galerkin semi-implicit formulation for the Navier-Stokes equations in nonhydrostatic mesoscale modeling
- A discontinuous \(hp\) finite element method for convection-diffusion problems
- A discontinuous \(hp\) finite element method for diffusion problems
- A discontinuous Galerkin method with weighted averages for advection-diffusion equations with locally small and anisotropic diffusivity
- A Flexible Inner-Outer Preconditioned GMRES Algorithm
- A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure
- 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
- A generic interface for parallel cell-based finite element operator application
- A parallel algebraic multigrid method for elliptic problems with highly discontinuous coefficients
- A review of algebraic multigrid
- A stable and high-order accurate discontinuous Galerkin based splitting method for the incompressible Navier-Stokes equations
- Algebraic multigrid for discontinuous Galerkin discretizations of heterogeneous elliptic problems
- Algebraic multigrid techniques for discontinuous Galerkin methods with varying polynomial order
- An Elliptic Collocation-Finite Element Method with Interior Penalties
- An Interior Penalty Finite Element Method with Discontinuous Elements
- An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier-Stokes equations
- An overlapping Schwarz method for spectral element solution of the incompressible Navier-Stokes equations
- Approximation of the velocity by coupling discontinuous Galerkin and mixed finite element methods for flow problems
- Efficient nonlinear solvers for nodal high-order finite elements in 3D
- Error estimation in preconditioned conjugate gradients
- Estimation of penalty parameters for symmetric interior penalty Galerkin methods
- Fast simplicial finite element algorithms using Bernstein polynomials
- Fast simplicial quadrature-based finite element operators using Bernstein polynomials
- From \(h\) to \(p\) efficiently: implementing finite and spectral/hp element methods to achieve optimal performance for low- and high-order discretisations
- GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
- scientific article; zbMATH DE number 3874483 (Why is no real title available?)
- scientific article; zbMATH DE number 2065097 (Why is no real title available?)
- Improved energy estimates for interior penalty, constrained and discontinuous Galerkin methods for elliptic problems. I
- Methods of conjugate gradients for solving linear systems
- Numerical recipes. The art of scientific computing.
- Numerical Solution of the Navier-Stokes Equations
- Parallel iterative methods for sparse linear systems
- Semi-automatic sparse preconditioners for high-order finite element methods on non-uniform meshes
- Spectral methods for problems in complex geometries
- Spectral/hp Element Methods for Computational Fluid Dynamics
- Sur l'approximation de la solution des équations de Navier-Stokes par la méthode des pas fractionnaires. II
- The \texttt{deal.II} library, version 8.4
- Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
Cited in
(24)- Advantages of static condensation in implicit compressible Navier-Stokes DGSEM solvers
- Multigrid preconditioners for the hybridised discontinuous Galerkin discretisation of the shallow water equations
- High-order matrix-free incompressible flow solvers with GPU acceleration and low-order refined preconditioners
- Hybrid multigrid methods for high-order discontinuous Galerkin discretizations
- The \textsc{Dune} framework: basic concepts and recent developments
- Nektar++: design and implementation of an implicit, spectral/\(hp\) element, compressible flow solver using a Jacobian-free Newton Krylov approach
- Scaling to the stars -- a linearly scaling elliptic solver for \(p\)-multigrid
- High-performance implementation of discontinuous Galerkin methods with application in fluid flow
- Fast tensor product Schwarz smoothers for high-order discontinuous Galerkin methods
- Efficient multigrid solution of elliptic interface problems using viscosity-upwinded local discontinuous Galerkin methods
- 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
- Manycore parallel computing for a hybridizable discontinuous Galerkin nested multigrid method
- Textbook efficiency: massively parallel matrix-free multigrid for the Stokes system
- Efficient low-order refined preconditioners for high-order matrix-free continuous and discontinuous Galerkin methods
- A Robust Multigrid Solver for Isogeometric Analysis Based on Multiplicative Schwarz Smoothers
- Low-Order Preconditioning for the High-Order Finite Element de Rham Complex
- A multigrid preconditioner for Jacobian-free Newton-Krylov methods
- A Jacobian-free multigrid preconditioner for discontinuous Galerkin methods applied to atmospheric flows
- A matrix-free stabilized solver for the incompressible Navier-Stokes equations
- Subspace and auxiliary space preconditioners for high-order interior penalty discretizations in \(H(\mathrm{div})\)
- An implementation of tensor product patch smoothers on GPUs
- Unified interface flux evaluation in a general discontinuous Galerkin spectral element framework
- Adaptive Multilevel Discontinuous Galerkin Methods on GPUs
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)
