—In this study, we report development of two innovatory parallel implementations of the Multidimensional Positive Definite Advection Transport Algorithm (MPDATA). MPDATA is the main module of the multiscale fluid model EULAG [Prusa et al., Computers & Fluids, vol. 37, 2008]. Recently, the dynamical core of EULAG has been implemented into COSMO (Consortium for Small-scale Modeling) weather prediction framework and is expected to be in operational use.
The original code of MPDATA is written in FORTRAN77 and has excellent efficiency and scalability on conventional supercomputer architectures. The new C++ implementations are designed and optimized under modern CPU and GPU based high-performance computing platforms. A number of innovatory solutions have been employed in these codes, including 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. The new techniques allow to take the full advantage of the new architectures and accelerate the code execution.
The correctness and accuracy of the new implementations are examined based on a standard three-dimensional solid body rotation test case. Additionally, we focus on testing computational efficiency and scalability. In most runs and especially in simulations with larger computational grids the new codes performed better than the traditional implementation.
—MPDATA, advection solver, parallel computing.
B. Rosa, A. A. Wyszogrodzki, 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: email@example.com).
L. Szustak, K. Rojek, and R. Wyrzykowski are with the Czestochowa University of Technology, Dabrowskiego 69 Street, 42-201 Czestochowa, Poland.
Cite: Bogdan Rosa, Lukasz Szustak, Andrzej A. Wyszogrodzki, Krzysztof Rojek, Damian K. Wójcik, and Roman Wyrzykowski, "Adaptation of Multidimensional Positive Definite Advection Transport Algorithm to Modern High-Performance Computing Platforms," International Journal of Modeling and Optimization vol. 5, no. 3, pp. 171-176, 2015.