GPBiCG ソルバー(2)

やっとバグがとれた。今度はちゃんと収束するようになったけど、bicgstabよりも誤差が大きいのが気になる。

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

引き継がれているコードが

研究室で開発されたコードを使って研究してるんだけど、かなり怪しい部分があるのに気づいた。
本でチェックしたら、何か違うんですけど。検証結果とか残ってないのかな？

正直、直すの面倒くさい。

また学会が増えた

中国とイギリスの学会にいけとのお達示が合って大変だった。まだ、まとまってないのに。
取り合えず、中国の方は審査を通った。

それよりも、早く論文投稿させてくれ。

GPBiCGソルバーを作ってみた

何か、かなり安定しないんですけど。BiCGStabで解ける行列も解けなかったりする。 追記: ソースコードを公開しました。http://a-daily-life-in-the-office.blogspot.jp/2012/04/blog-post_26.html

paraviewのへんな仕様見つけた

connectivity
二次元では0から始めないとダメ
三次元では1からはじめないとダメ

三次元の書いた奴Fortran厨だろ

what's the hell is this?

If you writing paraview vtk output file, be careful for your connectivity table
-in 2D, the smallest node number must be zero
-in 3D, the smallest node number must be one

The 3d module is apparently written by Fortran lover.

bicgstabをpthreadで無駄に並列化してみた

ソース乗っけとく

pthread_bicgstab

Fortranで文字列処理

やりたくないこと＃１だよ。マジ死ぬ