1 | function varargout=alldifferent(varargin) |
---|
2 | |
---|
3 | switch class(varargin{1}) |
---|
4 | |
---|
5 | case 'double' |
---|
6 | x = varargin{1}; |
---|
7 | x = sort(x(:)); |
---|
8 | varargout{1} = all(diff(x) > 0); |
---|
9 | |
---|
10 | case 'sdpvar' |
---|
11 | varargout{1} = set(yalmip('addextendedvariable',mfilename,varargin{:}) == 1); |
---|
12 | |
---|
13 | case 'char' |
---|
14 | x = varargin{3}; |
---|
15 | [nx,mx] = size(x); |
---|
16 | x = reshape(x,nx*mx,1); |
---|
17 | |
---|
18 | [M,m] = derivebounds(x); |
---|
19 | if 0 & all(M(1) == M) & all(m(1) == m) & (nx*mx == (M(1)-m(1) + 1)) |
---|
20 | % No benefit seen in this model in this special case so we |
---|
21 | % don't use it to (reduce possibility of introducing a new bug) |
---|
22 | d = binvar(nx*mx,nx*mx,'full'); |
---|
23 | F = set(d*(m(1):M(1))' == x) + set(sum(d,1) == 1) + set(sum(d,2) == 1); |
---|
24 | else |
---|
25 | % Add constraint |x(i)-x(j)| > 1 |
---|
26 | pairs = nchoosek(1:nx*mx,2); |
---|
27 | d = binvar(length(pairs),1); |
---|
28 | x1 = x(pairs(:,1)); |
---|
29 | x2 = x(pairs(:,2)); |
---|
30 | |
---|
31 | F = set(x1 - x2 > 1-(1+M(pairs(:,1))-m(pairs(:,2))).*d); |
---|
32 | F = F + set(x2 - x1 > 1-(1+M(pairs(:,2))-m(pairs(:,1))).*(1-d)); |
---|
33 | end |
---|
34 | varargout{1} = F; |
---|
35 | varargout{2} = struct('convexity','milp','monotonicity','milp','definiteness','milp','extra','marker'); |
---|
36 | varargout{3} = varargin{3}; |
---|
37 | end |
---|