0

我有一个图拉普拉斯算子,我需要找出最大的“k”特征值和特征向量。我正在使用这样的东西:-

#L= laplacian matrix.
eigVal,eigVectors = eigsh(L, k,  which='LA')

这给了我大致正确的结果,但是出了点问题,我得到的 eig 值略大于 1(比如 1.05)。在我的情况下,特征值的上限为 1。当使用 MATLAB 和其他平台时,我得到了想要的结果。

我在这里做错了什么?有什么方法可以并行计算特征向量和值?(我正在考虑使用 pyCuda。)

4

1 回答 1

1

你确定你的拉普拉斯算子的 Python 实现是正确的吗?您是否仔细检查过输入矩阵是否对称?

如果手头没有您的特定矩阵,很难说到底出了什么问题。您可以保存矩阵并将其放在互联网上的某个地方吗?

编辑:删除了对先前行为的提及eigs*——该例程在此之前没有eigsh名称,所以这里不是这种情况。

于 2012-06-17T15:49:01.427 回答