LAUSR

In this paper the implementation of an iterative solver based on the Generalized Minimum Residual Method (GMRES) with complex-valued coefficients is explored, with application to the Boundary Element Method (BEM). The solver is designed to be implemented in a GPU (Graphic Processing Unit) device, exploiting its massively parallel capabilities. The framework is in the context of linear water wave problems in the frequency domain. The main objective of the proposed solver is the direct acceleration of existing standard BEM codes, by transfering to the GPU the solver task. The CUDA programming language is used, exploiting the particular architecture of the GPU device for complex-valued systems. To explore the performances of the solver, two linear water wave problems have been tested: the frequency-dependent added mass and damping matrices of a 3D floating body, and the Helmholtz equation in a 2D domain. The parallelized GMRES solver is implemented in a NVidia GeForce GTX 970 graphic card, and shows drastic reductions in computing times when compared with its CPU implementation.