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.

Accelerating Algebraic Multigrid Solvers on NVIDIA GPUs

Development of Krylov and AMG linear solvers for large-scale sparse matrices on GPUs

An adaptive finite element multigrid solver using GPU acceleration

AmgX: A Library for GPU Accelerated Algebraic Multigrid and Preconditioned Iterative Methods

Parallel Smoothers for Matrix-Based Geometric Multigrid Methods on Locally Refined Meshes Using Multicore CPUs and GPUs

Multi-GPU aggregation-based AMG preconditioner for iterative linear solvers

Optimizing Finite Volume Method Solvers on Nvidia GPUs.

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

Parallel Triangular Solvers on GPU

Parallel Implementation of the Coordinates-partitioning Based Aggregation-type Algebraic Multigrid Preconditioners

Towards a Heterogeneous Architecture Solver for the Incompressible Navier–Stokes Equations

Exploiting mesh structure to improve multigrid performance for saddle point problems

Generalized GPU Acceleration for Applications Employing Finite-Volume Methods.

A Massively Parallel Algebraic Multigrid Preconditioner based on Aggregation for Elliptic Problems with Heterogeneous Coefficients

Extending DD-$α$AMG on heterogeneous machines

GPU-accelerated Linear Algebra for Coupled Solvers in Industrial CFD Applications with OpenFOAM

Research of the Algebraic Multigrid Method for Electron Optical Simulator

Optimization and acceleration of flow simulations for CFD on CPU/GPU architecture

Algebraic temporal blocking for sparse iterative solvers on multi-core CPUs

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

On Solving Groundwater Flow and Transport Models with Algebraic Multigrid Preconditioning