やっとバグがとれた。今度はちゃんと収束するようになったけど、bicgstabよりも誤差が大きいのが気になる。
残差のnorm、解のnormで収束判定するのかいつも悩む。大体は解の相対誤差で判定しているけれども、一つ余計にベクトルがひつようになる。残差と解のnormが厳密にリンクしてないスキームとかもあるらしいので、ソルバーの専門家出ないので怖くて手が出せない。
追記: ソースコードを公開しました。http://a-daily-life-in-the-office.blogspot.jp/2012/04/blog-post_26.html
0 件のコメント:
コメントを投稿