×

求解大型特征值问题的带内-外重启的块Chebyshev-Davidson方法。 (英语) Zbl 1203.65077号

摘要:针对大型对称/厄米特征值问题,我们提出了一种使用切比雪夫多项式滤波器的块Davidson型子空间迭代。该方法由三个基本组成部分组成。第一种是构造高效块切比雪夫多项式滤波器的自适应过程;第二种是Chebyshev-Davidson迭代中的内外重启技术,该技术降低了与使用大维子空间相关的计算成本;第三种是渐进滤波技术,如果可用,它可以充分利用大量良好的初始向量,而不需要使用大的块大小。基于密度泛函理论计算的几个哈密顿矩阵的数值实验表明了该方法的有效性和鲁棒性。

理学硕士:

2015财年65 矩阵特征值和特征向量的数值计算
15B57号 厄米特矩阵、斜厄米特阵和相关矩阵
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Parlett,B.N.,《对称特征值问题》。对称特征值问题,应用数学经典,第20卷(1998),SIAM:SIAM费城,宾夕法尼亚州·Zbl 0885.65039号
[2] Stewart,G.W.,《矩阵算法II:特征系统》(2001),SIAM:SIAM Philadelphia·Zbl 0984.65031号
[3] Saad,Y.,《大型特征值问题的数值方法》(1992),John Wiley:John Wiley纽约·Zbl 0991.65039号
[4] Golub,G.H。;Loan,C.F.V.,《矩阵计算》(1996),约翰霍普金斯大学出版社:约翰霍普金大学出版社,马里兰州巴尔的摩·Zbl 0865.65009号
[5] (Bai,Z.;Demmel,J.;Dongarra,J.,Ruhe,A.;van der Vorst,H.,《代数特征值问题的解决模板》(2000),SIAM:SIAM Philadelphia,PA)·Zbl 0965.65058号
[6] Lanczos,C.,解线性微分和积分算子特征值问题的迭代方法,J.Res.Nat.Bur。标准,45,255-282(1950)
[7] Sorensen,D.C.,多项式滤波器在步长Arnoldi方法中的隐含应用,SIAM J.矩阵分析。申请。,13, 357-385 (1992) ·Zbl 0763.65025号
[8] Wu,K。;Simon,H.,《大型对称特征值问题的厚重启动Lanczos方法》,SIAM J.矩阵分析。申请。,22602-616(2000年)·Zbl 0969.65030号
[9] Stewart,G.W.,大型特征问题的Krylov-Schur算法,SIAM J.矩阵分析。申请。,23, 601-614 (2001) ·Zbl 1003.65045号
[10] Davidson,E.R.,大型实对称矩阵的几个最低特征值和相应特征向量的迭代计算,J.Compute。物理。,17, 87-94 (1975) ·Zbl 0293.65022号
[11] 摩根R.B。;Scott,D.S.,计算稀疏对称矩阵特征值的Davidson方法的推广,SIAM J.Sci。统计计算。,7, 817-825 (1986) ·Zbl 0602.65020号
[12] Sleijpen,G.L.G。;van der Vorst,H.A.,线性特征值问题的Jacobi-Davidson迭代法,SIAM J.矩阵分析。申请。,17, 401-425 (1996) ·Zbl 0860.65023号
[13] 索伦森特区。;Yang,C.,大规模特征值计算的截断rq迭代,SIAM J.矩阵分析。申请。,19, 1045-1073 (1998) ·Zbl 0918.65023号
[14] Stathopoulos,A.,有限内存下厄米特特征问题的近似最优预处理方法。第一部分:寻找一个特征值,SIAM J.Sci。计算。,29, 2, 481-514 (2007) ·Zbl 1137.65019号
[15] Stathopoulos,A。;McCombs,J.R.,有限记忆条件下hermitian本征问题的近似最优预条件方法。第二部分:寻找许多特征值,SIAM J.Sci。计算。,29, 5, 2162-2188 (2007) ·Zbl 1151.65320号
[16] Brandt,A。;McCormick,S。;Ruge,J.,微分特征问题的多重网格方法,SIAM J.Sci。统计计算。,4, 244-260 (1983) ·Zbl 0517.65083号
[17] 曼德尔,J。;McCormick,S.,一种用于=λBu在复合网格上,J.Compute。物理。,80, 442-452 (1989) ·Zbl 0674.65013号
[18] Briggs,E.L。;沙利文,D.J。;Bernholc,J.,《具有多重网格加速度的大规模电子结构计算》,Phys。版本B,52,R5471-5474(1995)
[19] 科斯蒂纳,S。;Táasan,S.,《大规模特征值问题的自适应多重网格技术:二维和三维薛定谔问题的解决方案》,《物理学》。E版,51,3704-3717(1995)
[20] 丽芙,O。;Brandt,A.,《多尺度本征基计算:O(N)log(N)中的本征函数》,(Barth,T.;Chan,T.,Haimes,R.,多尺度和多分辨率方法,计算机科学工程讲义,第20卷(2002),Springer-Verlag),347-357·Zbl 0989.65129号
[21] Vömel,C。;托莫夫,S.Z。;马奎斯,O.A。;坎宁,A。;王立伟。;Dongarra,J.J.,用于大规模纳米系统电子结构计算的最新特征解算器,J.Compute。物理。,227, 15, 7113-7124 (2008) ·Zbl 1141.82346号
[22] 萨博,A。;Ostlund,N.S.,《现代量子化学:高级电子结构理论导论》(1996),多佛出版社:纽约多佛出版社
[23] 帕尔·R。;杨伟,《原子和分子密度泛函理论》(1989),牛津大学出版社
[24] Rutishauser,H.,F.L.Bauer同步迭代法的计算方面,数值。数学。,13, 4-13 (1969) ·Zbl 0182.21304号
[25] Rutishauser,H.,对称矩阵的同步迭代法,数值。数学。,16, 205-223 (1970) ·兹比尔0239.65037
[26] Rutishauser,H.,对称矩阵的同步迭代法,(Wilkinson,J.H.;Reinsh,C.,《自动计算手册(线性代数)》,第二卷(1971年),Springer-Verlag),284-302
[27] Ruhe,A.,计算大型稀疏对称矩阵特征值的带Lanczos算法的实现方面,数学。计算。,33, 680-687 (1979) ·Zbl 0443.65022号
[28] 格里姆斯,R.G。;刘易斯,J.G。;Simon,H.D.,解稀疏对称广义特征问题的移位块Lanczos算法,SIAM J.矩阵分析。申请。,15, 228-272 (1994) ·Zbl 0803.65044号
[29] Knyazev公司。;Neymeyr,K.,在局部最优块共轭梯度法中使用多重网格预条件有效解决对称特征值问题,Electron。事务处理。数字。分析。,15, 38-55 (2003) ·Zbl 1031.65126号
[30] 周,Y。;Saad,Y.,大型对称特征值问题的块Krylov-Schur方法,Numer。藻类。,47, 4, 341-359 (2008) ·Zbl 1153.65330号
[31] 周,Y。;Saad,Y.,大型对称特征值问题的Chebyshev-Davidson算法,SIAM J.矩阵分析。申请。,29, 3, 954-971 (2007) ·Zbl 1151.65321号
[32] Zhou,Y.,Jacobi-Davidson研究,瑞利商迭代,逆迭代广义Davidson和牛顿更新,Numer。线性代数应用。,13, 8, 621-642 (2006) ·Zbl 1174.65372号
[33] Chebyshev,P.L.,《公共图书馆》,梅姆。阿卡德。科学。佩特斯布,7539-568(1854)
[34] Rogonski,W.W.,《多项式的一些初等不等式》,数学。天然气。,39, 7-12 (1955) ·Zbl 0064.01802号
[35] Rivlin,T.J.,《函数逼近导论》(2003),多佛·Zbl 0189.06601号
[36] Lehoucq,R.B。;索伦森特区。;Yang,C.,《ARPACK用户指南:用隐式重启Arnoldi方法解决大规模特征值问题》(1998),SIAM:SIAM Philadelphia·Zbl 0901.65021号
[37] Saad,Y.,Chebyshev求解非对称特征值问题的加速技术,数学。计算。,42, 166, 567-588 (1984) ·Zbl 0539.65013号
[38] 周,Y。;萨阿德,Y。;蒂亚戈,M.L。;Chelikowsky,J.R.,使用Chebyshev滤波子空间迭代的自洽场计算,J.Compute。物理。,219, 1, 172-184 (2006) ·Zbl 1105.65111号
[39] 周,Y。;Li,R.-C,《大埃尔米特矩阵谱的定界》,线性代数应用。(2010)
[40] Notay,Y.,部分对称特征值问题的Jacobi-Davidson和共轭梯度组合,Numer。线性代数应用。,9, 1, 21-44 (2002) ·Zbl 1071.65516号
[41] 丹尼尔,J。;Gragg,W.B。;考夫曼,L。;Stewart,G.W.,更新Gram-Schmidt QR分解的重正交化和稳定算法,数学。计算。,30, 772-795 (1976) ·Zbl 0345.65021号
[42] Watkins,D.S。;Elsner,L.,特征值问题分解型算法的收敛性,线性代数应用。,41, 19-47 (1991) ·Zbl 0712.65025号
[43] Lehoucq,R.B.,《隐式重新启动Arnoldi方法和子空间迭代》,SIAM J.Matrix Ana。申请。,23, 551-562 (2001) ·Zbl 1004.65044号
[44] T.Davis,佛罗里达大学稀疏矩阵收集,技术报告,佛罗里达大学,2009年(《北美文摘》,第97卷,第23期,1997年6月7日)。;T.Davis,佛罗里达大学稀疏矩阵收集,技术报告,佛罗里达大学,2009年(《北美文摘》,第97卷,第23期,1997年6月7日)。
[45] 科恩,W。;Sham,L.J.,《包含交换和相关效应的自洽方程》,《物理学》。修订版,140,A1133-A1138(1965)
[46] Kohn,W.,诺贝尔讲座:物质波函数和密度泛函的电子结构,Rev.Mod。物理。,71, 5, 1253-1266 (1999)
[47] 切里科夫斯基,J.R。;Troullier,N。;Saad,Y.,有限差分赝势法:无基础的电子结构计算,Phys。修订稿。,72, 1240-1243 (1994)
[48] Chelikowsky,J.,应用于纳米结构的伪电位密度泛函方法,J.Phys。D: 申请。物理。,33,R33-R50(2000)
[49] Knyazev,A.V.,朝向最优预处理特征解算器:局部最优块预处理共轭梯度法,SIAM J.Sci。计算。,23, 2, 517-541 (2001) ·Zbl 0992.65028号
[50] Fokkema,D.R。;Sleijpen,G.L.G。;van der Vorst,H.A.,雅各比·戴维森风格二维码QZ公司矩阵铅笔简化算法,SIAM J.Sci。计算。,20, 1, 94-125 (1998) ·兹比尔0924.65027
[51] 巴格拉马,J。;Calvetti,D。;Reichel,L.,IRBL:用于大规模厄米特特征问题的隐式重新启动块-Lanczos方法,SIAM J.Sci。计算。,24, 1650-1677 (2003) ·Zbl 1044.65027号
[52] 巴格拉马,J。;Calvetti,D。;Reichel,L.,irbleigs:一个用于计算大型稀疏厄米矩阵的几个特征对的MATLAB程序,ACM-Trans。数学。软件,5337-348(2003)·Zbl 1070.65529号
[53] 萨阿德,Y。;Stathopoulos,A.公司。;J.切利科夫斯基。;Wu,K。;奥尤特,S.,电子结构计算中大特征值问题的解决,BIT,36,3,563-578(1996)·Zbl 0862.65059号
[54] Stathopoulos,A。;奥胡特,S。;萨阿德,Y。;J.切利科夫斯基。;Kim,H.,预测材料性能的并行方法和工具,计算。科学。工程,219-32(2000)
[55] Wu,K。;坎宁,A。;西蒙,H.D。;Wang,L.-W.,电子结构计算的厚重启动Lanczos方法,J.Compute。物理。,154, 156-173 (1999) ·Zbl 0963.82048号
[56] 阿本茨,P。;赫特马纽克,美国。;Lehoucq,R.B。;Tuminara,R.S.,《使用AMG预处理迭代方法进行大规模三维模型分析的特征解算器比较》,国际期刊编号。方法工程,64,2,204-236(2005)·Zbl 1093.74021号
[57] 周,Y。;萨阿德,Y。;蒂亚戈,M.L。;Chelikowsky,J.R.,使用Chebyshev滤波子空间加速度的并行自我场计算,Phys。版本E,74,6,066704-1-0667040-8(2006)
[58] 蒂亚戈,M.L。;周,Y。;Alemany,M。;萨阿德,Y。;Chelikowsky,J.R.,《铁从原子到体的磁性演化》,《物理学》。修订稿。,97, 147201 (2006)
[59] 切里科夫斯基,J.R。;萨阿德,Y。;Chan,T.-L。;蒂亚戈,M.L。;扎亚克,A.T。;周,Y.,《网格上的赝势:应用于硅纳米晶体的电子、光学和振动特性》,J.Compute。西奥。纳米科学。,6, 6, 1247-1261 (2009)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。