MATLAB で SDPT3 と SeDuMi

MATLAB Version 7.5.0.338で SDPT3(4.0beta) と SeDuMi (1.1)を使ってみよう。

1: SDPT3 4.0beta

> Installmex
> startup
> [blk,At,C,b] = read_sdpa('./sdplib/theta5.dat-s');
> [obj,X,y,Z] = sqlp(blk,At,C,b);

2: SeDuMi 1.1

SeDuMi の最新版は 64 bit版の Matlab では動作しない。理由は以下の通り。でも SDPT3 は 64bit 対応済み。昨年カナダで会ったときに作者に直接聞いたら SDPT3 の方は１日ぐらいで対応できたそうだ。

NOTE2: Matlab is now using a different interface for sparse matrices in 64 bit versions. They changed it without advance notice and there is no backward compatible workaround. Thus, SeDuMi will NOT work in 64 bit versions of Matlab under linux. We are testing the Windows versions right now.

ということで、まず Matlab を 32bit で起動する(SDPA フォーマットファイルの実行)

> mex -setup
> matlab -glnx86
> install_sedumi
> fname='./sdplib/mcp500-1.dat-s';
> [At,b,c,K]=fromsdpa(fname);
> [x,y,info]=sedumi(At,b,c,K);