看到这个函数(矩阵,向量积):OpenMP:'共享'的预定'共享'?
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
结果在相同的错误。)
我编译你的函数没有问题,使用'g ++ 4.6.3'。 – Massimiliano
我正在使用'i686-apple-darwin11-llvm-gcc-4.2'。这可能是这个编译器版本的问题。我正在尝试升级到gcc 4.7 – clstaudt
@ cls:我在开玩笑。我不使用OpenMP,所以从未见过这个GCC错误。所以你得到了一个赞成,只是为了给我一个新的GCC错误(你并不是真的)。 –