2012年7月20日金曜日

cholmod_solve は自動的にはorderingを調整しない

CHOLMOD のデモプログラムを見るかぎりだと、cholmod_solve は右辺ベクトルと解ベクトルにordering を自動的に適用するように見えるが、調べてみるとそうなっていない。
このあたりは、ユーザマニュアルには書かれていないので、実際に実行して数値的に確認する以外に方法はない。
おそらく、ユーザマニュアルが実際のプログラムをきちんと反映できていないのであろう。
ユーザマニュアルを整備するのはプログラムを書くのと同じぐらい大変なので、このあたりの事情はよく分かる。

とりあえず、ここのところは修正したが、なぜかステップ長が本来の長さよりも短く算出されてしまう。
6次元の行列程度の問題なのに80反復ぐらいかかってしまう。
このバグを取り除くのが次の仕事になりそうだ。

今日の作業内容:SDPA-C 5h
今日のランチ:サイゼリヤ みぞれハンバーグ&牛肉コロッケ
明日の予測作業時間:4h


0 件のコメント:

コメントを投稿