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 Multigrid Solver for PDE-Constrained Optimization with Uncertain Inputs

A multigrid solver for PDE-constrained optimization with uncertain inputs

Multigrid preconditioning of linear systems for semi-smooth Newton methods applied to optimization problems constrained by smoothing operators

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

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

Exploiting mesh structure to improve multigrid performance for saddle point problems

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

A multigrid reduction framework for domains with symmetries

A Symmetric Multigrid-Preconditioned Krylov Subspace Solver for Stokes Equations

Algebraic multigrid methods for metric-perturbed coupled problems

A geometric multigrid preconditioning strategy for DPG system matrices

Multigrid preconditioning of singularly perturbed convection-diffusion equations

Multigrid for two-sided fractional differential equations discretized by finite volume elements on graded meshes

Multigrid optimization methods for linear and bilinear elliptic optimal control problems

A reinforcement learning strategy to automate and accelerate h/p-multigrid solvers

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

Scheduling massively parallel multigrid for multilevel Monte Carlo methods

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

Scalable DPG Multigrid Solver for Helmholtz Problems: A Study on Convergence

A Note on the Convergence of Multigrid Methods for the Riesz–Space Equation and an Application to Image Deblurring

A Scalable Interior-Point Gauss-Newton Method for PDE-Constrained Optimization with Bound Constraints