A new acceleration technique for the computation of first order hydrodynamic coefficients for floating bodies in frequency domain and in deep water is proposed. It is based on the classical boundary element method (BEM) which requires solving a boundary integral equation for distributions of sources and/or dipoles and evaluating integrals of Kelvin’s Green function and its derivatives over panels. The Kelvin’s Green function includes two Rankine sources and a wave term. In present study, for the two Rankine sources, analytical integrations of strongly singular kernels are adopted for the linear density distributions. It is shown that these analytical integrations are more accurate and faster than numerical integrations. The wave term is obtained by solving Clément’s ordinary differential equations (ODEs) [1] and an adaptive numerical quadrature is performed for integrations over the panels. It is shown here that the computational time of the wave term by solving the ODEs is greatly reduced compared to the classical integration method [7].