twobody

• 
  slip の厳密解
  • 
    C 版を開始。
    • 
      maxima は遅くてかつ低次で止まるので、 数値的に lubrication の捻出が不可能なため。
    • 
      最初、 PARI/GP で組み直してみたが、速度は体感的には変わらず (コンパイラは試してない)。
    • 
      なので、結局 C 。それでも最初 GMP を考えたが、多項式を目の前に断念。
    • 
      ってことで、普通に浮動少数点。ついでに気になってた long double 版も書いてみた。
      • 
        FreeBSD は 何か変態なこと になっているようで、 fpsetprec で明示的に仮数部を 64 bit に指定する必要がある。 あとpowlが見つからず、もう面倒なので、自分で普通に書いた、 どうせ整数しか引き数とらないし。
    • 
      速度は、しかし期待した程、速くない…… メモリー不足で止まらないのが advantage か。
    • 
      簡単な cache も付けた。実計算では効果ないが、再計算が楽になる。
  • 
    lubrication
    • 
      何気なく Rev.Mod.Phys. 見てたら Goldenfeld の review を見つける。 彼が何で stokes flow? と思ったら、(その昔、近所で少し流行ってた)例の RG のテクニックを stokes flow に使ってみた、という話らしい。
    • 
      寝る前に横になって斜め読みしたところ、 stokes じゃなくて、 stokes 近似の有限 Re への展開の話。 例の Proudman-Peason 辺りの話。 technical には流体で一つの潮流をなす boundary layer とかの matched asymptotics の話 (なので、 lubrication theory には実はもろにかぶってくる話)。