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
- scientific article; zbMATH DE number 3874483 (Why is no real title available?)
- scientific article; zbMATH DE number 2065097 (Why is no real title available?)
- A Flexible Inner-Outer Preconditioned GMRES Algorithm
- A conservative discontinuous Galerkin semi-implicit formulation for the Navier-Stokes equations in nonhydrostatic mesoscale modeling
- A discontinuous Galerkin method with weighted averages for advection-diffusion equations with locally small and anisotropic diffusivity
- A discontinuous \(hp\) finite element method for convection-diffusion problems
- A discontinuous \(hp\) finite element method for diffusion problems
- 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
- Improved energy estimates for interior penalty, constrained and discontinuous Galerkin methods for elliptic problems. I
- Methods of conjugate gradients for solving linear systems
- Numerical Solution of the Navier-Stokes Equations
- Numerical recipes. The art of scientific computing.
- 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
(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
- Multigrid preconditioners for the hybridised discontinuous Galerkin discretisation of the shallow water equations
- Low-Order Preconditioning for the High-Order Finite Element de Rham Complex
- 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 Robust Multigrid Solver for Isogeometric Analysis Based on Multiplicative Schwarz Smoothers
- A multigrid preconditioner for Jacobian-free Newton-Krylov methods
- 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
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)
