Accelerating Lagrangian Particle Dispersion in the Atmosphere with OpenCL

Research Paper

Presented by Paul Harvey, University of Glasgow

FLEXPART is a popular simulator that models the transport and diffusion of air pollutants, based on the Lagrangian approach. It is capable of regional and global simulation and supports both forward and backwards runs. A complex model like this contains many calculations suitable for parallelisation. Recently, a GPU-accelerated version of the simulator (FLEXCPP) has been written in C++/CUDA. As CUDA is only supported on NVIDIA GPUs, such an implementation is tied to a single hardware vendor, and is not able to take advantage of other hardware acceleration platforms.

This paper presents an OpenCL/C++ version of FLEX-CPP, called FlexOcl. This simulator provides all the functionality of FLEXCPP, and has been extended to include modelling of the decay of radioactive particles. A performance comparison between the two simulators has been performed on GPU, and the performance of FlexOcl has also been evaluated on the Intel Xeon Phi. Our preliminary results show that the OpenCL code performs as well as the CUDA code on GPUs, and that the MIC outperforms the GPU platforms for this type of application.