julia - 特征值问题

Question

我正在尝试对 julia 中的矩阵进行对角化，但矩阵本身太大而无法计算。但我可以有效地计算矩阵向量积。我正在寻找一种解决方案，以便我可以将矩阵向量积用于 eigs 并使用 lanzcos 方法解决。是否有一个包可以以矩阵向量积的形式给我一个输出？

有人可以帮忙吗？

Answer 1

您可以使用 Julia 0.6.2 中包含的基于 ARPACK 的函数，也可以使用https://github.com/JuliaMath/IterativeSolvers.jleigs中相关方法的纯 Julia 实现。后者不是基于 Lanczos 的，但我想它并不重要，只要它给你特征值。