Paper detail

A PCG Implementation of an Elliptic Kernel in an Ocean Global Circulation Model Based on GPU Libraries

In this paper an inverse preconditioner for the numerical solution of an elliptic Laplace prob- lem of a global circulation ocean model is presented. The inverse preconditiong technique is adopted in order to efficiently compute the numerical solution of the elliptic kernel by using the Conjugate Gradient (CG) method. We show how the performance and the rate of conver- gence of the solver are linked to the discretized grid resolution and to the Laplace coefficients of the oceanic model. Finally, we describe an easy-to-implement version of the solver on the Graphical Processing Units (GPUs) by means of scientific computing libraries and we discuss its performance.