OpenMP:为“shared”预先定义的“shared”？

Question

查看此函数(矩阵-向量积)：

std::vector<double> times(std::vector<std::vector<double> > const& A, std::vector<double> const& b, int m, int n) {

    std::vector<double> c;
    c.resize(n);

    int i, j;
    double sum;

    #pragma omp parallel for default(none) private(i, j, sum) shared(m, n, A, b, c)
    for (i = 0; i < m; ++i) {
        sum = 0.0;
        for (j = 0; j < n; j++) {
            sum += A[i][j] * b[j];
        }
        c[i] = sum;
    }

    return c;
}

当尝试使用OpenMP编译此代码时，编译器会失败，并显示以下信息：

Invoking: GCC C++ Compiler
g++ -O0 -g3 -Wall -c -fmessage-length=0 -fopenmp -MMD -MP -MF"src/OpemMPTutorial.d" -MT"src/OpemMPTutorial.d" -o "src/OpemMPTutorial.o" "../src/OpemMPTutorial.cpp"
../src/OpemMPTutorial.cpp:127: warning: ignoring #pragma omp end
../src/OpemMPTutorial.cpp: In function 'std::vector<double, std::allocator<double> > times(std::vector<std::vector<double, std::allocator<double> >, std::allocator<std::vector<double, std::allocator<double> > > >&, std::vector<double, std::allocator<double> >&, int, int)':
../src/OpemMPTutorial.cpp:200: error: 'b' is predetermined 'shared' for 'shared'
../src/OpemMPTutorial.cpp:200: error: 'A' is predetermined 'shared' for 'shared'
make: *** [src/OpemMPTutorial.o] Error 1

这是怎么回事？

(请注意，简单地删除const会导致相同的错误。)

Answer 1

这是由于gcc-4.2中对OpenMP的支持不足造成的。使用gcc-4.7可以毫无问题地编译代码片段。

Answer 2

我有一个非常类似的问题，在我从OpenMP指令的shared部分删除了共享的const变量后，我体验到这样的程序可以用苹果的GCC 4.2编译。它们被预先确定为共享的，因为它们是常量，并且不需要为每个线程制作副本。而且编译器似乎不接受在知道已经…的情况下显式地告诉它

我还会删除default(none)规范(但请参阅comment below)。OpenMP旨在减少显式规范，所以让它来做它的工作吧。