Abstract:To address the problem of massive data computation in the process of gravity field calculation, we jointly complete a parallel acceleration algorithm for gravity field calculation based on the least squares method using MPI(massage passing interface) and CUDA(compute unified device architecture). MPI is used to complete the task allocation of complex processes and achieve global-level parallel acceleration. We write a parallel acceleration program for large-scale matrix multiplication based on CUDA and adapted for different types of matrices. At the same time, we use MPI to further subdivide the calculation process of the law matrix to compress the peak memory value within the sub-process. We complete the 30 and 120 order gravity field simulation calculation tasks on a single machine. The results show that when inverting the 30 order gravity field, the acceleration ratio can reach 180; when inverting the 120 order gravity field, parallel computing only takes 2 hours for a single iteration, while it cannot be calculated in serial mode.