#include <time.h>
#include <intrin.h>
#include <omp.h>
#include <iostream>

#define _CRT_SECURE_NO_WARNINGS

using namespace std;
#define SIZE_DARR 100000000
double arr[SIZE_DARR];

__int64 hz_cpu() { // замер частоты
	clock_t t_clock;
	__int64 t_tsc;
	t_clock = clock() + CLOCKS_PER_SEC;
	t_tsc = __rdtsc(); // взять TSC
	while (clock() < t_clock); // отсчет одной секунды
	return (__rdtsc() - t_tsc); // частота в герцах
}

int mhz_cpu() {
	clock_t clock_tick1, clock_tick2;
	__int64 cpu_tick1, cpu_tick2;
	int rdtsc_tick;
	double usec; // время в микросекундах
	clock_tick1 = clock();
	while ((clock_tick2 = clock()) == clock_tick1); // пропуск остатка текущего тика clock
	cpu_tick1 = __rdtsc(); // взять TSC
	while (clock() == clock_tick2); // отсчет одного тика clock
	cpu_tick2 = __rdtsc() - cpu_tick1; // сколько натикал счетчик TSC за один тик clock
	// вычисляем частоту в мегагерцах
	usec = 1000000.0 / CLOCKS_PER_SEC; // время одного тика clock в микросекундах
	return int(cpu_tick2 / usec);
}

void my_func(__int64 count) { // заполняет массив синусами
	for (__int64 i = 0; i < count; i++)
		arr[i] = asin((double)((i / count) % 100007));
}

int main() {
	setlocale(LC_CTYPE, "rus"); // кириллизация
	double a, b, result = 0, h;
	__int64 t_start, time = 0, freq, minTime = 100000000000000;

	cout << "Введите a: "; cin >> a;
	cout << "\nВведите b: "; cin >> b;

	my_func(SIZE_DARR);
	//h = (b - a) / n;
	omp_set_dynamic(0);
	cout << "proc = " << omp_get_num_procs() << "\n";
	cout << "threads = " << omp_get_num_threads() << "\n";

	for (int flow = 1; flow <= 8; flow++) {// количество потоков
		cout << "\n\nflow = " << flow << "\n_______________" << "\n_______________";

		for (int n = 100; n <= 1000000; n *= 10) { //гранулярность задачи 10-1000000
			h = (b - a) / n;

			for (int z = 0; z < 5; z++) {
				//t_start = __rdtsc();
#				pragma omp parallel reduction(+:result) num_threads(flow)
				{
					t_start = __rdtsc();
					for (int i = 1; i < n; i++) {
						result += asin((a + i * h) - h / 2);
					}
				}

				time = __rdtsc() - t_start;
				freq = 1000000.0 * mhz_cpu();
				time = 1000000000.0 * time / freq;
				if (time < minTime)
					minTime = time;
				result *= h;
			}

			cout << "\ncountN = " << n;
			cout << "\ntime = " << minTime << " мс.\n____________\n";
			minTime = 100000000000000;
		}

	}
	system("pause");
	return 0;
}