ScaLAPACK中为何无GEEV等效函数？非厄米矩阵特征值求解求助

一、为何无类似PZGEEV的通用复非厄米特征求解器？ #

ScaLAPACK的核心设计目标是分布式内存架构下的高效并行计算，复非厄米矩阵通用特征值求解（对应LAPACK的ZGEEV）的并行化实现存在极高门槛：

• 这类算法的QR迭代过程涉及大量全局数据依赖与同步通信，在分布式环境下很难实现可扩展的高效版本；
• ScaLAPACK优先覆盖了需求更广、并行化难度更低的模块（如厄米矩阵特征值、SVD等），对非厄米通用特征求解这类小众且实现成本极高的模块，未纳入核心库范畴。

二、现有工具结合数学技巧的替代方案 #

针对你这种由厄米矩阵扰动得到的非厄米系统，推荐两种可行方案：

方案1：构造扩充厄米矩阵求解 #

对于非厄米矩阵 ( A = H + V )（( H ) 为原厄米矩阵，( V ) 为扰动项），可构造2n阶厄米矩阵：
[
B = \begin{pmatrix}
0 & A \
A^H & 0
\end{pmatrix}
]
该矩阵的特征值为 ( \pm \lambda )，其中 ( \lambda ) 正是原非厄米矩阵 ( A ) 的特征值。你可以直接调用ScaLAPACK的PZHEEV求解 ( B ) 的特征值，再提取对应原矩阵的特征值（注意去重与符号匹配）。

• 优势：复用成熟的厄米特征求解器，并行效率有保障；
• 劣势：矩阵规模翻倍，内存与计算量开销相应增加。

方案2：使用并行Schur分解工具 #

ScaLAPACK提供了PZGESCHUR，专门用于复非厄米矩阵的并行Schur分解。Schur分解会将 ( A ) 转化为上三角矩阵 ( T = Q^H A Q )，( T ) 的对角元就是 ( A ) 的特征值，虚部信息完整保留。操作步骤如下：

• 调用PZGESCHUR完成Schur分解，得到上三角矩阵 ( T )；
• 提取 ( T ) 的对角元，即为原矩阵的特征值。

注意：PZGESVD计算的是奇异值（非负实数），与特征值属于完全不同的数值分解，这也是单独使用它无法满足需求的核心原因。

三、额外提示 #

对于扰动幅度较小的非厄米系统，上述两种方案的数值稳定性足以满足物理研究需求，特征值的虚部精度可以得到有效保障。

内容的提问来源于stack exchange，提问作者Stephen Gant