一种求解鞍点问题的预处理并行算法-应用数学和力学.PDFVIP

　 ， ３５ ９ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年９月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．９，Ｓｅｐ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）０９⁃１０１１⁃０９ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 一种求解鞍点问题的预处理并行算法∗ １ １ ２ 姜晓林 吕全义 谢公南 ，　 ，　 （１． 西北工业大学 应用数学系，西安 ７１０１２９； ２． 西北工业大学 机电学院；工程仿真与宇航计算技术联合实验室，西安 ７１００７２） （我刊编委谢公南来稿） 摘要：　 研究了一种求解鞍点问题的并行预处理变形共轭梯度算法．通过应用迭代法进行预处理 后，再采用变形共轭梯度求解的模式．首先构造系数矩阵近似逆的多项式表达式，以此作为预处理 矩阵的逆矩阵，对方程组进行预处理；然后采用变形共轭梯度法并行求解预处理后的线性方程组． 为减少运算量，采用迭代方式并行计算多项式与向量的乘法运算．通过调整迭代次数，即调整多项 式次数，检验各种次数的多项式进行预处理后的求解方程的效果．数值试验结果表明，该算法明显 优于未预处理的变形共轭梯度法，且当预处理迭代次数取４时效果最好． 关　 键　 词：　 鞍点问题；　 并行算法；　 变形共轭梯度法；　 预处理方法 中图分类号：　 Ｏ２４６　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．０９．００７ 引　 　 言 鞍点问题是一类非常重要的线性方程组的求解问题，常见于科学计算以及工程应用领域， 如Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程的离散、第二类椭圆形方程的混合离散、计算流体力学问题、电磁学问题、 最小二乘问题、约束优化问题和非线性规划问题等．它的应用已经渗透到应用数学、力学、电学 等众多领域．对于鞍点问题，虽然串行算法已经比较完善，但是随着计算规模的增大，求解鞍点 问题的存储要求和计算量快速增加，单台处理机往往难以承受．目前并行算法被认为是解决大 规模鞍点问题最有效的途径之一，因此研究求解鞍点问题的高效并行算法具有重要的理论意 义和实际意义． 一般来说，鞍点问题的求解方法主要分为两类：基于矩阵分裂的定常迭代法和基于Ｋｒｙｌｏｖ ［１］ 子空间的不定常迭代法．对于定常迭代法，其中最为著名的就是Ｕｚａｗａ⁃ｔｙｐｅ迭代法 和ＨＳＳ类 ［２］ ［３］ 迭代法 ．对于不定常迭代法，Ｋｒｙｌｏｖ子空间法 是求解鞍点问题行之有效的方法．为了提高 Ｋｒｙｌｏｖ子空间迭代法的收敛速度，往往需要对系数矩阵进行预处理．国内外很多学者做了大量 的工作，并给出了许多行之有效的预条件子，主要有块对角预条件子、约束预条件子、ＨＳＳ预条 件子等．共轭梯度法（ｃｏｎｊｕｇａｔｅ ｇｒａｄｉｅｎｔ ｍｅｔｈｏｄ，简称为ＣＧ法）作为Ｋｒｙｌｏｖ 子空间法的一种基 本方法更使得到了广泛的应用．ＣＧ法具有存储量少、计算量少和适合并行计算等优点．但是它 ∗ 收稿日期：　 ２０１４⁃０４⁃２５；修订日期：　 ２０１４⁃０６⁃３０ 基金项目：　 陕西省自然科学基金（２００９ＪＭ１００８）；国家自然科学基金 作者简介：　 姜晓林（１９８９—），男，山东人，硕士生（Ｅ⁃ｍａｉｌ：ｊｉａｎｇｘｉａｏｌｉｎ１２３４＠１２６．ｃｏｍ）； 吕全义（１９６３—），女，沈阳人，副教授（通讯作者． Ｅ⁃ｍａｉｌ：ｌｕｑｕａｎ＠ｎｗｐｕ．ｅｄｕ．ｃｎ）． １０１１