对称逐步超松驰预处理共轭梯度法(SSOR-PCG)是一种求解大型稀疏对称正定线性方程组的非常有效的迭代法。SSOR-PCG并行化的难点在于每步迭代都要求解2个三角方程组。采用一种改进的SSOR-PCG并行求解有限元方程组,并采用多色排序技术提高并行度。基于MPI模型开发了并行程序,通过测试,选择了有效的MPI通信函数。
Abstract
The method of symmetric successive over relaxation-preconditioned conjugate gradient (SSOR-PCG) is a very effective iterative method for solving large scale sparse symmetric positive-definite linear set of equations. The difficulty in the parallelization of the SSOR-PCG lies in solving two triangular equation systems in each iteration. In this research,an improved SSOR-PCG is applied to parallel solve finite element equations,and the multicolor ordering technique is used to increase the degree of parallelism. A MPI-based parallel program is coded and choice is made for efficient MPI communication routines by tests.
关键词
SSOR-PCG /
并行计算 /
多色排序 /
有限元方程组 /
MPI
Key words
SSOR-PCG /
parallel computing /
multicolor ordering /
finite element equations /
MPI (Message Passing Interface)
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}
参考文献
[1] 吕 涛,石济民,林振宝.区域分解算法-偏微分方程数值解法新技术[M].北京:科学出版社,1994.(LV Tao,SHIH T M,LIEM C B. Domain Decomposition Methods: New Numerical Techniques for Solving PDE[M]. Beijing: Science Press,1994.(in Chinese))
[2] ADAMS L,ORTEGA J M. A Multi-color SOR Method for Parallel Computation[C]∥International Association of Computing and Communication (IACC). Proceedings of the 1982 International Conference on Parallel Processing. Bellaire,Michigan,August 24-27,1982: 53-58.
[3] 张汝清.并行计算结构力学[M].重庆:重庆大学出版社,1993.(ZHANG Ru-qing. Parallel Computing in Structural Mechanics[M]. Chongqing: Chongqing University Press,1993.(in Chinese))
[4] 吴建平,王正华,李晓梅.稀疏线性方程组的高效求解与并行计算[M].长沙:湖南科学技术出版社,2004.(WU Jian-ping,WANG Zheng-hua,LI Xiao-mei. Efficient Solution and Parallel Computing of Sparse Linear Equations[M]. Changsha: Hunan Science & Technology Press,2004.(in Chinese))
[5] POOLE E L ,ORTEGA J M. Multicolor ICCG Methods for Vector Computers[J]. SIAM Journal on Numerical Analysis,1987,24(6): 1394-1418.
[6] 林绍忠.对称逐步超松驰预处理共扼梯度法的改进迭代格式[J].数值计算与计算机应用,1997,(4):266-270.(LIN Shao-zhong. Improved Iterative Format of Symmetric Successive over Relaxation - Preconditioned Conjugated Gradient Method[J]. Journal of Numerical Methods and Computer Applications,1997,(4):266-270.(in Chinese))
[7] 林绍忠,苏海东.大体积混凝土结构仿真应力分析快速算法及应用[J].长江科学院院报,2003,20(6):19-22.(LIN Shao-zhong,SU Hai-dong. Fast Algorithms for Stress Analysis Simulating Construction Process of Massive Concrete Structures and Applications[J]. Journal of Yangtze River Scientific Research Institute,2003,20(6):19-22.(in Chinese))
[8] 程 昭,陈胜宏,田胜利,等.大规模P型有限元方程组的修正SSOR-PCG解法[J].上海交通大学学报,2002,36(11):1608-1611.(CHENG Zhao,CHEN Sheng-hong,TIAN Sheng-li,et al. Modified SSOR-PCG Solver for Large Scale P-Version FEM Equations[J]. Journal of Shanghai Jiaotong University,2002,36(11): 1608-1611.(in Chinese))
[9] 李 静.陈健云,周 晶.改进的SSOR-PCG迭代法在接触问题研究中应用[J].大连理工大学学报,2006,46(4):533-537.(LI Jing,CHEN Jian-yun,ZHOU Jing. Application of an Improved SSOR-PCG Iteration Method to Contact Problems[J]. Journal of Dalian University of Technology,2006,46(4): 533-537.(in Chinese))
[10]陈国荣,李皇胜,李红健.大体积混凝土温度场高速求解方法的改进[J].河海大学学报,2009,37(4):396-399.(CHEN Guo-rong,LI Huang-sheng,LI Hong-jian. Improvement of Fast Solution Method for Mass Concrete Temperature Field[J]. Journal of Hohai University (Natural Sciences,2009,37(4): 396-399.(in Chinese))
[11]张 杨.高阶三维数值计算方法理论与应用研究[D].南京:河海大学,2009.(ZHANG Yang. Theory and Application for High Order Three Dimensional Numerical Method[D]. Nanjing: Hohai University,2009.(in Chinese))
[12]李 根.基于模拟的水岩耦合变形破坏过程及机理研究[D].大连:大连理工大学,2011.(LI Gen. Simulation-Based Research on Rock Deformation and Failure Process under Water-Rock Coupling[D]. Dalian: Dalian University of Technology,2011.(in Chinese))
[13]都志辉.高性能计算并行编程技术:MPI并行程序设计[M].北京: 清华大学出版社,2001.(DU Zhi-hui. Parallel Programming Techniques in High Performance Computing: MPI Parallel Programming[M]. Beijing: Tsinghua University Press,2001.(in Chinese))
[14]孙惠泉.图论及其应用[M].北京:科学出版社,2004.(SUN Hui-quan. Graph Theory and Its Application[M] Beijing: Science Press,2004.(in Chinese))
基金
中央级公益性科研院所基本科研业务费项目(CKSF2011016)