Abstract:The International Journal of High Performance Computing Applications, Ahead of Print. In recent years, solvers for finite-element discretizations of linear or linearized saddle-point problems, like the Stokes and Oseen equations, have become well established. There are two main classes of preconditioners for such systems: those based on a block-factorization approach and those based on monolithic multigrid. Both classes of preconditioners have several critical choices to be made in their composition, such as the selection of a suitable relaxation scheme for monolithic multigrid. From existing studies, some insight can be gained as to what options are preferable in low-performance computing settings, but there are very few fair comparisons of these approaches in the literature, particularly for modern architectures, such as GPUs. In this paper, we perform a comparison between a Block-Triangular preconditioner and monolithic multigrid methods with the three most common choices of relaxation scheme – Braess-Sarazin, Vanka, and Schur-Uzawa. We develop a performant Vanka relaxation algorithm for structured-grid discretizations, which takes advantage of memory efficiencies in this setting. We detail the behavior of the various CUDA kernels for the multigrid relaxation schemes and evaluate their individual arithmetic intensity, performance, and runtime. Running a preconditioned FGMRES solver for the Stokes equations with these preconditioners allows us to compare their efficiency in a practical setting. We show that monolithic multigrid can outperform Block-Triangular preconditioning, and that using Vanka or Braess-Sarazin relaxation is most efficient. Even though multigrid with Vanka relaxation exhibits reduced performance on the CPU (up to 100% slower than Braess-Sarazin), it is able to outperform Braess-Sarazin by more than 20% on the GPU, making it a competitive algorithm, especially given the high amount of algorithmic tuning needed for effective Braess-Sarazin relaxation.

Local Refinement and Multilevel Preconditioning: Implementation and Numerical Experiments

A Local Grid Refinement Scheme for B‐spline Material Point Method

Efficient additive Schwarz preconditioning for hypersingular integral equations on locally refined triangulations

Multilevel Preconditioner with Stable Coarse Grid Corrections for the Helmholtz Equation.

Novel Multilevel Preconditioners for the Systems Arising from Plane Wave Discretization of Helmholtz Equations with Large Wave Numbers

Goal-Oriented Adaptivity and Multilevel Preconditioning for the Poisson-Boltzmann Equation

An anisotropic multilevel preconditioner for solving the Helmholtz equation with unequal directional sampling intervals

Multilevel Tau preconditioners for symmetrized multilevel Toeplitz systems with applications to solving space fractional diffusion equations

Optimal preconditioners for nonsymmetric multilevel Toeplitz systems with application to solving non-local evolutionary partial differential equations

Hybrid Localized Spectral Decomposition for multiscale problems

A Non-Nested Multilevel Method for Meshless Solution of the Poisson Equation in Heat Transfer and Fluid Flow

Two-level hybrid Schwarz Preconditioners for The Helmholtz Equation with high wave number

Preconditioning Low Rank Generalized Minimal Residual Method (GMRES) for Implicit Discretizations of Matrix Differential Equations

On the Preconditioners for Fast Multipole Boundary Element Methods for 2D Multi-Domain Elastostatics

Scalable multigrid methods for immersed finite element methods and immersed isogeometric analysis

Monolithic Algebraic Multigrid Preconditioners for the Stokes Equations

Exploiting mesh structure to improve multigrid performance for saddle-point problems

CBS constants for multilevel splitting of graph-Laplacian and application to preconditioning of discontinuous Galerkin systems

Robust Preconditioners for DG-Discretizations with Arbitrary Polynomial Degrees

High level implementation of geometric multigrid solvers for finite element problems: Applications in atmospheric modelling

Coupling parallel adaptive mesh refinement with a nonoverlapping domain decomposition solver