root/applications/robust/main.cpp @ 1219

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

Line 
1
2/*!
3\file
4\brief Robust
5\author Vasek Smidl
6
7 */
8
9#include "robustlib.h"
10#include <vector>
11using namespace itpp;
12
13using namespace bdm;
14
15int main ( int argc, char* argv[] ) {
16       
17        emlig* emlig1 = new emlig(3);   
18
19        for(int i = 0;i<500;i++)
20        {
21                double condition[4];           
22
23                condition[0] = rand()/1000.0;
24                condition[1] = rand()/1000.0;
25                condition[2] = rand()/1000.0;
26                condition[3] = rand()/1000.0;
27
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
47        emlig1->add_condition(condition);
48
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
55        emlig1->step_me(2);
56
57        condition = "2.0 1.0 0.0";
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        }
66        */
67
68        return 0;
69}
Note: See TracBrowser for help on using the browser.