×

对称正定矩阵HSS预条件子的有效构造。 (英语) Zbl 07340708号

摘要:在求解具有稠密、病态、对称正定核矩阵的线性系统的迭代方法中,需要快速矩阵向量积和快速预处理操作。通过用(mathcal{H}^2)表示法或等效的快速多极方法表示法表示核矩阵,可以获得快速(线性缩放)矩阵向量积。本文研究了利用层次半可分(HSS)矩阵表示对此类矩阵进行预处理。在此之前,提出了一种算法来构造SPD核矩阵的HSS近似,该核矩阵被保证为SPD。然而,该算法具有二次代价,仅用于定义核矩阵的点的递归二进制划分。本文提出了一种构造SPD HSS近似的通用算法。重要的是,该算法使用SPD矩阵的(mathcal{H}^2)表示,将其计算复杂性从二次降低到拟线性。数值实验说明了这种SPD-HSS近似如何作为预处理器来求解由一系列核函数引起的线性系统。

MSC公司:

65F08个 迭代方法的前置条件
65层55 低阶矩阵逼近的数值方法;矩阵压缩
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] S.Ambikasaran和E.Darve,部分层次半可分矩阵的快速直接求解器,J.Sci。计算。,57(2013),第477-501页·兹比尔1292.65030
[2] S.Ambikasaran和E.Darve,反向快速多极方法,arXiv预印本,arXiv:1407.1772[math.NA],2014·兹比尔1292.65030
[3] S.Ambikasaran、M.O'Neil和K.R.Singh,层次矩阵的快速对称因式分解及其应用,arXiv预印本,arXiv:1405.0223[math.NA],2014年。
[4] S.Boörm、L.Grasedyck和W.Hackbusch,《层次矩阵》,讲座,注释21/2003,《自然科学院数学研究所》,2003年。
[5] L.Cambier、C.Chen、E.G.Boman、S.Rajamanickam、R.S.Tuminaro和E.Darve,使用低阶近似的代数稀疏嵌套剖分算法,SIAM J.矩阵分析。申请。,41(2020年),第715-746页·Zbl 1441.65048号
[6] S.Chandrasekaran、M.Gu和T.Pals,用于分层半可分表示的快速ULV分解求解器,SIAM J.矩阵分析。申请。,28(2006),第603-622页·Zbl 1120.65031号
[7] E.Chow和J.Skolnick,限制对细胞内大分子动力学模型的影响,Proc。国家。阿卡德。科学。美国,112(2015),第14846-14851页。
[8] P.Coulier、H.Pouransari和E.Darve,《逆快速多极子方法:使用快速近似直接解算器作为稠密线性系统的预条件》,SIAM J.Sci。计算。,39(2017年),第A761-A796页·Zbl 1365.65068号
[9] J.Feliu-Fabà,K.L.Ho,和L.Ying,椭圆偏微分方程的递归预处理分层插值因式分解,Commun。数学。科学。,18(2020年),第91-108页·Zbl 1437.65172号
[10] B.Fornberg和N.Flyer,用径向基函数求解偏微分方程,Acta Numer。,24(2015),第215页·Zbl 1316.65073号
[11] P.Ghysels、X.S.Li、F.-H.Rouet、S.Williams和A.Napov,使用随机抽样的新型HSS结构多额叶解算器的高效多核实现,SIAM J.Sci。计算。,38(2016),第S358-S384页·Zbl 1352.65092号
[12] L.Greengard和V.Rokhlin,粒子模拟的快速算法,J.Compute。物理。,73(1987),第325-348页·Zbl 0629.65005号
[13] L.Greengard和V.Rokhlin,三维拉普拉斯方程快速多极方法的新版本,Acta Numer。,6(1997年),第229-269页·Zbl 0889.65115号
[14] W.Hackbusch,基于矩阵的稀疏矩阵算法。第一部分:矩阵导论,计算,62(1999),第89-108页·Zbl 0927.65063号
[15] W.Hackbusch和S.Boërm,自适应矩阵的数据稀疏近似,计算,69(2002),第1-35页·Zbl 1012.65023号
[16] W.Hackbusch、B.Khoromskij和S.A.Sauter,《On(mathcal{H}^2)-矩阵》,摘自柏林施普林格大学应用数学讲座,2000年,第9-29页·Zbl 0963.65043号
[17] W.Hackbusch和B.N.Khoromskij,稀疏矩阵算法。第二部分:多维问题的应用,《计算》,64(2000),第21-47页·Zbl 0962.65029号
[18] N.Halko、P.Martinsson和J.Tropp,《寻找随机结构:构造近似矩阵分解的概率算法》,SIAM Rev.,53(2011),第217-288页·Zbl 1269.65043号
[19] M.J.Heaton、A.Datta、A.O.Finley、R.Furrer、J.Guinness、R.Guhaniyogi、F.Gerber、R.B.Gramacy、D.Hammerling、M.Katzfush、F.Lindgren、D.W.Nychka、F.Sun和Zammit-Mangion,大型空间数据分析方法之间的案例研究竞争,J.Agric。,生物与环境。《统计》,第24卷(2019年),第398-425页·Zbl 1426.62345号
[20] K.L.Ho和L.Greengard,通过递归骨架化快速直接求解结构化线性系统,SIAM J.Sci。计算。,34(2012年),第A2507-A2532页·Zbl 1259.65062号
[21] K.L.Ho和L.Ying,椭圆算子的层次插值因子分解:积分方程,Comm.Pure Appl。《数学》,69(2016),第1314-1353页·Zbl 1344.65123号
[22] H.Huang,X.Xing和E.Chow,H2Pack:使用代理点方法的内核矩阵的高性能矩阵包,ACM Trans。数学。软件,47(2020),3。
[23] L.Y.Kolotilina和A.Y.Yeremin,因子化稀疏近似逆预处理i.理论,SIAM J.矩阵分析。申请。,14(1993),第45-58页·Zbl 0767.65037号
[24] 李彦宏、杨浩、马丁、何国良和英,蝴蝶因式分解,多尺度模型。模拟。,13(2015),第714-732页·兹比尔1317.44004
[25] Y.Liu,X.Xing,H.Guo,E.Michielssen,P.Ghysels,and X.S.Li,通过随机矩阵向量乘法进行蝴蝶因子分解,arXiv预印本,arXiv:2002.03400[math.NA],2020年·Zbl 1462.65044号
[26] P.-G.Martinsson,用于计算矩阵的分层半可分表示的快速随机算法,SIAM J.matrix Anal。申请。,32(2011),第1251-1274页·Zbl 1237.65028号
[27] J.Rotne和S.Prager,聚合物中流体动力相互作用的变化处理,化学杂志。物理。,50(1969年),第4831-4837页。
[28] F.-H.Rouet、X.S.Li、P.Ghysels和A.Napov,使用随机化进行密集分层半可分矩阵计算的分布式内存包,ACM Trans。数学。《软件》,42(2016),第1-35页·Zbl 1369.65043号
[29] J.Xia,通用SPD矩阵的鲁棒有效eSIF预处理,arXiv预印本,arXiv:2007.03729[math.NA],2020。
[30] J.Xia、S.Chandrasekaran、M.Gu和X.S.Li,分层半可分矩阵的快速算法,Nume。线性代数应用。,17(2010年),第953-976页·Zbl 1240.65087号
[31] Xia和Z.Xin,通过结构不完全因式分解对一般SPD矩阵进行有效和稳健的预处理,SIAM J.矩阵分析。申请。,38(2017),第1298-1322页·Zbl 1386.15033号
[32] Z.Xin、J.Xia、S.Cauley和V.Balakrishnan,基于结构不完全因子分解的预处理技术的有效性和稳健性,Nume。线性代数应用。,27(2020),e2294·Zbl 07217203号
[33] X.Xing和E.Chow,一种用于核矩阵系统的块低秩逼近的有效方法,arXiv预印本,arXiv:11811.04134[math.NA],2018。
[34] X.Xing和E.Chow,在分层半可分矩阵逼近中保持正定性,SIAM J.矩阵分析。申请。,39(2018),第829-855页·Zbl 1391.65061号
[35] X.Xing和E.Chow,通过核矩阵的代理点进行插值分解,SIAM J.Matrix Ana。申请。,41(2020年),第221-243页·Zbl 1432.65052号
[36] H.Yamakawa,聚合物链在稀溶液中的传输特性:流体动力学相互作用,化学杂志。物理。,53(1970年),第436-443页。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。