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