Recovering numerical reproducibility in hydrodynamic simulations - LIRMM - Laboratoire d’Informatique, de Robotique et de Microélectronique de Montpellier Access content directly
Conference Papers Year : 2016

Recovering numerical reproducibility in hydrodynamic simulations


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⟩
275 View
456 Download



Gmail Mastodon Facebook X LinkedIn More