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.

Extrapolation-accelerated aggregation multigrid with block relaxations for GeneRank problems

Block-accelerated Aggregation Multigrid for Markov Chains with Application to PageRank Problems

A Cost-Effective Smoothed Multigrid with Modified Neighborhood-Based Aggregation for Markov Chains

Compact adaptive aggregation multigrid method for Markov chains

An economical aggregation algorithm for algebraic multigrid. (AMG)

An Aggregation-Based Two-Grid Method for Multilevel Block Toeplitz Linear Systems

Two Schemes to Improve the Implementation of the Aggregation Based Algebraic Multigrid Preconditioner

The Extrapolation-Accelerated Multilevel Aggregation Method in PageRank Computation

Analysis on aggregation and block smoothers in multigrid methods for block Toeplitz linear systems

Krylov subspace algorithms for computing GeneRank for the analysis of microarray data mining.

Fast algebraic multigrid for block-structured dense systems arising from nonlocal diffusion problems

Efficient Solution of 3D Elasticity Problems with Smoothed Aggregation Algebraic Multigrid and Block Arithmetics

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

Generalizing Reduction-Based Algebraic Multigrid

Two Novel Aggregation-Based Algebraic Multigrid Methods

Exploiting mesh structure to improve multigrid performance for saddle point problems

Light Weight Coarse Grid Aggregation for Smoothed Aggregation Algebraic Multigrid Solver

Understanding the GeneRank Model

A class of accelerated GADMM-based method for multi-block nonconvex optimization problems

Inference of Large-scale Time-delayed Gene Regulatory Network with Parallel MapReduce Cloud Platform

The finest level acceleration of multilevel aggregation for markov chains