标签: eigenvalue arpack
我正在使用Arpack来计算问题的特征值\ lambda Mx = Ax,采用具有复位移的反向移位方法。然后,我使用znaupd e zneupd。 一切似乎都没问题。我写了一个函数,它采用Arpack计算的特征值和特征向量,并将它们放在原始问题中: \ lambda Mx -Ax =残差。残差的范数应该几乎为零或非常小。无论如何,我得到了一个非常高的10 ^ -3的错误!你知道问题出在哪里吗? 我希望很清楚。
谢谢,
Nora