基于加权GMRES的多项式预处理广义极小残差法.pdfVIP

第31卷 第5期 工 程 数 学 学 报 v01．31No．5 2O14年10月 CHINESEJOURNALOFENGINEERING MATHEMATICS Oct．2014 doi：10．3969／j．issn．1005-3085．2014．05．007 文章编号：1005—3085(2014)05—0697-10 基于加权GMRES的多项式预处理广义极小残差法木 关朋燕 ，李春光 ， 景何仿 f1．北方民族大学信息与计算科学学院，银川 750021； 2一北方民族大学数值计算与工程应用研究所，银川 750021) 摘 要：Essai对求解非对称线性方程组 的广义极小残差法fGMRES)进行改进 ，提 出加 权GMRES方法fWGMRES1．该方法具有 良好 的收敛效果，但计算量较大 ．本文 利用加权GMRES本身构造出一种有效的多项式预处理并与加权GMRES结合，构造 了一种新算法．该算法能够显著地减少加权GMRES迭代次数，并能减少运算量和储 存量．数值算例表明，相对于加权GMRES方法来说，新算法减少了运算时间和迭代 次数． 关键词：多项式预处理；加权Arnoldi算法；加权GMRES算法；迭代法；平行板突扩管 分类号：AMS(2000)65F10；65F50 中图分类号：O241．6 文献标识码：A 1 引言 许多科学计算问题经过离散可归结为解线性方程组 Ax： b， (1) 其中A ∈CⅣ Ⅳ非奇异，b是列向量．解这类方程组的一类方法就是基于最小残差法 X =Xo+qn-1(A)ro， (2) 其中％一l(x)为多项式且deg‰一 n一1，使得IIII最小，或 x ∈Xo+％ (ro，)， r上 (r0，A)， r=b—Ax表示残差， (，A)=span{r0，Ar0，… ，一 r0)是Krylov子空间． 令P =b— ‰一()，此为残差多项式，于是r =P ()r0．在这类方法中应用最广 的是GMRES方法 [1]，GMRES是在Krylov子空间上进行迭代，该算法的收敛性 由残差 模的单调性得出，而且具有 良好的数值稳定性，但具有较长的迭代公式，因此每步迭代 的计算量与储存量都线性增长，同时难 以有效的并行运算 ．Saad采用循环 GMRES法， 但GMRES(m)2【_会失去超线性收敛 ，产生停滞 ．1998年，Essai对GMRES方法进行 收稿 日期：2012—12—31．作者简介：关朋燕 (1985年5月生)，女，硕士．研究方向：数值代数 基金项目：国家自然科学基金． t通讯作者：李春光 E—mail：cglizd@hotmail．com 698 工 程 数 学 学 报 第3l卷 改进 ，提出加权GMRES方法 3【】，其优 点在于具有 良好 的收敛性 ，但也增加 了计算 量．Nachtigal等 _【l1】讨论了GMRES和其他方法相结合策略．文献[12～17]是将加权 思想与其他方法相结合文章．本文采用多项式预处理加权GMRES(m)法及多项式预 处理加权GMRES(m)法．其主要问题是找到有效的低阶多项式P()，迭代解应用 到 (A)Ax=Pn()6中，这会产生一个Krylov子空间 K (p(A)A；ro)=span{ro，Pn(A)Aro，…， (A)A)m--1to}． 它是Krylov子空间 (+)(一)+(；r0)的子空间．本文直接利用加权GMRES方法本身 来构造预处理多项式，多项式的预处理和混合迭代[]极为相似，其思想就是尽可能减少迭 代次数．随着迭代次数的减少，存储量和运算量也会相应减少．