unsigned num_threads = std::thread::hardware_concurrency();

if(argc > 1) {

num_threads = std::atoi(argv[1]);

}

int rounds {5};

std::cout << std::setw(12) << "# blocks"

<< std::setw(12) << "OpenMP"

<< std::setw(12) << "TBB"

<< std::setw(12) << "Taskflow"

<< std::setw(12) << "speedup1"

<< std::setw(12) << "speedup2"

<< '\n';

for(int S=1; S<=4096; S += 128) {

M = N = S;

B = 8;

MB = (M/B) + (M%B>0);

NB = (N/B) + (N%B>0);

double omp_time {0.0};

double tbb_time {0.0};

double tf_time {0.0};

init_matrix();

for(int j=0; j<rounds; ++j) {

omp_time += measure_time_omp(num_threads).count();

tbb_time += measure_time_tbb(num_threads).count();

tf_time += measure_time_taskflow(num_threads).count();

}

destroy_matrix();

std::cout << std::setw(12) << MB*NB

<< std::setw(12) << omp_time / rounds / 1e3

<< std::setw(12) << tbb_time / rounds / 1e3

<< std::setw(12) << tf_time / rounds / 1e3

<< std::setw(12) << omp_time / tf_time

<< std::setw(12) << tbb_time / tf_time

<< std::endl;

}

return 0;