未使用collapse子句的OpenMP数据竞争问题排查

核心问题：共享变量col引发的竞争条件 #

你的代码中，col变量是在main函数的全局作用域（相对于并行区域）声明的，属于OpenMP并行区域的共享变量。当多个线程并行处理外层row循环时，每个线程在执行内层col循环时都会修改这个共享的col值——线程间的互相干扰会导致col的取值混乱，进而错误地访问数组元素，最终出现结果不一致的情况。

• 单线程运行时，不存在线程间干扰，col的修改是串行执行的，因此结果固定符合预期。
• 添加collapse(2)后，OpenMP会将两层循环合并为一个统一的迭代空间，row和col都会被自动私有化，每个迭代由独立线程处理，彻底避免了共享col带来的竞争，因此结果始终正确。

解决方案（无需使用collapse） #

你只需要确保col是每个线程的私有变量，有两种简单实现方式：

方式1：将col的声明移至内层循环内部

让col成为内层循环的局部变量，每个线程的内层循环都会拥有独立的col实例：

#include <stdio.h>
#include <stdlib.h>

int main()
{
    size_t nrows = 10, ncols = 10;
    unsigned int *cells;
    int row;

    cells = calloc(nrows * ncols, sizeof *cells);

#pragma omp parallel for
    for (row = 0; row < nrows; row++) {
        int col; // 将col移至此处，变为线程私有
        for (col = 0; col < ncols; col++)
            if (row * col % 10)
                cells[row * ncols + col] = 1; // 修正索引公式：row*ncols而非row*nrows，避免后续nrows/ncols不等时出错
    }

    for (row = 0; row < nrows; row++) {
        for (int col = 0; col < ncols; col++)
            printf("%d ", cells[row * ncols + col]);
        printf("\n");
    }
    free(cells); // 别忘了释放内存
}

方式2：使用private(col)显式声明私有变量

通过OpenMP的private子句，强制将col设为线程私有：

#pragma omp parallel for private(col)
for (row = 0; row < nrows; row++) {
    for (col = 0; col < ncols; col++)
        if (row * col % 10)
            cells[row * ncols + col] = 1;
}

补充说明 #

• 你猜测的伪共享并非问题根源，这里的错误是共享循环变量导致的逻辑竞争，而非缓存行冲突。因此添加atomic或flush无法解决问题。
• 修正数组索引公式为row * ncols + col，这是二维数组按行存储的正确索引方式（即使当前nrows=ncols，也建议修正以避免后续参数变化时出错）。

内容的提问来源于stack exchange，提问作者Huidae Cho