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.

An AMG saddle point preconditioner with application to mixed Poisson problems on adaptive quad/cube meshes

Monolithic Algebraic Multigrid Preconditioners for the Stokes Equations

An approximate block factorization preconditioner for mixed-dimensional beam-solid interaction

New preconditioner strategy for solving block four-by-four linear systems: An application to the saddle-point problem from 3D Stokes equation

A RELAXED HSS PRECONDITIONER FOR SADDLE POINT PROBLEMS FROM MESHFREE DISCRETIZATION

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

Exploiting mesh structure to improve multigrid performance for saddle point problems

Regularized DPSS Preconditioners for Generalized Saddle Point Linear Systems

A multigrid reduction framework for domains with symmetries

An Amg Preconditioner for Solving the Navier-Stokes Equations with A Moving Mesh Finite Element Method

Solving the Poisson equation on small aspect ratio domains using unstructured meshes

A Block Positive-Semidefinite Splitting Preconditioner for Generalized Saddle Point Linear Systems.

Double Saddle-Point Preconditioning for Krylov Methods in the Inexact Sequential Homotopy Method

Robust PRESB Preconditioning of a 3-Dimensional Space-Time Finite Element Method for Parabolic Problems

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

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

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

Permuted preconditioning for extended saddle point problem arising from Neumann boundary control

Fast Auxiliary Space Preconditioners for Linear Elasticity in Mixed Form

Block Preconditioning Methods for Asymptotic Preserving Scheme Arising in Anisotropic Elliptic Problems

Matrix-free GPU-accelerated saddle-point solvers for high-order problems in $H(\mathrm{div})$