Julia 特征值问题

我试图在julia中对一个矩阵进行对角化，但是这个矩阵本身太大，无法计算。但是我可以有效地计算矩阵向量积。我正在寻找一个解决方案，这样我可以将矩阵向量积应用到EIG中，并使用lanzcos方法进行求解。是否有一个软件包以矩阵向量积的形式给我输出

---

有人可以帮忙吗？

您可以使用Julia 0.6.2中包含的基于ARPACK的

eigs

函数，或者您可以使用中相关方法的纯Julia实现。后者不是基于Lanczos的，但我想这并不重要，只要它能给出特征值。

谢谢你的回答。我一直在使用迭代解算器包。我只是想知道如何创建一个操作员产品表示以作为输入。我一直在尝试使用LinearRoperator包。最简单的解决方案是将其作为

AbstractMatrix

的子类型。