使用 GMRES 解决密集线性系统的库

计算科学 正则 拉帕克 预处理 密集矩阵 格瑞斯
2021-12-24 15:13:34

我有一个 fortran 90 代码,我想用 GMRES 解决一个密集的线性系统。我更喜欢重新启动的带有预处理的 GMRES。你知道我可以使用一些图书馆吗?现在我正在使用LAPACKdgetrf和LU 解决系统问题。dgetrs

我还想补充一点,矩阵是正方形的,对角线占优且可逆。

谢谢!

1个回答

PETSC 可能是一个不错的选择。不是超级用户友好,但它有很多选择。