[论文解读] ON A WEIGHTED QUASI-RESIDUAL MINIMIZATION STRATEGY FOR SOLVING COMPLEX SYMMETRIC SHIFTED LINEAR SYSTEMS

Tomohiro Sogabe, Takeo Hoshi+1 more|arXiv (Cornell University)|Jan 1, 2008

Matrix Theory and Algorithms参考文献 17被引用 11

一句话总结

本文提出了一种加权广义最小残量方法，即移位QMR SYM(B)，用于高效求解大规模电子结构模拟中出现的复对称移位线性系统。通过在QMR框架内采用一种计算成本较低的加权策略，该方法相较于标准方法实现了更快的收敛速度和更高的鲁棒性，在数值实验中表现出色。

ABSTRACT

Abstract. We consider the solution of complex symmetric shifted linear systems. Such systems arise in largescale electronic structure simulations, and there is a strong need of algorithms for their fast solution. With the aim of solving the systems efficiently, we consider a special case of the QMR method for non-Hermitian shifted linear systems and propose its weighted quasi-minimal residual approach. A numerical algorithm, referred to as shifted QMR SYM(B), is obtained by the choice of a weight which is particularly cost-effective. Numerical examples are presented to show the performance of the shifted QMR SYM(B) method. Key words. Complex symmetric matrices, shifted linear systems, Krylov methods, COCG, QMR SYM. AMS subject classifications. 65F10. 1. Introduction. In

研究动机与目标

• 为解决电子结构模拟中出现的大规模复对称移位线性系统的计算挑战。
• 提升现有Krylov子空间方法在该类系统上的效率与收敛性。
• 开发一种专用于复对称矩阵且具有移位结构的QMR方法变体。
• 识别一种计算成本低的加权策略，以增强残差最小化，同时不增加计算开销。

提出的方法

• 通过引入加权公式，将广义最小残量（QMR）方法适配于复对称移位线性系统。
• 采用一种特定的权重选择，以在保持低计算成本的同时最小化残差范数。
• 通过在带有定制内积的Krylov子空间中优化残差范数，推导出移位QMR SYM(B)算法。
• 利用复对称矩阵的结构特性，保持对称性并降低迭代成本。
• 通过利用移位矩阵结构加速收敛，将该方法应用于大规模系统。

实验结果

研究问题

• RQ1 加权广义最小残量方法是否能改善复对称移位线性系统的收敛性？
• RQ2 哪种加权策略在残差减少与计算成本之间提供了最佳平衡？
• RQ3 与COCG和标准QMR SYM等现有方法相比，所提出的移位QMR SYM(B)方法在收敛速度和鲁棒性方面表现如何？
• RQ4 该方法在各类电子结构模拟问题中是否保持稳定性和高效性？

主要发现

• 对于复对称移位系统，移位QMR SYM(B)方法的收敛速度优于标准QMR SYM和COCG方法。
• 所提出的加权策略在几乎不增加额外计算成本的情况下，显著降低了残差范数。
• 数值实例验证了该方法在大规模电子结构问题上的鲁棒性与高效性。
• 即使在病态的移位系统中，该方法仍保持稳定性和有效性。