Single Scattering of Aspherical Particles in DDA Calculations on GPUs Using OpenCL
University of Bremen, Institute of Environmental Physics
Technical Report, Institute of Environmental Physics, University of Bremen, January 2011
@article{huntemann2011single,
title={Single Scattering of Aspherical Particles in DDA Calculations on GPUs Using OpenCL},
author={Huntemann, M. and Heygster, G. and Hong, G.},
year={2011}
}
The global distribution and climatology of ice clouds are among the main uncertainties in climate modelling and prediction. In order to retrieve ice cloud properties from remote sensing measurements, the scattering properties of all cloud ice particle types must be known. The Discrete Dipole Approximation (DDA) simulates scattering of radiation by arbitrarily shaped particles and is thus suitable for cloud ice crystals. The DDA models the particle as a collection of equal dipoles on a lattice, and is computationally much more expensive than approximations restricted to more regularly shaped particles. On a single computer the calculation for an ice particle of a specific size, for a given scattering plane at one specific wavelength can take several days. We have ported the core routines of the scattering suite "Amsterdam DDA" (ADDA) to the Open Computing Language (OpenCL), a framework for programming parallel devices like PC graphics cards (Graphics Processing Units, GPUs) or multicore CPUs. In a typical case we can achieve a speed-up on a GPU as compared to a CPU by a factor of 5 in double precision and a factor of 15 in single precision. Spreading the work load over multiple GPUs will allow calculating the scattering properties even of large cloud ice particles.
September 24, 2011 by hgpu
