Application of numerical integration and CUDA parallel in the calculation of ocean gravity gradient

Through Stokes kernel function and gravity anomaly, it is possible to calculate the gravity gradient disturbance on the geoid and its external space. For this Stokes’ integral expression, we apply Laguerre wavelet numerical integration to improve the accuracy of its computational results. Meanwhile, compute unified device architecture (CUDA) is used to implement parallel computing on the Graphic Processing Unit (GPU) for speeding up. The full tensors of gravity gradient in the experimental ocean area with 3°×2°are computed. Compared to serial computing, its computing acceleration ratio can be more than 10 times faster. The results of the vertical gravity gradient are compared and validated from the public model from the University of California San Diego.