前回、簡単な例で逐次代入法を考察した。koideforest.hatenadiary.com図の使い方のテストもかねて、その内容を図にしてみた。収束因子を使わない場合、は以下のようになる。 ただの長方形をグルグルまわるだけで一向に収束しないのが分かる。収束因子を使う…

自分の今の環境では、MKL11.3を使ってfortran90でプログラムを弄っている。違う環境では、ここで言うことは当てはまらないかも知れない。逆行列を計算したい。私の願いはそれだけであった。一般の行列に対して対応可能なsubroutineは、GETRF (LU分解)とGETRI…