Numerical resolution of conservation laws with OpenCL
Institut de Recherche Mathematique Avancee (IRMA), CNRS: UMR7501 – Universite de Strasbourg
hal-00759131 (30 November 2012)
@unpublished{crestetto:hal-00759131,
hal_id={hal-00759131},
url={http://hal.archives-ouvertes.fr/hal-00759131},
title={Numerical resolution of conservation laws with OpenCL},
author={Crestetto, Ana{"i}s and Helluy, Philippe and Jung, Jonathan},
keywords={Discontinuous Galerkin, Finite Volume, GPU, OpenCL},
language={Anglais},
affiliation={Institut de Recherche Math{‘e}matique Avanc{‘e}e – IRMA , CALVI – INRIA Nancy – Grand Est / IECN / LSIIT / IRMA},
year={2012},
pdf={http://hal.archives-ouvertes.fr/hal-00759131/PDF/esaim-opencl.pdf}
}
We present several numerical simulations of conservation laws on recent multicore processors, such as GPU’s, using the OpenCL programming framework. Depending on the chosen numerical method, different implementation strategies have to be considered, for achieving the best performance. We explain how to program efficiently three methods: a finite volume approach on a structured grid, a high order Discontinuous Galerkin (DG) method on an unstructured grid and a Particle-In-Cell (PIC) method. The three methods are respectively applied to a two-fluid computation, a Maxwell simulation and a Vlasov-Maxwell simulation.
December 2, 2012 by hgpu
