Block Iterative Methods and Recycling for Improved Scalability of Linear Solvers

Pierre Jolivet (1) , Pierre-Henri Tournier (2)

1 ENSEEIHT - Ecole Nationale Supérieure d'Electrotechnique, d'Electronique, d'Informatique, d'Hydraulique et de Télécommunications
2 ALPINES - Algorithms and parallel tools for integrated numerical simulations

Contemporary large-scale Partial Differential Equation (PDE) simulations usually require the solution of large and sparse linear systems. Moreover, it is often needed to solve these linear systems with different or multiple Right-Hand Sides (RHSs). In this paper, various strategies will be presented to extend the scalability of existing linear solvers using appropriate recycling strategies or block methods—i.e., by treating multiple right-hand sides simultaneously. The scalability of this work is assessed by performing simulations on up to 8,192 cores for solving linear systems arising from various physical phenomena modeled by Poisson's equation, the system of linear elasticity, or Maxwell's equation. This work is shipped as part of on open-source software, readily available and usable in any C, C++, or Python code. In particular, some simulations are performed on top of a well-established library, PETSc, and it is shown how our approaches can be used to decrease time to solution down by 30%.

Mots clés: Iterative methods; Distributed algorithms; Maxwell’s equation

Domaines: Equations aux dérivées partielles [math.AP]

Pierre Jolivet, Pierre-Henri Tournier. Block Iterative Methods and Recycling for Improved Scalability of Linear Solvers. SC16 - International Conference for High Performance Computing, Networking, Storage and Analysis, Nov 2016, Salt Lake City, Utah, United States. ⟨hal-01357998⟩