The use of block Krylov subspace methods for computing the solution to a sequence of shifted linear systems using subspace recycling was first proposed in [Soodhalter, SISC 2016], where a recycled shifted block GMRES algorithm (rsbGMRES) was proposed. Such methods use the equivalence of the shifted system to a Sylvester equation and exploit the shift invariance of the block Krylov subspace generated from the Sylvester operator. This avoids the need for initial residuals to span the same subspace and allows for a viable restarted Krylov subspace method with recycling for solving sequences of shifted systems. In this paper we propose to develop these types of methods using unprojected Krylov subspaces. In doing so we show how one can overcome the difficulties associated with developing methods based on projected Krylov subspaces such as rsbGMRES, while also allowing for practical methods to fit within a well known residual projection framework. In addition, unprojected methods are known to be advantageous when the projector is expensive to apply, making them of significant interest for High-Performance Computing applications. We develop an unprojected rsbFOM and unprojected rsbGMRES. We also develop a procedure for extracting shift dependent harmonic Ritz vectors over an augmented block Krylov subspace for shifted systems yielding an approach for selecting a new recycling subspace after each cycle of the algorithm. Numerical experiments demonstrate the effectiveness of our methods.


翻译:在[Soodhalter, SISC 2016] 中,首次提议使用Cret Krylov 子空间来计算转换线性系统序列的解决方案。在[Soodhalter, SISC 2016] 中,首次提议使用Crept Krylov 子空间来计算使用子空间回收利用的转换线性系统序列的解决方案。这些方法使用Sylvester 方程式的等同系统,并利用Sylvester 操作者产生的Krylov 子空间块的变换。这避免了初始残渣跨越同一子空间的需要,并允许采用可行的重新启用Krylov 子空间方法,用于解决转移系统的序列。在本文件中,我们提议利用未预测的Krylov 子空间来开发这些类型的方法。我们如何克服基于预测的Krylov 子空间(如 rsbGMRES) 的开发方法,同时允许实际方法适应众所周知的残余预测框架。此外,在投影器应用费用昂贵时,已知的Krylovallov 子空间再利用这些方法,从而对高精度的高级计算系统产生重大的兴趣分析程序。我们每个变压变压的系统后,我们再变换的系统将一个不动的系统。

0
下载
关闭预览

相关内容

专知会员服务
50+阅读 · 2020年12月14日
Linux导论,Introduction to Linux,96页ppt
专知会员服务
77+阅读 · 2020年7月26日
强化学习最新教程,17页pdf
专知会员服务
174+阅读 · 2019年10月11日
VCIP 2022 Call for Demos
CCF多媒体专委会
1+阅读 · 2022年6月6日
AIART 2022 Call for Papers
CCF多媒体专委会
1+阅读 · 2022年2月13日
【ICIG2021】Check out the hot new trailer of ICIG2021 Symposium9
中国图象图形学学会CSIG
0+阅读 · 2021年12月17日
【ICIG2021】Check out the hot new trailer of ICIG2021 Symposium8
中国图象图形学学会CSIG
0+阅读 · 2021年11月16日
【ICIG2021】Check out the hot new trailer of ICIG2021 Symposium1
中国图象图形学学会CSIG
0+阅读 · 2021年11月3日
【ICIG2021】Latest News & Announcements of the Industry Talk1
中国图象图形学学会CSIG
0+阅读 · 2021年7月28日
Hierarchically Structured Meta-learning
CreateAMind
26+阅读 · 2019年5月22日
Transferring Knowledge across Learning Processes
CreateAMind
27+阅读 · 2019年5月18日
Unsupervised Learning via Meta-Learning
CreateAMind
42+阅读 · 2019年1月3日
A Technical Overview of AI & ML in 2018 & Trends for 2019
待字闺中
16+阅读 · 2018年12月24日
国家自然科学基金
1+阅读 · 2015年12月31日
国家自然科学基金
0+阅读 · 2013年12月31日
国家自然科学基金
1+阅读 · 2012年12月31日
国家自然科学基金
1+阅读 · 2012年12月31日
国家自然科学基金
0+阅读 · 2012年12月31日
国家自然科学基金
0+阅读 · 2012年12月31日
国家自然科学基金
0+阅读 · 2011年12月31日
国家自然科学基金
0+阅读 · 2011年12月31日
Arxiv
0+阅读 · 2022年10月25日
Arxiv
0+阅读 · 2022年10月20日
VIP会员
相关资讯
VCIP 2022 Call for Demos
CCF多媒体专委会
1+阅读 · 2022年6月6日
AIART 2022 Call for Papers
CCF多媒体专委会
1+阅读 · 2022年2月13日
【ICIG2021】Check out the hot new trailer of ICIG2021 Symposium9
中国图象图形学学会CSIG
0+阅读 · 2021年12月17日
【ICIG2021】Check out the hot new trailer of ICIG2021 Symposium8
中国图象图形学学会CSIG
0+阅读 · 2021年11月16日
【ICIG2021】Check out the hot new trailer of ICIG2021 Symposium1
中国图象图形学学会CSIG
0+阅读 · 2021年11月3日
【ICIG2021】Latest News & Announcements of the Industry Talk1
中国图象图形学学会CSIG
0+阅读 · 2021年7月28日
Hierarchically Structured Meta-learning
CreateAMind
26+阅读 · 2019年5月22日
Transferring Knowledge across Learning Processes
CreateAMind
27+阅读 · 2019年5月18日
Unsupervised Learning via Meta-Learning
CreateAMind
42+阅读 · 2019年1月3日
A Technical Overview of AI & ML in 2018 & Trends for 2019
待字闺中
16+阅读 · 2018年12月24日
相关基金
国家自然科学基金
1+阅读 · 2015年12月31日
国家自然科学基金
0+阅读 · 2013年12月31日
国家自然科学基金
1+阅读 · 2012年12月31日
国家自然科学基金
1+阅读 · 2012年12月31日
国家自然科学基金
0+阅读 · 2012年12月31日
国家自然科学基金
0+阅读 · 2012年12月31日
国家自然科学基金
0+阅读 · 2011年12月31日
国家自然科学基金
0+阅读 · 2011年12月31日
Top
微信扫码咨询专知VIP会员