C ++ OpenMP переменная видимость в параллельных задачах

Question

Не понимаю, где я могу ошибаться.Если скомпилировано без поддержки openmp, код работает правильно.Но с переменными openmp, похоже, неправильная видимость.

У меня было следующее намерение.Каждый поток имеет свой собственный max_private, в котором он находит локальный максимум.Затем в критическом разделе обнаруживается глобальный максимум.

#include <iostream>
#include <vector>

typedef std::vector<long> Vector;

long max(const Vector& a, const Vector& b)
{
    long max = 0;
    #pragma omp parallel
    {
        long max_private = 0;

        #pragma omp single
        {
            for (   Vector::const_iterator a_it = a.begin();
                    a_it != a.end();
                    ++a_it)
            {
                #pragma omp task
                {
                    for (   Vector::const_iterator b_it = b.begin();
                            b_it != b.end();
                            ++b_it)
                    {
                        if (*a_it + *b_it > max_private) {
                            max_private = *a_it + *b_it;
                        }
                    }
                }
            }
        }

        #pragma omp critical
        {
            std::cout << max_private << std::endl;
            if (max_private > max) {
                max = max_private;
            }
        }
    }
    return max;
}

int main(int argc, char* argv[])
{
    Vector a(100000);
    Vector b(10000);
    for (long i = 0; i < a.size(); ++i) {
        a[i] = i;
    }
    for (long i = 0; i < b.size(); ++i) {
        b[i] = i * i;
    }

    std::cout << max(a, b) << std::endl;

    return 0;
}

Я не хочу использовать параллель для, поскольку в последнем случае я собираюсь использовать структуры данных, которые не поддерживают итераторы произвольного доступа.

Я использовал компилятор g ++ - 4.4.

Answer 1

Получил подробный ответ на форуме OpenMP. http://openmp.org/forum/viewtopic.php?f=3&t=912&start=0

Пришлось сделать max_private threadprivate.

Answer 2

Вам нужно уменьшить openmp, чтобы найти максимум. Это доступно в OpenMP 3.0

сделать сокращение на

long max(const Vector& a, const Vector& b)
{
    long max_val = 0;  // <- change the name for clarity
    #pragma omp parallel reduction(max : max_val)
    { 
      [...] 
      #pragma omp critical
      {
        std::cout << max_private << std::endl;
        if (max_private > max_val) {
            max_val = max_private;
        }
      }
    }
    return max_val; 
}

max_val является приватным в каждом потоке, и в конце цикла происходит уменьшение, которое возвращает максимум каждого потока open_mp.

Answer 3

Я думаю, что вам нужны некоторые атрибуты переменных. По умолчанию все переменные являются общими, и, так как вы хотите, чтобы max_private был для каждого потока, вы можете с радостью добавить предложение «private» в его объявление (тогда оно может пройти немного быстрее).

Для максимального вам нужно поместить его в блок omp и пометить его как «общий», если он все еще не работает должным образом (хотя он должен быть в порядке, если он у вас есть)

Итак, я думаю, что вы хотите #pragma omp parallel private(max_private) и поместить обе переменные вне основного блока omp.