×

PetRBF-一种并行的高斯径向基函数插值算法。 (英语) 兹比尔1231.65026

摘要:我们开发了一种用于径向基函数(RBF)插值的并行算法,该算法具有(O(N))复杂性,需要(O(N))存储,并且可扩展到上千个进程。该算法使用gmres迭代解算器和限制加性Schwarz方法(RASM)作为预条件和快速矩阵向量算法。以前的快速rbf方法是使用多二次基函数和多谐基函数开发的,最多可达到O(N(log N))的复杂度。相比之下,本方法使用的高斯函数与域的方差较小,但有足够的重叠。这是我们的主要目标应用程序流体模拟粒子方法中的常见选择。高斯基函数的快速衰减允许迭代解算器快速收敛,即使rasm中的子域非常小。同时表明插补精度可以达到机床精度。目前的方法是使用petsc库(开发人员版本)并行实现的。数值实验证明了它在1024个Blue Gene/L处理器(700 MHz PowerPC处理器)上以106 s(19次迭代,误差容限为\(10^{-15}\))的时间处理具有超过5000万个数据点的rbf插值问题中的能力。并行代码在开源模型中是免费的。

MSC公司:

65D05型 数值插值
2005年5月 并行数值计算
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 南卡罗莱纳州阿特卢里。;Zhu,T.,计算力学中的一种新的无网格局部Petrov-Galerkin(MLPG)方法,计算。机械。,22, 117-127 (1998) ·兹比尔0932.76067
[2] 巴布斯卡,I。;Melenk,J.M.,单位分割法,国际数字杂志。方法工程,40,727-758(1997)·Zbl 0949.65117号
[3] 巴莱,S。;Buschelman,K。;埃伊霍特,V。;格罗普,W.D。;考希克,D。;Knepley,M。;Curfman-McInnes,L。;B.F.史密斯。;Zhang,H.,(PETSc用户手册,技术报告ANL-95/11-3.0.0版(2008),阿贡国家实验室)
[4] Barba,L.A.,《计算高雷诺数旋涡流:采用完全无网格公式的高精度方法》,(并行计算流体动力学-多学科应用(2005),Elsevier B.V),305-312
[5] 巴巴,洛杉矶。;Leonard,A.,《网络循环初始条件下三角涡的出现和演化》,《物理学》。流体,19,1,017101(2007)·Zbl 1146.76319号
[6] 巴巴,洛杉矶。;伦纳德(A.Leonard)。;Allen,C.B.,《粘性涡方法的进展——基于径向基函数插值的无网格空间自适应》,国际期刊《数值方法》。流体,47,5,387-421(2005)·Zbl 1085.76052号
[7] 巴巴,洛杉矶。;Rossi,L.F.,《粒子方法的全局场插值》,J.Comp。物理。,229, 4, 1292-1310 (2010) ·Zbl 1329.76282号
[8] Beatson,R.K。;Cherrie,J.B。;Mouat,C.T.,径向基函数的快速拟合:基于预处理GMRES迭代的方法,Adv.Comp。数学。,11、2-3、253-270(1999年11月)
[9] Beatson,R.K。;轻型,W.A。;Billings,S.,径向基函数插值方程的快速解:区域分解方法,SIAM J.Sci。计算。,22, 5, 1717-1740 (2000) ·Zbl 0982.65015号
[10] Beatson,R.K。;Newsam,G.N.,径向基函数的快速评估:I,Comp。数学。应用。,24, 12, 7-19 (1992) ·Zbl 0765.65021号
[11] Beatson,R.K。;Powell,M.J.D.,《采用拉格朗日函数近似的薄板样条插值迭代方法》(Griffiths,D.F.;Watson,G.A.,《数值分析》(1993),Longman Sci。技术),17-39·兹伯利0795.65004
[12] Belytschko,T。;吕义勇。;Gu,L.,无元素伽辽金方法,国际数值杂志。方法工程,37,2,229-256(1994)·Zbl 0796.73077号
[13] 比林斯,S.D。;Beatson,R.K。;Newsam,G.N.,使用连续全球表面的地球物理数据插值,地球物理学,67,61810-1822(2002)
[14] Buhmann,M.D.,径向基函数。《理论与实现》(2003),剑桥大学出版社·Zbl 1038.41001号
[15] 蔡,X.C。;Sarkis,M.,一般稀疏线性系统的限制加性Schwarz预条件,SIAM J.Sci。计算。,21, 792-797 (1997) ·Zbl 0944.65031号
[16] 卡尔,J.C。;弗里特·W·R。;Beatson,R.K.,医学成像用径向基函数曲面插值,IEEE Trans。医学图像。,16, 1, 96-107 (1997)
[17] Chen,W.,新RBF配置方案和核RBF及其应用,Lect。注释计算。科学。工程,26,75-86(2002)·Zbl 1016.65094号
[18] Cherrie,J.B。;Beatson,R.K。;Newsam,G.N.,《径向基函数的快速计算:(R^N)中广义多二次曲面的方法》,SIAM J.Sci。计算。,23, 5, 1549-1571 (2002) ·Zbl 1009.65007号
[19] Duarte,C.A。;Tinsley Oden,J.,使用云的h-p自适应方法,计算。方法应用。机械。工程,139,1-4,237-262(1996)·Zbl 0918.73328号
[20] Efstathiou,大肠杆菌。;Gander,M.J.,《为什么限制加法Schwarz比加法Schvarz收敛得更快》,BIT-Numer。数学。,43, 945-959 (2003) ·Zbl 1045.65027号
[21] Fasshauer,G.E.,《用径向基函数配置求解偏微分方程》,(Le Mehaute,A.;Rabut,C.;Schumaker,L.L.,《曲面拟合和多分辨率方法》(1997),范德比尔特大学出版社),131-138·Zbl 0938.65140号
[22] Fasshauer,G.E.,《利用MATLAB的无网格近似方法》(2007年),《世界科学》·Zbl 1123.65001号
[23] Faul,A.C.公司。;Goodsell,G。;Powell,M.J.D.,《多维多重二次插值的Krylov子空间算法》,IMA J.Numer。分析。,2005年1月25日至24日·Zbl 1070.65006号
[24] Franke,R.,《分散数据插值:一些方法的测试》,数学。公司。,38, 157, 181-200 (1982) ·Zbl 0476.65005号
[25] Greengard,L。;Strain,J.,快速高斯变换,SIAM J.Sci。统计师。计算。,12、1、79-94(1991年1月)
[26] Gumerov,N.A。;Duraiswami,R.,双调和方程的快速多极方法,(技术报告UMIACS-TR-2005-29,马里兰大学计算机科学系技术报告(2005))·Zbl 1103.65122号
[27] Gumerov,N.A。;Duraiswami,R.,通过预处理Krylov迭代进行快速径向基函数插值,SIAM J.Sci。计算。,29, 5, 1876-1899 (2007) ·Zbl 1154.65303号
[28] Hardy,R.L.,多重二次双调和方法的理论和应用,计算。数学。申请。,19, 8/9, 163-208 (1990) ·Zbl 0692.65003号
[29] 休斯·T·J·R。;Cottrell,J.A。;Bazilevs,Y.,等几何分析:CAD,有限元,NURBS,精确几何和网格细化,计算。方法应用。机械。工程,194,39-41,4135-4195(2005)·Zbl 1151.74419号
[30] Ingber,M.S。;Chen,C.S。;Tanski,J.A.,《使用径向基函数和并行区域分解求解三维扩散方程的无网格方法》,《国际数值》。方法工程,60,2183-2201(2004)·Zbl 1178.76276号
[31] Kansa,E.J.,《多元二次曲面——一种应用于计算流体动力学的离散数据近似方案》,I.《曲面近似和偏导数估计》,《计算》。数学。申请。,19, 8/9, 127-145 (1990) ·Zbl 0692.76003号
[32] Kansa,E.J.,《多元二次曲面——一种离散数据近似方案及其在计算流体动力学中的应用》,第二卷。抛物型、双曲型和椭圆型偏微分方程的解,计算。数学。申请。,19, 8/9, 147-161 (1990) ·Zbl 0850.76048号
[33] Kansa,E.J。;Hon,Y.C.,用多二次径向基函数规避病态调节问题:椭圆偏微分方程的应用,计算。数学。申请。,39123-137(2000年)·Zbl 0955.65086号
[34] Leonard,A.,《流动模拟的旋涡方法》,J.Comp。物理。,37, 289-335 (1980) ·Zbl 0438.76009号
[35] 李,J。;陈,Y。;Pepper,D.,一维和二维地下水污染物运移建模的径向基函数法,计算。机械。,23, 10-15 (2003) ·Zbl 1045.76548号
[36] 李,J。;Hon,Y.C.,径向基无网格方法的区域分解,数值。方法部分差异。Equ.、。,20, 3, 450-462 (2004) ·Zbl 1048.65124号
[37] Ling,L。;Kansa,E.J.,用域分解方法对径向基函数进行预处理,数学。计算。型号。,40, 1413-1427 (2004) ·Zbl 1077.41008号
[38] Ling,L。;Kansa,E.J.,径向基函数配置方法的最小二乘预处理,高级计算。数学。,23, 31-54 (2005) ·Zbl 1067.65136号
[39] 刘伟凯。;S·6月。;李,S。;Adee,J。;Belytschko,T.,《结构动力学的再生核粒子方法》,国际数值杂志。方法工程师,38,101165-1679(1995)·Zbl 0840.73078号
[40] 马迪奇,W.R。;Nelson,S.A.,多元插值和条件正定函数:II,数学。公司。,54, 189, 211-230 (1990) ·Zbl 0859.41004号
[41] Mai-Duy,N。;Tran-Cong,T.,微分方程数值解的薄板样条间接RBFN方法,计算。模型。工程师科学。,4, 1, 85-102 (2003) ·Zbl 1148.76351号
[42] Michelli,C.A.,《散乱数据的插值:距离矩阵和条件正定函数》,Constr。约211-22(1986年)·Zbl 0625.41005号
[43] Powell,M.J.D.,薄板样条快速评估的截断Laurent展开,数值。算法,599-120(1993)·Zbl 0790.65010号
[44] Roussos,G。;Baxter,B.J.C.,径向基函数的快速评估,J.Compute。申请。数学。,180, 1, 51-70 (2005) ·Zbl 1069.65008号
[45] Saad,Y.,《稀疏线性系统的迭代方法》(2003),SIAM·Zbl 1002.65042号
[46] B.F.史密斯。;比约斯塔德,体育。;Gropp,W.D.,区域分解,椭圆偏微分方程的并行多层方法(1996),剑桥大学出版社·兹比尔0857.65126
[47] 托雷斯,C.E。;Barba,L.A.,通过定位和迭代使用高斯函数进行快速径向基函数插值,J.Comp。物理。,228, 4976-4999 (2009) ·Zbl 1168.65318号
[48] Wong,S.M。;尊敬的Y.C。;Li,T.S。;Chung,S.L.,使用多二次格式模拟时间相关问题的多区域分解,计算。数学。申请。,37, 23-43 (1999) ·Zbl 0951.76066号
[49] 横田,R。;Sheel,T.K。;Obi,S.,使用纯拉格朗日涡方法计算各向同性湍流,J.Comp。物理。,226, 1589-1606 (2007) ·Zbl 1121.76046号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。