×

GMRES完全停滞。 (英语) Zbl 1025.65022号

在他们的引言中,作者指出,“我们研究了一个奇怪的问题:广义最小残差(GMRES)算法在以初始猜测(x^{(0)}=0)开始并使用精确算法时,计算(m)迭代(x^}=cdots=x^{(m)}=0\)而没有任何进展部分m级停滞。如果\(m=n-1\),我们称之为完全停滞GMRES。在这种情况下,GMRES将在迭代时计算精确解。”
假设待解矩阵方程为(Ax=b),矩阵(A)可对角化且具有谱分解(A=V\Lambda V^{-1}。写入\(y=V^{-1}b\)并用(y)表示其条目为\(y)分量的对角矩阵。进一步,用(Z_m)表示由(A)的第一个(m)特征值构造的Vandermonde矩阵。最后,表示\(e_1=[1,0,\ldots,0]^T\)。这个停滞,停滞方程式\[Z^H_{m+1}\bar{Y}(V^HV)Y=e_1\]被证明完全表征了(m)-停滞。该证明基于Krylov矩阵的一个众所周知的因式分解。
在(A)的完全停滞和不同特征值的情况下,停滞方程成为两个矩阵(G(lambda)=Z_n)的等式^{-H}e_1\)和(F_V(y)=\bar{y}(V^HV)y\)。这种特征值和特征向量影响的分离允许对停滞方程进行几何解释。
作者继续讨论了几个特殊情况:(1)正规矩阵,如果(G(y))有复实项或负实项,则完全停滞是不可能的;(2) 不可能完全停滞的厄米矩阵和实对称矩阵;当且仅当特征值的相位角均匀分布在单位球面上时,酉矩阵可能停滞。一些讨论给出了已知结果的新证明,而一些结果本身是新的。一些讨论集中于从理论结果中导出的特定数值示例。

MSC公司:

65层10 线性系统的迭代数值方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Brown,P.N.,《Arnoldi和GMRES算法的理论比较》,SIAM J.Sci。统计计算。,12, 58-78 (1991) ·Zbl 0719.65022号
[2] de Sturler,E.,最优Krylov子空间方法的截断策略,SIAM J.Numer。分析。,36, 864-889 (1999) ·Zbl 0960.65031号
[4] 欧·恩斯特。;艾尔曼,M.,《克里洛夫子空间方法理论的几何方面》,《数值学报》。,10, 251-312 (2001) ·Zbl 1105.65328号
[6] 格林鲍姆,A。;普塔克,V。;Strakoš,Z.,对于GMRES,任何非增量收敛曲线都是可能的,SIAM J.Matrix Anal。申请。,17, 465-469 (1996) ·兹比尔0857.65029
[8] Higham,N.,《数值算法的准确性和稳定性》(1996),SIAM:SIAM Philadelphia·Zbl 0847.65010号
[9] Ipsen,I.C.F.,GMRES残差的表达式和界限,BIT,40,524-533(2000)·Zbl 0962.65030号
[10] 兰卡斯特,P。;Tismenetsky,M.,《矩阵理论》(1985),学术出版社:纽约学术出版社·Zbl 0516.15018号
[11] 马歇尔,A.W。;Olkin,I.,矩阵缩放以实现指定的行和列总和,Numer。数学。,12, 83-90 (1968) ·Zbl 0165.17401号
[12] Morgan,A.P.,《利用工程和科学问题的连续性解决多项式系统》(1987),普伦蒂斯·霍尔:普伦蒂斯霍尔·恩格尔伍德·克利夫斯,新泽西州·Zbl 0733.65031号
[13] 摩根会计师事务所。;Sommese,A.J.,系数参数多项式延拓,应用。数学。计算。,29, 123-160 (1989) ·Zbl 0664.65049号
[14] 新墨西哥州纳希蒂加尔。;南卡罗来纳州雷迪。;Trefethen,L.N.,《非对称矩阵迭代有多快》,SIAM J.《矩阵分析》。申请。,13, 778-795 (1992) ·Zbl 0754.65036号
[15] 佩奇,C.C。;桑德斯,M.A.,稀疏不定线性方程组的求解,SIAM J.Numer。分析。,12, 617-629 (1975) ·Zbl 0319.65025号
[16] 萨阿德,Y。;Shultz,M.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 856-869 (1986) ·Zbl 0599.65018号
[17] 萨阿德,Y。;Wu,K.,DQGMRES:基于不完全正交化的准最小残差算法,Numer。线性代数应用。,3, 329-343 (1996) ·Zbl 0906.65033号
[18] Simoncini,V.,重新启动的GMRES的新变体,线性代数应用。,第6页,第61-77页(1999年)·Zbl 0983.65033号
[19] Sosonkina,M。;Watson,L.T。;Walker,H.F。;Kapania,R.K.,一种新的自适应GMRES算法,用于实现高精度,Numer。线性代数应用。,5, 275-297 (1998) ·Zbl 0937.65037号
[20] 斯图尔特,G.W。;Sun,J.-G.,矩阵微扰理论(1990),学术出版社:纽约学术出版社·Zbl 0706.65013号
[21] Vorst,H.V.D。;Vuik,C.,GMRESR:一系列嵌套的GMRES方法,Numer。线性代数应用。,1, 369-386 (1994) ·兹比尔083965040
[22] 怀斯,S.M。;Sommese,A.J。;Watson,L.T.,算法801:POLSYS_PLP:用于求解多项式方程组的分区线性乘积同伦代码,ACM Trans。数学。软件,26176-200(2000)·Zbl 1137.65353号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。