とっても大きな SDP その４

mDIM の値自体は以前解いていた SDP よりも小さいが、各行列のサイズが以下のように非常に大きいので、総合的に見ればこの SDP の方がサイズが大きいと言うことができる。

76554 (= mDIM)
22 (= nBLOCK)
18 18 18 18 153 153 324 153 153 324 648 324 324 816 2754 2754 816 8604 8604 2754 2754 -694 (= bLOCKsTRUCT)

現在は以下の指標まで進んでいる。あと数反復で終わりそうである。
38 1.5e-07 3.8e-11 4.9e-10 -9.10e+01 -9.10e+01 7.7e-01 5.8e-01 1.00e-01
39 7.0e-08 3.8e-11 2.0e-10 -9.10e+01 -9.10e+01 8.1e-01 5.4e-01 1.00e-01
40 3.6e-08 3.8e-11 8.9e-11 -9.10e+01 -9.10e+01 9.3e-01 6.2e-01 1.00e-01
41 1.6e-08 3.8e-11 3.6e-11 -9.10e+01 -9.10e+01 9.4e-01 7.2e-01 1.00e-01
42 5.6e-09 3.8e-11 1.1e-11 -9.10e+01 -9.10e+01 8.0e-01 6.3e-01 1.00e-01
43 2.4e-09 3.8e-11 3.8e-12 -9.10e+01 -9.10e+01 7.1e-01 5.4e-01 1.00e-01

メモリ的にはそろそろ限界に近づいているが、SDPARA を完全 64bit 化すれば、もっと大きな問題を扱うこともできるだろう。今は Schur complement 行列(mDIM x mDIM) を分割しているのだが、int = 32bit の制約から、この分割数を多くする必要があるので、メモリを多く消費する傾向にある。