algorithm - 高效的 EigenSolver 实现

Question

我正在寻找一种高效的特征求解器(语言并不重要，尽管我会使用 C# 进行编程)，它利用了现代 CPU 中的多核功能。能够直接使用 pardiso solver是一大优点。我的矩阵大多是sparse matrix ，因此理想的求解器应该能够利用这一事实并大大提高内存使用率和性能。

到目前为止我只找到了LAPACK和 ARPACK .拉帕克，as implemented in Intel MKL , 是一个很好的候选者，因为它提供了多核优化。但似乎 drivers LAPACK 内部不直接使用 pardiso 求解器，此外，它们似乎没有利用稀疏矩阵(但我不确定这一点)。

另一方面，ARPACK 似乎很难在 Windows 环境中设置，而并行版本 PARPACK , doesn't work so well .奖励点是它可以与 pardiso solver 一起使用.

最好的是具有多核加速的 Intel MKL + ARPACK。不确定是否有任何现有的实现已经完成了我想做的事情？

Answer 1

我正在处理一个与您所说的需求非常相似的问题。我正在考虑盛宴: http://www.ecs.umass.edu/~polizzi/feast/index.htm我现在正在努力让它发挥作用，但它看起来很完美。如果您使用它，我很想听听您的体验。干杯内德