×

参数偏微分方程的交替最小二乘TT-交叉混合算法。 (英语) Zbl 1418.65170号

小结:我们考虑使用低秩张量列(TT)分解的参数偏微分方程的近似解。例如,此类参数PDE出现在工程应用中的不确定性量化问题中。我们提出了一种交替最小二乘法和TT交叉法的混合算法。它计算整个解的TT近似值,这在寻求多个感兴趣的量时是有益的。例如,通过最大熵方法计算概率密度函数时可能需要这样做[M.卡维拉德M.约瑟夫,“数字通信系统性能评估中的最大熵和矩方法”,IEEE Trans。Comm.34,No.12,1183–1189(1986;doi:10.1109/TCOM.1986.1096484)]. 新算法利用并保留了随机配置方案中离散化算子的块对角结构。这就解开了TT表示中空间自由度和参数自由度的计算。特别是,它只需要在几个参数值下求解独立的偏微分方程,从而允许使用现有的高性能偏微分方程求解器。在我们的数值实验中,我们将新算法应用于随机扩散方程,并将其与TT格式的预处理最速下降法、(多级)拟蒙特卡罗和维数自适应稀疏网格方法进行了比较。对于足够平滑的随机场,新方法的速度要快几个数量级。

MSC公司:

65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65楼30 其他矩阵算法(MSC2010)
15A69号 多线性代数,张量微积分
65立方米 随机微分和积分方程的数值解
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] I.Babuška,F.Nobile,and R.Tempone,{带随机输入数据的椭圆偏微分方程的随机配置方法},SIAM J.Numer。分析。,45(2007),第1005-1034页·Zbl 1151.65008号
[2] M.Bachmayr、A.Cohen和W.Dahmen,{参数偏微分方程:稀疏或低阶近似?},IMA J.Numer。分析。,38(2018),第1661-1708页·Zbl 1477.65189号
[3] J.Ballani和L.Grasedyck,参数相关PDEs输出量的层次张量近似,SIAM/ASA J.Uncertain。数量。,3(2015),第852-872页·Zbl 1327.65010号
[4] J.Ballani和D.Kressner,《简化基方法:从低秩矩阵到低秩张量》,SIAM J.Sci。计算。,38(2016),第A2045-A2067页·Zbl 1382.15036号
[5] M.Barrault,Y.Maday,N.C.Nguyen和A.T.Patera,{一种“经验插值”方法:应用于偏微分方程的有效降基离散化},C.R.Math。阿卡德。科学。巴黎。,339(2004),第667-672页·Zbl 1061.65118号
[6] M.Bebendorf,{边界元矩阵的近似},数值。数学。,86(2000),第565-589页·Zbl 0966.65094号
[7] M.Bebendorf,{多元函数的自适应交叉逼近},Const。约,34(2011),第149-179页·Zbl 1248.41049号
[8] M.Bebendorf和C.Kuske,{函数生成高阶张量的变量分离},科学杂志。计算。,61(2014),第145-165页·兹比尔1307.65052
[9] M.Bieri和C.Schwab,{椭圆sPDEs}的稀疏高阶有限元,计算。方法应用。机械。工程,198(2009),第1149-1170页·Zbl 1157.65481号
[10] C.Bierig和A.Chernov,{用多级蒙特卡罗最大熵方法逼近概率密度函数},J.Compute。物理。,314(2016),第661-681页·Zbl 1349.65007号
[11] H.-J.Bungartz和M.Griebel,《稀疏网格》,《数值学报》。,13(2004),第147-269页·Zbl 1118.65388号
[12] J.Charrier和A.Debussche,具有对数正态系数的椭圆偏微分方程的弱截断误差估计,Stoch。部分差异。埃克。分析。计算。,1(2013年),第63-93页·Zbl 1273.65166号
[13] A.Cohen,R.DeVore,and C.Schwab,{一类椭圆sPDEs}的最佳N项Galerkin逼近的收敛速度。计算。数学。,10(2010),第615-646页·Zbl 1206.60064号
[14] J.Dick、F.Y.Kuo、Q.T.Le Gia、D.Nuyens和C.Schwab,{带随机场输入的仿射参数算子方程的高阶QMC Petrov-Galerkin离散化},SIAM J.Numer。分析。,52(2014),第2676-2702页·Zbl 1326.65013号
[15] S.Dolgov、K.Anaya-Izquierdo、C.Fox和R.Scheichl,《张量列分解中多元概率分布的近似和抽样》,预印本,2018年·兹比尔1436.62192
[16] S.Dolgov、B.N.Khoromskij、A.Litvinenko和H.G.Matthies,{随机系数的多项式混沌展开和张量列格式的随机偏微分方程解},SIAM/ASA J.不确定性。数量。,3(2015),第1109-1135页·Zbl 1329.65271号
[17] S.V.Dolgov和D.V.Savostyanov,{高维线性系统的交替最小能量方法},SIAM J.Sci。计算。,36(2014),第A2248-A2271页·Zbl 1307.65035号
[18] M.Eigel、M.Pfeffer和R.Schneider,{具有层次张量表示的自适应随机Galerkin FEM},Numer。数学。,136(2017),第765-803页·Zbl 1397.65262号
[19] M.Espig、L.Grasedyck和W.Hackbusch,{使用光纤交叉的黑盒低张量-库近似},Constr。约,30(2009),第557-597页·Zbl 1181.65064号
[20] T.Gerstner和M.Griebel,{尺寸自适应张量积求积},《计算》,71(2003),第65-87页·Zbl 1030.65015号
[21] S.A.Goreinov、I.V.Oseledets、D.V.Savostyanov、E.E.Tyrtyshnikov和N.L.Zamarashkin,{\it How to find A good submatrix},摘自《矩阵方法:理论、算法和应用》,V.Olshevsky和E.Tyrityshniko编辑,《世界科学》,新泽西州哈肯萨克,2010年,第247-256页·Zbl 1215.65078号
[22] I.Graham,F.Kuo,D.Nuyens,R.Scheichl,and I.Sloan,{随机系数椭圆偏微分方程的拟蒙特卡罗方法及其应用},J.Compute。物理。,230(2011),第3668-3694页·Zbl 1218.65009号
[23] L.Grasedyck,{张量的层次奇异值分解},SIAM J.矩阵分析。申请。,31(2010),第2029-2054页·Zbl 1210.65090号
[24] W.Hackbusch,《椭圆微分方程:理论和数值处理》,Springer-Verlag,柏林,1992年·Zbl 0755.35021号
[25] W.Hackbusch,{张量空间和数值张量微积分},Springer-Verlag,柏林,2012·Zbl 1244.65061号
[26] W.Hackbusch和S.Kuöhn,《张量表示的一种新格式》,J.Fourier Ana。申请。,15(2009),第706-722页·Zbl 1188.15022号
[27] S.Holtz、T.Rohwedder和R.Schneider,《张量列格式中张量优化的交替线性方案》,SIAM J.Sci。计算。,34(2012),第A683-A713页·Zbl 1252.15031号
[28] M.Kavehrad和M.Joseph,{数字通信系统性能评估中的最大熵和矩方法},IEEE Trans。Comm.,34(1986),第1183-1189页。
[29] B.N.Khoromskij,{多维偏微分方程的张量数值方法:理论分析和初步应用},载于《CEMRACS 2013–复杂系统建模与仿真:随机和确定性方法》,ESAIM Proc。调查48,EDP科学。,Les Ulis,2015年,第1-28页·Zbl 1382.65461号
[30] B.N.Khoromskij和I.V.Oseledets,参数相关和随机椭圆PDEs的{数量-TT配置近似},计算。方法应用。数学。,10(2010),第376-394页·Zbl 1283.65039号
[31] B.N.Khoromskij和C.Schwab,参数和随机椭圆偏微分方程的{张量结构Galerkin近似},SIAM J.Sci。计算。,33(2011),第1-25页·Zbl 1243.65009号
[32] A.Klimke,《稀疏网格插值工具箱:用户指南》,IANS技术报告2007/017,斯图加特大学,2007年。
[33] A.Klimke和B.Wohlmuth,{it Algorithm 847:spinterp:MATLAB}中的分段多线性分层稀疏网格插值,ACM Trans。数学。《软件》,31(2005),第561-579页·Zbl 1136.65308号
[34] T.G.Kolda和B.W.Bader,{张量分解和应用},SIAM Rev.,51(2009),第455-500页·Zbl 1173.65029号
[35] F.Kuo,{\it格规则生成向量}。
[36] F.Kuo,R.Scheichl,C.Schwab,I.Sloan和E.Ullmann,{对数正态扩散问题的多级拟蒙特卡罗方法},数学。公司。,86(2017),第2827-2860页·Zbl 1368.65005号
[37] F.Y.Kuo,C.Schwab,I.H.Sloan,{发现了一类随机系数椭圆偏微分方程}的多层拟蒙特卡罗有限元方法。公司。数学。,15(2015),第411-449页·Zbl 1318.65006号
[38] H.Matthies和A.Keese,{线性和非线性椭圆随机偏微分方程的Galerkin方法},计算。方法应用。机械。工程,194(2005),第1295-1331页·Zbl 1088.65002号
[39] H.Niederreiter,{\it准蒙特卡罗方法和伪随机数},Bull。阿默尔。数学。《社会学杂志》,84(1978),第957-1041页·Zbl 0404.65003号
[40] F.Nobile,L.Tamellini,F.Tesei,and R.Tempone,{\it具有对数正态扩散系数的椭圆偏微分方程的自适应稀疏网格算法},摘自《稀疏网格与应用》,斯图加特,2014年,施普林格,查姆,2016年,第191-220页·Zbl 1339.65016号
[41] F.Nobile和F.Tesei,{带控制变量的对数正态系数椭圆偏微分方程的多层蒙特卡罗方法},Stoch。部分差异。埃克。分析。计算。,3(2015),第398-444页·Zbl 1334.60133号
[42] I.V.Oseledets,{张量-应变分解},SIAM J.科学。计算。,33(2011),第2295-2317页·Zbl 1232.15018号
[43] I.V.Oseledets,{低阶张量格式函数的构造表示},Constr。约,37(2013),第1-18页·Zbl 1282.15021号
[44] I.V.Oseledets、S.Dolgov、V.Kazeev、D.Savostyanov、O.Lebedeva、P.Zhlobich、T.Mach和L.Song,{it TT-Toolbox},2011年。
[45] I.V.Oseledets和S.V.Dolgov,《线性系统的解和TT格式中的矩阵反演》,SIAM J.Sci。计算。,34(2012),第A2718-A2739页·Zbl 1259.65071号
[46] I.V.Oseledets和E.E.Tyrtyshnikov,{\it多维数组的TT交叉近似},线性代数应用。,432(2010),第70-88页·Zbl 1183.65040号
[47] U.Schollwo¨ck,{\it密度矩阵重整化群},修订版。物理。,77(2005),第259-315页·Zbl 1205.82073号
[48] S.R.White,量子重整化群的密度矩阵算法},Phys。B版,48(1993),第10345-10356页。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。