市來健吾の日記

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

ryuon


  • 高速な implicit algorithm の実装、 phase 2 - top level への組み込み:

    • 現状を整理。 Brownian dynamics 関連の parameters は、 "BD-scheme" で integrator を、 "BD-nl-solver" で implicit scheme の solver を指定している。

    • なので、 BD-scheme に "SI-connector" を、 BD-nl-solver に "fastSI" を追加することにしよう。

    • コードにどう手を入れていくか、見ていく。

      • 今回書いた新しいコードは struct BD_imp を受けて動くようになっているので、 問題無し(というか、ここしばらくは、この一般化に時間を費やしていた、と)。

      • evolver として bd-imp.c の BD_imp_ode_evolve から fastSI 系のコードを呼べばよい。

      • そのルーチン(wrapper, 名前は BD_evolve_fastSI かな) を書けばよい、と。

    • 今日はここまで。

    • 8/11/2008: 実装、完了。もうちょっと。