Efficient algebraic multigrid preconditioners on clusters of GPUs

Date published

2019-05-10

Free to read from

Supervisor/s

Journal Title

Journal ISSN

Volume Title

Publisher

World Scientific Publishing

Department

Type

Article

ISSN

0129-6264

Format

Citation

Abdullahi Hassan A, Cardellini V, D’Ambra P, et al., (2019) Efficient algebraic multigrid preconditioners on clusters of GPUs. Parallel Processing Letters, Volume 29, Issue 1, 2019, Article Number 1950001

Abstract

Many scientific applications require the solution of large and sparse linear systems of equations using Krylov subspace methods; in this case, the choice of an effective preconditioner may be crucial for the convergence of the Krylov solver. Algebraic MultiGrid (AMG) methods are widely used as preconditioners, because of their optimal computational cost and their algorithmic scalability. The wide availability of GPUs, now found in many of the fastest supercomputers, poses the problem of implementing efficiently these methods on high-throughput processors. In this work we focus on the application phase of AMG preconditioners, and in particular on the choice and implementation of smoothers and coarsest-level solvers capable of exploiting the computational power of clusters of GPUs. We consider block-Jacobi smoothers using sparse approximate inverses in the solve phase associated with the local blocks. The choice of approximate inverses instead of sparse matrix factorizations is driven by the large amount of parallelism exposed by the matrix-vector product as compared to the solution of large triangular systems on GPUs. The selected smoothers and solvers are implemented within the AMG preconditioning framework provided by the MLD2P4 library, using suitable sparse matrix data structures from the PSBLAS library. Their behaviour is illustrated in terms of execution speed and scalability, on a test case concerning groundwater modelling, provided by the Jülich Supercomputing Center within the Horizon 2020 Project EoCoE.

Description

Software Description

Software Language

Github

Keywords

Clusters of GPUs, algebraic multigrid, block-Jacobi smoothers, sparse approximate inverses

DOI

Rights

Attribution-NonCommercial 4.0 International

Relationships

Relationships

Supplements

Funder/s