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.

A robust all-at-once multigrid method for the Stokes control problem

Efficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems

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

A staggered-grid multilevel incomplete LU for steady incompressible flows

A Symmetric Multigrid-Preconditioned Krylov Subspace Solver for Stokes Equations

Multigrid methods for the Stokes problem on GPU systems

Symbol based Convergence Analysis in Block Multigrid Methods with Applications for Stokes Problems

Monolithic Algebraic Multigrid Preconditioners for the Stokes Equations

Solution Techniques for the Stokes System: A priori and a posteriori modifications, resilient algorithms

Coarse-grid operator optimization in multigrid reduction in time for time-dependent Stokes and Oseen problems

RBF collocation and hybrid-LHI methods for Stokes systems and its application to controllability problems

A Multigrid Discretization of Discontinuous Galerkin Method for the Stokes Eigenvalue Problem

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

A two-grid Adaptive Finite Element Method for the Dirichlet Boundary Control Problem Governed by Stokes Equation

An augmented Lagrangian preconditioner for the control of the Navier--Stokes equations

Scalable robust solvers for unstructured FE geodynamic modeling applications: Solving the Stokes equation for models with large localized viscosity contrasts

Optimal Control of the Multiphase Stefan Problem

Robust finite element solvers for distributed hyperbolic optimal control problems

Novel Approach for solving the discrete Stokes problems based on Augmented Lagrangian and Global Techniques: Application to Saddle-Point Linear Systems from Incompressible flow

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

Rotated block diagonal preconditioners for Navier-Stokes control problems