我正在使用 Arpack 计算问题的特征值$\lambda Mx = Ax$具有复杂移位的反向移位方法。$A$和$M$是真实的，$M$是对称的。然后，我使用 znaupd e zneupd。我使用superLU来解决系统。一切似乎都很好。我编写了一个函数，该函数采用 Arpack 计算的特征值和特征向量并将它们放入原始问题中：$\lambda Mx -Ax = \text{residual}$. 残差的范数应该几乎为零或非常小。无论如何，我得到一个错误$10^{-3}$这是相当高的！你知道问题出在哪里吗？我希望这很清楚。

谢谢，

诺拉