A CUBLAS-CUDA Implementation of PCG Method of an Ocean Circulation Model
R. Farina, Salvatore Cuomo, Pasquale De Michele, Theodore E. Simos, George Psihoyios, Ch. Tsitouras, Zacharias A. Anastassi · AIP conference proceedings · 2011
A numerical model of an ocean global circulation is presented. It consists of the discretization of Laplace’s problem by means of finite differences scheme of second order that gives a linear system solved by Preconditioned Conjugate Gradient Method (PCG). In this work, we observe that the performance and the accuracy of PCG solver depend on the grid resolution and the ratio of Laplace’s coefficients. Moreover, a case study of an implementation of PCG solver with diagonal preconditioner on multi‐core GPU architecture based on CUBLAS library is proposed.