并行求解多维递归方程组的三种Krylov 子空间迭代方法.pdfVIP

第２９卷第１１期　　　 计算机应用与软件 Ｖｏｌ２９Ｎｏ．１１ ２０１２年１１月　　 ＣｏｍｐｕｔｅｒＡｐｐｌｉｃａｔｉｏｎｓａｎｄＳｏｆｔｗａｒｅ Ｎｏｖ．２０１２ 并行求解多维递归方程组的三种Ｋｒｙｌｏｖ子空间迭代方法 李　芳　尹万旺　刘　鑫　陆林生 （江南计算技术研究所　江苏无锡２１４０８３） 摘　要　　多维递归方程组在并行求解时存在串并行不一致问题，提供三种Ｋｒｙｌｏｖ子空间迭代求解方法———ＰＣＧ／ＡＴＣＧ和ＧＭＲＥＳ 来解决这一问题，并采用典型算例对这三种Ｋｒｙｌｏｖ子空间迭代方法进行正确性验证和加速比测试。试验表明这三种Ｋｒｙｌｏｖ子空间 迭代法在并行规模较大的情况下，均能够正确求解多维递归方程组，并且加速特性良好。 关键词　　并行　多维　递归方程组　Ｋｒｙｌｏｖ子空间迭代法　ＰＣＧ　ＡＴＣＧ　ＧＭＲＥＳ 中图分类号　ＴＮ４０１　　　　文献标识码　Ａ ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１０００３８６ｘ．２０１２．１１．０２２ ＴＨＲＥＥＫＲＹＬＯＶＳＵＢＳＰＡＣＥＩＴＥＲＡＴＩＯＮＭＥＴＨＯＤＳＦＯＲＰＡＲＡＬＬＥＬＣＯＭＰＵＴＡＴＩＯＮ ＯＦＭＵＬＴＩＤＩＭＥＮＳＩＯＮＲＥＣＵＲＳＩＶＥＥＱＵＡＴＩＯＮＳ ＬｉＦａｎｇ　ＹｉｎＷａｎｗａｎｇ　ＬｉｕＸｉｎ　ＬｕＬｉｎｓｈｅｎｇ （ＪｉａｎｇｎａｎＩｎｓｔｉｔｕｔｅｏｆＣｏｍｐｕｔｉｎｇＴｅｃｈｎｏｌｏｇｙ，Ｗｕｘｉ２１４０８３，Ｊｉａｎｇｓｕ，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ　　Ｔｈｅｐａｒａｌｌｅｌｓｏｌｖｉｎｇｏｆｍｕｌｔｉｄｉｍｅｎｓｉｏｎｒｅｃｕｒｓｉｖｅｅｑｕａｔｉｏｎｓｈａｓｉｎｃｏｎｆｏｒｍｉｔｙｂｅｔｗｅｅｎｔｈｅｓｅｒｉａｌｃｏｍｐｕｔａｔｉｏｎａｎｄｔｈｅｐａｒａｌｌｅｌ ｃｏｍｐｕｔａｔｉｏｎ．Ｔｏｓｏｌｖｅｔｈｉｓｐｒｏｂｌｅｍ，ｉｎｔｈｅｄｉｓｓｅｒｔａｔｉｏｎｗｅｐｒｅｓｅｎｔｔｈｒｅｅＫｒｙｌｏｖｓｕｂｓｐａｃｅｉｔｅｒａｔｉｏｎｍｅｔｈｏｄｓ，ＰＣＧ，ＡＴＣＧａｎｄＧＭＲＥＳ．Ｓｏｍｅ ｔｙｐｉｃａｌｎｕｍｅｒｉｃａｌｅｘａｍｐｌｅｓａｒｅｕｓｅｄｔｏｖｅｒｉｆｙｔｈｅｃｏｒｒｅｃｔｎｅｓｓａｎｄｃｏｎｄｕｃｔｔｈｅａｃｃｅｌｅｒａｔｉｏｎｒａｔｉｏｔｅｓｔｆｏｒｔｈｅｓｅｔｈｒｅｅｍｅｔｈｏｄｓ．Ｔｅｓｔｒｅｓｕｌｔｓｓｈｏｗ ｔｈａｔｔｈｅｓｅｔｈｒｅｅＫｒｙｌｏｖｓｕｂｓｐａｃｅｉｔｅｒａｔｉｏｎｍｅｔｈｏｄｓａｒｅａｂｌｅｔｏｃｏｒｒｅｃｔｌｙｒｅｓｏｌｖｅｔｈｅｍｕｌｔｉｄｉｍｅｎｓｉｏｎｒｅｃｕｒｓｉｖｅｅｑｕａｔｉｏｎｓｗｉｔｈｓａｔｉｓｆｉｅｄ ａｃｃｅｌｅｒａｔｉｏｎｐｒｏｐｅｒｔｙｅｖｅｎｉｎｔｈｅｃｏｎｄｉｔｉｏｎｏｆｂｉｇｇｅｒｐａｒａｌｌｅｌｓｃａｌｅ． Ｋｅｙｗｏｒｄｓ　　Ｐａｒａｌｌｅｌ　Ｍｕｌｔｉｄｉｍｅｎｓｉｏｎ　Ｒｅｃｕｒｓｉｖｅｅｑｕａｔｉｏｎｓ　Ｋｒｙｌｏｖｓｕｂｓｐａｃｅｉｔｅｒａｔｉｏｎｍｅｔｈｏｄ　ＰＣＧ　ＡＴＣＧ　ＧＭＲＥＳ 算及通信过程，往往可以获得相当好的并行加速比，然而对于块 ０　引　言 数较多的多维递归方程组，尽管可以很好地解决串并行一致性 问题，但要达到高效的流水线并行尤为困难。 在科学与工程计算领域，如计算流体力学（ＣＦＤ）、结构与非 并行迭代法的每个迭代步均通信，从而使块边界面上的值 结构问题的有限元分析、石油地震数据处理、数值天气预报等， 不断更新，这样通过多次迭代的方法将相邻区域的信息逐渐传 数值模拟均是将模拟方程离散成线性代数方程组求解，很多复 过来，并逐步逼近串行强隐式的计算结果。Ｋｒｙｌｏｖ子空间迭代 杂问题都归结为求解多维递归方程组。随着科学技术的发展， 法是当前求解大型稀疏矩阵线性代数方程组的有效方法，已经 求解问题的规模越来越大，复杂性越来越高，应用问题的多块化