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.

Block-circulant with circulant-block preconditioners for two-dimensional spatial fractional diffusion equations

R. Chan's circulant-based approximate inverse preconditioning iterative method for solving second-order space fractional advection–dispersion equations with variable coefficients

Circulant-based approximate inverse preconditioners for a class of fractional diffusion equations

An efficient multigrid method with preconditioned smoother for two-dimensional anisotropic space-fractional diffusion equations

Strang-type Preconditioners for Solving Fractional Diffusion Equations by Boundary Value Methods

Diagonal and normal with Toeplitz-block splitting iteration method for space fractional coupled nonlinear Schrödinger equations with repulsive nonlinearities

A -Preconditioner for Space Fractional Diffusion Equation with Non-separable Variable Coefficients

Spectral analysis for preconditioning of multi-dimensional Riesz fractional diffusion equations

A novel fourth-order scheme for two-dimensional Riesz space fractional nonlinear reaction-diffusion equations and its optimal preconditioned solver

A block $α$-circulant based preconditioned MINRES method for wave equations

Mixed Precision Block-Jacobi Preconditioner: Algorithms, Performance Evaluation and Feature Analysis

Block [math]-Circulant Preconditioners for Parabolic Optimal Control Problems

Practical Restrictively Preconditioned Conjugate Gradient Methods for a Class of Block Two-by-Two Linear Systems

On $τ$-preconditioners for a quasi-compact difference scheme to Riesz fractional diffusion equations with variable coefficients

A reduced conjugate gradient basis method for fractional diffusion

An IRK‐QCD Scheme for the Space Fractional Diffusion Equation and Block‐Lower‐Triangle Preconditioners for the Corresponding Linear Systems

On Τ Matrix-Based Approximate Inverse Preconditioning Technique for Diagonal-Plus-toeplitz Linear Systems from Spatial Fractional Diffusion Equations

Uniform block-diagonal preconditioners for divergence-conforming HDG Methods for the generalized Stokes equations and the linear elasticity equations

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

Effective Block Preconditioners for Fluid Dynamics Coupled to Reduced Models of a Non-Local Nature

A τ Matrix Based Approximate Inverse Preconditioning for Tempered Fractional Diffusion Equations