torous の計算 on bue9。 一応 N=1000 も作ったが、まずは N=200 から。 今回は dt を 0.1 に刻んでおく。

Brownian dynamics 検証 run - 粒子編。 bue9 で週末走らせておいた結果： JGdP だと dt=5e-5 は問題なし。 dt=0.0005 だと途中で止まってた。 bue9(SUSE-Linux) の RYUON の整備。 libstokes の pgcc 版ができた。 が、結果は期待はずれ。 ということで、 b…