root/applications/robust/main.cpp @ 1267

/*!
\file
\brief Robust
\author Vasek Smidl 

 */

#include "robustlib.h"
#include <vector>
using namespace itpp;

using namespace bdm;

int main ( int argc, char* argv[] ) {
        
        emlig* emlig1 = new emlig(4);

        /*
        emlig1->set_correction_factors(4);

        for(int j = 0;j<emlig1->correction_factors.size();j++)
        {
                for(set<my_ivec>::iterator vec_ref = emlig1->correction_factors[j].begin();vec_ref!=emlig1->correction_factors[j].end();vec_ref++)
                {
                        for(int i=0;i<(*vec_ref).size();i++)
                        {
                                cout << (*vec_ref)[i];
                        }

                        cout << endl;
                }
        }
        */

        //emlig1->step_me(0);

        for(int i = 0;i<500;i++)
        {
                cout << "Step:" << i << endl;

                double condition[4];            

                condition[0] = rand()/1000.0;
                condition[1] = rand()/1000.0;
                condition[2] = rand()/1000.0;
                condition[3] = rand()/1000.0;
                /*condition[4] = rand()/1000.0;
                condition[5] = rand()/1000.0;
                condition[6] = rand()/1000.0;
                condition[7] = rand()/1000.0;
                condition[8] = rand()/1000.0;*/
                

                vec* condition_vec = new vec(condition,4);
                emlig1->add_condition(*condition_vec);

                emlig1->step_me(i);
                
                vector<int> sizevector;
                for(int s = 0;s<=emlig1->number_of_parameters;s++)
                {
                        sizevector.push_back(emlig1->statistic_rowsize(s));
                }
        }


        
        /*
        emlig1->step_me(1);

        vec condition = "2.0 0.0 1.0";  

        emlig1->add_condition(condition);

        vector<int> sizevector;
        for(int s = 0;s<=emlig1->number_of_parameters;s++)
        {
                sizevector.push_back(emlig1->statistic_rowsize(s));
        }

        emlig1->step_me(2);

        condition = "2.0 1.0 0.0";

        emlig1->add_condition(condition);

        sizevector.clear();
        for(int s = 0;s<=emlig1->number_of_parameters;s++)
        {
                sizevector.push_back(emlig1->statistic_rowsize(s));
        }
        */

        return 0;
}