我想获取特定实对称矩阵的所有特征值和特征向量。 numpy.linalg.eigh
显然可以做到这一点,但是,该矩阵具有特殊的稀疏结构,可以对向量进行线性缩放点积。因此,我想使用scipy.sparse.linalg.eigsh
(它允许使用LinearOperator
代替输入数组),并使用隐式重启的Lanczos方法。
我的问题是scipy.sparse.linalg.eigsh
不允许计算所有特征值和特征向量(即k=n
),并且输入矩阵的秩通常等于n
。有什么办法可以解决这个问题,或者其他功能是否允许类似的功能?