在python的numpy和Scipy中找到大矩阵的特征向量的问题

Question

我正在使用Python，并尝试使用numpy.linalg.eig()（对于稠密矩阵）或scipy.sparse.linalg.eigsh（对于稀疏矩阵）来计算巨大（20.000x20.000和更大）实对称矩阵的所有特征向量。 这适用于最多20.000列和行的矩阵。 但是对于更大的矩阵，计算出的特征向量的条目全为零。

这里的python有什么问题，如何在不减少可计算特征向量的情况下解决它？

使用这两种方法，计算出的特征向量的所有条目均为0。