Scipy arpack eigs与eigsh特征值的数量

Question

在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)

Answer 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

这会是一个奇怪的错误......