市來健吾の日記

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

ryuon


  • libstokes の高速化:

    • 行列版で LAPACK を使っている部分に無駄があった、という話。 (cf. 11/25/2007

    • まずは、該当部分の修正が本当に高速化されるか確認:

      • test/check-lapack-solve-lin.c を書く。

      • 結果:かなり高速化されることが分かる。 100x100 の行列で約 2.6 倍、 1000x1000 なら 6.8 倍だ!

    • この修正が適用可能な部分を探す。

      • 結局、使えるのは2ヶ所のみだった。 つまり、 F version と FT version の res 問題および mob-lub 問題の行列版。