×

移位线性系统的一种更简单的GMRES及其自适应变体。 (英语) Zbl 1413.65063号

作者考虑了两个线性系统的解(Ax=b)和(Ax=b),其中种子矩阵(A)和加法矩阵(A=A-σI)是非奇异的。众所周知,由移位矩阵(A)生成的Krylov空间(K_K(A,b))与由(A)产生的K_K空间(b)相同,前提是我们从相同的向量(b)开始。通过生成单个搜索空间,可以获得两个系统的近似解。然而,如果在GMRES方法中使用重新启动,则种子系统的剩余向量\(rm\)和加法系统的剩余矢量\。要求\(\hat r_m\)和\(r_m\)的共线,A.弗罗默U.Glässner公司[SIAM J.Sci.Compute.19,No.1,15-26(1998;Zbl 0912.65023号)]已经为移位线性系统开发了重新启动的GMRES方法的变体,该方法允许一个人以仅生成一个搜索空间为代价同时求解多个线性系统(此处表示为GMRES-Sh)。作者在本文中提出了一种基于由H.F.沃克L.Zhou(L.周)【数值线性代数应用1,第6期,571-581(1994;Zbl 0838.65030号)](生成的方法称为SGMRES-Sh)。然而,较简单的GMRES的主要问题是,由于计算基础的不稳定,它本质上是不稳定的。因此,作者提出了一种基于Krylov子空间基自适应选择的SGMRES-Sh稳定变种,使其条件数保持在合理水平。在数值实验部分,他们给出了一个示例,说明自适应SGMRES-Sh方法比普通SGMRES-Sh方法表现得好得多。通过对GMRES-Sh、SGMRES-Sh和自适应SGMRES-Sh的矩阵向量乘积数量和所用时间进行数值比较,证明了该方法的有效性。

MSC公司:

65层10 线性系统的迭代数值方法
65克50 舍入误差
65英尺35英寸 矩阵范数、条件、缩放的数值计算
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 大型Sylvester-like观测器矩阵方程的Datta,Arnoldi方法和部分谱分配的相关算法,线性代数应用154-156 pp 447–(1991)·Zbl 0734.65037号 ·doi:10.1016/0024-3795(91)90378-A
[2] Sweet,循环约简算法的并行和矢量变体,SIAM J Stat Sci Compute 9 pp 89–(1988)·Zbl 0651.65019号 ·doi:10.1137/0909050
[3] Darnell,晶格QCD迭代方法的特征值收缩,Nucl Phys B Proc Suppl 129 pp 856–(2004)·doi:10.1016/S0920-5632(03)02734-8
[4] Simoncini,线性系统Krylov子空间方法的最新计算发展,数字线性代数应用14第1页–(2007)·Zbl 1199.65112号 ·doi:10.1002/nla.499
[5] Walker,《一个更简单的GMRES》,数字线性代数应用1第571页–(1994年)·Zbl 0838.65030号 ·doi:10.1002/nla.1680010605
[6] Saad,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J Stat Sci Compute 7 pp 856–(1986)·Zbl 0599.65018号 ·doi:10.1137/0907058
[7] 稀疏线性系统的Saad Y迭代方法2003
[8] Boojhawon,用调和Ritz向量扩充的重新启动更简单的GMRES,Future Gener Compute Syst 20 pp 389–(2004)·Zbl 1062.65039号 ·doi:10.1016/j.future.2003.07.004
[9] Lin,Simpler GMRES with deflated restarting,数学计算模拟82 pp 2238–(2012)·Zbl 1253.65044号 ·doi:10.1016/j.matcom.2012.05.019
[10] Liu,具有多个右侧的非对称系统的Simpler block GMRES,《电子传递数值分析》30第1页–(2008)·Zbl 1188.65034号
[11] Zhong,《具有多个右侧线性系统的带放气重启的灵活和自适应简单块GMRES》,《计算应用数学杂志》282第139页–(2015)·Zbl 1309.65038号 ·doi:10.1016/j.cam.2014.12.040
[12] Carpentieri,解非对称线性系统的BiCOR和CORS迭代算法,SIAM J Sci Compute 33 pp 3020–(2011)·Zbl 1251.65045号 ·数字对象标识代码:10.1137/100794031
[13] Gu,移位线性系统族的BiCR型方法,计算数学应用68第746页–(2014)·Zbl 1362.65039号 ·doi:10.1016/j.camwa.2014.07.029
[14] Jing,关于求解线性系统的新迭代方法和比较结果,J Comput Appl Math 220 pp 74–(2008)·Zbl 1159.65037号 ·doi:10.1016/j.cam.2007.07.035
[15] Jing,移位线性系统的重新启动加权全正交化方法,计算数学应用57 pp 1583–(2009)·Zbl 1186.65061号 ·doi:10.1016/j.camwa.2008.10.088
[16] Jing,复杂非对称线性系统COCR方法的Lanczos型变体,《计算物理杂志》228,第6376页–(2009)·Zbl 1173.65316号 ·doi:10.1016/j.jcp.2009.05.022
[17] Sogabe,COCR方法对求解具有复杂对称矩阵的移位线性系统的扩展,东亚应用数学杂志第1期第97页–(2011)·Zbl 1287.65023号 ·doi:10.4208/eajam.260410.240510a
[18] Sogabe,带复对称矩阵的广义移位线性系统的解,《计算物理杂志》231 pp 5669–(2012)·Zbl 1277.65021号 ·doi:10.1016/j.jcp.2012.04.046
[19] Imakura,解移位线性系统的重新启动移位GMRES方法的有效变体,《数学研究应用33》第127页–(2013)·兹比尔1289.65057
[20] Frommer,《移位线性系统的重新启动GMRES》,SIAM科学计算杂志,19页,第15页–(1998年)·Zbl 0912.65023号 ·doi:10.1137/S1064827596304563
[21] Jiránek,简单GMRES的自适应版本,《数值算法》53第93页–(2010)·Zbl 1188.65033号 ·doi:10.1007/s11075-009-9311-2
[22] Jiránek,《如何使更简单的GMRES和GCR更稳定》,SIAM科学计算杂志30页1483–(2008)·兹比尔1176.65036
[23] Liesen,最小二乘残差和最小残差法,SIAM科学计算杂志23页1503–(2002)·Zbl 1012.65037号 ·doi:10.1137/S1064827500377988
[24] Li,求解多右手边移位系统的一种新种子投影方法,国际数学计算物理量子工程杂志7(4),第386页–(2013)
[25] Agullo,Block GMRES方法,具有不精确的故障和放气重启,SIAM J Matrix Ana Appl 35(4)第1625页–(2014)·Zbl 1355.65051号 ·数字对象标识代码:10.1137/140961912
[26] Morgan,Restarted block-GMRES with deflacting of attributive values,《应用数学》15 pp 222–(2005)·Zbl 1074.65043号 ·doi:10.1016/j.apnum.2004.09.028
[27] Davis,佛罗里达大学稀疏矩阵集合,ACM数学软件翻译38页1:1–(2011)·Zbl 1365.65123号 ·doi:10.1145/2049662.2049663
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。