市來健吾の日記

プログラマ、(元)物理屋(ナノテク、流体)

ryuon


  • 高速化の試み(前処理)。

    • 前処理付き逐次法の実装。

      • LU 分解を使った線形方程式の解法ルーチン inv_ILU を書く。

      • まずは CG の前処理版を書く。

      • libiter の orthomin を前処理化。

        • まずは今のコードを ALTAS 化する(ほとんど効果はなかったけど)。

        • 実装した。一応動いている様子 (ただし、今は LU 分解は完全なので一発で終わってて、 本当のチェックにはなっていない)。

      • 対角成分の反転による前処理も実装(inv_diag)。

        • これでもうまく走っているので、よしとする。

      • libiter の gmres.c にも手を入れる。

        • 実装済み。 とりあえず inv_diag, inv_ILU とも、まともに動いている模様。