Volume 4 Number 3 (May 2014)
Home > Archive > 2014 > Volume 4 Number 3 (May 2014) >
IJAPM 2014 Vol.4(3): 188-195 ISSN: 2010-362X
DOI: 10.7763/IJAPM.2014.V4.281

Porting Multiscale Fluid Model EULAG to Modern Heterogeneous Architectures

Bogdan Rosa, Miłosz Ciżnicki, Krzysztof A. Rojek, Damian K. Wójcik, Piotr K. Smolarkiewicz, and Roman Wyrzykowski
Abstract—The goal of this study is to adapt the multiscale fluid solver EULAG [Prusa et al., Computers & Fluids, vol 37, 2008] to the future GPU-based high-performance computing platforms. The EULAG model has a proven record of successful applications in a range of environmental fluid dynamics, and excellent efficiency and scalability on conventional supercomputer architectures. Currently, the model is being implemented as a new dynamical core of COSMO (Consortium for Small-scale Modeling) weather prediction framework. The EULAG code combines features of a stencil and point wise computations. Its communication scheme consists of both halo exchange subroutines and global reduction functions. Within the project, two main modules of EULAG, namely the multidimensional positive definite advection transport algorithm, MPDATA, and the variational generalized conjugate residual, GCR, elliptic pressure solver are analyzed and optimized. Relevant techniques have been chosen and applied to accelerate code execution on modern GPU architectures: stencil decomposition, block decomposition (with weighting analysis between computation and communication), reduction of intercache communication by partitioning of cores into independent teams, cache reusing and vectorization.
Testing and validation of the new GPU implementation have been carried out based on modeling decaying turbulence of a homogeneous incompressible fluid in a triply-periodic cube. Simulations performed using the standard version of EULAG and its new GPU implementation give similar solutions. Preliminary results of the parallel performance of the new implementation show a promising increase in terms of computational efficiency.

Index Terms—Anelastic model, EULAG, hybrid architectures, parallel computing.

B. Rosa and D. K. Wójcik are with the Institute of Meteorology and Water Management - National Research Institute, Podleśna 61 Street, 01-673 Warsaw, Poland (e-mail: bogdan.rosa@imgw.pl).
P. K. Smolarkiewicz is with the European Centre for Medium-Range Weather Forecasts, Reading, RG2 9AX, UK.
M. Ciżnicki is with the Poznan Supercomputing and Networking Center, Noskowskiego 10 Street, 61-704 Poznań, Poland.
K. A. Rojek and R. Wyrzykowski are with the Czestochowa University of Technology, Dabrowskiego 69 Street, 42-201 Czestochowa, Poland.

Cite: Bogdan Rosa, Miłosz Ciżnicki, Krzysztof A. Rojek, Damian K. Wójcik, Piotr K. Smolarkiewicz, and Roman Wyrzykowski, "Porting Multiscale Fluid Model EULAG to Modern Heterogeneous Architectures," International Journal of Applied Physics and Mathematics vol. 4, no. 3, pp. 188-195, 2014.

General Information

ISSN: 2010-362X (Online)
Abbreviated Title: Int. J. Appl. Phys. Math.
Frequency: Quarterly
DOI: 10.17706/IJAPM
Editor-in-Chief: Prof. Haydar Akca 
Abstracting/ Indexing:  Index Copernicus, EI (INSPEC, IET), Chemical Abstracts Services (CAS), Nanowerk Database, Google Scholar, EBSCO, etc.
E-mail: ijapm@iap.org
  • May 14, 2019 News!

    IJAPM Vol 7, No 4-Vol 8, No 3 have been indexed by EI (Inspec)   [Click]

  • Jun 17, 2020 News!

    Vol 10, No 3 has been published with online version   [Click]

  • Apr 02, 2020 News!

    The paper published in Vol 10, No 2 has received dois from Crossref

  • Mar 25, 2020 News!

    Vol 10, No 2 has been published with online version     [Click]

  • Dec 13, 2019 News!

    The paper published in Vol 10, No 1 has received dois from Crossref

  • Read more>>