解大型稀疏线性方程组的整体松弛并行多分裂法.pdfVIP

第一章引言 3 cao【12，16]等对二级多分裂也进行了深入的研究； 的牛顿一并行矩阵多分裂算法，建立了收敛性定理，估计了收敛速度．总之，关于多分裂 的研究及结果如雨后春笋般出现．本文研究的是求解线性代数方程组Az=6的整体松弛 (非定常)并行多分裂(多参数)迭代法，通过选用多个松弛因子，我们的方法覆盖了许多 文献中已见的并行多分裂迭代法，具有很强的普遍性．本文详细地比较了并行多分裂之间 的敛散速度，并在一台MPP并行机上进行多次试验．大量的数值试验表明：我们的方法 确实对文献的方法进行了有效的改进． 第二章 整体松弛并行多分裂ToR法及敛散速度的比较 §2．1整体松弛soR方法 许多大型科学与工程计算问题最终归结为求解线性代数方程组 Az=6 (21．1) 其中4∈丑”“是非奇异矩阵，￡∈』妒为待求解向量，且设∥∈，∥为(21．1)的解． 定义2．1．1令地：肌，风均是Ⅳ×Ⅳ矩阵，☆=l，2，…，Ⅱ，若满足 o)‘4=尬一^k，n行1存在． ∞∑取=』(J为Ⅳ×Ⅳ单位矩阵)，岛是非负对角矩阵． 则称三元组(慨，肌，西)(≈=l，2，．．，o)为矩阵A的一个多分裂 方程组(2．11)可以写成 z=峰1％。+螈16，≈=1，2，…，＆ 用权矩阵鼠组合这n个方程，得到 — Z 酞Z Z+ Mb 盯Z+∞ 帆 玩 一 悔 k 。∑M 。∑m 。∑M M G= M ‰ ‰ 研 —k 。∑㈦ 。∑㈦ 其迭 中代 壮黻 1 z(研+1)=．H茁《巩)+Gbm=0 称为多分裂迭代格式．口称为多分裂迭代法的迭代矩阵，其算法如下 算法1(并行多分裂方法)． 任取初始近似z(o)∈冗”． 对m=o，1，…重复下面的i)ji)两步，直到收敛 i)对☆=1，2，…，n，(并行)求解玑． i)对☆=l，2，…，n．(并行)求解鲰． A靠№=Ⅳ≈z(m)+b (2l2) 第二章整体松弛并行多分裂TOR法及敛散速度的比较 5 ii)计算 。(”+1)=F最靴． (2．1．3) k=l 我们看到，在(2．1．2)中每个玑的计算是相互独立的．因此，若有一台具有一个主机 的某些对角元素为零，则瓠的相应分量无需计算，从而大大节省工作量．由此可见＆也 起到调配各处理机上工作量的作用．我们应尽量选择B使得各处理机间大致做到负载平 衡，从而减少同步等待的开销．可见，算法1具有自然的并行性． 定义2．1．2 』4=(血“)∈．RⅣ。Ⅳ，zⅣ。Ⅳ={A∈R。v。ⅣIaiJ≤o，b乍≠J)． 若A—B≥0，则表示为A≥B； (3)若A∈zⅣ。Ⅳ且A一1≥o，称A为M矩阵； (4)若似)为M矩阵，称A为日矩阵． 引理2．1．3【18]若A，B∈月Ⅳ。Ⅳ，D=d{n9(A)． (n)若A是M矩阵，则D≥o，且D非奇异；