在 Fortran 90 中实现直接求解器?

计算科学 线性求解器 稀疏矩阵
2021-12-10 01:30:05

我的问题可能是基本的,但当我感到困惑时,它非常重要。在这里,我应该解决以下等式:

Ax=B

据我了解,我可以选择使用直接求解器或迭代求解器。我曾尝试使用迭代求解器 GMRES(在 Fortran 90 上的实现非常简单),但现在由于某些原因我想使用直接求解器,并且从资源和评论中,我将其中三个入围:

  1. UMFPACK
  2. 超级路
  3. 腮腺炎

我认为前两个是使用 C 开发的,可以从 Fortran 90 中调用;第三个是在 Fortran 90 中。

现在,我的问题是如何从我的 Fortran 90 代码中实现这个或调用函数,即。是否有可靠的程序来安装或调用库?

似乎有大量文献超出了我的理解范围。

为了简化,我有一个稀疏的对称矩阵A, 以天际线格式和 rhs 矢量存储B, 我如何计算x使用 Fortran 90 中的直接稀疏求解器之一?

0个回答
没有发现任何回复~