[37] | 1 | function m = monolistcoeff(n,d,summax) |
---|
| 2 | %MONOLISTCOEFF Internal function used in SOS programs |
---|
| 3 | |
---|
| 4 | % Author Johan Löfberg |
---|
| 5 | % $Id: monolistcoeff.m,v 1.5 2005/02/20 17:04:28 johanl Exp $ |
---|
| 6 | |
---|
| 7 | if length(d)==1 |
---|
| 8 | d = repmat(full(d),1,full(n)); |
---|
| 9 | end |
---|
| 10 | |
---|
| 11 | if nargin==2 |
---|
| 12 | summax = (sum(d)); |
---|
| 13 | end |
---|
| 14 | |
---|
| 15 | if 0 |
---|
| 16 | dmax = sum(d);%max(full(d));%sum(d); |
---|
| 17 | base = eye(n); |
---|
| 18 | m = base; |
---|
| 19 | for i = 1:dmax-1 |
---|
| 20 | temp=[]; |
---|
| 21 | for k = 1:n |
---|
| 22 | temp = [temp;m+repmat(base(k,:),size(m,1),1)]; |
---|
| 23 | end |
---|
| 24 | ii=find(~any(temp-repmat(d,size(temp,1),1)>0,2)); |
---|
| 25 | temp=temp(ii,:); |
---|
| 26 | ii=find(~any(sum(temp,2)-summax>0,2)); |
---|
| 27 | temp=temp(ii,:); |
---|
| 28 | |
---|
| 29 | m = [m;temp]; |
---|
| 30 | [ii,jj,kk]=uniquesafe(m,'rows'); |
---|
| 31 | m = m(jj,:); |
---|
| 32 | end |
---|
| 33 | % Add constant |
---|
| 34 | m = [zeros(1,n);m]; |
---|
| 35 | else |
---|
| 36 | m = (0:d(1))'; |
---|
| 37 | for i = 2:n |
---|
| 38 | z = (0:d(i))'; |
---|
| 39 | m = [kron(m,ones(size(z,1),1)) kron(ones(size(m,1),1),z)]; |
---|
| 40 | [ii,jj,kk]=uniquesafe(m,'rows'); |
---|
| 41 | m = m(jj,:); |
---|
| 42 | m = m(sum(m,2)<=summax,:); |
---|
| 43 | end |
---|
| 44 | end |
---|
| 45 | m = full(m); |
---|