multithreading - Clang を使用した OpenMP の C++11 代替

Question

Clang は (まだ) OpenMP をサポートしていませんが、C++11 で「並列 for」を実装することは可能ですか?

Answer 1

OpenMP バージョン:

// parallelfor_gcc.cpp
// g++ -O2 -Wall -std=c++11 -fopenmp parallelfor_gcc.cpp
#include <cmath>
#include <vector>
int main() {
  unsigned int size = 1e8;
  std::vector<double> vect(size);
#pragma omp parallel for
  for (unsigned int i=0; i<size; i++) {
    vect[i] = sin(2*M_PI*i/(double)size);
  }
  return 0;
}

C++11 バージョン:

// parallelfor_clang.cpp
// clang++ -O4 -Wall -std=c++11 -lpthread parallelfor_clang.cpp
#include <cmath>
#include <thread>
#include <vector>
void parallelFor(const unsigned int size, 
                 std::function<void(const unsigned int)> func) {
  const unsigned int nbThreads = std::thread::hardware_concurrency();
  std::vector < std::thread > threads;
  for (unsigned int idThread = 0; idThread < nbThreads; idThread++) {
    auto threadFunc = [=, &threads]() {
      for (unsigned int i=idThread; i<size; i+=nbThreads) {
        func(i);
      }
    };
    threads.push_back(std::thread(threadFunc));
  }
  for (auto & t : threads) t.join();
}
int main() {
  unsigned int size = 1e8;
  std::vector<double> vect(size);
  auto myFunc = [=, &vect](unsigned int i){
    vect[i] = sin(2*M_PI*i/(double)size);
  };
  parallelFor(size, myFunc);
  return 0;
}

OpenMP 句 (firstprivate...) も同じ方法で実装できますが、(少し) 手間がかかります...