Revision 1219, 1.2 kB
(checked in by sindj, 14 years ago)
|
Dodelani deleni - uz funguje. Je treba jeste dodelat predavani vrcholu a triangulaci a vypocet integralu v jednotlivych polyhedronech. Pak uz slucovani. JS
|
Rev | Line | |
---|
[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 | |
---|
| 15 | int main ( int argc, char* argv[] ) { |
---|
| 16 | |
---|
[1219] | 17 | emlig* emlig1 = new emlig(3); |
---|
[1186] | 18 | |
---|
[1219] | 19 | for(int i = 0;i<500;i++) |
---|
| 20 | { |
---|
| 21 | double condition[4]; |
---|
[1208] | 22 | |
---|
[1219] | 23 | condition[0] = rand()/1000.0; |
---|
| 24 | condition[1] = rand()/1000.0; |
---|
| 25 | condition[2] = rand()/1000.0; |
---|
| 26 | condition[3] = rand()/1000.0; |
---|
[1216] | 27 | |
---|
[1219] | 28 | vec* condition_vec = new vec(condition,4); |
---|
| 29 | emlig1->add_condition(*condition_vec); |
---|
| 30 | |
---|
| 31 | // emlig1->step_me(i); |
---|
| 32 | |
---|
| 33 | vector<int> sizevector; |
---|
| 34 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
| 35 | { |
---|
| 36 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
| 37 | } |
---|
| 38 | } |
---|
| 39 | |
---|
| 40 | |
---|
| 41 | |
---|
| 42 | /* |
---|
| 43 | emlig1->step_me(1); |
---|
| 44 | |
---|
| 45 | vec condition = "2.0 0.0 1.0"; |
---|
| 46 | |
---|
[1208] | 47 | emlig1->add_condition(condition); |
---|
| 48 | |
---|
[1216] | 49 | vector<int> sizevector; |
---|
| 50 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
| 51 | { |
---|
| 52 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
| 53 | } |
---|
| 54 | |
---|
[1219] | 55 | emlig1->step_me(2); |
---|
[1216] | 56 | |
---|
[1219] | 57 | condition = "2.0 1.0 0.0"; |
---|
[1216] | 58 | |
---|
| 59 | emlig1->add_condition(condition); |
---|
| 60 | |
---|
| 61 | sizevector.clear(); |
---|
| 62 | for(int s = 0;s<=emlig1->number_of_parameters;s++) |
---|
| 63 | { |
---|
| 64 | sizevector.push_back(emlig1->statistic_rowsize(s)); |
---|
| 65 | } |
---|
[1219] | 66 | */ |
---|
[1216] | 67 | |
---|
[976] | 68 | return 0; |
---|
| 69 | } |
---|