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 問題の行列版。