我有一个图拉普拉斯,我需要找出最大的'k'特征值和特征向量。我正在使用这样的东西: -
#L= laplacian matrix.
eigVal,eigVectors = eigsh(L, k, which='LA')
这给了我大致正确的结果,但出了点问题,我的eig值略大于1(比如说1.05)。在我的例子中,当使用MATLAB和其他平台时,特征值的上限为1.我得到了理想的结果。
我在这做错了什么?有没有办法可以并行化特征向量和值的计算? (我正在考虑pyCuda。)
答案 0 :(得分:1)
您确定您的Laplacian的Python实现是否正确?你仔细检查了吗输入矩阵是对称的吗?
如果没有您手头的特定矩阵,很难说究竟出了什么问题。你能保存矩阵并把它放在互联网上吗?
编辑:删除提及eigs*
之前的行为 - 例程之前没有eigsh
名称,所以情况并非如此。