ryuon

• 
  Brownian dynamics
  • 
    cholesky と chebyshev 及び matrix と atimes の benchmark をしておく。
    • 
      時間積分は semi-implicit, Pe = 1 で flow は、なし。
    • 
      結果： 確かに N>21 で chebyshev-atimes が一番速くなっているが、 スケーリングは N³ みたい。
    • 
      挙動が理解できないので、 RYUON/stokes/src/bench3 の結果を一応確認。 結果はどちらも N² だ。
    • 
      で、少し考えて分かった。 F の no-lub の mob は、 線形方程式を解かないので、そもそも N² なんだ。
    • 
      つまり、非線形方程式の解法で N¹ を食っているということらしい。