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.

Performance comparison of interpolation operators of multigrid methods based on distance weight and area (volume) intersection approaches Part I: Finite volume discretization

A quasi-three dimensional multi-level Green's function interpolation method for multilayered structures

A Directional Interpolation Pattern for Multilevel Green's Function Interpolation Method

The Implementation of Multilevel Green's Function Interpolation Method for Full-Wave Electromagnetic Problems

A Novel Interpolation Scheme for Full-Wave Electromagnetic Simulations

A comparison of interpolation techniques for non-conformal high-order discontinuous Galerkin methods

Fast, high-order numerical evaluation of volume potentials via polynomial density interpolation

Efficient Interpolation Scheme for Multilevel Fast Multipole Algorithm

Operator-difference approximations on two-dimensional merged Voronoi-Delaunay grids

Interpolation Conditions for Linear Operators and Applications to Performance Estimation Problems

Order Preserving Interpolation for Summation-by-Parts Operators at Non-Conforming Grid Interfaces

A Novel Variational Multiscale Interpolating Element-Free Galerkin Method for Generalized Oseen Problems

A high-efficiency sharp-interface immersed boundary method based on multi-linear interpolation

Comparison of Interpolation Functions in Control Volume Finite Element Method and Numerical Analysis

A variational multiscale interpolating element-free Galerkin method for convection-diffusion and Stokes problems

An Efficient Quadratic Interpolation Scheme for a Third-Order Cell-Centered Finite-Volume Method on Tetrahedral Grids

Comparative Evaluation of Meshfree Data Transfer Methods in Hydromechanical Problems

A comparison of element agglomeration algorithms for unstructured geometric multigrid

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

Unified geometric multigrid algorithm for hybridized high-order finite element methods