×
奥米德·卡达尼;安德烈·利亚明;克里斯蒂安·卡本霍夫特

GPU加速混合预处理共轭梯度方法在计算地质力学中的大型三维问题中的应用。 （英文） Zbl 1443.65045号

计算。数学。申请。 69，第10期，1114-1131（2015）.
摘要：本文提出了一种共轭梯度法的混合预处理技术，并讨论了它在图形处理单元（GPU）上的并行实现用于在计算地质力学的非线性有限元极限分析（FELA）的背景下，解决因将内点方法应用于圆锥优化问题而产生的大型稀疏线性系统。对于大型3D问题，由于内存需求和计算时间呈指数级增长，通常使用直接求解器的成本非常昂贵。此外，使用优化框架产生的所谓鞍点系统也不例外。另一方面，尽管预处理迭代方法具有适度的存储需求，因此可以应用于比直接方法大得多的问题，但它们通常会表现出大量迭代以达到收敛。在本文中，我们证明了使用所提出的混合预处理器以及GPU上的详细实现可以有效地解决这个问题。此外，数值结果验证了该方法的鲁棒性和有效性。

引用于三文件

MSC公司：

65层10 线性系统的迭代数值方法
65F08个 迭代方法的前置条件
65千5 数值数学规划方法
65K10码 数值优化和变分技术
65日元10 特定类别体系结构的数值算法
2008年6月 地球物理问题的计算方法

关键词：

近似逆预处理器;不完全Cholesky因子分解;极限分析;预处理共轭梯度法;圆锥规划;图形处理单元
PDF格式BibTeX公司 XML格式引用
\textit{O.Kardani}等人，《计算》。数学。申请。69，第10号，1114--1131（2015；Zbl 1443.65045）
全文： DOI程序

参考文献：

