YALMIP-R20160930.zip

function varargout = yalmip(varargin) %YALMIP Returns various information about YALMIP % % YALMIP can be used to check version numbers and find the SDPVAR and SET % objects available in workspace % % EXAMPLES % V = YALMIP('version') % Returns version % YALMIP('nvars') % Returns total number of declared variables % % If you want information on how to use YALMIP % http://yalmip.github.io % % See also YALMIPTEST persistent prefered_solver internal_sdpvarstate internal_setstate if nargin==0 help yalmip return end if isempty(internal_sdpvarstate) if exist('OCTAVE_VERSION', 'builtin') more off end internal_sdpvarstate.monomtable = spalloc(0,0,0); % Polynomial powers table internal_sdpvarstate.hashedmonomtable = []; % Hashed polynomial powers table internal_sdpvarstate.hash = []; internal_sdpvarstate.boundlist = []; internal_sdpvarstate.variabletype = spalloc(0,0,0); % Pre-calc linear/quadratic/polynomial/sigmonial internal_sdpvarstate.intVariables = []; % ID of integer variables internal_sdpvarstate.binVariables = []; % ID of binary variables internal_sdpvarstate.semicontVariables = []; internal_sdpvarstate.uncVariables = []; % ID of uncertain variables (not used) internal_sdpvarstate.parVariables = []; % ID of parametric variables (not used) internal_sdpvarstate.extVariables = []; % ID of extended variables (for max,min,norm,sin, etc) internal_sdpvarstate.auxVariables = []; % ID of auxilliary variables (introduced when modelling extended variables) internal_sdpvarstate.auxVariablesW = []; % ID of uncertain auxilliary variables (introduced when modelling uncertain extended variables) internal_sdpvarstate.logicVariables = []; % ID of extended logic variables (for or, nnz, alldifferent etc) internal_sdpvarstate.complexpair = []; internal_sdpvarstate.internalconstraints = []; internal_sdpvarstate.ExtendedMap = []; internal_sdpvarstate.ExtendedMapHashes = []; internal_sdpvarstate.DependencyMap = sparse(0); internal_sdpvarstate.DependencyMapUser = sparse(0); internal_sdpvarstate.sosid = 0; internal_sdpvarstate.sos_index = []; internal_sdpvarstate.sos_data = []; internal_sdpvarstate.sos_ParV = []; internal_sdpvarstate.sos_Q = []; internal_sdpvarstate.sos_v = []; internal_sdpvarstate.optSolution{1}.info = 'Initialized by YALMIP'; internal_sdpvarstate.optSolution{1}.variables = []; internal_sdpvarstate.optSolution{1}.optvar =[]; internal_sdpvarstate.optSolution{1}.values =[]; internal_sdpvarstate.activeSolution = 1; internal_sdpvarstate.nonCommutingTable = []; internal_sdpvarstate.nonHermitiannonCommutingTable = []; internal_sdpvarstate.containsSemivar = false; try warning off Octave:possible-matlab-short-circuit-operator catch end end if isempty(internal_setstate) internal_setstate.LMIid = 0; internal_setstate.duals_index = []; internal_setstate.duals_data = []; internal_setstate.duals_associated_index = []; internal_setstate.duals_associated_data = []; end switch varargin{1} case 'clearsolution' internal_sdpvarstate.optSolution{1}.variables = []; internal_sdpvarstate.optSolution{1}.optvar =[]; internal_sdpvarstate.optSolution{1}.values =[]; internal_sdpvarstate.activeSolution = 1; case 'monomtable' varargout{1} = internal_sdpvarstate.monomtable; n = size(internal_sdpvarstate.monomtable,1); if size(internal_sdpvarstate.monomtable,2) < n % Normalize the monomtalbe. Some external functions presume the % table is square varargout{1}(n,n) = 0; internal_sdpvarstate.monomtable = varargout{1}; need_new = size(internal_sdpvarstate.monomtable,1) - length(internal_sdpvarstate.hash); internal_sdpvarstate.hash = [internal_sdpvarstate.hash ; 3*gen_rand_hash(size(internal_sdpvarstate.monomtable,1),need_new,1)]; internal_sdpvarstate.hashedmonomtable = internal_sdpvarstate.monomtable*internal_sdpvarstate.hash; end if nargout == 2 varargout{2} = internal_sdpvarstate.variabletype; elseif nargout == 4 varargout{2} = internal_sdpvarstate.variabletype; varargout{3} = internal_sdpvarstate.hashedmonomtable; varargout{4} = internal_sdpvarstate.hash; end case 'setmonomtable' % New monom table internal_sdpvarstate.monomtable = varargin{2}; if nargin>=4 % User has up-dated the hash tables him self. internal_sdpvarstate.hashedmonomtable=varargin{4}; internal_sdpvarstate.hash = varargin{5}; end if size(internal_sdpvarstate.monomtable,2)>length(internal_sdpvarstate.hash) need_new = size(internal_sdpvarstate.monomtable,1) - length(internal_sdpvarstate.hash); internal_sdpvarstate.hash = [internal_sdpvarstate.hash ; 3*gen_rand_hash(size(internal_sdpvarstate.monomtable,1),need_new,1)]; end if size(internal_sdpvarstate.monomtable,1)>size(internal_sdpvarstate.hashedmonomtable,1) % Need to add some hash values need_new = size(internal_sdpvarstate.monomtable,1) - size(internal_sdpvarstate.hashedmonomtable,1); temp = internal_sdpvarstate.monomtable(end-need_new+1:end,:); internal_sdpvarstate.hashedmonomtable = [internal_sdpvarstate.hashedmonomtable;temp*internal_sdpvarstate.hash]; end if nargin >= 3 && ~isempty(varargin{3}) internal_sdpvarstate.variabletype = varargin{3}; if length(internal_sdpvarstate.variabletype) ~=size(internal_sdpvarstate.monomtable,1) error('ASSERT') end else internal_sdpvarstate.variabletype = zeros(size(internal_sdpvarstate.monomtable,1),1)'; nonlinear = ~(sum(internal_sdpvarstate.monomtable,2)==1 & sum(internal_sdpvarstate.monomtable~=0,2)==1); if ~isempty(nonlinear) %mt = internal_sdpvarstate.monomtable; internal_sdpvarstate.variabletype(nonlinear) = 3; quadratic = sum(internal_sdpvarstate.monomtable,2)==2; internal_sdpvarstate.variabletype(quadratic) = 2; bilinear = max(internal_sdpvarstate.monomtable,[],2)<=1; internal_sdpvarstate.variabletype(bilinear & quadratic) = 1; sigmonial = any(0>internal_sdpvarstate.monomtable,2) | any(internal_sdpvarstate.monomtable-fix(internal_sdpvarstate.monomtable),2); internal_sdpvarstate.variabletype(sigmonial) = 4; end end case 'variabletype' varargout{1} = internal_sdpvarstate.variabletype; case {'addextendedvariable','addEvalVariable'} varargin{2} disp('Obsolete use of the terms addextendedvariable and addEvalVariable'); error('Obsolete use of the terms addextendedvariable and addEvalVariable'); case 'defineVectorizedUnitary' varargin{2} = strrep(varargin{2},'sdpvar/',''); % Clean due to different behaviour of the function mfilename in ML 5,6 and 7 % Is this operator variable already defined correct_operator = []; if ~isempty(internal_sdpvarstate.ExtendedMap) OperatorName = varargin{2}; Arguments = {varargin{3:end}}; this_hash = create_trivial_hash(firstSDPVAR(Arguments)); correct_operator = find([internal_sdpvarstate.ExtendedMapHashes == this_hash]); if ~isempty(correct_operator) correct_operator = correct_operator(strcmp(OperatorName,{internal_sdpvarstate.ExtendedMap(correct_operator).fcn})); end