×

一种使用低阶近似的代数稀疏嵌套剖分算法。 (英语) Zbl 1441.65048号

作者感兴趣的是求解大型对称正定稀疏线性系统(Ax=b\),(A\in\mathbb{R}^{N\次N})并提出了一种基于嵌套分割、稀疏化和低秩压缩的新算法。
该方法基于下面描述的分层插值因子分解思想K.L.Ho先生李英【公共纯应用数学69,第8期,1415-1451(2016;Zbl 1353.35142号)]. 然而,存在一些差异、改进和新功能。例如:该算法是完全通用的,可以应用于任何对称正定矩阵。唯一需要的输入是稀疏矩阵本身。如果一些几何信息可用,则可以使用这些信息来提高排序和聚类的质量;在算法中加入了一个额外的步骤来缩放对角线块,以很小的额外成本大大提高了预处理器的精度;使用正交(而不是插值)变换,提高稳定性,并确保当(A)为SPD时,预条件保持SPD。作者在一些大型问题上对该算法进行了评估,结果表明该算法具有近似线性的标度。分解时间大致为(O(N)),迭代次数随(N)缓慢增长。

MSC公司:

65层50 稀疏矩阵的计算方法
65F08个 迭代方法的前置条件
65层10 线性系统的迭代数值方法
65年20月 数值算法的复杂性和性能
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] S.Ambikasaran,《稠密数值线性代数的快速算法及其应用》,博士论文,斯坦福大学,加利福尼亚州斯坦福市,2013年。
[2] P.Amestoy、C.Ashcraft、O.Boiteau、A.Buttari、J.-Y.L'Excellent和C.Weisbecker,通过块低阶表示改进多额叶方法,SIAM J.Sci。计算。,37(2015),第A1451-A1474页,https://doi.org/10.1137/120903476。 ·Zbl 1314.05111号
[3] P.R.Amestoy,I.S.Duff,J.-Y.L'Excellent,和J.Koster,MUMPS:通用分布式内存稀疏解算器,《应用并行计算》,施普林格,柏林,海德堡,2000年,第121-130页。
[4] E.Anderson、Z.Bai、C.Bischof、L.S.Blackford、J.Demmel、J.Dongarra、J.Du Croz、A.Greenbaum、S.Hammarling、A.McKenney和D.Sorensen,《LAPACK用户指南》,SIAM,费城,1999年,https://doi.org/10.1137/1.9780898719604。 ·Zbl 0934.65030号
[5] J.Barnes和P.Hut,一种分层的力计算算法,《自然》,324(1986),第446-449页。
[6] M.Bebendorf,具有非光滑系数的一般二阶椭圆算子Galerkin矩阵的有效反演,数学。公司。,74(2005),第1179-1199页·Zbl 1330.65173号
[7] M.Bebendorf和W.Hackbusch,系数为(L_)-的椭圆算子的逆FE-矩阵的(mathcal{H})-矩阵逼近的存在性,Numer。数学。,95(2003),第1-28页·Zbl 1033.65100号
[8] M.J.Berger和S.H.Bokhari,多处理器上非均匀问题的分区策略,IEEE Trans。计算。,36(1987年),第570-580页。
[9] J.H.Bramble,《多重网格方法》,第294卷,CRC出版社,1993年·Zbl 0786.65094号
[10] A.Brandt、S.McCormick和J.Ruge,《自动多重网格解决方案的代数多重网格(AMG)及其在大地测量计算中的应用》,计算研究所报告,科罗拉多州柯林斯堡,1982年。
[11] 陈振华,揭示秩的QR分解,线性代数应用。,88(1987),第67-82页·Zbl 0624.65025号
[12] S.Chandrasekaran、P.Dewilde、M.Gu、T.Pals、X.Sun、A.-J.van der Veen和D.White,《序列半可分表示的一些快速算法》,SIAM J.矩阵分析。申请。,27(2005),第341-364页,https://doi.org/10.1137/S0895479802405884。 ·Zbl 1091.65063号
[13] S.Chandrasekaran,P.Dewilde,M.Gu和N.Somasunderam,关于离散椭圆偏微分方程Schur补的非对角块的数值秩,SIAM J.矩阵分析。申请。,31(2010),第2261-2290页,https://doi.org/10.1137/090775932。 ·Zbl 1209.65032号
[14] S.Chandrasekaran、M.Gu和W.Lyons,《分层半可分表示的快速稳定自适应解算器》,未出版,2004年。
[15] S.Chandrasekaran、M.Gu和T.Pals,用于分层半可分表示的快速ULV分解求解器,SIAM J.矩阵分析。申请。,28(2006),第603-622页,https://doi.org/10.1137/S0895479803436652。 ·兹比尔1120.65031
[16] C.Chen、L.Cambier、E.G.Boman、S.Rajamanickam、R.S.Tuminaro和E.Darve,《适用于冰盖建模的病态系统的鲁棒分层求解器》,J.Compute。物理。,396(2019),第819-836页·Zbl 1452.65079号
[17] Y.Chen、T.A.Davis、W.W.Hager和S.Rajamanickam,算法887:CHOLMOD,超节点稀疏Cholesky因子分解和更新/停机,ACM Trans。数学。软件,35(2008),22。
[18] H.Cheng、Z.Gimbutas、P.-G.Martinsson和V.Rokhlin,关于低秩矩阵的压缩,SIAM J.Sci。计算。,26(2005),第1389-1404页,https://doi.org/10.1137/030602678。 ·Zbl 1083.65042号
[19] M.Christie和M.Blunt,《第十个SPE比较解决方案项目:升级技术的比较》,石油工程师学会,2001年。
[20] T.A.Davis和Y.Hu,佛罗里达大学稀疏矩阵收集,ACM Trans。数学。软件,38(2011),1·Zbl 1365.65123号
[21] M.Faverge、G.Pichon、P.Ramet和J.Roman,《关于在PaStiX中使用(mathcal{H})矩阵算法的初步研究》,法国图卢兹快速直接求解研讨会,2015年;可在线访问https://hal.inia.fr/hal-01187882。
[22] R.P.Fedorenko,解椭圆差分方程的松弛方法,计算。数学。数学。物理。,1(1962年),第1092-1096页·Zbl 0163.39303号
[23] J.Feliu-Fabà,K.L.Ho,和L.Ying,椭圆偏微分方程的递归预处理层次插值因子分解,预印本,https://arxiv.org/abs/1808.01364, 2018. ·Zbl 1437.65172号
[24] W.Fong和E.Darve,黑盒快速多极子方法,J.Compute。物理。,228(2009),第8712-8725页·Zbl 1177.65009号
[25] A.George,规则有限元网格的嵌套剖分,SIAM J.Numer。分析。,10(1973年),第345-363页,https://doi.org/10.1137/0710032。 ·Zbl 0259.65087号
[26] P.Ghysels,X.S.Li,F.-H.Rouet,S.Williams和A.Napov,一种使用随机采样的新型HSS结构多前沿解算器的高效多核实现,SIAM J.Sci。计算。,38(2016),第S358-S384页,https://doi.org/10.1137/15M1010117。 ·Zbl 1352.65092号
[27] G.H.Golub和C.F.Van Loan,《矩阵计算》,第4卷,约翰霍普金斯大学出版社,马里兰州巴尔的摩,2013年·Zbl 1268.65037号
[28] L.Grasedyck、R.Kriemann和S.Le Borne,基于区域分解的LU预处理,数值。数学。,112(2009),第565-600页·Zbl 1178.65140号
[29] L.Greengard和V.Rokhlin,粒子模拟的快速算法,J.Comput。物理。,73(1987),第325-348页,https://doi.org/10.1016/0021-9991(87)90140-9. ·Zbl 0629.65005号
[30] M.Gu和S.C.Eisenstat,计算强秩揭示QR因式分解的高效算法,SIAM J.Sci。计算。,17(1996),第848-869页,https://doi.org/10.1137/0917055。 ·Zbl 0858.65044号
[31] W.Hackbusch,基于矩阵的稀疏矩阵算法。第一部分:矩阵导论,计算,62(1999),第89-108页·Zbl 0927.65063号
[32] W.Hackbusch,多网格方法和应用,Springer Ser。计算。数学。4,施普林格,柏林,海德堡,2013年·Zbl 0595.65106号
[33] W.Hackbusch和S.Boírm,自适应矩阵的数据解析近似,计算,69(2002),第1-35页·Zbl 1012.65023号
[34] W.Hackbusch和S.Bo¨rm,(mathcal{H}^2)-积分算子的插值矩阵逼近,应用。数字。数学。,43(2002),第129-143页·Zbl 1019.65103号
[35] P.Heínon,P.Ramet,and J.Roman,PaStiX:稀疏对称正定系统的高性能并行直接求解器,parallel Compute。,28(2002),第301-321页·Zbl 0984.68208号
[36] M.R.Hestenes和E.Stiefel,《求解线性系统的共轭梯度方法》,《自然科学研究杂志》。《标准》,49(1952),第409-436页·Zbl 0048.09901号
[37] K.L.Ho和L.Ying,椭圆算子的层次插值因式分解:微分方程,通信纯应用。数学。,69(2016),第1415-1451页·Zbl 1353.35142号
[38] G.Karypis和V.Kumar,用于划分不规则图的快速高质量多级方案,SIAM J.Sci。计算。,20(1998)第359-392页,https://doi.org/10.1137/S1064827595287997。 ·Zbl 0915.68129号
[39] E.Liberty、F.Woolfe、P.-G.Martinsson、V.Rokhlin和M.Tygert,矩阵低阶近似的随机算法,Proc。国家。阿卡德。科学。美国,104(2007),第20167-20172页·Zbl 1215.65080号
[40] R.J.Lipton、D.J.Rose和R.E.Tarjan,广义嵌套解剖,SIAM J.Numer。分析。,16(1979年),第346-358页,https://doi.org/10.1137/0716027。 ·Zbl 0435.65021号
[41] L.Miranian和M.Gu,揭示LU分解的强秩,线性代数应用。,367(2003),第1-16页·兹比尔1020.65016
[42] C.C.Paige和M.A.Saunders,稀疏线性方程组的解,SIAM J.Numer。分析。,12(1975年),第617-629页,https://doi.org/10.1137/0712047。 ·Zbl 0319.65025号
[43] H.Pouransari、P.Coulier和E.Darve,《使用扩展稀疏化和低秩近似的稀疏矩阵快速分层求解器》,SIAM J.Sci。计算。,39(2017年),第A797-A830页,https://doi.org/10.1137/15M1046939。 ·Zbl 1365.65072号
[44] A.Prokopenko、C.M.Siefert、J.J.Hu、M.Hoemmen和A.Klinvex,《Ifpack用户指南1.0》,技术报告SAND2016-5338,新墨西哥州阿尔伯克基Sandia国家实验室,2016年。
[45] Y.Saad,ILUT:双阈值不完全LU因式分解,Numer。线性代数应用。,1(1994年),第387-402页·Zbl 0838.65026号
[46] Y.Saad,《稀疏线性系统的迭代方法》,第2版,SIAM,费城,2003年,https://doi.org/10.1137/1.9780898718003。 ·Zbl 1031.65046号
[47] Y.Saad和M.H.Schultz,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计师。计算。,7(1986),第856-869页,https://doi.org/10.1137/0907058。 ·Zbl 0599.65018号
[48] Y.Saad和J.Zhang,BILUM:一般稀疏线性系统的多重消除和多级ILU预处理程序的块版本,SIAM J.Sci。计算。,20(1999),第2103-2121页,https://doi.org/10.1137/S16482759732753X。 ·Zbl 0956.65026号
[49] P.G.Schmitz和L.Ying,二维一般网格上椭圆问题的快速直接求解器,J.Compute。物理。,231(2012),第1314-1338页·Zbl 1408.65022号
[50] K.Stuében,《代数多重网格综述》,J.Compute。申请。数学。,128(2001),第281-309页·Zbl 0979.65111号
[51] D.A.Sushnikova和I.V.Oseledets,对称正定稀疏矩阵的“压缩和消除”解算器,SIAM J.Sci。计算。,40(2018),第A1742-A1762页,https://doi.org/10.1137/16M1068487。 ·兹比尔1391.65052
[52] I.K.Tezaur、M.Perego、A.G.Salinger、R.S.Tuminaro和S.F.Price,Albany/FELIX:一个用于高级分析的并行、可扩展和稳健的有限元一阶斯托克斯近似冰盖解算器,Geo。模型开发(在线),8(2015)。
[53] H.A.van der Vorst,Bi-CGSTAB:非对称线性系统解的Bi-CG的一个快速且平滑收敛的变体,SIAM J.Sci。统计师。计算。,13(1992年),第631-644页,https://doi.org/10.1137/0913035。 ·Zbl 0761.65023号
[54] Xia,一般大型稀疏矩阵的高效结构化多前沿因式分解,SIAM J.Sci。计算。,35(2013年),第A832-A860页,https://doi.org/10.1137/120867032。 ·Zbl 1266.15022号
[55] J.Xia,随机稀疏直接解算器,SIAM J.矩阵分析。申请。,34(2013),第197-227页,https://doi.org/10.1137/12087116X。 ·Zbl 1269.65029号
[56] J.Xia、S.Chandrasekaran、M.Gu和X.S.Li,大型结构线性方程组的超快速多面方法,SIAM J.矩阵分析。申请。,31(2009),第1382-1411页,https://doi.org/10.1137/09074543X。 ·Zbl 1195.65031号
[57] J.Xia、S.Chandrasekaran、M.Gu和X.S.Li,分层半可分矩阵的快速算法,数值。线性代数应用。,17(2010年),第953-976页·Zbl 1240.65087号
[58] Xia和M.Gu,秩结构对称正定矩阵的鲁棒近似Cholesky分解,SIAM J.矩阵分析。申请。,31(2010),第2899-2920页,https://doi.org/10.1137/090750500。 ·Zbl 1217.65061号
[59] J.Xia和Z.Xin,通过结构化不完全因子分解对一般SPD矩阵进行有效和稳健的预处理,SIAM J.Matrix Anal。申请。,38(2017),第1298-1322页,https://doi.org/10.1137/17M1124152。 ·兹伯利1386.15033
[60] K.Yang、H.Pouransari和E.Darve,保证收敛的稀疏层次解算器,国际。J.数字。方法工程,120(2019),第964-986页。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。