Evaluation of the Coarray Fortran Programming Model on the Example of a Lattice Boltzmann Code Klaus Sembritzki Georg Hager Bettina Krammer Friedrich-Alexander University Friedrich-Alexander University Université de Versailles Erlangen-Nuremberg Erlangen-Nuremberg St-Quentin-en-Yvelines Erlangen Regional Computing Erlangen Regional Computing Exascale Computing Center (RRZE) Center (RRZE) Research (ECR) Martensstrasse 1 Martensstrasse 1 45 Avenue des Etats-Unis 91058 Erlangen, Germany 91058 Erlangen, Germany 78000 Versailles, France
[email protected] [email protected] [email protected] Jan Treibig Gerhard Wellein Friedrich-Alexander University Friedrich-Alexander University Erlangen-Nuremberg Erlangen-Nuremberg Erlangen Regional Computing Erlangen Regional Computing Center (RRZE) Center (RRZE) Martensstrasse 1 Martensstrasse 1 91058 Erlangen, Germany 91058 Erlangen, Germany
[email protected] [email protected] ABSTRACT easier to learn. The Lattice Boltzmann method is an explicit time-stepping sche- With the Fortran 2008 standard, coarrays became a native Fort- me for the numerical simulation of fluids. In recent years it has ran feature [2]. The “hybrid” nature of modern massively parallel gained popularity since it is straightforward to parallelize and well systems, which comprise shared-memory nodes built from multico- suited for modeling complex boundary conditions and multiphase re processor chips, can be addressed by the compiler and runtime flows. Starting from an MPI/OpenMP-parallel 3D prototype im- system by using shared memory for intra-node access to codimen- plementation of the algorithm in Fortran90, we construct several sions. Inter-node access may either by conducted using an existing coarray-based versions and compare their performance and requi- messaging layer such as MPI or SHMEM, or carried out natively red programming effort to the original code, demonstrating the per- using the primitives provided by the network infrastructure.