TY - JOUR
T1 - PHIST
T2 - A Pipelined, Hybrid-Parallel Iterative Solver Toolkit
AU - Thies, Jonas
AU - Röhrig-Zöllner, Melven
AU - Overmars, Nigel
AU - Basermann, Achim
AU - Ernst, Dominik
AU - Hager, Georg
AU - Wellein, Gerhard
PY - 2020/11
Y1 - 2020/11
N2 - The increasing complexity of hardware and software environments in high-performance computing poses big challenges on the development of sustainable and hardware-efficient numerical software. This article addresses these challenges in the context of sparse solvers. Existing solutions typically target sustainability, flexibility, or performance, but rarely all of them. Our new library PHIST provides implementations of solvers for sparse linear systems and eigenvalue problems. It is a productivity platform for performance-aware developers of algorithms and application software with abstractions that do not obscure the view on hardware-software interaction. The PHIST software architecture and the PHIST development process were designed to overcome shortcomings of existing packages. An interface layer for basic sparse linear algebra functionality that can be provided by multiple backends ensures sustainability, and PHIST supports common techniques for improving scalability and performance of algorithms such as blocking and kernel fusion. We showcase these concepts using the PHIST implementation of a block Jacobi-Davidson solver for non-Hermitian and generalized eigenproblems. We study its performance on a multi-core CPU, a GPU, and a large-scale many-core system. Furthermore, we show how an existing implementation of a block Krylov-Schur method in the Trilinos package Anasazi can benefit from the performance engineering techniques used in PHIST.
AB - The increasing complexity of hardware and software environments in high-performance computing poses big challenges on the development of sustainable and hardware-efficient numerical software. This article addresses these challenges in the context of sparse solvers. Existing solutions typically target sustainability, flexibility, or performance, but rarely all of them. Our new library PHIST provides implementations of solvers for sparse linear systems and eigenvalue problems. It is a productivity platform for performance-aware developers of algorithms and application software with abstractions that do not obscure the view on hardware-software interaction. The PHIST software architecture and the PHIST development process were designed to overcome shortcomings of existing packages. An interface layer for basic sparse linear algebra functionality that can be provided by multiple backends ensures sustainability, and PHIST supports common techniques for improving scalability and performance of algorithms such as blocking and kernel fusion. We showcase these concepts using the PHIST implementation of a block Jacobi-Davidson solver for non-Hermitian and generalized eigenproblems. We study its performance on a multi-core CPU, a GPU, and a large-scale many-core system. Furthermore, we show how an existing implementation of a block Krylov-Schur method in the Trilinos package Anasazi can benefit from the performance engineering techniques used in PHIST.
KW - GPU computing
KW - hybrid parallel computing
KW - Jacobi-Davidson
KW - performance engineering
KW - Sparse eigenvalue solvers
UR - http://www.scopus.com/inward/record.url?scp=85096837577&partnerID=8YFLogxK
U2 - 10.1145/3402227
DO - 10.1145/3402227
M3 - Article
AN - SCOPUS:85096837577
SN - 0098-3500
VL - 46
JO - ACM Transactions on Mathematical Software
JF - ACM Transactions on Mathematical Software
IS - 4
M1 - 3402227
ER -