×

混合精度对称矩阵的迭代对角化及其在电子结构计算中的应用。 (英语) Zbl 1264.65056号

摘要:大矩阵的对角化是电子结构计算等许多应用中的计算瓶颈。我们表明,在保持64位精度的同时,利用32位浮点运算可以实现30%以上的加速。此外,在我们的实现中,大多数计算成本高昂的操作都是由三级BLAS/LAPACK例程执行的,从而在大多数平台上实现最佳性能。我们还讨论了考虑非对角线元素的特定问题预条件的有效性。

MSC公司:

2015财年65 矩阵特征值和特征向量的数值计算
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Martin,R.M.,《电子结构:基本理论与实践方法》(2004),剑桥大学出版社·Zbl 1152.74303号
[2] 马克思,D。;Hutter,J.,《Ab Initio分子动力学基本理论和高级方法》(2009),剑桥大学出版社
[3] 萨阿德,Y。;切里科夫斯基,J.R。;Shontz,S.M.,材料电子结构计算的数值方法,SIAM Rev.,52,3-54(2010)·邮编:1185.82004
[4] Bai,Z。;德梅尔,J。;Dongarra,J。;Ruhe,A。;van der Vorst,H.,《代数特征值问题求解模板:实用指南》(2000),SIAM·Zbl 0965.65058号
[5] 巴伯林,M。;Buttari,A。;Dongarra,J。;Kurzak,J。;Langou,J。;Langou,J。;Luszczek,P。;Tomov,S.,用混合精度算法加速科学计算,计算。物理学。Comm.,180,2526-2533(2009年)·兹比尔1197.65240
[6] Papadrakakis,M。;Bitoulas,N.,大型和病态问题预处理共轭梯度算法的准确性和有效性,计算。方法应用。机械。工程,109219-232(1993)·Zbl 0845.73072号
[7] Dongarra,J.J。;莫勒,C.B。;Wilkinson,J.H.,《提高计算特征值和特征向量的准确性》,SIAM J.Numer。分析。,20, 23-45 (1983) ·Zbl 0523.65021号
[8] Drygalla,V.,利用混合精度计算对称矩阵和奇异值的特征值,PAMM,8,10809-10810(2008)·Zbl 1392.65029号
[9] Choe,Y.-K。;Tsuchida,E。;Ikeshoji,T。;山中社。;Hyodo,S.,《聚合物电解质膜中质子动力学的性质》,nafion:第一原理分子动力学研究,Phys。化学。化学。物理。,11, 3892-3899 (2009)
[10] Choe,Y.-K。;津田,E。;Ikeshoji,T。;奥希拉,A。;Kidena,K.,《模拟水合聚合物电解质膜磺化聚醚砜(SPES)的从头算模拟研究》,J.Phys。化学。B、 1142411-2421(2010)
[11] Sameh,A.H。;Wisniewski,J.A.,广义特征值问题的轨迹最小化算法,SIAM J.Numer。分析。,19, 1243-1259 (1982) ·Zbl 0493.65017号
[12] Stich,I。;R车。;Parrinello,M。;Baroni,S.,《能量泛函的共轭梯度最小化:电子结构计算的新方法》,Phys。B版,39,4997-5004(1989)
[13] 佩恩,M.C。;特特,M.P。;阿伦·D.C。;阿里亚斯,T.A。;Joannopoulos,J.D.,从头算总能量计算的迭代最小化技术:分子动力学和共轭梯度,修订版。物理。,64, 1045-1097 (1992)
[14] Edelman,A。;Smith,S.T.,关于类特征问题的共轭梯度方法,BIT,36494-508(1996)·兹比尔0856.65037
[15] 出版社,W.H。;Teukolsky,S.A。;维特林,W.T。;Flannery,B.P.,《Fortran中的数字配方》(1992),剑桥大学出版社·Zbl 0778.65002号
[16] King-Smith,R.D。;Vanderbilt,D.,钙钛矿化合物铁电性的第一原理研究,Phys。版本B,49,5828-5844(1994)
[17] Annett,J.F.,Kohn-Sham密度泛函理论算法的效率,计算。马特。科学。,4, 23-42 (1995)
[18] 北马尔扎里。;范德比尔特,D。;Payne,M.C.,集成密度泛函理论从头算金属和有限温度绝缘体的分子动力学,物理学。修订稿。,791337-1340(1997年)
[19] 萨米赫,A。;Tong,Z.,对称广义特征值问题的迹最小化方法,J.Compute。申请。数学。,123, 155-175 (2000) ·Zbl 0964.65038号
[20] Goto,K。;van de Geijn,R.,《3级BLAS的高性能实现》,ACM Trans。数学。软件,35,1-14(2008)
[21] Kosloff,R。;Tal-Ezer,H.,网格上计算薛定谔方程本征函数和本征值的直接松弛方法,化学。物理学。莱特。,127, 223-230 (1986)
[22] 阿里亚斯,T.A。;佩恩,M.C。;Joannopoulos,J.D.,《从头算分子动力学:分析连续能量泛函和迭代解洞察》,《物理学》。修订稿。,69, 1077-1080 (1992)
[23] Sawamura,A。;Kohyama,M。;Keishi,T.,基于平面波的电子结构计算的有效预处理方案,计算。马特。科学。,第14页,第4-7页(1999年)
[24] 韦伯,V。;VandeVondele,J。;Hutter,J。;Niklasson,A.M.N.,正交约束下的直接能量函数最小化,J.Chem。物理。,128, 084113 (2008)
[25] 法特伯特,J.-L。;Bernholc,J.,《走向基于网格的密度泛函理论方法:优化非正交轨道和多重网格加速》,Phys。版本B,621713-1722(2000)
[26] Torsti,T.,《电子结构计算中的三种实空间离散化技术》,Phys。统计解决方案。(b) ,2431016-1053(2006)
[27] Tsuchida,E.,液体甲酰胺的从头算分子动力学研究,化学杂志。物理。,1214740-4746(2004年)
[28] M.Guidon。;希夫曼,F。;Hutter,J。;VandeVondele,J.,《使用混合密度泛函的从头算分子动力学》,J.Chem。物理。,128, 214104 (2008)
[29] Davidson,E.R.,大型实对称矩阵的几个最低特征值和相应特征向量的迭代计算,J.Compute。物理。,17, 87-94 (1975) ·Zbl 0293.65022号
[30] Hutter,J。;吕蒂,H.P。;Parrinello,M.,通过迭代子空间中的直接反演进行基于平面波的密度函数计算中的电子结构优化,计算。马特。科学。,2, 244-248 (1994)
[31] Sleijpen,G.L.G。;Van der Vorst,H.A.,线性特征值问题的Jacobi-Davidson迭代方法,SIAM Rev.,4267-293(2000)·Zbl 0949.65028号
[32] Ozaki,T.,《利用局部化基函数进行大规模电子结构计算的高效低阶标度方法》,Phys。B版,82,075131(2010)
[33] 鲍勒,D.R。;宫崎骏,T.,《用密度泛函理论计算数百万个原子:线性标度显示其潜力》,J.Phys.:康登斯。Matter,22074207(2010)
[34] 津田,E。;Tsukada,M.,基于自适应有限元方法的大尺度电子结构计算,J.Phys。Soc.Jpn.公司。,67, 3844-3858 (1998)
[35] 刘博士。;Nocedal,J.,《关于大规模优化的有限内存BFGS方法》,数学。程序。,45, 503-528 (1989) ·Zbl 0696.90048号
[36] Tsuchida,E.,《电子结构计算的有效算法》,J.Phys。Soc.Jpn.公司。,71, 197-203 (2002)
[37] 吉尔,体育。;Leonard,M.W.,用于大规模无约束优化的有限记忆约简-海西方法,SIAM J.Optim。,14, 380-401 (2003) ·Zbl 1046.65045号
[38] Gonze,X。;Lee,C.,《动力学矩阵》,玻恩有效电荷,介电常数张量,以及密度泛函微扰理论中的原子间力常数,物理学。B版,55,10355-10368(1997)
[39] Putrino,A。;塞巴斯蒂亚尼,D。;Parrinello,M.,广义变分密度泛函微扰理论,J.Chem。物理。,113, 7102-7109 (2000)
[40] 巴罗尼,S。;de Gironcoli,S。;Dal Corso,A。;Giannozzi,P.,密度泛函微扰理论中的声子和相关晶体特性,Rev.Mod。物理。,73, 515-562 (2001)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。