やっぱり自分で書いたライブラリーが一番使いやすい。今までヤコビとかで頑張っていたんだけど、compressed row storaged ( CRS ) Bi Conjugate Gradient Stable ソルバーつっくったよ。
めっさ早い。
256x256の７重バントマトリックスかでも50スッテプ以内で収束する。さっきやっとこさOpen MPで並列化が完了したさ。
前処理も時間があったらしてみたい。
09dec2011追記
ソースをホームページにのけってみた。cでかいてあります。makeすれば動くはず
0 件のコメント:
コメントを投稿