ryuon
libiter の整備。
昨日の変更の house keeping など。
test でこけるものの調査。 どうも Bi-CGSTAB は実装に sensitive みたい。 厄介なのは一回発散してから再び(見掛け上)収束しているもの。 (暫定的に)収束だけでなく上限チェックを入れておく。
前処理付き CG をテスト環境から libiter に持ってきておく。
高速化の試み(前処理)。
LU 分解のテストコード(test-ilu.c)の整理。
これもベタァーっと書いていたものを、パーツごとに切り出し。
前処理チェック関数群も libiter/test 並にしておく。
Newton-GMRES の勉強