[论文解读] Real dqds for the nonsymmetric tridiagonal eigenvalue problem

Carla Ferreira, Beresford Ν. Parlett|arXiv (Cornell University)|Jan 24, 2012

Matrix Theory and Algorithms参考文献 18被引用 3

一句话总结

本文提出了一种用于计算实非对称三对角矩阵特征值的实数 dqds 算法，仅使用实数算术运算，通过三步变换实现高效的复数位移。其性能与复数 dqds 代码相当，且比 Ehrlich-Aberth 方法快达 100 倍，同时具备稳健的元素增长控制和动态变换拒绝机制，确保了稳定性和速度。

ABSTRACT

We present a new transform, triple dqds, to help to compute the eigenvalues of a real tridiagonal matrix C using real arithmetic. The algorithm uses the real dqds transform to shift by a real number and triple dqds to shift by a complex conjugate pair. We present what seems to be a new criteria for splitting the current pair L,U. The algorithm rejects any transform which suffers from excessive element growth and then tries a new transform. Our numerical tests show that the algorithm is about 100 times faster than the Ehrlich-Aberth method of D. A. Bini, L. Gemignani and F. Tisseur. Our code is comparable in performance to a complex dqds code and is sometimes 3 times faster.

研究动机与目标

• 开发一种仅使用实数算术运算的快速且精确的实非对称三对角矩阵特征值求解器。
• 在避免复数算术运算的同时，通过共轭对实现复数位移，保持数值稳定性和效率。
• 设计一种鲁棒的变换拒绝机制，用于检测并丢弃导致过度元素增长的步骤。
• 在保持实数算术优势的同时，实现与复数 dqds 算法相当的性能。
• 为非对称三对角特征值问题提供一种实用且高性能的替代方案，优于现有方法如 Ehrlich-Aberth 和复数 QR/dqds。

提出的方法

• 该算法采用一种称为 tridqds 的三步变换，结合实数 dqds 与通过共轭对实现的复数位移，以实数算术模拟双位移。
• 通过一系列相似变换，使用矩阵 $ Z_i $, $ ilde{L}_i^{-1} $, $ ilde{L}_i $, $ Y_i^{-1} $, 和 $ Y_i $，对状态变量 $ x_r, y_r, z_r, x_l, y_l $ 进行递归更新。
• 通过比率 $ \hat{u}_i / \text{max}(u_i, l_i) $ 监控元素增长，若增长超过阈值则拒绝该变换。
• 基于当前的 $ L,U $ 对采用动态位移策略，选择能最大化收敛速度的位移，同时避免不稳定性。
• 该方法在 MATLAB 中实现，并包含降维与分裂准则，用于隔离已收敛的特征值。
• 核心计算通过递归更新 $ x, y, z $ 状态向量完成，每一步均进行归一化与缩放以保持数值稳定性。

实验结果

研究问题

• RQ1 能否设计一种仅使用实数算术的算法，使其在非对称三对角特征值问题上的性能与复数 dqds 相当？
• RQ2 如何在不引入数值不稳定性的情况下，有效在实数算术中应用复数位移？
• RQ3 何种标准可可靠检测并拒绝导致过度元素增长的变换，以确保算法鲁棒性？
• RQ4 与 Ehrlich-Aberth 和复数 dqds 等现有方法相比，三步 tridqds 变换在速度和精度方面表现如何？
• RQ5 该算法能否在避免复数算术计算开销的同时，保持高相对精度？

主要发现

• 所提出的实数 tridqds 算法比 Bini、Gemignani 和 Tisseur 的 Ehrlich-Aberth 方法快约 100 倍。
• 该算法的性能与复数 dqds 代码相当，某些情况下甚至快达三倍。
• 对过度元素增长变换的动态拒绝机制显著提升了数值稳定性与收敛可靠性。
• 该算法仅使用实数算术即可成功计算实非对称三对角矩阵的所有特征值，即使特征值为共轭复数对亦成立。
• 分裂与降维策略能有效隔离已收敛的特征值，从而高效计算完整特征谱。
• 该方法保持了高相对精度，与已知的 dqds 在奇异值计算中的优势一致，并将其拓展至非对称特征值问题。