openmp/C++ 简单并行区域返回不一致的结果答案

【问题标题】：openmp/C++ simple parallel region returns incosistent resultsopenmp/C++ 简单并行区域返回不一致的结果
【发布时间】：2016-06-04 20:29:40
【问题描述】：

所以我正在尝试学习 OpenMP API，我已经掌握了基础知识，但这部分让我感到困惑：（可测试的代码！）

#include <iostream>
#include <cstdlib>
#include <string>
#include <omp.h>
#include <unistd.h>

using namespace std;
const int col = 10, row = 10;
int c[][11] = {{-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
               {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1}  };

int main(int argc, char** argv) 
{
int temp[3] = {-1, -1, -1};
int id;
for (unsigned short i = 0; i < 10; i++){
    #pragma omp parallel shared(c) firstprivate(i) private(temp)  
    {
        #pragma ivdep
        #pragma omp for schedule(static, 1) //private(id)//,m_i, ind)
        for(unsigned short j = 0; j < 10; j++){
            // calculate c
            /*  old
            temp[0] = c[i-1][j-1]+3;
            temp[1] = c[i-1][j]-4;
            temp[2] = c[i][j-1]-5;
            c[i][j] = temp[0]; 
            c[i][j] = temp[1];
            c[i][j] = temp[2];
            */
            temp[0] = c[i][j-1]+3;
            c[i][j] = temp[0]; 
        }
    }
}
// print c
for (int i = 0; i< row + 1; i++)
{
    for (int j = 0; j< col + 1; j++)
    {
        cout << c[i][j] << "/";
    }
    cout << endl;
}
 } // end main

我相信代码是自我解释的，我只想补充一点，我只想并行化内部循环，因为我希望外部循环与 c 数组中的依赖项一起使用。

还可以使用row 和col 变量来控制数组的维度。

我使用这个 bash 脚本来运行整个程序

#!/bin/bash
g++ -O0 -o xx par_test.cpp -fopenmp
export OMP_NUM_THREADS=1 # or 2 or 4
./xx

但是当我使用多个 1 线程时，输出会有所不同。 OMP_NUM_THREADS=1 的输出（一致）

OMP_NUM_THREADS=2 的输出（每次执行都会返回不同的内容）

OMP_NUM_THREADS=4 的输出（每次执行都返回不同的东西）

我也尝试过的事情：

带 & 不带矢量化（即带 & 不带 #pragma ivdep）
threadprivate(temp) 这里也没有运气。

由于单线程执行总是具有相同的输出，我相当确定我的变量共享存在问题，尽管它让我无法理解..

提前感谢您的宝贵时间。

【问题讨论】：

您正在使用 openmp 和 simd 指令。这可能会出错。尝试删除#pragma ivdep
@hr0m 根据我的理解，#pragma ivdep 向编译器建议下一个循环可以向量化，不是吗？因此，如果我将其删除，则不会发生矢量化（即已经这样做了）。无论如何，正如我上面提到的，我已经尝试过使用和不使用这个特定的#pragma，但行为保持不变 - 顺便说一句，代码是可测试的。
omfg 我真是瞎了眼！！！您在循环本身中有依赖关系！！！我不知道你为什么使用临时工。如果您正在更改 c[i][j]，它可能已经在上一步中更改。记住你是平行的！
并行化外部循环，而不是内部循环，因为循环携带的依赖关系。

标签： c++ multithreading parallel-processing openmp

【解决方案1】：

您正在此处对j 进行并行循环：

    #pragma omp for schedule(static, 1)
    for(unsigned short j = 0; j < 10; j++){
        // calculate c
        temp[0] = c[i][j-1]+3;
        c[i][j] = temp[0]; 
    }

您引用了单元格c[i][j-1]，但该值可能正在另一个线程中计算。你怎么能确定这已经发生了？

【讨论】：

哦，现在我明白了。所以基本上，计算将是正确的，只要我不使用在同一个 j 循环中计算的c 单元格（即，如果我使用 @ 987654325@ 而不是 c[i][j-1] 都可以吗？）
是的，这是正确的，因为循环中没有改变。
或者即使我使用了c 的固定单元格（例如c[0][0] ），对于计算，它仍然应该正确计算吗？那么，我只需要调节对单元格的写入而不是读取？
如果您使用这些值进行计算，您必须考虑在循环中更改了哪些值。如果您的一个线程更改c[0][0] 而其他线程之一使用它进行计算，那么您就有麻烦了。