ryuon

• 
  高速な implicit algorithm の実装：
  • 
    前回のテストの続き：
    • 
      N=100 で NITSOL がこける問題、解決。
      • 
        問題は、またもや config の設定 （else if すべき所が一ヶ所 if だけになっていた、というおバカなもの）。
    • 
      以上で no-HI 版のテストは終了とする。
  • 
    HI 化：
    • 
      ノートに戻って、定式化を確認。
      • 
        ふむ、 right-hand side vector のみの変更で対応できる。
    • 
      しかるべく、実装。
      • 
        できた。 straight chain の F の決め打ちで、しかも non periodic だけ。
      • 
        ついでなので corrector step の初期と update の統合もしておく。 （詳しい理由は知らないけど、微妙に違うんだよな。 この違いについては、後でちょっと調べておこう。）
      • 
        これまでのチェックで通ることを確認ののち、古い（非統合）コードは破棄。
    • 
      これで、一応 HI 化はできたはす。
      • 
        不十分な点がもう一つ、 connector から position を作るのに、 COM を手で与えている。
    • 
      テスト：
      • 
        既存のテストに flag_noHI を追加。 一通り走るようになった！！
      • 
        と思っていたら、 何かの加減で check_fastSI_solve_cubic が文句を言うようになった。
        • 
          見ると、いろいろ下らない bug や、思い違いがあった。
        • 
          それらを修正後も、まだどこかで nan になっている模様……。
      • 
        7/3/2008: テストの続き。
    • 
      今月はここまで。