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

RevLine 
[976]1
2/*!
3\file
4\brief Robust
5\author Vasek Smidl
6
7 */
8
9#include "robustlib.h"
[1216]10#include <vector>
[1208]11using namespace itpp;
[976]12
[1213]13using namespace bdm;
[976]14
15int 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}
Note: See TracBrowser for help on using the browser.