main.cpp 2.19 KB
Newer Older
1
#include "pls/pls.h"
2 3 4

#include "benchmark_runner.h"
#include "benchmark_base/fft.h"
5

6
using namespace comparison_benchmarks::base;
7

8
void pls_conquer(fft::complex_vector::iterator data, fft::complex_vector::iterator swap_array, int n) {
9
  if (n < 2) {
10
    return;
11 12
  }

13
  fft::divide(data, swap_array, n);
14
  if (n <= fft::RECURSIVE_CUTOFF) {
15 16
    fft::conquer(data, swap_array, n / 2);
    fft::conquer(data + n / 2, swap_array + n / 2, n / 2);
17
  } else {
18
    pls::spawn([data, n, swap_array]() {
19
      pls_conquer(data, swap_array, n / 2);
20
    });
21
    pls::spawn([data, n, swap_array]() {
22
      pls_conquer(data + n / 2, swap_array + n / 2, n / 2);
23
    });
24
    pls::sync();
25
  }
26 27

  fft::combine(data, n);
28 29
}

30 31
constexpr int MAX_NUM_TASKS = 16;
constexpr int MAX_STACK_SIZE = 4096 * 1;
32

33
int main(int argc, char **argv) {
34
  auto settings = benchmark_runner::parse_parameters(argc, argv);
35

36 37 38 39 40 41
  fft::complex_vector data(settings.size_);
  fft::complex_vector swap_array(settings.size_);
  fft::fill_input(data);

  string test_name = to_string(settings.num_threads_) + ".csv";
  string full_directory = settings.output_directory_ + "/PLS_v3/";
42
  benchmark_runner runner{full_directory, test_name};
43

44
  pls::scheduler scheduler{(unsigned) settings.num_threads_, MAX_NUM_TASKS, MAX_STACK_SIZE};
45

46 47 48 49
  if (settings.type_ == benchmark_runner::benchmark_settings::ISOLATED) {
    printf("Running isolated measurement...\n");
    runner.enable_memory_stats();
    runner.pre_allocate_stats();
50

51 52 53 54 55 56
    runner.run_iterations(settings.iterations_, [&]() {
      scheduler.perform_work([&]() {
        pls_conquer(data.begin(), swap_array.begin(), settings.size_);;
      });
    }, [&]() {
      fft::fill_input(data); // Reset data before each run
57
    });
58 59 60 61 62
    runner.commit_results(true);
  } else {
    printf("Running periodic measurement...\n");
    runner.enable_wall_time_stats();
    runner.pre_allocate_stats();
63

64 65 66 67 68 69 70 71
    runner.run_periodic(settings.iterations_, settings.interval_period_, settings.interval_deadline_, [&]() {
      scheduler.perform_work([&]() {
        pls_conquer(data.begin(), swap_array.begin(), settings.size_);;
      });
      fft::fill_input(data); // Reset data before each run
    });
    runner.commit_results(true);
  }
72

73
  return 0;
74
}