Parallel 5 point SOR for solving the Convection Diffusion equation using graphics processing units
Department of Informatics and Telecommunications, University of Athens, Panepistimiopolis, 15784, Athens, Greece
18th International Conference, Euro-Par 2012, Volume: 7484, 2012
@article{cotronis2012parallel,
title={Parallel 5 point SOR for solving the Convection Diffusion equation using graphics processing units},
author={Cotronis, Y and Konstantinidis, E and Louka, MA and Missirlis, NM},
year={2012}
}
In this paper we study a parallel form of the SOR method for the numerical solution of the Convection Diffusion equation suitable for GPUs using CUDA. To exploit the parallelism offered by GPUs we consider the fine grain parallelism model. This is achieved by considering the local relaxation version of SOR. More specifically, we use SOR with red black ordering with two sets of parameters omega_ij and omega’_ij . The parameter omega_ij is associated with each red (i+j even) grid point (ij), whereas the parameter omega’_ij is associated with each black (i+j odd) grid point (ij). The use of a parameter for each grid point avoids the global communication required in the adaptive determination of the best value of omega and also increases the convergence rate of the SOR method [3]. We present our strategy and the results of our effort to exploit the computational capabilities of GPUs under the CUDA environment. Additionally, a program for the CPU was developed as a performance reference. Significant performance improvement was achieved with the three developed GPU kernel variations which proved to have different pros and cons.
May 20, 2013 by hgpu
