Fortran中的迭代求解器库

我编写了一个代码，它使用一些离散化方法来求解大型偏微分方程组，基本上涉及到在每个时间步多次求解大型稀疏系统Ax=b

我目前使用PARDISO解算器（来自英特尔MKL库），它是a的直接LU分解，用于解算系统。我想将此方法与使用迭代解算器进行比较（如果雅可比矩阵变化不大，使用预条件器可能会更好，因为我可以在多个时间步上使用相同的预条件器）

那么，我的问题是，您建议在fortran中为稀疏迭代求解器使用什么库？我找到了一本1993年写的书（SLATEC），所以我想知道是不是最近写的更有表现力

谢谢：）

谢谢你的评论，似乎正是我想要的，现在只需要学习如何将C调用链接到fortran：）

我还要补充：

---

哦，好吧。。。因为它太旧了，并不意味着它就不好看！：]其中一些可能会引起兴趣：如果你不介意其他语言，你一定要去看看。IMO，它正是你想要的。PETSC支持C、C++和FORTRAN.PESC提供了一个合适的FORTRAN API。你应该看看这些例子。