1 | |
---|
2 | /*! |
---|
3 | \file |
---|
4 | \brief Robust |
---|
5 | \author Vasek Smidl |
---|
6 | |
---|
7 | */ |
---|
8 | |
---|
9 | #include "robustlib.h" |
---|
10 | #include <vector> |
---|
11 | using namespace itpp; |
---|
12 | |
---|
13 | //using namespace bdm; |
---|
14 | |
---|
15 | const int emlig_size = 2; |
---|
16 | |
---|
17 | |
---|
18 | int main ( int argc, char* argv[] ) { |
---|
19 | |
---|
20 | |
---|
21 | emlig* emlig1 = new emlig(emlig_size); |
---|
22 | |
---|
23 | /* |
---|
24 | emlig1->set_correction_factors(4); |
---|
25 | |
---|
26 | for(int j = 0;j<emlig1->correction_factors.size();j++) |
---|
27 | { |
---|
28 | for(set<my_ivec>::iterator vec_ref = emlig1->correction_factors[j].begin();vec_ref!=emlig1->correction_factors[j].end();vec_ref++) |
---|
29 | { |
---|
30 | cout << j << " "; |
---|
31 | |
---|
32 | for(int i=0;i<(*vec_ref).size();i++) |
---|
33 | { |
---|
34 | cout << (*vec_ref)[i]; |
---|
35 | } |
---|
36 | |
---|
37 | cout << endl; |
---|
38 | } |
---|
39 | }*/ |
---|
40 | |
---|
41 | vec condition1 = "0.3 1.0 1.0"; |
---|
42 | emlig1->add_condition(condition1); |
---|
43 | |
---|
44 | vec condition2 = "-0.5 -1.0 1.0"; |
---|
45 | emlig1->add_condition(condition2); |
---|
46 | |
---|
47 | //vec condition3 = "-0.3 0.5 0.5"; |
---|
48 | //emlig1->add_condition(condition3); |
---|
49 | |
---|
50 | // emlig1->step_me(0); |
---|
51 | |
---|
52 | //emlig1->remove_condition(condition1); |
---|
53 | |
---|
54 | |
---|
55 | |
---|
56 | |
---|
57 | |
---|
58 | /* |
---|
59 | for(int i = 0;i<500;i++) |
---|
60 | { |
---|
61 | cout << endl << "Step:" << i << endl; |
---|
62 | |
---|
63 | |
---|
64 | |
---|
65 | double condition[emlig_size+1]; |
---|
66 | |
---|
67 | for(int k = 0;k<=emlig_size;k++) |
---|
68 | { |
---|
69 | condition[k] = (rand()-RAND_MAX/2)/1000.0; |
---|
70 | } |
---|
71 | |
---|
72 | |
---|
73 | vec* condition_vec = new vec(condition,emlig_size+1); |
---|
74 | emlig1->add_condition(*condition_vec); |
---|
75 | |
---|
76 | /* |
---|
77 | for(polyhedron* toprow_ref = emlig1->statistic.rows[emlig_size]; toprow_ref != emlig1->statistic.end_poly; toprow_ref = toprow_ref->next_poly) |
---|
78 | { |
---|
79 | cout << ((toprow*)toprow_ref)->probability << endl; |
---|
80 | } |
---|
81 | */ |
---|
82 | /* |
---|
83 | cout << emlig1->statistic_rowsize(emlig_size) << endl << endl; |
---|
84 | |
---|
85 | /* |
---|
86 | if(i-emlig1->number_of_parameters >= 0) |
---|
87 | { |
---|
88 | pause(30); |
---|
89 | } |
---|
90 | */ |
---|
91 | |
---|
92 | // emlig1->step_me(i); |
---|
93 | |
---|
94 | /* |
---|
95 | vector<int> sizevector; |
---|
96 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
97 | { |
---|
98 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
99 | } |
---|
100 | */ |
---|
101 | //} |
---|
102 | |
---|
103 | |
---|
104 | |
---|
105 | |
---|
106 | /* |
---|
107 | emlig1->step_me(1); |
---|
108 | |
---|
109 | vec condition = "2.0 0.0 1.0"; |
---|
110 | |
---|
111 | emlig1->add_condition(condition); |
---|
112 | |
---|
113 | vector<int> sizevector; |
---|
114 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
115 | { |
---|
116 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
117 | } |
---|
118 | |
---|
119 | emlig1->step_me(2); |
---|
120 | |
---|
121 | condition = "2.0 1.0 0.0"; |
---|
122 | |
---|
123 | emlig1->add_condition(condition); |
---|
124 | |
---|
125 | sizevector.clear(); |
---|
126 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
127 | { |
---|
128 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
129 | } |
---|
130 | */ |
---|
131 | |
---|
132 | return 0; |
---|
133 | } |
---|