OpenMP并行for循环数组赋值出现随机结果问题求助

核心问题：数据竞争与写入顺序混乱 #

你的代码中，每个迭代i会同时修改d[i]、d[i+1]、d[i+2]三个位置，结合schedule(static,64)的块划分（4个线程各负责64个迭代），必然出现不同线程交叉写入同一数组元素的情况：

• 负责i=0~63的线程0，当i=62时会修改d[64]；负责i=64~127的线程1，当i=64时也会修改d[64]。两个线程的写入操作并行执行，最终d[64]的值取决于哪个线程先完成写入，导致结果随机。
• 同理，d[128]会被线程1的i=126（写入c[126]）和线程2的i=128（写入a[128]）交叉修改，出现错误。

你对static调度的理解有误：它只是给线程分配固定的迭代块，线程之间仍然是并行执行，而非按顺序串行。即使添加#pragma omp critical，如果只是包裹单个赋值操作，无法解决写入顺序的问题；如果包裹整个循环体，会直接让循环串行执行，完全失去并行意义。

解决方案 #

方案1：修正代码逻辑，消除冗余写入 #

观察串行执行的结果，每个d[j]最终会被i=j的迭代写入a[j]，之前的d[j] = b[j-1]和d[j] = c[j-2]都会被覆盖。如果这是你的真实需求，可直接简化代码：

#pragma omp parallel num_threads(4)
{
    #pragma omp for schedule(static, 64)
    for(int i = 0; i < 256; i++)
    {
        d[i] = a[i];
    }
}

方案2：调整迭代划分，避免跨块写入 #

如果必须保留原有的多位置写入逻辑，需要确保每个d元素只被一个线程修改。可以手动划分迭代范围，让每个线程负责的区域包含所有会被它修改的位置：

#pragma omp parallel num_threads(4)
{
    int tid = omp_get_thread_num();
    // 每个线程处理的起始i比块起始早2，避免跨块写入后续线程的区域
    int start = tid * 64;
    int end = (tid + 1) * 64;
    // 最后一个线程处理到255
    if(tid == 3) end = 256;
    // 起始位置不能小于0
    if(start > 2) start -= 2;

    for(int i = start; i < end; i++)
    {
        if(i >=0 && i <256) d[i] = a[i];
        if( i + 1 < 256 ) d[i+1] = b[i];
        if( i + 2 < 256 ) d[i+2] = c[i];
    }
}

这种方式下，每个线程会提前处理前2个迭代，确保后续线程负责的区域不会被其他线程修改，同时避免越界。

方案3：使用ordered指令保证执行顺序 #

如果需要严格按照串行的迭代顺序执行写入操作，可添加ordered指令，但这会牺牲大部分并行性能：

#pragma omp parallel num_threads(4)
{
    #pragma omp for schedule(static, 64) ordered
    for(int i = 0; i < 256; i++)
    {
        #pragma omp ordered
        {
            d[i] = a[i];
            if( i + 1 < 256 )
                d[i+1] = b[i];
            if( i + 2 < 256 )
                d[i+2] = c[i];
        }
    }
}

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