Recovering numerical reproducibility in hydrodynamic simulations - Archive ouverte HAL Access content directly
Conference Papers Year : 2016

Recovering numerical reproducibility in hydrodynamic simulations

(1) , (1) , (2)


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.
Fichier principal
Vignette du fichier
LaNhDe_2016 (1).pdf (1.72 Mo) Télécharger le fichier
Origin : Files produced by the author(s)

Dates and versions

lirmm-01274671 , version 1 (16-02-2016)



Philippe Langlois, Rafife Nheili, Christophe Denis. Recovering numerical reproducibility in hydrodynamic simulations. ARITH: Computer Arithmetic, Jul 2016, Silicon Valley, Santa Clara, CA, United States. pp.63-70, ⟨10.1109/ARITH.2016.27⟩. ⟨lirmm-01274671⟩
266 View
401 Download



Gmail Facebook Twitter LinkedIn More