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: 実装、完了。もうちょっと。