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

Не понимайте, где я понимаю это превратно. Если скомпилировано без поддержки 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 компилятора.