我写了一段代码，使用一些离散化方法求解大型偏微分方程系统，基本上包括在每个时间步长多次求解大型稀疏系统Ax=b。
我目前使用PARDISO求解器（来自英特尔MKL库），它是a的直接LU因子分解来求解系统。我想将这种方法与迭代求解器的使用进行比较（使用预处理器后，迭代求解器可能会表现得更好，因为如果我的雅可比矩阵变化不大，我可以在许多时间步长上使用相同的预处理器）。
那么，我的问题是，你建议fortran中的稀疏迭代求解器使用什么库？我发现了一个（SLATEC），它写于1993年，所以我想知道是否有更具表演性的东西是最近写的？
感谢：）
我还想添加：
LIS
AGMG
哦，好吧。。。线性代数解算器的完整列表
感谢您的评论，PETSc似乎正是我想要的，现在只需要学习如何将C调用链接到fortran中：）