readsedumi.m

半定规划

%%*******************************************************************
%%  Convert problem from SeDuMi format to DSDP format
%%
%%  [AC,b] = readsedumi(At,b,C,K)
%%
%%  Input: At, b, C, K = Data in SeDuMi format.
%%
%% 
%% DSDP: version 5.0
%% Copyright (c) 2004 by
%% S. Benson Y. Ye
%% Last modified: 2 Jan 04
%%******************************************************************

  function [AC,b]=readsedumi(At,bb,c,K);
%%
%%  First, load the matlab file containing At, c, b, and K
%%
  if (size(c,1) == 1), c = c'; end;
  if (size(bb,1) == 1), bb = bb'; end;
  [nn,mm] = size(At); if (max(size(c)) == 1); c = c*ones(nn,1); end; 

  if ~isfield(K,'l'); K.l = 0; end   
  if ~isfield(K,'q'); K.q = 0; end
  if ~isfield(K,'s'); K.s = 0; end
  if K.l == 0 | isempty(K.l); K.l = 0; end;
  if sum(K.q) == 0 | isempty(K.q); K.q = 0; end
  if sum(K.s) == 0 | isempty(K.s); K.s = 0; end 
%%
%%
%%
  AC=cell(1,3);
block=1;
top=1;

   m = length(bb);
   b = -bb;

   if (K.l > 0) 
      AC{block,1} = 'LP';
      AC{block,2} = K.l;
      A = At(top:top+n-1,:);
      CC = c(top:top+n-1,:);
      AC{block,3} = [-A(:,1:m) CC(:,1)];
      block = block+1;
      top = top+K.l;
   end
   if (K.q > 0) 
      error(' Cannot accept SOCP cones\n');
   end

   if (K.s > 0) 
      for i = 1:length(K.s)
         n = K.s(i);
         AC{block,1} = 'SDP';
         AC{block,2} = n;
         A = At(top:top+n^2-1,:);
         CC = c(top:top+n^2-1,:);
         indicies = triu(reshape(1:n^2,n,n));
         indicies = indicies(find(indicies));
         AC{block,3} = [-A(indicies,1:m) CC(indicies,1)];
         block = block+1;
         top = top+n*n;
       end
   end  
 
%%