shift-invert 方法是否存在可逆性问题?

计算科学 稀疏矩阵 特征值 scipy
2021-12-12 12:06:58

请注意,我对数值方法的背景几乎为零。

我了解SciPy 教程中描述的移位反转方法

上述链接的主要论点如下。假设我们要解决一个广义特征值问题Ax=λMx. 然后,移位反转方法将问题转换为

(AσM)1Mx=1λσx.

我的问题:如果AσM不可逆?特别是我对这个案子感兴趣M=Iσ=0. 然后,上式简化为A1x=λ1x. 但是,我知道A有零作为特征值,所以上述问题确实发生了!

提问的动机如下。使用eigshSciPy 中的函数,我得到具有最小量级的特征值是有序的104,但它应该是零!

我用了σ=0eigsh函数参数中,我的矩阵是稀疏的 Hermitian,大小10000×10000.

0个回答
没有发现任何回复~