ryuon

• 
  stokes3 の polydisperse 対応
  • 
    まず非周期で mono と poly のチェックを済ませておく (Durlofsky et al. 1987 で)。 two-body-res (50項) と twobody (100+alpha 項) の精度の差が出て来るが、問題無し。 (表の数値の一つの最後の一桁が四捨五入で一つ繰り上がったが、 多分、こっちの値が正しいものだろう。)
  • 
    stokes3 のテスト(SC N=8) に戻る。スピードは許容範囲になっている。 mono と poly の対応は、上の２体の厳密解の精度由来の小さな差のみ。
  • 
    周期系のチェック。 Brady et al. (1988) の表１で比較。確認済み。