high performance computing on graphics processing units: hgpu.org

AIMS: We discuss an implementation of our 3D radiative transfer (3DRT) framework with the OpenCL paradigm for general GPU computing. METHODS: We implemented the kernel for solving the 3DRT problem in Cartesian coordinates with periodic boundary conditions in the horizontal (x,y) plane, including the construction of the nearest neighbor ^* and the operator splitting step. RESULTS: We present the results of both a small and a large test case and compare the timing of the 3DRT calculations for serial CPUs and various GPUs. CONCLUSIONS: The latest available GPUs can lead to significant speedups for both small and large grids compared to serial (single core) computations.