データ整理とSOCP

今日は、この前の論文の続きで、SOCP で組んだプログラムで、どの程度の精度が出せるかを調べることにした。
特化したアルゴリズムを作るのもいいが、まずは SOCP でどのくらいできるのか判断したほうがいいかと考え直した。

データについては、インターネット経由で探したところ、9個中8個までは揃った。あとひとつは、論文までは特定できるがその論文のどのデータを使ったかを特定するのが大変なので、元論文の人にメールで聞いてみようと思う。

ただ、データについて整理して気がついたが、元論文の表はあちこちに違いがあるようだ。たとえば、130個しかないデータに対して、140個のデータとして計算しているようなところがある。単純に TeX に書くときに間違えたことも考えられるし、元論文の後にデータが更新されたこともありえる。
いずれにしても、結果をストレートに比較できない点には後で注意が必要そうだ。

プログラムについては、Xeon で走らせているので、あとはコンピュータに任せておいて、明日は別の仕事を進めることにしようと思う。
途中までの結果を見ると、ペナルティパラメータを大きくすると、計算時間がかかるようになるだけでなく、誤差も大きくなっている。
現在のプログラムでは、終了条件を目的関数値の変動としているが（前の反復と比較して1e-3など）、これは相対値に変更したほうがいいかもしれない。
いずれにしても、現在のプログラムの結果を見てから判断しよう。

今日の内容：SOCP のプログラム
今日のBGM: FF8 OST [1-4]
今日のランチ：つかさ　まかない丼
明日の予測作業時間：6h