[1] Chen，W.F.，极限分析和土壤塑性（1975），Elsevier·Zbl 0354.73033号
[2] Chuang，P.H.，用线性规划进行地质力学稳定性分析。一： 公式，J.Geotech。工程，1181696-1715（1992）
[3] Chuang，P.H.，用线性规划进行地质力学稳定性分析。二： 应用，J.Geotech。工程师，1181716-1726（1992）
[4] Sloan，S.W.，《使用有限元和线性规划的上限分析》，《国际数值杂志》。分析。方法地质技术。，13, 263-282 (1989) ·Zbl 0687.73040号
[5] Krabbenhöft，K。；Lyamin，A.V.，使用二阶锥编程计算凸轮粘土塑性，计算。方法应用。机械。工程，209239-249（2012）·兹比尔1243.74015
[6] Krabbenhøft，K。；Lyamin，A.V。；Sloan，S.W.，一些塑性问题的圆锥规划公式和求解，国际。《固体结构杂志》，44，1533-1549（2007）·Zbl 1119.74035号
[7] Krabbenhöft，K。；Lyamin，A.V。；黄，J.，使用数学规划方法的颗粒接触动力学，计算。方法应用。机械。工程，43，165-176（2012）
[8] 安德森，E.D。；Roos，C。；Terlaky，T.，关于实现二次曲线二次优化的原对偶内点方法，数学。程序。，95, 249-277 (2003) ·Zbl 1030.90137号
[9] Wright，S.J.，原始对偶内点方法（1997），SIAM:SIAM费城·Zbl 0863.65031号
[10] O.卡达尼。；Lyamin，A.V。；Krabbenhöft，K.，《有限元极限分析中大型稀疏系统预处理技术的比较研究》，IAENG Int.J.Appl。数学。，43, 195-203 (2013) ·Zbl 1512.74094号
[11] Krabbenhöft，K。；Lyamin，A.V。；Sloan，S.W.，《使用半定规划的三维Mohr-Coulomb极限分析》，Comm.Numer。方法工程，241107-1119（2008）·Zbl 1154.74050号
[12] 卢斯蒂格，I.J。；Marsten，R.E。；Shanno，D.F.，《线性规划的内点方法：计算现状》，ORSA J.Compute。，6, 1-14 (1994) ·Zbl 0798.90100号
[13] Mehrotra，S.，关于原对偶内点方法的实现，SIAM J.Optim。，2, 575-601 (1992) ·Zbl 0773.90047号
[14] Gondzio，J.，线性规划原对偶方法中的多重中心性校正，计算。最佳方案。申请。，6, 137-156 (1996) ·Zbl 0860.90084号
[15] Saad，Y.，《稀疏线性系统的迭代方法》（2003），SIAM:SIAM Philadelphia·Zbl 1002.65042号
[16] Benzi，M.，《大型线性系统的预处理技术：综述》，J.Comput。物理。，182, 418-477 (2002) ·Zbl 1015.65018号
[17] 奥利维拉，A.R.L。；Sorensen，D.C.，线性规划内点方法中的一类新的大规模线性系统预处理器，线性代数应用。，394, 1-24 (2005) ·Zbl 1071.65088号
[18] Bocanegra，S。；坎波斯，E.F。；Oliveira，A.R.L.，使用混合预处理程序求解由内点方法产生的大规模线性系统，计算。最佳方案。申请。，36149-164（2007年）·Zbl 1148.90350号
[19] Fontova，M.I.V。；奥利维拉·A.R.L.D。；Campos，F.F.，内点方法混合预条件实现的启发式，Pesqui。操作。，31, 3, 579-591 (2011)
[20] Bocanegra，S。；卡斯特罗，J。；Oliveira，A.R.，《用混合预条件器改进大型块角问题的内点方法》，欧洲期刊Oper。第231号决议，第263-273号决议（2013年）·Zbl 1317.90195号
[21] 小岛，M。；瑞穗，S。；Yoshise，A.，线性规划的原对偶内点算法，（Megiddo，N.，《数学规划进展：内点和相关方法》（1989），Springer:Springer New York），29-47·Zbl 0708.90049号
[22] 内斯特罗夫，Y。；Nemirovsky，A.，（凸规划中的内点多项式方法。凸规划中内点多项式的方法，SIAM应用数学和数值数学研究（1994），SIAM:SIAM Philadelphia）·Zbl 0824.90112号
[23] Field，M.R.，《提高分解稀疏近似逆预处理器的性能》（技术报告HDL-TR-98-199）（1998年），日立都柏林实验室：爱尔兰都柏林日立都林实验室
[24] Benzi，M。；Tuma，M.，稀疏近似逆预条件的比较研究，应用。数字。数学。，30, 305-340 (1999) ·Zbl 0949.65043号
[25] Benzi，M。；库卢姆，J.K。；Tuma，M.，共轭梯度法的稳健近似逆预处理，SIAM J.Sci。计算。，22, 1318-1332 (2000) ·Zbl 0985.65035号
[26] Bollhofer，M。；Saad，Y.，带枢轴的因子化近似逆预条件，SIAM J.矩阵分析。申请。，23, 692-705 (2002) ·Zbl 0999.65027号
[27] Huckle，T.，预处理和平滑的分解稀疏近似逆，Selçuk J.Appl。数学。，1, 63 (2000) ·兹比尔1011.65018
[29] 哈尔琴科，S.A。；Kolotilina，L.Y。；Nikishin，A.A。；Yeremin，A.Y.，一种稳健的AINV型方法，用于构造因子形式的稀疏近似逆预条件，Numer。线性代数应用。，8, 165-179 (2001) ·兹比尔1051.65057
[30] Benzi，M。；Tuma，M.，非对称线性系统的稀疏近似逆预条件，SIAM J.Sci。计算。，19, 968-994 (1998) ·Zbl 0930.65027号
[31] O.卡达尼。；Lyamin，A.V。；Krabbenhöft，K.，计算地质力学中出现的大型稀疏和高度病态系统的并行预处理共轭梯度法，国际计算杂志。科学。工程（2014）
[32] Chen，K.，《矩阵预处理技术与应用》（2005），剑桥大学出版社·兹比尔1079.65057
[33] Meijerink，J.A。；van der Vorst，H.A.，系数矩阵为对称M矩阵的线性系统的迭代求解方法，数学。计算。，31, 148-162 (1977) ·Zbl 0349.65020号
[34] Tismenetsky，M.，求解大型稀疏线性系统的新预处理技术，线性代数应用。，154, 331-353 (1991) ·Zbl 0732.65030号
[35] 卡菲里，S。；达帕佐，M。；马里诺，M。；Mucherino，A。；Toraldo，G.，带约束的大规模二次规划问题的内点求解器，J.Optim。理论应用。，129, 55-75 (2006) ·Zbl 1139.90021号
[36] 卡菲里，S。；达帕佐，M。；德西蒙，V。；di Serafino，D.，关于大型二次型问题势约简软件中KKT系统的迭代解，计算。最佳方案。申请。，38, 27-45 (2007) ·Zbl 1171.90548号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。