大规模矩阵的MPI并行求逆算法设计与分析.pdfVIP

第34卷第5期 大 地 测 量 与 地 球 动 力 学 Vo1．34 Nc】．5 2014年 10月 J0URNALOFGEODESYANDGEODYNAMICS ()(·1．．20I4 文章编号 ：1671—5942(2014)05—0120-05 大规模矩阵的MPI并行求逆算法设计与分析 周 浩 罗志才 钟 波 ， ，1)武汉大学测绘学院，武汉 430079 、 ＼2)地球空间环境与大地测量教育部蓖点实验室，武汉 430079／ 摘 要 针对最小二乘解算中的大规模矩阵求逆问题，基于MPI实现了高阶稠密对称 定矩 的( 一Jtrdm 并行求逆算法，减少了计算耗时；通过优化矩阵读写、存储等方式降低了单个计算节点的内存耗用量，拓展了算法 的可移植性。通过并行读写效率、单节点进程数、加速比和相对效率等凶素评价算法的计箅效率，计箅结果表明： 通过引入并行读写、减小单个计算节点的负荷等方式，并行求逆的相对效率峰值可达60％。以卫星嚷力场反浈 为 例，采用曙光集群 的8个计算节点分别恢复截断阶次为 120、240的地球重力场模 ，求逆耗时为229S、7395s， 单个节点的内存耗用峰值为205MB、1．57GB，反演精度可达 10 量级 ，表明该算法能够快速稳定地获取最小二乘 问题的最优估值。 关键词 MPI；对称正定矩阵；矩阵求逆；Gauss—Jordan算法；最小二乘法 中图分类号 ：P209 文献标识码 ：A DESIGN AND ANALYSIS OF MPIPARALLEL ALGoRITHM FoR LARGE．SCALE MATRIX INVERSIoN ZhouHao ．LuoZhieai’ andZhongBo， ／1)Schoolo厂GeodessandGeomatics，WuhanUniversity，Wuhan 430079 ＼ ＼2)KeyLabofGeospaceEnvironmentandGeodesy，Ministryo／’Education．z￡rznUniversity． ￡ “， 430079J Abstract Inordertosolvetheprobleminlargescalemau-ixinversionwithleastslaresolvingtl1e}】aIa·11 nI— ， gorithm toinversethehigh—orderdensesymmetricalpositivedefinematrixwithGauss—Jordanmethodonthebasisof MP1wasproposed．Consideringthehighmemorydemandonsingleprocessorforthisalgorithm ，theMPII／0 inler— facewasintrodueed，which can reducecommunicationaltimebetween each nodesimultane(}usly ． 1ntrodu(．i