#include <iostream>
#include <omp.h>
#include <time.h>
#include <chrono>
#include <iomanip> 

int main()
{
	double accuracy;
	long numberOfLoops = 10;
	std::cout << "Set accuracy: \n";
	std::cin >> accuracy;
	std::cout << " \n";
	double pi = 0;
	long minTime = INFINITY;
	long maxTime = 0;
	long sq_time = 0;
	long mt_time = 0;

	auto start = std::chrono::high_resolution_clock::now();
	auto finish = std::chrono::high_resolution_clock::now();


	//Sequenial
	std::cout << "Sequential:\n";
	for(long k = 0; k<numberOfLoops; k++)
	{	
		pi = 0;
		start = std::chrono::high_resolution_clock::now();
		////////////////////////
		for (double i = 1; i <= accuracy; i++)
		{
			pi += (1 / (i*i));
		}
		////////////////////////
		finish = std::chrono::high_resolution_clock::now();
		long currentTime = std::chrono::duration_cast<std::chrono::nanoseconds>(finish - start).count();
		sq_time += currentTime;
		if (currentTime > maxTime) maxTime = currentTime;
		if (currentTime < minTime) minTime = currentTime;
	}

	//std::cout << "PI number =	" << sqrt(pi * 6) << "\n";

	std::cout << "Avg time  = " << sq_time/numberOfLoops << "ns\n";
	std::cout << "Pi  =	" << std::setprecision(30) << sqrt(pi*6) << "ns\n\n";
	pi = 0;


	minTime = INFINITY;
	maxTime = 0;

	//Parallel
	std::cout << "Parallel:\n";
	for (long k = 0; k<numberOfLoops; k++)
	{
		pi = 0;
		start = std::chrono::high_resolution_clock::now();
		////////////////////////
#pragma omp parallel for
		for (double i = 1; i <= accuracy; i++)
		{
			pi += (1 / (i*i));
		}
		////////////////////////


		finish = std::chrono::high_resolution_clock::now();
		long currentTime = std::chrono::duration_cast<std::chrono::nanoseconds>(finish - start).count();
		mt_time += currentTime;
		if (currentTime > maxTime) maxTime = currentTime;
		if (currentTime < minTime) minTime = currentTime;
	}
	//std::cout << "PI number =	" << sqrt(pi * 6) << "\n";

	std::cout << "Avg time  = " << mt_time / numberOfLoops << "ns\n";
	std::cout << "Pi  =	" << std::setprecision(30) << sqrt(pi * 6) << "ns\n\n";
	pi = 0;
	

	std::cin.ignore(2);
	return 0;
}