我有一个50,000 x 50,000的密集矩阵或更大的矩阵。如果我使用numpy或scipy-package,则所有本征向量的条目均为0。如果我使用scipy.sparse仅计算1000-8000个本征向量,则将得到正确的本征向量。但我都需要。
出现此问题的原因是什么? 我可以使用LAPACK或ARPACK来计算正确的特征向量吗?
请注意,我的矩阵是networkx图的表示,因此是稀疏矩阵。我使用numpy.linalg
将它们转换为密集矩阵,否则我将使用scipy.sparse.linalg
。