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.

Analysis of a New Space-Time Parallel Multigrid Algorithm for Parabolic Problems

Parallel Conjugate Gradient Method Based on Spline Difference Method for the One-Dimensional Heat Equation

An optimized Space-Time Multigrid algorithm for parabolic PDEs

A parallel space-time multigrid method for the eddy-current equation

A Space-Time Multigrid Method for Space-Time Finite Element Discretizations of Parabolic and Hyperbolic PDEs

Efficient Multigrid Reduction-in-Time for Method-of-Lines Discretizations of Linear Advection

The research of discontinuous galerkin p-multigrid solver

Multigrid reduction-in-time convergence for advection problems: A Fourier analysis perspective

Parallel Smoothers for Matrix-Based Geometric Multigrid Methods on Locally Refined Meshes Using Multicore CPUs and GPUs

Parallel-in-time Multilevel Krylov Methods: A Prototype

Matrix-free multigrid block-preconditioners for higher order Discontinuous Galerkin discretisations

Multigrid Reduction‐In‐Time Convergence for Advection Problems: A Fourier Analysis Perspective

A new ParaDiag time-parallel time integration method

Parallel discontinuous Galerkin finite element method for computing hyperbolic conservation law on unstructured meshes

A Fully Parallelizable Space-Time Multilevel Monte Carlo Method for Stochastic Differential Equations with Additive Noise

Analysis of the Parareal Approach Based on Discontinuous Galerkin Method for Time‐dependent Stokes Equations

Efficient Parallel Algorithms for Parabolic Problems.

Parallel-in-time solution of hyperbolic PDE systems via characteristic-variable block preconditioning

An energy-efficient GMRES–multigrid solver for space-time finite element computation of dynamic poroelasticity

A Two-Grid Discontinuous Galerkin Method for a Kind of Nonlinear Parabolic Problems

Analysis of Parallel Multigrid Methods in Real-Time Fluid Simulation