2008-06-24 ryuon phys 高速な implicit algorithm の実装: さて、今度こそ本当に実装するぞ! 作業場にて実装開始。 まずは3次方程式の解から……。えっと、千里の道も一歩から、と。 6/26/2008: バグっている。 丸一日かかって(家に持ち帰って)、何とか書き上げた。 とりあえず straight chain の決め打ち また、である……後で一般化するのが、また大変なのは分かっているが……。 あと、とりあえず noHI だけ。 ただし force law は ILC 以外6種類、一通り組みました。 対応する非線形解法アプローチも GSL multiroot 版と NITSOL 版とも組んだ。 今のところ、全て未チェック。 6/25/2008: チェック開始。