the-mastercoder
diff --git a/‎benchmarks/parallel_pipeline/main.cpp‎
Lines changed: 78 additions & 0 deletions b/‎benchmarks/parallel_pipeline/main.cpp‎
Lines changed: 78 additions & 0 deletions
diff --git a/‎benchmarks/parallel_pipeline/omp.cpp‎
Lines changed: 49 additions & 0 deletions b/‎benchmarks/parallel_pipeline/omp.cpp‎
Lines changed: 49 additions & 0 deletions
diff --git a/‎benchmarks/parallel_pipeline/parallel_pipeline.hpp‎
Lines changed: 8 additions & 0 deletions b/‎benchmarks/parallel_pipeline/parallel_pipeline.hpp‎
Lines changed: 8 additions & 0 deletions
diff --git a/‎benchmarks/parallel_pipeline/seq.cpp‎
Lines changed: 16 additions & 0 deletions b/‎benchmarks/parallel_pipeline/seq.cpp‎
Lines changed: 16 additions & 0 deletions
@@ -0,0 +1,78 @@
+#include "parallel_pipeline.hpp"
+#include <CLI11.hpp> 
+
+int main(int argc, char* argv[]) {
+
+  CLI::App app{"Parallel Pipeline"};
+
+  unsigned num_threads {1}; 
+  app.add_option("-t,--num_threads", num_threads, "number of threads (default=1)");
+
+  unsigned num_rounds {1};  
+  app.add_option("-r,--num_rounds", num_rounds, "number of rounds (default=1)");
+
+  std::string model = "tf";
+  app.add_option("-m,--model", model, "model name tbb|omp|tf (default=tf)")
+     ->check([] (const std::string& m) {
+        if(m != "tbb" && m != "omp" && m != "tf") {
+          return "model name should be \"tbb\", \"omp\", or \"tf\"";
+        }
+        return "";
+     });
+
+  unsigned num_lines {8};
+  app.add_option("-l,--num_lines", num_lines, "num of lines (default=8)");
+
+  std::string pipes = "ssssssss";
+  app.add_option("-p,--pipes", pipes, "the chain of pipes (default=ssssssss)")
+      ->check([pipes] (const std::string& p) {
+        if (p[0] == 'p') {
+          return "the first pipe should be \"s\" (serial)";
+        }
+        else if (pipes.size() > 10) {
+          return "no more than 10 pipes";
+        }
+        else if (pipes.size() == 0) {
+          return "at least one pipe is required";
+        }
+        return "";
+      });
+
+  CLI11_PARSE(app, argc, argv);
+
+  std::cout << "model="       << model       << ' '
+            << "num_threads=" << num_threads << ' '
+            << "num_rounds="  << num_rounds  << ' '
+            << "num_lines="   << num_lines   << ' '
+            << "pipes="       << pipes       << ' '
+            << std::endl;
+
+  std::cout << std::setw(12) << "size"
+            << std::setw(12) << "Runtime"
+            << '\n';
+
+  size_t  log_length = 23;
+  
+  for(size_t i = 1; i <= log_length; ++i) {
+
+    size_t L = 1 << i;
+    
+    double runtime {0.0};
+    
+    for(unsigned j = 0; j < num_rounds; ++j) {
+      if(model == "tf") {
+        runtime += measure_time_taskflow(pipes, num_lines, num_threads, L).count();
+      }
+      else if(model == "tbb") {
+        runtime += measure_time_tbb(pipes, num_lines, num_threads, L).count();
+      }
+      else if(model == "omp") {
+        runtime += measure_time_omp(pipes, num_lines, num_threads, L).count();
+      }
+    }
+
+    std::cout << std::setw(12) << L 
+              << std::setw(12) << runtime / num_rounds / 1e3
+              << std::endl;
+  }
+}
@@ -0,0 +1,49 @@
+#include "parallel_pipeline.hpp"
+#include <omp.h>
+
+// linear_chain_omp
+void linear_chain_omp(size_t length, unsigned num_threads) {
+
+  size_t counter = 0;
+  size_t *D = new size_t [length]; 
+
+  #pragma omp parallel num_threads(num_threads)
+  {
+    #pragma omp single
+    {
+      for(size_t i=0; i<length; ++i) {
+
+        if(i==0) {
+          #pragma omp task firstprivate(i) depend(out:D[i])
+          {
+            ++counter;
+          }
+        }
+        else {
+          #pragma omp task firstprivate(i) depend(out:D[i]) depend(in:D[i-1])
+          {
+            ++counter;
+          }
+        }
+      }
+    }
+  }
+
+  assert(counter == length);
+
+  delete [] D;
+}
+
+std::chrono::microseconds measure_time_omp(
+  std::string pipes,
+  unsigned length,
+  unsigned num_threads,
+  size_t size
+) {
+  auto beg = std::chrono::high_resolution_clock::now();
+  //linear_chain_omp(length, num_threads);
+  auto end = std::chrono::high_resolution_clock::now();
+  return std::chrono::duration_cast<std::chrono::microseconds>(end - beg);
+}
+
+
@@ -0,0 +1,8 @@
+#include <chrono>
+#include <string>
+#include <cassert>
+
+std::chrono::microseconds measure_time_taskflow(std::string, unsigned, unsigned, size_t);
+std::chrono::microseconds measure_time_tbb(std::string, unsigned, unsigned, size_t);
+std::chrono::microseconds measure_time_omp(std::string, unsigned, unsigned, size_t);
+
@@ -0,0 +1,16 @@
+#include <iostream>
+#include <chrono>
+
+#include "levelgraph.hpp"
+
+std::chrono::microseconds measure_time_seq(LevelGraph& graph){
+  auto beg = std::chrono::high_resolution_clock::now();
+  for(size_t l=0; l<graph.level(); l++){ 
+    for(size_t i=0; i<graph.length(); i++){
+      Node& n = graph.node_at(l, i);
+      n.mark();
+    }
+  }
+  auto end = std::chrono::high_resolution_clock::now();
+  return std::chrono::duration_cast<std::chrono::microseconds>(end - beg);
+}