ryuon

• 
  高速な implicit algorithm の実装：
  • 
    昨日の続き、テストのテスト。
    • 
      6/27/2008: さらに続き。
  • 
    全然テストになっていなかった （つまり、勝率５割なんて大嘘だったということ、やれやれ）。
  • 
    一般的な非線形解法の debug は、終わったかな。
  • 
    一方、高速 scheme は全然ダメであることが判明：
    • 
      調べてみたら、３次方程式が正しく解けていない。
    • 
      自分でコードを書いているのだが（6/24/2008）、そいつがうまく動いていない模様。参ったね。
    • 
      結局、これで一日潰す。しかもまだ動いていない……。
    • 
      帰宅後、 GSL を覗くとあったよ。 gsl_poly_solve_cubic. こいつは（当然）きちんと解けているな。