#include "../performance/clock.h"
#include <vector>
#include <boost/random.hpp>


//#include rand
using namespace std;

double const D2I_CONVERSION_FACTOR = 1.5 * 4503599627370496.0;

inline int d2i(double val)
{
    union {
        int i[2];
        double d;
    } dlong;
    dlong.d = val + D2I_CONVERSION_FACTOR;
    return dlong.i[0];
}

int main()
{
	vector <double> data;
	int const nbData = 100;
	int const nbIter = 10000000;
	boost::mt19937 generatorEngine;
	boost::uniform_real<double> generator(std::numeric_limits<int>::min()+10, std::numeric_limits<int>::max()-10);
//	boost::uniform_real<double> generator(0, std::numeric_limits<int>::max()-10);

	cout << generator.min() << " -> " << generator.max() << endl;
	for(int i = 0 ; i<nbData ; ++i)
	{
		data.push_back(generator(generatorEngine));
	}
	std::vector<int> basicConversion(nbData);
	std::vector<int> convolutedConversion(nbData);
	std::vector<int>::iterator itInt;
	std::vector<double>::iterator itDouble;
	std::vector<double>::iterator const endIt = data.end();

	Clock clk;
	clk.reset();
	for (int i=0 ; i<nbIter ; ++i)
	{
		itInt = basicConversion.begin();
		itDouble = data.begin();
		while (itDouble != endIt)
		{
			*(itInt++) = *(itDouble++);
		}
	}
	double basicTime = clk.run();
	cout << "Basic: " << basicTime << endl;

	clk.reset();
	for (int i=0 ; i<nbIter ; ++i)
	{
		itInt = convolutedConversion.begin();
		itDouble = data.begin();
		while (itDouble != endIt)
		{
			//*(itInt++) = d2i(*(itDouble++)-0.5);
			// Décommenter pour nombres négatifs
			*itInt = d2i(*itDouble + ((*itDouble > 0) ? -0.5 : 0.5));
			++itInt;
			++itDouble;
		}
	}
	double convolutedTime = clk.run();
	cout << "Convoluted: " << convolutedTime << endl;

	itDouble = data.begin();
	itInt = basicConversion.begin();
	std::vector<int>::iterator itInt2 = convolutedConversion.begin();
	int nbError = 0;
	while (itInt != basicConversion.end())
	{
		if (*itInt != *itInt2)
		{
			cout << "Error: " << *itDouble << " -> " << *itInt << " ; " << *itInt2 << endl;
			nbError++;
		}
		++itInt;
		++itDouble;
		++itInt2;
	}
	cout << "Rapport :\nTaux d'erreur : " << double(nbError)/double(nbData) << 
		"\nConvoluted/Basic : " << convolutedTime/basicTime << "\nBasic/Convoluted : "
		<< basicTime/convolutedTime << endl;
}