在scipy的ARPACK绑定中，无法计算矩阵的所有特征值 . 但是，我发现eigsh能够计算n - 1个特征值，而eigs只能计算n - 2个特征值 . 任何人都可以验证这实际上是ARPACK的基本限制而不是scipy中的错误吗？
这是示例代码：
import scipy.sparse, scipy.sparse.linalg
t = scipy.sparse.eye(3,3).tocsr()
l,v = scipy.sparse.linalg.arpack.eigs(t,k=2)
l,v = scipy.sparse.linalg.arpack.eigsh(t,k=2)
1 回答
这是ARPACK限制：
http://forge.scilab.org/index.php/p/arpack-ng/source/tree/master/SRC/dnaupd.f
http://forge.scilab.org/index.php/p/arpack-ng/source/tree/master/SRC/dsaupd.f
弄错了将是一个奇怪的错误......