×
马蒂亚斯·博尔霍夫;奥拉夫·申克;法比奥·韦波西奥

一种高性能的块近似LU因子分解预处理器算法。 （英语） Zbl 1460.65046号

申请。数字。数学。 162, 265-282 (2021).
摘要：导致求解线性系统的许多应用问题都使用预处理的Krylov子空间解算器来计算其解。在最流行的预处理方法中，不完全因式分解方法可以是单层方法，也可以是多级框架中的方法。我们将提出一种基于巧妙地在初始和整个因式分解过程中阻塞系统的块不完全因次分解。我们的目标是通过Krylov子空间方法开发代数块预条件，以有效地求解此类系统。我们将演示我们的逐层近似算法如何在现代体系结构上以数量级的方式优于单层逐块标量方法，从而为其在各种多级不完全因子分解方法或其他核心部分依赖于高效不完全因子化算法的应用中的未来应用铺平了道路。

MSC公司：

65层50 稀疏矩阵的计算方法
65F08个 迭代方法的前置条件
65层10 线性系统的迭代数值方法

关键词：

稀疏矩阵;不完全LU分解;块结构方法;稠密矩阵核;块ILU方法

软件：

YSMP公司;VBARMS公司;MUMPS公司;UMFPACK公司;超级LU;伊波特;手臂;业务流程工具包;帕迪索
PDF格式BibTeX公司 XML格式引用
\textit{M.Bollhöfer}等人，应用。数字。数学。162265--282（2021；Zbl 1460.65046）
全文： 内政部

参考文献：

