A parallel numerical algorithm by combining MPI and OpenMP programming models with applications in gravity field recovery

被引:2
作者
Lu, Biao [1 ,2 ]
Luo, Zhicai [3 ,4 ]
Zhong, Bo [5 ]
Zhou, Hao [3 ,4 ]
机构
[1] Sez Geodesia & Geomatica, Dipt Civile & Ambientale Ingn, Milan, Italy
[2] Univ Twente, Fac Geoinformat Sci & Earth Observat ITC, Dept Appl Earth Sci, Enschede, Netherlands
[3] Huazhong Univ Sci & Technol, PGMF & Sch Phys, MOE Key Lab Fundamental Phys Quant Measurement, Hubei Key Lab Gravitat & Quantum Phys, Wuhan, Peoples R China
[4] Huazhong Univ Sci & Technol, Inst Geophys, PGMF, Wuhan, Peoples R China
[5] Wuhan Univ, Sch Geodesy & Geomatics, MOE Key Lab Geospace Environm & Geodesy, Wuhan, Peoples R China
基金
荷兰研究理事会;
关键词
gravity field recovery; GOCE; GRACE; polar gravity data; MPI; OpenMP; parallel numerical algorithm; LINEAR ALGEBRA SUBPROGRAMS; CRUSTAL THICKNESS; ANTARCTICA; STRATEGIES; GRAVIMETRY; STRESS; SET;
D O I
10.3389/feart.2023.1080879
中图分类号
P [天文学、地球科学];
学科分类号
07 ;
摘要
Satellite gravimetry missions have enabled the calculation of high-accuracy and high-resolution Earth gravity field models from satellite-to-satellite tracking data and gravitational gradients. However, calculating high maximum degree/order (e.g., 240 or even higher) gravity field models using the least squares method is time-consuming due to the vast amount of gravimetry observations. To improve calculation efficiency, a parallel algorithm has been developed by combining Message Passing Interface (MPI) and Open Multi-Processing (OpenMP) programming models to calculate and invert normal equations for the Earth gravity field recovery. The symmetrical feature of normal equations has been implemented to speed up the calculation progress and reduce computation time. For example, the computation time to generate the normal equation of an IGGT_R1 test version of degree/order 240 was reduced from 88 h to 27 h by considering the symmetrical feature. Here, the calculation was based on the high-performance computing cluster with 108 cores in the School of Geodesy and Geomatics, at Wuhan University. Additionally, the MPI parallel Gaussian-Jordan elimination method was modified to invert normal equation matrices and scaled up to 100 processor cores in this study while the traditional method was limited in a certain number of processors. Furthermore, the Cholesky decomposition from the ScaLAPACK library was used to compare with the parallel Gauss-Jordan elimination method. The numerical algorithm has effectively reduced the amount of calculation and sped up the calculation progress, and has been successfully implemented in applications such as building the gravity field models IGGT_R1 and IGGT_R1C.
引用
收藏
页数:10
相关论文
共 67 条
[1]  
[Anonymous], 2012, GFZ GRACE LEVEL 2 PR, DOI DOI 10.2312/GFZ.B103-12020
[2]   Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition [J].
Baur, Oliver ;
Austen, Gerrit ;
Kusche, Juergen .
JOURNAL OF GEODESY, 2008, 82 (4-5) :207-221
[3]   Tailored least-squares solvers implementation for high-performance gravity field research [J].
Baur, Oliver .
COMPUTERS & GEOSCIENCES, 2009, 35 (03) :548-556
[4]  
Blackford L. S., 1997, SCALAPACK USERS GUID
[5]   An updated set of Basic Linear Algebra Subprograms (BLAS) [J].
Blackford, LS ;
Demmel, J ;
Dongarra, J ;
Duff, I ;
Hammarling, S ;
Henry, G ;
Heroux, M ;
Kaufman, L ;
Lumsdaine, A ;
Petitet, A ;
Pozo, R ;
Remington, K ;
Whaley, RC .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2002, 28 (02) :135-151
[6]   Antarctic crustal thickness from satellite gravity: Implications for the Transantarctic and Gamburtsev Subglacial Mountains [J].
Block, Adrienne E. ;
Bell, Robin E. ;
Studinger, Michael .
EARTH AND PLANETARY SCIENCE LETTERS, 2009, 288 (1-2) :194-203
[7]   Use of High Performance Computing for the Rigorous Estimation of Very High Degree Spherical Harmonic Gravity Field Models [J].
Brockmann, Jan Martin ;
Roese-Koerner, Lutz ;
Schuh, Wolf-Dieter .
GRAVITY, GEOID AND HEIGHT SYSTEMS, 2014, 141 :27-33
[8]   A concept for the estimation of high-degree gravity field models in a high performance computing environment [J].
Brockmann, Jan Martin ;
Roese-Koerner, Lutz ;
Schuh, Wolf-Dieter .
STUDIA GEOPHYSICA ET GEODAETICA, 2014, 58 (04) :571-594
[9]  
Bruinsma S., 2010, 2010 ESA LIVING PLAN
[10]   ESA's satellite-only gravity field model via the direct approach based on all GOCE data [J].
Bruinsma, Sean L. ;
Foerste, Christoph ;
Abrikosov, Oleg ;
Lemoine, Jean-Michel ;
Marty, Jean-Charles ;
Mulet, Sandrine ;
Rio, Marie-Helene ;
Bonvalot, Sylvain .
GEOPHYSICAL RESEARCH LETTERS, 2014, 41 (21) :7508-7514