A 3D radiative transfer framework. VIII. OpenCL implementation
Abstract
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.
- Publication:
-
Astronomy and Astrophysics
- Pub Date:
- September 2011
- DOI:
- 10.1051/0004-6361/201117051
- Bibcode:
- 2011A&A...533A.127H
- Keywords:
-
- radiative transfer;
- methods: numerical