Line | |
---|
1 | function massive |
---|
2 | ops{1} = sdpsettings('sos.cong',0,'sos.model',1,'verbose',0); |
---|
3 | ops{2} = sdpsettings('sos.cong',1,'sos.model',2,'verbose',0); |
---|
4 | ops{3} = sdpsettings('sos.cong',0,'sos.newton',0,'verbose',0,'sos.extlp',0); |
---|
5 | |
---|
6 | x = sdpvar(1,1); |
---|
7 | y = sdpvar(1,1); |
---|
8 | t = sdpvar(1,1); |
---|
9 | s = sdpvar(1,1); |
---|
10 | F = set(sos(1+y^2-t))+set(sos(1+x^2-s))+set(t>0)+set(s>0.5); |
---|
11 | obj = -s-t; |
---|
12 | for i = 1:length(ops) |
---|
13 | i |
---|
14 | fail = regresstest(F,obj,ops{i}); |
---|
15 | mbg_asserttolequal(fail,0); |
---|
16 | end |
---|
17 | |
---|
18 | |
---|
19 | function fail = regresstest(F,obj,ops,pv); |
---|
20 | |
---|
21 | if nargin==3 |
---|
22 | pv = []; |
---|
23 | end |
---|
24 | |
---|
25 | ops.sos.model = 1; |
---|
26 | solvesos(F,obj,ops,pv); |
---|
27 | obj1 = double(obj); |
---|
28 | p1s = checkset(F(find(is(F,'sos')))); |
---|
29 | p1e = checkset(F(find(~is(F,'sos')))); |
---|
30 | |
---|
31 | ops.sos.model = 2; |
---|
32 | solvesos(F,obj,ops,pv); |
---|
33 | obj2 = double(obj); |
---|
34 | p2s = checkset(F(find(is(F,'sos')))); |
---|
35 | p2e = checkset(F(find(~is(F,'sos')))); |
---|
36 | |
---|
37 | fail = 0; |
---|
38 | |
---|
39 | if abs(obj1-obj2) > 1e-4 |
---|
40 | fail = 1; |
---|
41 | end |
---|
42 | |
---|
43 | if any(p1s>1e-4) |
---|
44 | fail = 2; |
---|
45 | p1s |
---|
46 | end |
---|
47 | if any(p2s>1e-4) |
---|
48 | fail = 2; |
---|
49 | p2s |
---|
50 | end |
---|
51 | if any(p1e<-1e-4) |
---|
52 | fail = 2; |
---|
53 | p1e |
---|
54 | end |
---|
55 | if any(p2e<-1e-4) |
---|
56 | fail = 2; |
---|
57 | p2e |
---|
58 | end |
---|
59 | if fail==0 |
---|
60 | disp('Correct solution'); |
---|
61 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.