[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 | |
---|
[976] | 17 | int main ( int argc, char* argv[] ) { |
---|
| 18 | |
---|
[1186] | 19 | |
---|
[1268] | 20 | emlig* emlig1 = new emlig(emlig_size); |
---|
| 21 | |
---|
[1267] | 22 | /* |
---|
| 23 | emlig1->set_correction_factors(4); |
---|
[1266] | 24 | |
---|
[1267] | 25 | for(int j = 0;j<emlig1->correction_factors.size();j++) |
---|
| 26 | { |
---|
| 27 | for(set<my_ivec>::iterator vec_ref = emlig1->correction_factors[j].begin();vec_ref!=emlig1->correction_factors[j].end();vec_ref++) |
---|
| 28 | { |
---|
[1268] | 29 | cout << j << " "; |
---|
| 30 | |
---|
[1267] | 31 | for(int i=0;i<(*vec_ref).size();i++) |
---|
| 32 | { |
---|
| 33 | cout << (*vec_ref)[i]; |
---|
| 34 | } |
---|
| 35 | |
---|
| 36 | cout << endl; |
---|
| 37 | } |
---|
[1268] | 38 | }*/ |
---|
| 39 | |
---|
[1267] | 40 | |
---|
[1234] | 41 | //emlig1->step_me(0); |
---|
| 42 | |
---|
[1219] | 43 | for(int i = 0;i<500;i++) |
---|
| 44 | { |
---|
[1271] | 45 | cout << endl << "Step:" << i << endl; |
---|
| 46 | |
---|
| 47 | |
---|
[1208] | 48 | |
---|
[1268] | 49 | double condition[emlig_size+1]; |
---|
[1220] | 50 | |
---|
[1268] | 51 | for(int k = 0;k<=emlig_size;k++) |
---|
| 52 | { |
---|
| 53 | condition[k] = rand()/1000.0; |
---|
| 54 | } |
---|
| 55 | |
---|
[1216] | 56 | |
---|
[1268] | 57 | vec* condition_vec = new vec(condition,emlig_size+1); |
---|
[1219] | 58 | emlig1->add_condition(*condition_vec); |
---|
[1271] | 59 | |
---|
| 60 | cout << emlig1->statistic_rowsize(emlig_size) << endl << endl; |
---|
[1268] | 61 | |
---|
[1271] | 62 | if(i-emlig1->number_of_parameters >= 0) |
---|
| 63 | { |
---|
| 64 | pause(30); |
---|
| 65 | } |
---|
[1219] | 66 | |
---|
[1271] | 67 | // emlig1->step_me(i); |
---|
[1219] | 68 | |
---|
| 69 | vector<int> sizevector; |
---|
| 70 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
| 71 | { |
---|
| 72 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
| 73 | } |
---|
| 74 | } |
---|
| 75 | |
---|
| 76 | |
---|
| 77 | |
---|
| 78 | /* |
---|
| 79 | emlig1->step_me(1); |
---|
| 80 | |
---|
| 81 | vec condition = "2.0 0.0 1.0"; |
---|
| 82 | |
---|
[1208] | 83 | emlig1->add_condition(condition); |
---|
| 84 | |
---|
[1216] | 85 | vector<int> sizevector; |
---|
| 86 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
| 87 | { |
---|
| 88 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
| 89 | } |
---|
| 90 | |
---|
[1219] | 91 | emlig1->step_me(2); |
---|
[1216] | 92 | |
---|
[1219] | 93 | condition = "2.0 1.0 0.0"; |
---|
[1216] | 94 | |
---|
| 95 | emlig1->add_condition(condition); |
---|
| 96 | |
---|
| 97 | sizevector.clear(); |
---|
| 98 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
| 99 | { |
---|
| 100 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
| 101 | } |
---|
[1219] | 102 | */ |
---|
[1216] | 103 | |
---|
[976] | 104 | return 0; |
---|
| 105 | } |
---|