我编写了一个代码，它使用一些离散化方法解决了大型PDE系统，这些方法主要涉及在每个时间步长多次解决大型稀疏系统Ax = b。
我目前使用PARDISO求解器（来自intel MKL库），它是A的直接LU分解来解决系统问题。我想比较这种方法与迭代求解器的使用（使用预处理器，可能表现得更好，因为如果我的雅可比矩阵变化不太大，我可以在很多时间步长上使用相同的预处理器。）
那么我的问题是，你对fortran中的稀疏迭代求解器建议使用哪个库？我发现了一个（SLATEC）是在1993年写的，所以我想知道是否有一些更新的性能是最近写的？
谢谢：）
答案 0 :(得分：1)
哦，好吧...... complete list of Linear Algebra Solvers
答案 1 :(得分：0)
感谢您的评论，PETSc似乎正是我所寻找的，只需要学习如何将C调用链接到fortran：）