[1] Ajiz，医学硕士。；Jennings，A.，一种稳健的不完全Choleski共轭梯度算法，国际期刊Numer。方法工程，20，5，949-966（1984）·Zbl 0541.65019号
[2] 阿梅斯托，P。；Davis，T.A。；达夫，I.S.，《近似最小度排序算法》，SIAM J.矩阵分析。申请。，17, 4, 886-905 (1996) ·兹比尔0861.65021
[3] Amestoy，P.R.（艾姆斯泰，P.R.）。；达夫，I.S。；科斯特，J。；L'Excellence，J.-Y.，使用分布式动态调度的完全异步多前沿解算器，SIAM J.Matrix Ana。申请。，23, 1, 15-41 (2001) ·Zbl 0992.65018号
[4] Amestoy，P.R.（艾姆斯泰，P.R.）。；达夫，I.S。；L'Excellent，J.-Y。；罗伯特·Y。；Rouet，F.-H。；Uçar，B.，《关于在非核心环境中计算稀疏矩阵的逆项》，SIAM J.Sci。统计计算。，34、4、A1975-A1999（2012）·Zbl 1252.05108号
[5] Benzi，M.，《大型线性系统的预处理技术：综述》，J.Compute。物理。，182, 418-477 (2002) ·Zbl 1015.65018号
[6] Benzi，M。；霍斯，J.C。；Tůma，M.，预处理高度不定和非对称矩阵，SIAM J.Sci。计算。，22, 4, 1333-1353 (2000) ·Zbl 0985.65036号
[7] 博勒费尔，M。；Saad，Y.，基于逆ILU构建的多层预处理子，SIAM J.Sci。计算。，27, 5, 1627-1650 (2006) ·Zbl 1104.65037号
[8] 博勒费尔，M。；格罗特，M.J。；Schenk，O.，异质介质中亥姆霍兹方程的代数多层预条件，SIAM J.Sci。计算。，31, 5, 3781-3805 (2009) ·Zbl 1203.65273号
[9] Carpentieri，B。；廖，J。；Sosonkina，M.，VBARMS：稀疏线性系统的可变块代数递归多级解算器，J.Compute。申请。数学。，259, 164-173 (2014) ·Zbl 1291.65095号
[10] 查普曼，A。；萨阿德，Y。；Wigton，L.，CFD问题的高阶ILU预条件，国际期刊数值。方法流体，33，6677-788（2000）·兹比尔0959.76077
[11] 周，E。；Heroux，M.A.，块预处理的面向对象框架，ACM Trans。数学。软质。，24, 2, 159-183 (1998) ·兹伯利0930.65052
[12] 柯蒂斯，F.E。；胡贝尔，J。；申克，O。；Wächter，A.，关于用不精确步长计算实现非线性优化的内点算法的注释，数学。程序。，136, 1, 209-227 (2012) ·Zbl 1255.49043号
[13] Davis，T.A.，算法832:UMFPACK V4.3——非对称模式多波前方法，ACM Trans。数学。软质。，30, 2, 196-199 (2004) ·Zbl 1072.65037号
[14] D'Azevedo，E.F。；福赛斯，宾夕法尼亚州。；Tang，W.-P.，《寻找一种具有高水平填充的成本效益的ILU预处理器》，BIT编号。数学。，32, 3, 442-463 (1992) ·Zbl 0761.65017号
[15] Dolan，D.E。；Moré，J.J.，《性能曲线基准优化软件》，数学。程序。，91, 2, 201-213 (2002) ·邮编：1049.90004
[16] 达夫，I.S。；Koster，J.，《将大型条目排列到稀疏矩阵对角线的算法的设计和使用》，SIAM J.Matrix Anal。申请。，20, 4, 889-901 (1999) ·Zbl 0947.65048号
[17] 达夫，I.S。；Pralet，S.，稀疏对称不定问题的缩放和旋转策略，SIAM J.矩阵分析。申请。，27, 2, 313-340 (2005) ·Zbl 1092.65037号
[18] 达夫，I.S。；Pralet，S.，《稀疏对称不定系统序列和并行解的稳定混合旋转策略》，SIAM J.Matrix Ana。申请。，29, 3, 1007-1024 (2007) ·Zbl 1145.65021号
[19] 达夫，I.S。；埃里斯曼，A.M。；Reid，J.K.，《稀疏矩阵的直接方法》（2017），牛津大学出版社·Zbl 1364.65067号
[20] 艾森斯塔特，S.C。；古尔斯基，M.C。；舒尔茨，M.H。；谢尔曼，A.H.，耶鲁稀疏矩阵包I：对称码，国际数学家杂志。方法工程，18，8，1145-1151（1982）·Zbl 0492.65012号
[21] 新泽西州古尔德。；斯科特·J·A。；Hu，Y.，用于求解大型稀疏对称线性方程组的稀疏直接求解器的数值评估，ACM Trans。数学。软质。，33, 10:1-10:31 (2007) ·兹比尔1365.65129
[22] 古普塔，A。；George，T.，《改进不完全因子分解预处理性能的自适应技术》，SIAM J.Sci。计算。，32, 1, 84-110 (2010) ·Zbl 1209.65036号
[23] M.哈格曼。；Schenk，O.，对称不定线性系统预处理的加权匹配，SIAM J.Sci。计算。，28, 2, 403-420 (2006) ·Zbl 1111.65042号
[24] Hénon，P。；Ramet，P。；Roman，J.，关于寻找有效ILU（k）因式分解的近似超节点，并行计算。，34, 345-362 (2008)
[25] Hysom，D。；Pothen，A.，不完全因子预处理的可扩展并行算法，SIAM J.Sci。计算。，2194-2215年6月22日（2001年）·Zbl 0986.65048号
[26] Hysom，D。；Pothen，A.，《基于水平的不完全LU分解：图形模型和算法》（2002年11月），技术代表UCRL-JC-150789
[27] 琼斯，M.T。；Plassmann，P.E.，改进的不完全Cholesky因子分解，ACM Trans。数学。软质。，21, 1, 5-17 (1995) ·Zbl 0886.65024号
[28] Karypis，G。；Kumar，V.，《划分不规则图的快速高质量多级方案》，SIAM J.Sci。计算。，20, 1, 359-392 (1998) ·兹比尔0915.68129
[29] Kershaw，D.S.，线性方程组迭代解的不完全Cholesky共轭梯度法，J.Compute。物理。，26, 43-65 (1978) ·Zbl 0367.65018号
[30] Kurzak，J。；盖茨，M。；亚尔汗，A。；山崎，I。；吴，P。；Luszczek，P。；Finney，J。；Dongarra，J.，并行BLAS性能报告（2018年），ICL-UT-18-01（2018年4月）
[31] 拉萨尔，D。；Karypis，G.，多线程图划分（2013），明尼苏达大学计算机科学与工程系，明尼阿波利斯明尼苏打大学计算机科学和工程系，技术代表。
[32] 李，N。；萨阿德，Y。；Chow，E.，通用稀疏矩阵ILU的Crout版本，SIAM J.Sci。计算。，25, 2, 716-728 (2004) ·Zbl 1042.65025号
[33] Li，X.S.，《超级逻辑单元概述：算法、实现和用户界面》，ACM Trans。数学。软质。，31, 3, 302-325 (2005) ·Zbl 1136.65312号
[34] Li，X.S。；Shao，M.，部分枢轴不完全LU因式分解的超节点方法，ACM-Trans。数学。软质。，37, 4, 43:1-43:20 (2011) ·Zbl 1365.65078号
[35] 林，C.-J。；Moré，J.J.，记忆有限的不完全Cholesky因子分解，SIAM J.Sci。计算。，21, 1, 24-45 (1999) ·兹伯利0941.65033
[36] Maurer，H。；Mittelmann，H.，解决具有控制和状态约束的椭圆控制问题的优化技术：第1部分。边界控制，计算。最佳方案。申请。，16, 1, 29-55 (2000) ·Zbl 0961.93026号
[37] Ng，E.G。；佩顿，B.W。；Raghavan，P.，层次记忆计算机的分块不完全Cholesky预条件，（科学计算中的迭代方法IV，科学计算的迭代方法Ⅳ，计算和应用数学中的IMACS系列（1999），IMACS），211-221
[38] Olschowka，M。；Neumaier，A.，《高斯消去的新枢轴策略》，线性代数应用。，240, 131-151 (1996) ·Zbl 0852.65021号
[39] Saad，Y.，《寻找ILU预处理的精确和近似嵌段结构》，SIAM J.Sci。计算。，24, 4, 1107-1123 (2003) ·Zbl 1034.65020号
[40] Saad，Y.，《稀疏线性系统的迭代方法》（2003），SIAM出版物：费城SIAM出版物·Zbl 1002.65042号
[41] 萨阿德，Y。；Schultz，M.H.，GMRES:求解非对称线性系统的广义最小残差算法，SIAM J.Sci。统计计算。，7, 856-869 (1986) ·Zbl 0599.65018号
[42] 萨阿德，Y。；Suchomel，B.J.，ARMS：一般稀疏线性系统的代数递归多级解算器，Numer。线性代数应用。，9, 359-378 (2002) ·Zbl 1071.65001号
[43] 申克，O。；Gärtner，K.，用PARDISO求解非对称稀疏线性方程组，J.Futur。通用。计算。系统。，20, 3, 475-487 (2004)
[44] 申克，O。；Gärtner，K.，关于对称不定系统的快速因式分解旋转方法，Electron。事务处理。数字。分析。，23, 1, 158-179 (2006) ·Zbl 1112.65022号
[45] 申克，O。；Röllin，S。；Gupta，A.，半导体器件和电路模拟中非对称矩阵排列和缩放的影响，IEEE Trans。计算-辅助设计。集成。电路系统。，23, 3, 400-411 (2004)
[46] 申克，O。；瓦希特，A。；Hagemann，M.，基于Matching的预处理算法，用于解决大规模非凸内部点优化中的鞍点问题，计算。最佳方案。申请。，36, 2, 321-341 (2007) ·Zbl 1146.90055号
[47] 申克，O。；博勒费尔，M。；Römer，R.A.，SIGEST：关于Anderson本地化模型的大规模对角化技术，SIAM Rev.，50，91-112（2008）·Zbl 1136.65044号
[48] 申克，O。；瓦希特，A。；Weiser，M.，《大规模非凸约束优化的惯性揭示预处理》，SIAM J.Sci。计算。，31, 2, 939-960 (2008) ·Zbl 1194.35029号
[49] 申克，O。；瓦希特，A。；Weiser，M.，《大规模非凸约束优化的惯性释放预处理》，SIAM J.Sci。计算。，1993年2月31日至960年（2009年）·Zbl 1194.35029号
[50] 斯科特·J·A。；Tůma，M.，结构在不完全因子分解预条件中的重要性，BIT-Numer。数学。，51, 2, 385-404 (2011) ·Zbl 1221.65090号
[51] 瓦希特，A。；Biegler，L.T.，《关于大规模非线性规划的原始-对偶内点过滤器线搜索算法的实现》，数学。程序。，106，1，25-57（2006年）·Zbl 1134.90542号
[52] Watts，J.W.，油藏模拟压力方程迭代解的共轭梯度截断直接法，Soc.Pet。《工程师杂志》，21，345-353（1981）
[53] Xi，Y。；李，R。；Saad，Y.，《稀疏对称矩阵的低秩修正代数多层预条件器》，SIAM J.Matrix Anal。申请。，37, 1, 235-259 (2016) ·Zbl 1376.65036号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。