A portable OpenCL implementation of generic particle-mesh and mesh-particle interpolation in 2D and 3D

Publikation: Beitrag in FachzeitschriftForschungsartikelBeigetragenBegutachtung

Beitragende

  • Ferit Büyükkeçeci - , ETH Zurich (Autor:in)
  • Omar Awile - , ETH Zurich (Autor:in)
  • Ivo F. Sbalzarini - , ETH Zurich (Autor:in)

Abstract

Hybrid particle-mesh methods provide a versatile framework for simulating discrete and continuous systems. A key component is the forward and backward interpolation of particle data to mesh nodes. These interpolations typically account for a significant portion of the computational cost of a simulation. Due to its regular compute structure, interpolation admits SIMD parallelism, and several GPU-accelerated implementations have been presented in the literature. We build on these works to develop a streaming-parallel algorithm for interpolation in hybrid particle-mesh methods that works in both 2D and 3D and is free of assumptions about the particle density, the number of particle properties to be interpolated, and the particle indexing scheme. We provide a portable OpenCL implementation of the algorithm and benchmark its accuracy and performance. We show that with such a generic algorithm speedups of up to 15× over an 8-core multi-thread CPU implementation are possible if the data are already available on the GPU. The maximum speedup reduces to about 7× if the data first have to be transferred to the GPU. The benchmarks also expose several limitations of GPU acceleration, in particular for low-order and 2D interpolation schemes. The present algorithm is integrated and available in the open-source Parallel Particle Mesh (PPM) library as a hybrid MPI-OpenCL implementation.

Details

OriginalspracheEnglisch
Seiten (von - bis)94-111
Seitenumfang18
FachzeitschriftParallel Computing
Jahrgang39
Ausgabenummer2
PublikationsstatusVeröffentlicht - 2013
Peer-Review-StatusJa
Extern publiziertJa

Externe IDs

ORCID /0000-0003-4414-4340/work/159608292

Schlagworte

Schlagwörter

  • GPGPU, Interpolation, OpenCL, Particle-mesh method, PIC method, PPM library