Abstract: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 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 a monolithic multigrid method 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 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.

Analysis of a multigrid preconditioner for Crouzeix-Raviart discretization of elliptic PDE with jump coefficient

Convergence Analysis on Multigrid Methods for Elliptic Problems with Large Jumps in Coefficients

Multigrid method and multilevel additive preconditioner for mixed element method for non-self-adjoint and indefinite problems

Uniform convergence of multigrid V-cycle on adaptively refined finite element meshes for second order elliptic problems

Uniform Convergence of a Multigrid Method for Elliptic Equations with Anisotropic Coefficients

Multigrid preconditioning for discontinuous Galerkin discretizations of an elliptic optimal control problem with a convection-dominated state equation

The effect of approximate coarsest-level solves on the convergence of multigrid V-cycle methods

Rounding-Error Analysis of Multigrid V-Cycles

Two Schemes to Improve the Implementation of the Aggregation Based Algebraic Multigrid Preconditioner

Uniform Convergence of Adaptive Multigrid Methods for Elliptic Problems and Maxwell's Equations

A note on the convergence of multigrid methods for the Riesz-space equation and an application to image deblurring

Monolithic Algebraic Multigrid Preconditioners for the Stokes Equations

Multigrid Methods for Hellan–Herrmann–Johnson Mixed Method of Kirchhoff Plate Bending Problems

Multigrid method for symmetric Toeplitz tridiagonal matrix: Convergence analysis and application to fractional Feynman-Kac equation

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

Achieving $h$- and $p$-robust monolithic multigrid solvers for the Stokes equations

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

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

Exploiting mesh structure to improve multigrid performance for saddle point problems

Numerical Stability of the C Method and a Perturbative Preconditioning Technique to Improve Convergence.

Preconditioners for Elliptic Problems Via Non-Uniform Meshes.