特征值探针

Question

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

有人可以帮忙吗？

Answer 1

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