2018年1月4日木曜日

Matlab 2017bでYALMIPとSeDuMiを使えるようにした

Matlab 2017bでYALMIP, SeDuMiを使えるようにしたメモ.

環境

  • macOS: 10.11.6-x86_64
  • Xcode: 8.2.1

YALMIP

サイト(https://yalmip.github.io/download/)からダウンロードして解凍する.適当な場所に置く(~/MATLAB/YALMIP-master など).
Matlabでパスを通す.
addpath(genpath('/Users/hoge/Documents/MATLAB/YALMIP-master'))
savepath

SeDuMi

サイト(https://github.com/sqlp/sedumi)からダウンロード.macOS 64bit用のバイナリ(**.mexmaci64)を適当な場所に配置してパスを通す.
テストしてみるとSDPが解けていない.
yalmiptest
|                    SDP|        N/A|   Unknown problem in solver (Turn on 'debug' in sdpsettings) (関数または変数 'y_s' が未定義です。)|

SeDuMiバイナリの作成

どうやらsedumiがインストールできていないよう.
sedumi([1; 1], 1, [1 1])
無効な MEX ファイル '/Users/hoge/Documents/MATLAB/sedumi-master/whichcpx.mexmaci64':
dlopen(/Users/hoge/Documents/MATLAB/sedumi-master/whichcpx.mexmaci64, 6): Library not loaded:
@loader_path/libmwlapack.dylib
  Referenced from: /Users/hoge/Documents/MATLAB/sedumi-master/whichcpx.mexmaci64
  Reason: image not found。

エラー: pretransfo (line 211)
cpx = whichcpx(K);

エラー: sedumi (line 261)
[A,b,c,K,prep,origcoeff] = pretransfo(A,b,c,K,pars);
mexファイルがおかしいと言われるのでコンパイルしてみる.github上のソースファイルから修正が必要.ここを参考に:
Found the link for the fix.
-#define FORT(x) FORTRAN_WRAPPER(x)
+#define FORT(x) x
at line at line 45 in blksdp.h
修正したらSeDuMiのファイル直下でinstall_sedumiをオプション付きで実行する.
install_sedumi -rebuild
できた.
yalmiptest
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|                   Test|   Solution|                     Solver message|
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|   Core functionalities|        N/A|       Successfully solved (YALMIP)|
|                     LP|    Correct|      Successfully solved (LINPROG)|
|                     LP|    Correct|      Successfully solved (LINPROG)|
|                     QP|    Correct|     Successfully solved (QUADPROG)|
|                     QP|    Correct|     Successfully solved (QUADPROG)|
|                   SOCP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                   SOCP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                   SOCP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                    SDP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                    SDP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                    SDP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                    SDP|    Correct|   Successfully solved (SeDuMi-1.3)|
|                 MAXDET|    Correct|   Successfully solved (SeDuMi-1.3)|
|                 MAXDET|    Correct|   Successfully solved (SeDuMi-1.3)|
|          Infeasible LP|        N/A|       Infeasible problem (LINPROG)|
|          Infeasible QP|        N/A|      Infeasible problem (QUADPROG)|
|         Infeasible SDP|        N/A|    Infeasible problem (SeDuMi-1.3)|
|      Moment relaxation|    Correct|   Successfully solved (SeDuMi-1.3)|
|         Sum-of-squares|    Correct|   Successfully solved (SeDuMi-1.3)|
|           Bilinear SDP|        N/A|                 No suitable solver|
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
設定
Note that YALMIP uses historical default settings for SeDuMi. According to the user guide for SeDuMi 1.1, the YALMIP defaults should be updated as below:
S = sdpsettings('solver', 'sedumi', 'sedumi.eps', 1e-8, ...
                'sedumi.cg.qprec', 1, 'sedumi.cg.maxiter', 49, ...
                'sedumi.stepdif', 2);