[37] | 1 | function YESNO = is(X,property,additional) |
---|
| 2 | %IS Check property of variable. |
---|
| 3 | % d = IS(x,property) returns 1 if 'property' holds |
---|
| 4 | % |
---|
| 5 | % Properties possible to test are: 'real', 'symmetric', 'hermitian', |
---|
| 6 | % 'scalar', 'linear', 'bilinear','quadratic','sigmonial', 'homogeneous', 'integer', 'binary' |
---|
| 7 | |
---|
| 8 | % Author Johan Löfberg |
---|
| 9 | % $Id: is.m,v 1.1 2006/08/10 18:00:20 joloef Exp $ |
---|
| 10 | |
---|
| 11 | switch property |
---|
| 12 | case 'logic' |
---|
| 13 | YESNO = X.typeflag==12; |
---|
| 14 | case 'binary' |
---|
| 15 | YESNO = any(ismember(depends(X),yalmip('binvariables'))); |
---|
| 16 | case 'integer' |
---|
| 17 | YESNO = any(ismember(depends(X),yalmip('intvariables'))); |
---|
| 18 | case 'real' |
---|
| 19 | YESNO = isreal(X); |
---|
| 20 | case 'complex' |
---|
| 21 | YESNO = ~isreal(X.basis); |
---|
| 22 | case 'symmetric' |
---|
| 23 | YESNO = issymmetric(X); |
---|
| 24 | case 'hermitian' |
---|
| 25 | YESNO = ishermitian(X); |
---|
| 26 | case 'scalar' |
---|
| 27 | YESNO = prod(X.dim)==1; |
---|
| 28 | case 'compound' |
---|
| 29 | YESNO = any(ismember(getvariables(X),yalmip('extvariables'))); |
---|
| 30 | case 'linear' |
---|
| 31 | variabletype = yalmip('variabletype'); |
---|
| 32 | variabletype = variabletype(X.lmi_variables); |
---|
| 33 | YESNO = ~any(variabletype); |
---|
| 34 | case 'bilinear' |
---|
| 35 | variabletype = yalmip('variabletype'); |
---|
| 36 | variabletype = variabletype(X.lmi_variables); |
---|
| 37 | YESNO = all(variabletype<=1); |
---|
| 38 | case 'quadratic' |
---|
| 39 | variabletype = yalmip('variabletype'); |
---|
| 40 | variabletype = variabletype(X.lmi_variables); |
---|
| 41 | YESNO = all(variabletype<=2); |
---|
| 42 | |
---|
| 43 | case 'LBQS' |
---|
| 44 | % Fast code for use in display etc. |
---|
| 45 | % Checks linearity, bilinearity etc in one call. |
---|
| 46 | quadratic = 0; |
---|
| 47 | bilinear = 0; |
---|
| 48 | linear = 0; |
---|
| 49 | sigmonial = 0; |
---|
| 50 | |
---|
| 51 | variabletype = yalmip('variabletype'); |
---|
| 52 | variabletype = variabletype(X.lmi_variables); |
---|
| 53 | |
---|
| 54 | linear = all(variabletype==0); |
---|
| 55 | bilinear = all(variabletype<=1) & ~linear; |
---|
| 56 | quadratic = all(variabletype<=2) & ~bilinear; |
---|
| 57 | sigmonial = any(variabletype==4); |
---|
| 58 | |
---|
| 59 | YESNO = full([linear bilinear quadratic sigmonial]); |
---|
| 60 | |
---|
| 61 | |
---|
| 62 | |
---|
| 63 | case 'lpcone' |
---|
| 64 | base = X.basis; |
---|
| 65 | YESNO = all(base(:,1)==0); % No constant |
---|
| 66 | base = base(:,2:end); |
---|
| 67 | YESNO = YESNO & all(sum(base,2)==1); |
---|
| 68 | YESNO = YESNO & all(sum(base,1)==1); |
---|
| 69 | YESNO = YESNO & all(sum(base~=0,2)==1); |
---|
| 70 | YESNO = YESNO & all(sum(base~=0,1)==1); |
---|
| 71 | YESNO = full(YESNO); |
---|
| 72 | |
---|
| 73 | case 'shiftlpcone' |
---|
| 74 | base = X.basis; |
---|
| 75 | base = base(:,2:end); |
---|
| 76 | YESNO = all(sum(base,2)==1); |
---|
| 77 | YESNO = YESNO & all(sum(base,1)==1); |
---|
| 78 | YESNO = YESNO & all(sum(base~=0,2)==1); |
---|
| 79 | YESNO = YESNO & all(sum(base~=0,1)==1); |
---|
| 80 | YESNO = full(YESNO); |
---|
| 81 | |
---|
| 82 | case 'sdpcone' |
---|
| 83 | if isequal(X.conicinfo,[1 0]) |
---|
| 84 | YESNO = 1; |
---|
| 85 | return |
---|
| 86 | end |
---|
| 87 | base = X.basis; |
---|
| 88 | n = X.dim(1); |
---|
| 89 | YESNO = full(issymmetric(X) & nnz(base)==n*n & all(sum(base,2)==1) & all(base(:,1)==0)) & length(X.lmi_variables)==n*(n+1)/2 & isreal(base); |
---|
| 90 | |
---|
| 91 | case 'shiftsdpcone' |
---|
| 92 | |
---|
| 93 | if isequal(X.conicinfo,[1 0]) |
---|
| 94 | YESNO = 1; |
---|
| 95 | return |
---|
| 96 | elseif isequal(X.conicinfo,[1 1]) |
---|
| 97 | YESNO = 1; |
---|
| 98 | return |
---|
| 99 | end |
---|
| 100 | |
---|
| 101 | |
---|
| 102 | base = X.basis; |
---|
| 103 | n = X.dim(1); |
---|
| 104 | base(:,1)=0; |
---|
| 105 | YESNO = full(issymmetric(X) & nnz(base)==n*n & all(sum(base,2)==1)) & length(X.lmi_variables)==n*(n+1)/2 & isreal(X); |
---|
| 106 | if YESNO |
---|
| 107 | % Possible case |
---|
| 108 | % FIX : Stupidly slow and complex |
---|
| 109 | [i,j,k] = find(base'); |
---|
| 110 | Y = reshape(1:n^2,n,n); |
---|
| 111 | Y = tril(Y); |
---|
| 112 | Y = (Y+Y')-diag(sparse(diag(Y))); |
---|
| 113 | [uu,oo,pp] = unique(Y(:)); |
---|
| 114 | YESNO = isequal(i,pp+1); |
---|
| 115 | % YESNO = isequal(base,getbase(sdpvar(n,n))); |
---|
| 116 | end |
---|
| 117 | |
---|
| 118 | case 'socone' |
---|
| 119 | base = X.basis; |
---|
| 120 | n = X.dim(1); |
---|
| 121 | YESNO = X.dim(1)>1 & X.dim(2)==1 & length(X.lmi_variables)==n; |
---|
| 122 | if YESNO |
---|
| 123 | cb = base(:,1); |
---|
| 124 | vb = base(:,2:end); |
---|
| 125 | YESNO = YESNO & (nnz(cb)==0) & (nnz(vb-speye(n))==0); |
---|
| 126 | end |
---|
| 127 | |
---|
| 128 | case 'sigmonial' |
---|
| 129 | monomtable = yalmip('monomtable'); |
---|
| 130 | monomtable = monomtable(getvariables(X),:); |
---|
| 131 | YESNO = any(find(any(0>monomtable,2) | any(monomtable-fix(monomtable),2))); |
---|
| 132 | |
---|
| 133 | case 'general' |
---|
| 134 | evalvariables = yalmip('evalVariables'); |
---|
| 135 | YESNO = ~isempty(intersect(getvariables(X),evalvariables)); |
---|
| 136 | |
---|
| 137 | case 'nonlinear' |
---|
| 138 | YESNO = ~islinear(X); |
---|
| 139 | |
---|
| 140 | case 'homogeneous' |
---|
| 141 | [sqrList,CompressedList] = yalmip('nonlinearvariables'); |
---|
| 142 | [LinearTerms,NonlinearVariables] = getvariables(X,'both'); |
---|
| 143 | if isempty(NonlinearVariables) |
---|
| 144 | YESNO = nnz(getbasematrix(X,0))==0; |
---|
| 145 | else |
---|
| 146 | % No linear terms |
---|
| 147 | YESNO = isempty(LinearTerms); |
---|
| 148 | % No constant terms |
---|
| 149 | YESNO = YESNO & (nnz(getbasematrix(X,0))==0); |
---|
| 150 | % Largest degree+1 |
---|
| 151 | maxdegree = sum(any(CompressedList(find(ismember(CompressedList,NonlinearVariables)),:),1)); |
---|
| 152 | % All same degree |
---|
| 153 | YESNO = YESNO & all(all(CompressedList(find(ismember(CompressedList,NonlinearVariables)),2:maxdegree)>0)); |
---|
| 154 | end |
---|
| 155 | |
---|
| 156 | case 'sos' |
---|
| 157 | YESNO = (X.typeflag==11); |
---|
| 158 | |
---|
| 159 | otherwise |
---|
| 160 | error('Wrong input argument.'); |
---|
| 161 | end |
---|
| 162 | |
---|
| 163 | YESNO = full(YESNO); |
---|
| 164 | |
---|
| 165 | |
---|