如何将 netlib 中的 lanczos 代码用于大型稀疏对称矩阵?

计算科学 本征系统 正则 拉帕克 带状矩阵
2021-12-22 09:12:48

我想使用 lanczos 方法来计算大型稀疏对称矩阵(~50k x ~50k)的几个最低特征值和特征向量。http://www.netlib.org/lanczos/index.html我发现代码 leval 和 levec 适合我的情况。但是,该指令太复杂,找不到如何使用它,我不知道如何设置输入和输出文件。谁用过,请给点建议。谢谢!

1个回答

这并不能直接回答您的问题,但也许您想使用 ARPACK 之类的东西,例如scipy.sparse.linalg.eigsh在 python 中。