Exposing fine-grained parallelism in algebraic multigrid methods

N Bell, S Dalton, LN Olson - SIAM Journal on Scientific Computing, 2012 - SIAM
N Bell, S Dalton, LN Olson
SIAM Journal on Scientific Computing, 2012SIAM
Algebraic multigrid methods for large, sparse linear systems are a necessity in many
computational simulations, yet parallel algorithms for such solvers are generally
decomposed into coarse-grained tasks suitable for distributed computers with traditional
processing cores. However, accelerating multigrid methods on massively parallel
throughput-oriented processors, such as graphics processing units, demands algorithms
with abundant fine-grained parallelism. In this paper, we develop a parallel algebraic …
Algebraic multigrid methods for large, sparse linear systems are a necessity in many computational simulations, yet parallel algorithms for such solvers are generally decomposed into coarse-grained tasks suitable for distributed computers with traditional processing cores. However, accelerating multigrid methods on massively parallel throughput-oriented processors, such as graphics processing units, demands algorithms with abundant fine-grained parallelism. In this paper, we develop a parallel algebraic multigrid method which exposes substantial fine-grained parallelism in both the construction of the multigrid hierarchy as well as the cycling or solve stage. Our algorithms are expressed in terms of scalable parallel primitives that are efficiently implemented on the GPU. The resulting solver achieves an average speedup of in the setup phase and in the cycling phase when compared to a representative CPU implementation.
Society for Industrial and Applied Mathematics