2013-01-12 31 views
2

在scipy的ARPACK綁定中,無法計算矩陣的所有特徵值。但是,我發現eigsh能夠計算n-1個特徵值,而eigs只能計算n-2個特徵值。任何人都可以驗證這實際上是ARPACK的根本限制,而不是scipy中的錯誤?Scipy arpack eigs與eigsh的特徵值數量

下面是示例代碼:

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) 

回答

相關問題