-
Notifications
You must be signed in to change notification settings - Fork 0
/
testKMWalt.m
117 lines (91 loc) · 2.65 KB
/
testKMWalt.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
function testKMWalt(userFunction)
tol=1e-6;
n=7;
eps=.17;
testA = ...
[1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0;
-1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0;
0.17 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0;
0.17 -1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0;
0 0.17 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0;
0 0.17 -1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0;
0 0 0.17 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0;
0 0 0.17 -1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0;
0 0 0 0.17 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0;
0 0 0 0.17 -1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0;
0 0 0 0 0.17 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0;
0 0 0 0 0.17 -1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0;
0 0 0 0 0 0.17 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0;
0 0 0 0 0 0.17 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 1];
testb = [1; -0.17; 1; 0; 1; 0; 1; 0; 1; 0; 1; 0; 1; 0];
testB = [1 2 3 4 5 6 7 8 10 12 14 16 18 20];
testc = [0; 0; 0; 0; 0; 0; -1; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; ...
0];
[A,b,c,B]=userFunction(n,eps);
resA=norm(A-testA);
resb=norm(b-testb);
resc=norm(c-testc);
resB=numel(setdiff(testB,B));
fprintf('\n Teste mit n=%d und eps= %g \n',n,eps);
fprintf('------------------------------------\n');
if resA<tol
fprintf('Matrix: OK \n');
else
fprintf('Matrix: falsch, Fehler: %d \n', resA);
end
if resb<tol
fprintf('rechte Seite: OK \n');
else
fprintf('rechte Seite: falsch, Fehler: %d \n', resb);
end
if resc<tol
fprintf('Zielfunktion: OK \n');
else
fprintf('Zielfunktion: falsch, Fehler: %d \n', resc);
end
if resB==0
fprintf('Startbasis: OK \n');
else
fprintf('Startbasis: falsch, Fehler: %d Einträge \n', resB);
end
n=4;
eps=.42;
testA = ...
[1 0 0 0 1 0 0 0 0 0 0 0;
-1 0 0 0 0 1 0 0 0 0 0 0;
0.42 1 0 0 0 0 1 0 0 0 0 0;
0.42 -1 0 0 0 0 0 1 0 0 0 0;
0 0.42 1 0 0 0 0 0 1 0 0 0;
0 0.42 -1 0 0 0 0 0 0 1 0 0;
0 0 0.42 1 0 0 0 0 0 0 1 0;
0 0 0.42 -1 0 0 0 0 0 0 0 1];
testb = [1; -0.42; 1; 0; 1; 0; 1; 0];
testB = [1 2 3 4 5 7 9 11];
testc = [0; 0; 0; -1; 0; 0; 0; 0; 0; 0; 0; 0];
[A,b,c,B]=userFunction(n,eps);
resA=norm(A-testA);
resb=norm(b-testb);
resc=norm(c-testc);
resB=numel(setdiff(testB,B));
fprintf('\n Teste mit n=%d und eps= %g \n',n,eps);
fprintf('------------------------------------\n');
if resA<tol
fprintf('Matrix: OK \n');
else
fprintf('Matrix: falsch, Fehler: %d \n', resA);
end
if resb<tol
fprintf('rechte Seite: OK \n');
else
fprintf('rechte Seite: falsch, Fehler: %d \n', resb);
end
if resc<tol
fprintf('Zielfunktion: OK \n');
else
fprintf('Zielfunktion: falsch, Fehler: %d \n', resc);
end
if resB==0
fprintf('Startbasis: OK \n');
else
fprintf('Startbasis: falsch, Fehler: %d Einträge \n', resB);
end