ryuon
高速化の試み(前処理)。
前処理付き逐次法の実装。
LU 分解を使った線形方程式の解法ルーチン inv_ILU を書く。
まずは CG の前処理版を書く。
libiter の orthomin を前処理化。
まずは今のコードを ALTAS 化する(ほとんど効果はなかったけど)。
実装した。一応動いている様子 (ただし、今は LU 分解は完全なので一発で終わってて、 本当のチェックにはなっていない)。
対角成分の反転による前処理も実装(inv_diag)。
これでもうまく走っているので、よしとする。
libiter の gmres.c にも手を入れる。
実装済み。 とりあえず inv_diag, inv_ILU とも、まともに動いている模様。