fortran - Fortran 95 で線形方程式を解く

Question

私はプログラミングの初心者で、現在、Fortran 95 を使用する必要があるプロジェクトに取り組んでいます。線形方程式を解くためのサブルーチンはありますか。たとえば、2*2 行列の行列xを見つけるA*x=bなどです。A

この問題を解決するのに役立つ便利なリンクを教えていただければ幸いです。お手伝いありがとうございます

Answer 1

2x2 連立方程式の場合、det(A) の式はかなり単純なので (3x3 の場合も同様)、Cramer の規則をコーディングする必要があります。 https://en.wikipedia.org/wiki/Cramer%27s_rule

Answer 2

使用できるライブラリは多数あります。事実上の標準は、より多くのアルゴリズムを選択できるLAPACKパッケージです。OpenBLAS、Intel MKL、Sun Performance Library など、同じサブルーチンの無料の商用実装が多数あります。

システムが非常に大きい場合は、反復ソルバーを探します。それらの多くがあります。GMRES、BiCGSTAB、または同様のメソッドとその実装を検索してください (たとえば、http://people.sc.fsu.edu/~%20jburkardt/f_src/templates/templates.html )。