×

CARP-CG:用于强对流占优PDE的线性系统的鲁棒高效并行求解器。 (英语) Zbl 1195.65062号

小结:CARP-CG是CARP的共轭梯度(CG)加速,作者将其作为稀疏线性系统的鲁棒块并行方案介绍。CARP在块内执行Kaczmarz(KACZ)行投影,通过平均每个组件来自不同块的更新值,合并来自单独块的结果。平均运算等价于某些超空间中的某些KACZ行投影序列(“平均投影”),因此CARP等价于该超空间中具有循环松弛参数的KACZ。CARP的CG-加速是基于以下(顺序)CGMN算法的推广A.比约克T.Elfving公司[BIT,Nord.Tidskr.Inf.-behandl.19,145-1-63(1979;Zbl 0409.65022号)]它通过使用固定的松弛参数对线性系统的方程进行双重扫描来加速KACZ。CGMN被推广为允许循环松弛参数,因此所得到的方法,称为CGMNC,可以应用于超空间。超空间中的平均投影可以按任意顺序进行,因此可以在规则空间中使用CARP在双扫描中实现超空间中CGMNC。由此产生的算法CARP-CG与CARP一样健壮,但收敛速度明显更快。在由对流占优的椭圆偏微分方程导出的一些刚性线性系统上,将CARP-CG与重新启动的GMRES、Bi-CGSTAB和CGS进行了比较。结果表明,CARP-CG具有很强的鲁棒性,其运行时间与其他方法相比具有很强竞争力。还测试了CGNR的缩放版本,它与CARP-CG一样健壮,但速度较慢。

MSC公司:

