×

边界元矩阵的线性时间CUR近似。 (英语) Zbl 1431.65022号

小结:本文针对从边界元矩阵的层次形式获得的容许矩阵,提出了线性时间CUR近似算法。我们提出了一种称为几何抽样的新方法,利用问题所在域的信息获取最重要行和列的索引。我们的策略是根据边界元方法(BEM)定制的,因为它直接明确地使用包含问题几何信息的聚类树。我们的CUR算法的精度与使用带列旋转的截断QR分解(QRCP)和带全旋转的自适应交叉近似(ACA)创建的低阶近似相当,这是二次成本方法。与著名的部分旋转线性时间算法ACA相比,我们的算法总体上改善了收敛误差,并克服了ACA失败的一些情况。我们提供了几何采样创建的CUR近似的一般相对误差界。最后,我们评估了我们的算法在不同几何形状上定义的传统边界元问题上的性能。

MSC公司:

65平方英尺 线性系统和矩阵反演的直接数值方法
65层25 数值线性代数中的正交化
65层55 低阶矩阵逼近的数值方法;矩阵压缩
65号38 偏微分方程边值问题的边界元方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部 哈尔

参考文献:

[1] 格林加德。;Rokhlin,V.,《粒子模拟的快速算法》,J.Compute。物理。,73, 2, 325-348 (1987) ·Zbl 0629.65005号
[2] Koumoutsakos,P.,《三维N体问题的快速多极方法》,377-390(1995),美国宇航局埃姆斯/斯坦福大学湍流研究中心。
[3] 方,W。;Darve,E.,《黑盒快速多极子方法》,J.Compute。物理。,228, 8712-8725 (2009) ·Zbl 1177.65009号
[4] Martinsson,P.G。;Rokhlin,V.,《一维加速核相关快速多极子方法》,SIAM J.Sci。计算。,29, 3, 1160-1178 (2007) ·Zbl 1154.65318号
[5] Bebendorf,M.,《层次矩阵》(Hierarchical matrices,2008),施普林格:施普林格-莱比锡,德国)·Zbl 1151.65090号
[6] Börm,S.,《非局部算子的有效数值方法:H-2矩阵压缩、算法和分析》(2010),欧洲数学学会·Zbl 1208.65037号
[7] Hackbusch,W.,《层次矩阵:算法和分析》,(矩阵计算(2015),计算数学中的Springer系列:巴尔的摩计算数学的Springr系列)·Zbl 1336.65041号
[8] 3月,西。;Biros,G.,《高维快速核求和方法的远场压缩》,应用。计算。哈蒙。分析。,43, 1, 39-75 (2017) ·Zbl 1365.65267号
[9] 马奥尼,M。;Drineas,P.,用于改进数据分析的CUR矩阵分解,Proc。国家。阿卡德。科学。,106, 3, 697-702 (2009) ·Zbl 1202.68480号
[10] 邱,J。;Demanet,L.,骨架分解的次线性随机算法,SIAM J.矩阵分析。申请。,34, 3, 1361-1383 (2013) ·兹比尔1277.68296
[11] Ozdemir,N.A。;Lee,J.-F.,体积积分方程中矩阵压缩的低阶IE-QR算法,IEEE Trans。马格纳。,40, 2, 1017-1020 (2004)
[12] 卡普尔,S。;Long,D.,N体问题:IES3:高效静电和电磁模拟,IEEE Computat。科学。工程师,5,4,60-67(1998)
[13] C.Boutsidis,M.Mahoney,P.Drineas,《列子集选择问题的改进近似算法》,载于《第二十届ACM-SIAM离散算法研讨会论文集》,2009年,第968-977页·兹比尔1420.68235
[14] 沃罗宁,S。;Martinsson,P.G.,《CUR和插值矩阵分解的高效算法》,《高级计算》。数学。,43, 3, 495-516 (2017) ·Zbl 1369.65049号
[15] 顾,M。;Eisenstat,S.,计算强秩揭示QR因式分解的高效算法,SIAM J.矩阵分析。申请。,17, 4, 848-869 (1996) ·Zbl 0858.65044号
[16] 库马尔,N。;G.,S.,矩阵低阶逼近的文献综述,线性多线性代数,65,11,2212-2244(2017)·Zbl 1387.65039号
[17] 霍恩,R。;Johnson,C.,矩阵分析主题,(矩阵计算(1991),剑桥大学出版社:美国纽约剑桥大学出版社)·Zbl 0729.15001号
[18] Mirsky,L.,《对称规范函数和酉不变范数》,夸特。数学杂志。牛津大学。,11, 2, 50-59 (1960) ·Zbl 0105.01101号
[19] 德梅尔,J。;格里戈里,L。;顾,M。;Xiang,H.,利用柱旋转避免秩揭示QR因式分解的通信,SIAM J.矩阵分析。申请。(2015) ·Zbl 1327.65078号
[20] 潘,C.-T。;Tang,P.,QR因式分解揭示的奇异值的界,BIT-Numer。数学。,39940-756(1999年)·Zbl 0944.65042号
[21] Golub,G。;Van Loan,C.,矩阵计算,(矩阵计算(1996),约翰霍普金斯大学出版社:巴尔的摩约翰霍普金斯大学出版社)·Zbl 0865.65009号
[22] 阿亚拉,A。;克莱斯,X。;Grigori,L.,ALORA:仿射低阶近似,科学杂志。计算。,79, 2, 1135-1160 (2019) ·Zbl 1417.65120号
[23] S.Goreinov,I.Oseledets,D.Savostyanov,E.Tyrtyshnikov,N.Zamarashkin,《如何找到一个好的子矩阵》,研究报告08-10,ICM HKBU,香港九龙塘,2008·Zbl 1215.65078号
[24] 戈雷诺夫,S。;Tyrtyshnikov,E.,低秩矩阵逼近中的最大体积概念,Contemp。数学。,280, 47-51 (2001) ·Zbl 1003.15025号
[25] 戈雷诺夫,S。;Tyrtyshnikov,E.,Chebyshev范数中矩阵骨架近似的拟最优性,Dokl。数学。,83, 3, 374-375 (2011) ·Zbl 1252.65078号
[26] 乔伊夫里尔,A。;Magdon-Ismail,M.,选择最大体积子矩阵的指数不可逼近性,算法,65,1,159-176(2013)·Zbl 1259.68086号
[27] Bebendorf,M.,边界元矩阵近似,数值。数学。,86, 4, 565-589 (2000) ·Zbl 0966.65094号
[28] 格里戈里,L。;Cayrols,S。;Demmel,J.,基于LU因式分解的稀疏矩阵的低秩近似,列和行锦标赛枢轴,SIAM J.Sci。计算。,40、2、C181-C209(2018)·Zbl 1453.65090号
[29] Higham,N.,《数值算法的准确性和稳定性》(2002),工业和应用数学学会·Zbl 1011.65010号
[30] Sommerville,D.,《n维几何导论》(1958),多佛:纽约多佛·Zbl 0086.35804号
[31] 奥辛斯基,A。;Zamarashkin,N.,具有更好精度估计的伪骨架近似,线性代数应用。,537, 221-249 (2018) ·Zbl 1376.65073号
[32] Deshpande,A。;拉德马赫,L。;Vempala,S。;Wang,G.,通过体积采样的矩阵近似和投影聚类,理论计算。,2, 225-247 (2006) ·Zbl 1213.68702号
[33] Steinbach,O.,《椭圆边值问题的数值逼近方法:有限元和边界元》,(椭圆边值方程的数值逼近法:有限元与边界元(2008),Springer:Springer New York,USA)·Zbl 1153.65302号
[34] S.Rjasanow,稠密矩阵的自适应交叉逼近,收录于:国际边界元方法协会会议,IABEM,2002年。
[35] 德米尔,J。;格里戈里,L。;霍姆曼,M。;Langou,J.,通信最优并行和序列QR和LU分解(2008)
[36] 德梅尔,J。;格里戈里,L。;霍姆曼,M。;Langou,J.,通信优化并行和序列QR和LU分解,SIAM J.科学。计算。,34、1、A206-A239(2012)·Zbl 1241.65028号
[37] 格雷,A。;A.,M.,《统计学习中的N体问题》,《高级神经信息处理》。系统。,521-527 (2001)
[38] S.Dasgupta,Y.Freund,《随机投影树和低维流形》,载于《第四十届ACM计算理论研讨会论文集》,STOC,2008年,第537-546页·Zbl 1231.68114号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不声称其完整性或完全匹配。