/*!
\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;
}