65英尺50英寸 稀疏矩阵的计算方法
35J99型 椭圆方程和椭圆系统
65层20 超定系统伪逆的数值解
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Arioli,M。;达夫,I.S。;Noailles,J。;Ruiz,D.:稀疏矩阵的块投影方法,SIAM科学与统计计算杂志13,47-70(1992)·Zbl 0761.65015号
[2] Arioli,M。;达夫,I.S。;鲁伊斯,D。;Sadkane,M.:Block Lanczos加速Block cimmino方法的技术,SIAM科学与统计计算期刊161478-1111(1995)·Zbl 0839.65037号 ·doi:10.1137/0916086
[3] R.Barrett,M.Berry,T.F.Chan,J.Demmel,J.Donato,J.Dongarra,V.Eijkhout,R.Pozo,C.Romine,H.A.van der Vorst,《线性系统解的模板:迭代方法的构建块》,第二版,SIAM,宾夕法尼亚州费城,1994年。可从以下网址获得:<http://netlib2.cs.utk.edu/linalg/html_templates/templates.html>. ·Zbl 0814.65030号
[4] &焦虑;比约克。;Elfving,T.:计算线性方程组伪逆解的加速投影方法,Bit 19,145-163(1979)·Zbl 0409.65022号
[5] BLAS–基本线性代数子程序。可从以下网址获得:<http://www.netlib.org/blas>.
[6] R.Bramley,H.-C.Chen,U.Meier,A.Sameh,《关于一些并行预条件CG方案》,载:O.Axelsson,L.Yu。Kolotilina(编辑),《预条件共轭梯度方法国际会议论文集》,荷兰奈梅亨,1989年6月,《数学讲义》,第1457卷,施普林格,柏林,1990年,第17-27页·兹比尔0713.65028
[7] 布拉姆利,R。;Sameh,A.:平行行投影算法的区域分解,应用数值数学8303-315(1991)·Zbl 0738.65022号 ·doi:10.1016/0168-9274(91)90071-7
[8] 布拉姆利,R。;Sameh,A.:大型非对称线性系统的行投影方法,SIAM科学与统计计算杂志13,168-193(1992)·Zbl 0752.65024号 ·doi:10.1137/0913010
[9] Cimmino,G.:Calcolo approssimato per le soluzioni dei sistemi di equazioni lineari,《科学学报》第十六卷第二辑,anno IX 1,326-333(1938)·Zbl 0018.41802号
[10] 艾格蒙特,P.P.B。;Herman,G.T。;Lent,A.:大型分块线性系统的迭代算法及其在图像重建中的应用,线性代数及其应用40,37-67(1981)·Zbl 0466.65021号 ·doi:10.1016/0024-3795(81)90139-7
[11] Elfving,T.:一致和不一致线性方程的块迭代方法,数值数学35,1-12(1980)·Zbl 0416.65031号 ·doi:10.1007/BF01396365
[12] H.Elman,大型稀疏非对称线性方程组的迭代方法。技术报告229,耶鲁大学计算机科学系,美国康涅狄格州纽黑文,1982年。
[13] Elman,H。;Golub,G.:循环约化非自伴线性系统的迭代方法,《计算数学》54,第190期,第671-700页(1990年)·Zbl 0699.65021号 ·doi:10.2307/2008506
[14] 费尔南德斯,J.-J。;戈登,D。;Gordon,R.:《电子层析成像迭代重建算法的高效并行实现》,《并行与分布式计算杂志》68,第5期,626-640(2008)
[15] Ferziger,J.H。;Perić,M.:流体动力学的计算方法,(2002)·Zbl 0998.76001号
[16] Fletcher,R.:不定系统的共轭梯度法,数学讲义50673-89(1976)·Zbl 0326.65033号
[17] 弗洛伊德,R。;Nachtigal,N.:QMR:非厄米线性系统的准最小残差法,《数值数学》60,315-339(1991)·Zbl 0754.65034号 ·doi:10.1007/BF01385726
[18] 戈登,D。;Gordon,R.:《分量平均行投影:稀疏线性系统的稳健块并行方案》,SIAM科学计算杂志27,1092-1117(2005)·Zbl 1093.65033号 ·数字对象标识代码:10.1137/040609458
[19] 戈登,D。;Gordon,R.:CGMN重访:从椭圆偏微分方程导出的刚性线性系统的稳健和有效解,ACM数学软件事务35,No.3,18:1-18:27(2008)
[20] 戈登,D。;Gordon,R.:具有大非对角元素和不连续系数的线性系统的求解方法,《工程与科学中的计算机建模》53,第1期,23-45(2009)·Zbl 1231.65065号 ·doi:10.3970/cmes.2009.053.023
[21] D.Gordon,R.Gordon。几何尺度:一些不连续系数线性系统的简单有效预条件,计算与应用杂志;《应用数学》,出版。doi:10.1016/j.cam.2010.05.021·Zbl 1196.65066号 ·doi:10.1016/j.cam.2010.05.021
[22] Gresho,P。;Sani,R.:不可压缩流和有限元方法:平流扩散和等温层流,(1998)·Zbl 0941.76002号
[23] Hestenes,M.R。;Stiefel,E.:求解线性系统的共轭梯度方法,国家标准局研究杂志49,409-436(1952)·Zbl 0048.09901号
[24] Kaczmarz,S.:Angenäherte auflösung von systemen linearer gleichungen,《波兰科学公报》35,355-357(1937)·兹标0017.31703
[25] 卡马特,C。;Sameh,A.:在多处理器上求解非对称线性系统的投影方法,并行计算9,第3期,291-312(1989)·Zbl 0668.65028号 ·doi:10.1016/0167-8191(89)90114-2
[26] 卡马特,C。;Weeratunga,S.:在共享内存多处理器上实现两种投影方法:DEC VAX 6240,并行计算16,375-382(1990)·Zbl 0713.65025号 ·doi:10.1016/0167-8191(90)90075-K
[27] 卡马特,C。;Weeratunga,S.:非自伴椭圆偏微分方程数值解的投影方法,偏微分方程的数值方法8,59-76(1992)·Zbl 0752.65070号 ·doi:10.1002/num.1690080104
[28] Kincaid,D。;Young,D.:将为对称系统开发的迭代算法应用于非对称系统,椭圆问题求解器,353-359(1981)·Zbl 0467.65015号
[29] Kuck博士。;戴维森,E。;Lawrie博士。;Sameh,A.:《今日并行超级计算与雪松方法》,《科学》231967-974(1986)
[30] Lanczos,C.:通过最小化迭代求解线性方程组,国家标准局研究期刊49,33-53(1952)
[31] Meurant,G.:大型线性系统的计算机解,(1999)·Zbl 0934.65032号
[32] Saad,Y.:共轭梯度法多项式预处理的实际应用,SIAM科学与统计计算杂志6,865-881(1985)·Zbl 0601.65019号 ·doi:10.1137/0906059
[33] Y.Saad,SPARSKIT:稀疏矩阵计算的基本工具包。技术报告RIACS-90-20,高级计算机科学研究所,NASA艾姆斯研究中心,加利福尼亚州莫菲特菲尔德,1990年。
[34] Saad,Y.:稀疏线性系统的迭代方法,(2003)·兹比尔1031.65046
[35] 萨阿德,Y。;Schultz,M.H.:GMRES:求解非对称线性系统的广义最小残差算法,SIAM科学与统计计算杂志7856-869(1986)·Zbl 0599.65018号 ·doi:10.1137/0907058
[36] Sonneveld,P.:CGS:非对称线性系统的快速Lanczos型解算器,SIAM科学与统计计算杂志10,36-52(1989)·Zbl 0666.65029号 ·数字对象标识代码:10.1137/0910004
[37] R.S.Tuminaro、M.A.Heroux、S.A.Hutchinson、J.N.Shadid、AZTEC用户指南。技术报告SAND99-8801J,新墨西哥州阿尔伯克基市桑迪亚国家实验室,1999年。
[38] Van Der Vorst,H.A.:Bi-CGSTAB:非对称线性系统解的Bi-CG快速平滑收敛变体,SIAM科学与统计计算杂志13,631-644(1992)·Zbl 0761.65023号 ·doi:10.1137/0913035
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。