利用CPU和GPU混合并行方法快速构建海洋扰动重力梯度基准图  

Rapid Construction of Marine Disturbing Gravity Gradient Reference Map by CPU and GPU Hybrid Parallel Method

作　　者：黄炎 李姗姗[2] 吕明昊 范雕 谭勖立 冯进凯 HUANG Yan;LI Shanshan;LU Minghao;FAN Diao;TAN Xuli;FENG Jinkai(National Innovation Institute of Defense Technology,Academy of Military Science,Beijing 100071,China;Institute of Geographical Spatial Information,Information Engineering University,Zhengzhou 450001,China;396941 Troops,Beijing 100080,China)

机构地区：[1]军事科学院国防科技创新研究院,北京100071 [2]信息工程大学地理空间信息学院,河南郑州450001 [3]96941部队,北京100080

出　　处：《武汉大学学报(信息科学版)》2025年第3期515-527,共13页Geomatics and Information Science of Wuhan University

基　　金：国家自然科学基金(42174007,42174013)。

摘　　要：高精度、高分辨率的海洋扰动重力梯度基准图是将水下扰动重力梯度辅助惯性导航付诸于实践的关键技术之一,可依据边值问题理论,基于移去-恢复技术计算求得。基于传统串行算法存在计算效率低、耗时长等问题,为满足大范围乃至全球海洋扰动重力梯度基准图的快速构建需求,利用中央处理器(central processing unit,CPU)与图形处理器(graph processing unit,GPU)混合并行(CPU+GPU)编程平台,设计了一种高效的混合并行计算方案。首先,利用数组收缩膨胀方法,有效解决了CPU进行勒让德函数递推计算过程中内存读写冲突问题;然后,引入Hilbert空间填充曲线,将二维格网重力异常数据降维成为具有高度聚簇性的一维数组,通过其在GPU显存中的快速索引,实现了Stokes积分的高效计算。在CPU型号为Intel Xeon(R)Gold 6130、GPU型号为Tesla V100的计算机上进行实验,利用该混合并行方案计算6°×6°范围、1′分辨率的全张量海洋扰动重力梯度基准图仅需15.84 s,而传统串行方法需要35 min以上;计算2°×2°范围、30″分辨率基准图仅需22.7 s,串行方法则需要52 min;计算2°×2°范围、10″分辨率基准图需要26.7 min,串行方法则需要两天以上。在保证并行计算绝对误差小于等于1×10^(-6) E、T_(xx)+T_(yy)+T_(zz)绝对值的最大值小于0.1 E、均方根误差小于0.01 E的前提下,实现了全球全张量海洋扰动重力梯度基准图模型的快速构建。Objectives:High precision and high resolution marine disturbing gravity gradient reference map is one of the key technologies to put underwater disturbing gravity gradient assisted inertial navigation into practice.At present,gravity anomaly is mainly measured by satellite and calculated based on the boundary value problem theory and the removal recovery technique.Methods:The traditional serial algorithm has some problems such as low efficiency and long time consuming.In order to meet the needs of rapid construction of marine disturbing gravity gradient reference maps on a large scale and even around the world,this paper designs an efficient hybrid parallel computing scheme using central processing unit(CPU)+graph processing unit(GPU)hybrid parallel programming platform.First,the shrink and expansion method of array is proposed to solve the problem of memory conflict between read and write in Legendre recursive com putation.Then Hilbert space filling curve is used to reduce the dimension of two-dimensional grid gravity anomaly data into a one-dimensional array with high clustering.The efficient calculation of Stokes integral is realized by the fast index in GPU video memory.Results:The experiments are carried out on computers with Intel Xeon(R)Gold 6130 CPU and Tesla V100 GPU.The hybrid parallel scheme is used to calculate 6°×6°range and 1′resolution full tensor marine disturbing gravity gradient reference map and only takes 15.84 seconds,while the traditional serial method takes more than 35 minutes.To calculate 2°×2°range and 30′′resolution map,the proposed method only takes 22.7 seconds,while the serial method takes 52 minutes.To calculate 2°×2°range and 10′′resolution map,the proposed method only takes 26.7 minutes,while the serial method takes more than two days.Conclusions:On the premise of ensuring that the absolute error of parallel computing is less than or equal to 1×10^(-6) E,the maximum absolute value of T_(xx)+T_(yy)+T_(zz) is less than 0.1 E and the root mean square error is less

关 键 词：扰动重力梯度 边值问题 CPU GPU 混合并行 

分 类 号：P223[天文地球—大地测量学与测量工程]