У меня есть цикл, который должен быть хорошо распараллелен, вставив одну прагму openmp:

  boost::normal_distribution<double> ddist(0, pow(retention, i - 1));
  boost::variate_generator<gen &, BOOST_TYPEOF(ddist)> dgen(rng, ddist);
  // Diamond                                                                
  const std::uint_fast32_t dno = 1 << i - 1;
// #pragma omp parallel for
  for (std::uint_fast32_t x = 0; x < dno; x++)
    for (std::uint_fast32_t y = 0; y < dno; y++)
      {
        const std::uint_fast32_t diff = size/dno;
        const std::uint_fast32_t x1 = x*diff, x2 = (x + 1)*diff;
        const std::uint_fast32_t y1 = y*diff, y2 = (y + 1)*diff;
        double avg =
          (arr[x1][y1] + arr[x1][y2] + arr[x2][y1] + arr[x2][y2])/4;
        arr[(x1 + x2)/2][(y1 + y2)/2] = avg + dgen();
      }

(если я не сделаю ошибку, каждое выполнение не зависит вообще от других. Извините, что не весь код вставлен).

Тем не менее, мой вопрос - является ли поток RNG потокобезопасным? Кажется, они ссылаются на код gcc для gcc, поэтому, даже если код gcc является поточно-ориентированным, это может быть не так для других платформ.