|
Revision 1219, 1.2 kB
(checked in by sindj, 15 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 | } |
|---|