/* 
Compilation command :
g++ Vavilov03.cpp -I `root-config --incdir` `root-config --libs` -lMathMore -fopenmp

To run it with a single thread, set
export OMP_NUM_THREADS=1
To run with multiple threads, set
export OMP_NUM_THREADS=4
*/

#include "iostream"
#include "Math/VavilovAccurate.h"
#include "TRandom.h"
#include <omp.h>
#include "time.h"
using namespace std;

int main()
{

	TRandom r(0);
	double ran;
	double a;
	cout<<"Start "<<clock()/CLOCKS_PER_SEC<<endl;
	#pragma omp parallel for private(ran,a)
	for(int i=0;i<10000000;i++){
		double x = 5;
		ran = r.Uniform();
	
		//#pragma omp critical 
		a = 1-ROOT::Math::vavilov_accurate_cdf((-1-x)/3,1,ran)-ROOT::Math::vavilov_accurate_cdf_c((1-x)/3,1,ran);
		//cout<<a<<endl;
		//cout<<"this is thread #"<<omp_get_thread_num()<<endl;

	}
	cout<<"Stop "<<clock()/CLOCKS_PER_SEC<<endl;
	return 0;

}