[论文解读] Subspace Iteration with Approximate Spectral Projection

Ping Tang, Eric Polizzi|arXiv (Cornell University)|Feb 2, 2013

Matrix Theory and Algorithms被引用 10

一句话总结

本文提供了FEAST特征值算法的首个严格收敛性分析，表明其作为使用有理函数逼近谱投影算子的加速子空间迭代算法运行。该分析证明了FEAST在舍入误差存在的情况下仍能收敛，并为其鲁棒性及向非厄米特问题的扩展建立了理论基础。

ABSTRACT

The calculation of a segment of eigenvalues and their corresponding eigenvectors of a Hermitian matrix or matrix pencil has many applications. A new density-matrix-based algorithm has been proposed recently and a software package FEAST has been developed. The density-matrix approach allows FEAST's implementation to exploit a key strength of modern computer architectures, namely, multiple levels of parallelism. Consequently, the software package has been well received, especially in the electronic structure community. Nevertheless, theoretical analysis of FEAST has lagged. For instance, the FEAST algorithm has not been proven to converge. This paper offers a detailed numerical analysis of FEAST. In particular, we show that the FEAST algorithm can be understood as an accelerated subspace iteration algorithm in conjunction with the Rayleigh-Ritz procedure. The novelty of FEAST lies in its accelerator which is a rational matrix function that approximates the spectral projector onto the eigenspace in question. Analysis of the numerical nature of this approximate spectral projector and the resulting subspaces generated in the FEAST algorithm establishes the algorithm's convergence. This paper shows that FEAST is resilient against rounding errors and establishes properties that can be leveraged to enhance the algorithm's robustness. Finally, we propose an extension of FEAST to handle non-Hermitian problems and suggest some future research directions.

研究动机与目标

• 解决尽管FEAST特征值算法在电子结构计算中被广泛应用，但其缺乏理论收敛性分析的问题。
• 正式确立FEAST作为带有瑞利-里茨投影的加速子空间迭代算法的运行机制。
• 分析FEAST中使用的有理逼近的数值特性，以理解其稳定性和收敛行为。
• 证明FEAST对舍入误差具有鲁棒性，确保其在高性能计算环境中的实际可靠性。
• 将FEAST的理论框架扩展至非厄米特特征值问题，从而拓展其应用领域。

提出的方法

• 将FEAST重新表述为一种由逼近目标特征子空间上谱投影算子的有理矩阵函数加速的子空间迭代。
• 利用瑞利-里茨方法从每次迭代生成的子空间中提取里茨对。
• 通过研究近似谱投影算子的谱性质，分析子空间迭代的收敛性。
• 基于有理逼近对真实谱投影算子的精度，建立残差和收敛速率的界。
• 利用有理矩阵函数的结构，证明其在有限精度算术中对浮点误差的鲁棒性。
• 通过将有理投影适配至复围线积分，提出FEAST框架在非厄米特矩阵上的推广。

实验结果

研究问题

• RQ1 FEAST算法在有限精度算术下是否收敛，其收敛性的理论基础是什么？
• RQ2 谱投影算子的有理逼近如何影响FEAST的收敛速率和稳定性？
• RQ3 近似谱投影算子的哪些特性可确保其对舍入误差的数值鲁棒性？
• RQ4 在保持收敛性和效率的前提下，FEAST框架能否扩展至非厄米特特征值问题？
• RQ5 将FEAST解释为带有有理投影的加速子空间迭代，可获得哪些理论洞见？

主要发现

• FEAST被正式证明为一种由谱投影算子的有理逼近加速的子空间迭代算法。
• FEAST的收敛速率由有理逼近对谱投影算子的精度决定，逼近越精确，收敛越快。
• 由于有理矩阵函数的结构和瑞利-里茨过程，FEAST对舍入误差具有鲁棒性，确保了数值稳定性。
• 在理想条件下，算法保持二次收敛行为，且在有限精度算术中收敛性依然得以保持。
• 理论框架支持通过复围线积分和有理投影将FEAST推广至非厄米特问题。
• 该分析为通过改进有理逼近策略和预条件技术来提升FEAST的鲁棒性提供了理论基础。