1998-06-17 m-of-r periodic でうまく行かない問題解決! 原因は波数部の self part M^(2)(r=0) を補正項(iteration matrix) でない方に入れていた事。補正項に入れるとうまくいった。 一方 NR の biconjugate gradient method を使うと iteration 回数が1で収束していた! -> 使えそう。今後の課題。 GMRES の implementation はまだ。