三对角系统并行算法的研究概况.docVIP

　　三对角系统并行算法的研究概况 【摘 要】在 科学 和工程 计算 中，许多 问题 往往归结为三对角线性方程组的求解，其并行算法的 研究 具有重要意义。文章全面 总结 了当前求解三对角线性方程组的两类并行算法：直接解法和迭代解法，并介绍了其特点。 　　 【关键词】三对角线性方程组；分治策略；并行算法；算法可扩展性 一、概述 　　三对角线性方程组的求解是许多科学和工程计算中最重要也是最基本的问题之一。在核物理、流体力学、油藏工程、石油地震数据处理及数值天气预报等许多领域的大规模科学工程和数值处理中都会遇到三对角系统的求解问题。很多三对角线性方程组的算法可以直接推广到求解块三对角及带状线性方程组。由于在 理论 和实际 应用 上的重要性，近20年来三对角方程组的并行算法研究十分活跃。 　　大规模科学计算需要高性能的并行计算机。随着软硬件技术的 发展 ，高性能的并行计算机日新月异。现今，SMP可构成每秒几十亿次运算的系统，PVP和COPP和DSM可构成每秒万亿次运算或更高的系统。 　　高性能并行计算机只是给大型科学计算提供了计算工具。如何发挥并行计算机的潜在性能和对三对角系统进行有效求解，其关键在于抓住并行计算的特点进行并行算法的研究和程序的设计与实现。另外，对处理机个数较多的并行计算系统，在设计并行算法时必须解决算法的可扩展性，并对可扩展性进行研究和 分析 。 　　二、问题的提出 　　设三对角线性方程组为 AX=Y （1） 式中：A∈Rn×n非奇异，αij=0， 。X=(x1，x2，…xn)T Y=(y1，y2，…yn）T。 　　此系统在许多算法中被提出，因此研究其高性能并行算法是很有理论和实际意义的。 　　三、并行求解三对角系统的直接解法 　　关于三对角线性方程组的直接求解已经有大量并行算法，其中3090VF向量机的特点提出了局部循环约化法；P.Amodio针对分布式系统的特点改进了循环约化方法；最近针对此方法又提出对三对角方程组进行更大约化步的交替迭代策略。 　　（三）基于矩阵乘分解算法 　　将系数矩阵A分解成A=FT，方程Ax=b化为Fy=b和Tx=y两个方程组的并行求解。这种算法又可以分为两类： 　　1.重叠分解。如ichielse和A.Van der Vorst改变ichielse和A.Van der Vorst算法中的通信模式从单向串行改为双向并行，提出DPP算法，是 目前 最好的三对角方程组分布式算法之一。2000年骆志刚等中依据DPP算法，利用计算与通信重叠技术，减少处理机空闲时间取得了更好的并行效果。此类算法要求解P-1阶缩减系统。 　　2.不重叠分解。例如Lap; Sameh算法、Johsoon算法、Baron算法、Chaan Morrison公式将和的逆化为子矩阵逆的和。按矩阵分解方法，这种算法又可分为两类： 　　1.重叠分解。这类算法首先由Mehrmann在1990年提出，通过选择好的分解在计算过程中保持原方程组系数矩阵的结构特性，具有好的数值稳定性，需要求解P-1阶缩减系统。 　　2.不重叠分解。Sun等在1992年提出的并行划分LU算法PPT算法和并行对角占优算法PDD算法均属于这一类。需要求解2P-2阶缩减系统。其中PDD算法的通讯时间不随处理机的变化而变化，具有很好的可扩展性。X.H.Sun和W.Zhang在2002年提出了两层混合并行方法PTH ，其基本思想是在PDD中嵌入一个内层三对角解法以形成一个两层的并行，基本算法是PDD，三对角系统首先基于PDD分解。PTH算法也具有很好的可扩展性。 　　四、并行求解三对角系统的迭代解法 　　当稀疏线性方程组的系数矩阵不规则时，直接法在求解过程中会带来大量非零元素，增加了 计算 量、通信量和存储量，并且直接法不易并行，不能满足求解大规模 问题 的需要。因此通常使用迭代法来求解一般系数线性方程组和含零元素较多三对角线性方程组。迭代法包括古典迭代法和Krylov子空间迭代法。 　　古典迭代法包括Jacobi、Gauss-SEidel、SOR、SSOR等 方法 。通常采用红黑排序、多色排序和多分裂等技术进行并行计算。由于古典迭代法有收敛速度慢、并行效果不好等缺点， 目前 已较少用于直接求解大型稀疏线性方程组，而是作为预条件子和其它方法（如Krylov子空间方法）相结合使用。 　　Krylov子空间方法具有存储量小，计算量小且易于并行等优点，非常适合于并行求解大型稀疏线性方程组。结合预条件子的Krylov子空间迭代法是目前并行求解大型稀疏线性方程组的最主要方法。 　　给定初值X0，求解稀疏线性方程组AX=Y。设Km为维子空间，一般投影方法是从m维仿射子空间X0+Km中寻找近似解Xm使之满足Petrov-Galerkin条件