我正在研究一个流动物理问题（2D），最终归结为解决特征值问题。即使对于 60 x 60 的网格大小，特征值问题中的矩阵最终也会有 12.6 亿个元素（这个矩阵必须迭代求解！）。有时矩阵是稀疏的，有时它不取决于要解决的问题的性质。

我使用用 Fortran 编写的ARPACK来解决这个问题。令人惊讶的是，MATLAB 解决相同特征值问题的速度至少比 Fortran 解决方案快 5 倍。（我想知道怎么做！）

我的问题是：在计算社区中是否有任何已知的技术可以加速 Eigen 求解器？

并行化 ARPACK 是我尝试过的一种解决方案，但仍然达不到 MATLAB 的水平