Abstract : HPC simulations suffer from failures of numerical reproducibility because of floating-point arithmetic peculiarities. Different computing distributions of a parallel computation may yield different numerical results. We are interested in a finite element computation of hydrodynamic simulations within the openTelemac software where parallelism is provided by domain decomposition. One main task in a finite element simulation consists in building one large linear system and to solve it. Here the building step relies on element-by-element storage mode and the solving step applies the conjugated gradient algorithm. The subdomain parallelism is merged within these steps. We study why reproducibility fails in this process and which operations have to be corrected. We detail how to use compensation techniques to compute a numerically reproducible resolution. We illustrate this approach presenting the reproducible version of hydrodynamic simulations for one test cases provided with the openTelemac software suite.
https://hal-lirmm.ccsd.cnrs.fr/lirmm-01274671
Contributor : Rafife Nheili
<>
Submitted on : Tuesday, February 16, 2016 - 10:07:22 AM
Last modification on : Thursday, May 24, 2018 - 3:59:23 PM
Long-term archiving on : Tuesday, May 17, 2016 - 5:23:25 PM
Philippe Langlois, Rafife Nheili, Christophe Denis. Recovering numerical reproducibility in hydrodynamic simulations. ARITH: Computer Arithmetic, Jul 2016, Silicon Valley, Santa Clara, CA, United States. ⟨lirmm-01274671⟩
