×

广义极分解的稳定高效计算。 (英语) Zbl 1492.65116号

总结:我们提出了基于动态加权哈雷迭代计算矩阵广义极分解的方法。该方法在计算标准极分解时得到了很好的确定。一种稳定的实现是可用的,其中避免了矩阵求逆,而是使用QR分解。我们建立了该方法的自然推广,用于计算签名矩阵的广义极分解。同样,可以通过使用称为双曲线QR分解的广义QR分解来避免逆问题。然而,这种分解并没有表现出与正交分解相同的良好稳定性。我们通过推广CholeskyQR2方法克服了数值上的困难。该方法通过两个连续的Cholesky分解以稳定的方式计算标准QR分解。通过使用置换图基,可以获得更好的数值稳定性,即使对于条件恶劣的矩阵,也可以得到阶数为(10^{-14})的残差,而其他方法都失败了。

MSC公司:

65英尺60英寸 矩阵指数和相似矩阵函数的数值计算
65层55 低秩矩阵近似的数值方法;矩阵压缩
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] E.Anderson,不连续平面旋转和对称特征值问题,LAPACK工作注释150,2000,http://www.netlib.org/lapack/lawnspdf/lawn150.pdf。
[2] C.Ashcraft、R.G.Grimes和J.G.Lewis,精确对称不定线性方程求解器,SIAM J.矩阵分析。申请。,20(1999),第513-561页,https://doi.org/10.1137/S0895479896921。 ·Zbl 0923.65010号
[3] Z.Bai和J.W.Demmel,《并行非对称特征例程工具箱的设计》,第一部分,技术报告UCB/CSD-92-718,加州大学伯克利分校EECS系,1993年,http://www2.eecs.berkeley.edu/Pubs/TechRpts/1993/6014.html。
[4] P.Benner,哈密顿矩阵的辛平衡,SIAM J.Sci。计算。,22(2001),第1885-1904页,https://doi.org/10.1137/S1064827500367993。 ·Zbl 0985.65025号
[5] P.Benner、V.Khoromskia和B.N.Khoroskij,使用低秩张量因子分解计算Bethe-Salpeter激发能的简化基方法,摩尔物理学。,114(2016),第1148-1161页,https://doi.org/10.1080/00268976.2016.1149241。
[6] P.Benner和C.Penke,求解晶体系统Bethe-Salpeter特征值问题的高效准确算法,J.Compute。申请。数学。,400 (2022), 113650. ·兹比尔1473.65044
[7] P.Benner和C.Penke,GR分解及其与类Cholesky因子分解的关系,Proc。申请。数学。机械。,20(2021),e20200065,https://doi.org/10.1002/pamm.202000065。
[8] P.Benner和E.S.Quintana Ortí,用矩阵符号函数求解稳定的广义Lyapunov方程,Numer。算法,20(1999),第75-100页,https://doi.org/10.1023/A:1019191431273。 ·Zbl 0940.65035号
[9] X.Blase、I.Duchemin、D.Jacquemin和P.-F.Loos,《Bethe-Salpeter方程形式主义:从物理到化学》,J.Phys。化学。莱特。,11(2020年),第7371-7382页,https://doi.org/10.1021/acs.jpclett.0c01875。
[10] A.Bojanczyk、N.J.Higham和H.Patel,通过双曲QR因式分解求解不定最小二乘问题,SIAM J.矩阵分析。申请。,24(2003),第914-931页,https://doi.org/10.1137/S0895479802401497。 ·Zbl 1036.65035号
[11] Y.Bolshakov和B.Reichstein,不定标量积中的幺正等价:奇异值分解的类似物,线性代数应用。,222(1995),第155-226页,https://doi.org/10.1016/0024-3795(93)00295-B·Zbl 0828.15009号
[12] Y.Bolshakov、C.V.M.van der Mee、A.C.M.Ran、B.Reichstein和L.Rodman,有限维不定标量积空间中的极分解:一般理论,线性代数应用。,261(1997),第91-141页,https://doi.org/10.1016/S0024-3795(96)00317-5. ·Zbl 0881.15013号
[13] I.Borg和P.J.F.Groenen,《现代多维标度:理论和应用》,Springer-Verlag出版社,柏林,2005年,https://doi.org/10.1007/0-387-28981-X。 ·Zbl 1085.62079号
[14] J.R.Bunch和L.Kaufman,计算惯性和求解对称线性系统的一些稳定方法,数学。公司。,31(1977年),第163-179页,https://doi.org/10.1090/S0025-5718-1977-0428694-0。 ·Zbl 0355.65023号
[15] J.R.Bunch、L.Kaufman和B.Parlett,对称矩阵的分解,数值。数学。,27(1976年),第95-109页,https://doi.org/10.1007/BF01399088。 ·兹比尔0342.65026
[16] J.R.Bunch和B.N.Parlett,求解线性方程组对称不定方程组的直接方法,SIAM J.Numer。分析。,8(1971年),第639-655页,https://doi.org/10.1137/0708060。 ·Zbl 0199.49802号
[17] W.Bunse和A.Bunse-Gerstner,《数值线性代数》,斯图加特特伯纳出版社,1985年·Zbl 0547.65018号
[18] R.Byers,用矩阵符号函数求解代数Riccati方程,线性代数应用。,85(1987),第267-279页·Zbl 0611.65027号
[19] R.Byers和H.Xu,牛顿迭代极分解及其后向稳定性的新标度,SIAM J.矩阵分析。申请。,30(2008),第822-843页,https://doi.org/10.1137/070699995。 ·Zbl 1170.65019号
[20] M.Casida,分子的时间依赖密度函数响应理论,《密度函数方法的最新进展》,《世界科学》,新泽西州恩格尔伍德克利夫斯,1995年,第155-192页,https://doi.org/10.1142/9789812830586_0005。
[21] S.Chandrasekaran和A.H.Sayed,稳定广义Schur算法,SIAM J.矩阵分析。申请。,17(1996),第950-983页,https://doi.org/10.1137/S0895479895287419。 ·Zbl 0860.65016号
[22] J.J.Dongarra、J.R.Gabriel、D.D.Koelling和J.H.Wilkinson,具有时间反转对称性的厄米特矩阵的特征值问题,线性代数应用。,60(1984年),第27-42页,https://doi.org/10.1016/0024-3795(84)90068-5. ·兹伯利0559.65019
[23] I.S.Duff,MA57-A稀疏对称定和不定系统解的代码,ACM Trans。数学。《软件》,30(2004),第118-144页,https://doi.org/10.1145/992200.992202。 ·Zbl 1070.65525号
[24] A.Edelman和S.Jeong,《53矩阵因子分解:系统方法》,https://arxiv.org/abs/1204.08669, 2021.
[25] G.H.Golub和C.F.Van Loan,《矩阵计算》,约翰·霍普金斯数学研究所。科学。,第四版,约翰·霍普金斯大学出版社,巴尔的摩,2013年·Zbl 1268.65037号
[26] D.Henrion和P.Hippe,多项式矩阵J谱因式分解的双曲QR因式分解,第42届IEEE国际决策与控制会议论文集,2003年第4卷,第3479-3484页,https://doi.org/10.109/CDC.2003.1271685。
[27] N.J.Higham,(J)-正交矩阵:性质和生成,SIAM Rev.,45(2003),第504-519页,https://doi.org/10.1137/S0036144502414930。 ·Zbl 1034.65026号
[28] N.J.Higham,《矩阵的函数:理论与计算》,SIAM,费城,2008年,https://doi.org/10.1137/1.9780898717778。 ·Zbl 1167.15001号
[29] N.J.Higham、D.Mackey、N.Mackey和F.Tisseur,《保持矩阵群和矩阵平方根迭代的函数》,SIAM J.matrix Anal。申请。,26(2005),第849-877页,https://doi.org/10.1137/S0895479804442218。 ·兹比尔1079.65053
[30] N.J.Higham、C.Mehl和F.Tisseur,正则广义极分解,SIAM J.矩阵分析。申请。,31(2010),第2163-2180页,https://doi.org/10.1137/090765018。 ·Zbl 1222.15018号
[31] I.Houtzager,JQR/JRQ/JQL/JLQ分解,MATLAB中央文件交换,https://www.mathworks.com/matlabcentral/fileexchange/50329-jqr-jrqjql-jlq因子分解,2015。
[32] C.Kenney和A.J.Laub,矩阵符号函数的有理迭代方法,SIAM J.矩阵分析。申请。,12(1991),第273-291页,https://doi.org/10.1137/0612020。 ·Zbl 0725.65048号
[33] C.Kenney和A.J.Laub,《关于极分解的尺度牛顿法和矩阵符号函数》,SIAM J.matrix Ana。申请。,13(1992),第688-706页,https://doi.org/10.1137/0613044。 ·Zbl 0754.65043号
[34] U.Kintzel,有限维不定标量积空间中的Procrustes问题,线性代数应用。,402(2005),第1-28页,https://doi.org/10.1016/j.laa.2005.01.004。 ·Zbl 1072.15035号
[35] H.Li,H.Yang和H.Shao,双曲QR因子分解的扰动分析,计算。数学。申请。,63(2012),第1607-1620页,https://doi.org/10.1016/j.camwa.2012.03.036。 ·Zbl 1252.65080号
[36] H.Ltaief、D.Sukkari、A.Esposito、Y.Nakatsukasa和D.Keyes,分布式记忆系统上的大规模并行极分解,ACM Trans。并行计算。,6 (2019), https://doi.org/10.1145/3328723。 ·Zbl 1471.65023号
[37] D.S.Mackey、N.Mackeys和F.Tisseur,标量积空间中的结构化因式分解,SIAM J.矩阵分析。申请。,27(2005),第821-850页,https://doi.org/10.1137/040619363。 ·邮编1098.15005
[38] C.Mehl、V.Mehrmann和H.Xu,关于线性响应理论中出现的双重结构矩阵和铅笔,线性代数应用。,380(2004),第3-51页,https://doi.org/10.1016/S0024-3795(02)00455-X·Zbl 1083.65038号
[39] C.Mehl、A.C.M.Ran和L.Rodman,不定内积空间中正规算子的极分解,收录于Krein空间中的算子理论和非线性特征值问题,Oper。理论高级应用。162,Birkha¨user,巴塞尔,2006年,第277-292页,https://doi.org/10.1007/3-7643-7453-5_15。 ·Zbl 1108.47034号
[40] V.Mehrmann和F.Poloni,具有置换拉格朗日图基的加倍算法,SIAM J.矩阵分析。申请。,33(2012),第780-805页,https://doi.org/10.1137/10850773。 ·Zbl 1260.65059号
[41] Y.Nakatsukasa,Z.Bai和F.Gygi,优化哈雷迭代计算矩阵极分解,SIAM J.矩阵分析。申请。,31(2010),第2700-2720页,https://doi.org/10.1137/090774999。 ·Zbl 1215.65081号
[42] Y.Nakatsukasa和R.W.Freund,通过矩阵符号函数在两次迭代中精确计算基本矩阵分解:Zolotarev函数的幂,SIAM Rev.,58(2016),第461-493页,https://doi.org/10.1137/10990334。 ·Zbl 1383.15012号
[43] Y.Nakatsukasa和N.J.Higham,计算极分解迭代的向后稳定性,SIAM J.矩阵分析。申请。,33(2012),第460-479页,https://doi.org/10.1137/10857544。 ·Zbl 1252.65083号
[44] G.Onida、L.Reining和A.Rubio,《电子激发:密度-功能与多体格林函数方法》,修订版。物理。,74(2002),第601-659页,https://doi.org/10.1103/RevModPhys.74.601。
[45] C.Penke、A.Marek、C.Vorwerk、C.Draxl和P.Benner,偏对称特征值问题的高性能解及其在求解Bethe-Salpeter特征值问题中的应用,并行计算。,96 (2020), 102639, https://doi.org/10.1016/j.parco.2020.102639。
[46] F.Poloni,置换图矩阵及其应用,Springer,Cham,2015年,第107-129页,https://doi.org/10.1007/978-3-319-15260-8_5。 ·Zbl 1327.65083号
[47] 罗伯茨,线性模型简化和利用符号函数求解代数Riccati方程,国际。《控制杂志》,32(1980),第677-687页,https://doi.org/10.1080/00207178008922881。 ·Zbl 0463.93050号
[48] M.Shao、F.H.da Jornada、C.Yang、J.Deslippe和S.G.Louie,解决Bethe-Salpeter特征值问题的保结构并行算法,线性代数应用。,488(2016),第148-167页,https://doi.org/10.1016/j.laa.2015.09.036。 ·Zbl 1330.65059号
[49] S.Singer,不确定QR因子分解,BIT,46(2006),第141-161页,https://doi.org/10.1007/s10543-006-0044-5。 ·Zbl 1093.65029号
[50] S.Singer和S.Singer,不定因子分解的舍入误差和扰动界,线性代数应用。,309(2000),第103-119页,https://doi.org/10.1016/S0024-3795(99)00156-1. ·Zbl 0948.65042号
[51] D.Sukkari、H.Ltaief、A.Esposito和D.Keyes,分布式内存多核系统上基于QDWH的SVD软件框架,ACM Trans。数学。软件,45(2019),18,https://doi.org/10.1145/3309548。 ·Zbl 1471.65023号
[52] X.Sun和E.S.Quintana Ortí,块广义Schur分解的谱划分方法,数学。公司。,73(2004),第1827-1847页,https://doi.org/10.1090/S0025-5718-04-01667-9。 ·Zbl 1054.65034号
[53] K.Veselicí,线性系统的阻尼振动,数学课堂讲稿。,施普林格·弗拉格,柏林,2011年,https://doi.org/10.1007/978-3-642-21335-9。 ·兹比尔1232.37004
[54] D.Watkins,矩阵特征值问题,SIAM,费城,2007,https://doi.org/10.1137/1.9780898717808。 ·Zbl 1142.65038号
[55] Y.Yamamoto、Y.Nakatsukasa、Y.Yanagisawa和T.Fukaya,Cholesky QR2算法的舍入误差分析,Electron。事务处理。数字。分析。,44(2015),第306-326页·兹比尔1330.65049
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。