ALF project

• 
  Stokesilet in a spherical vessel, その３
  • 
    programming、 仕上げる
  • 
    テスト開始
    • 
      Ji 非対称だが、いいのか?
    • 
      球面上で速度を計算、 zero になってない、 parallel 部はok, 問題は perpendicular 部, (お粗末な) bug 発見
  • 
    テスト終了 (Rev 1.3)
  • 
    ついでに論文の typo を確認
    • 
      (5) 式の perpendicular 部(δ_km-e_ke _m) 第一項の " s " はいらない
    • 
      まぁ最初から typo だと思って coding してたが
  • 
    例外処理
    • 
      y = 0 の時、厳密を返す
    • 
      y // x の時、 漸近挙動…詰まった、 とりあえず TINY を噛ます
    • 
      Rev 1.4
  • 
    テスト、 y = 0 の時は ok, y // x の時、 ときどきおかしい…