私有和共享子句的 C++ OpenMP 计算错误

Question

我有一个for 循环要与 OpenMP 并行化，但有多个计算错误，可能是由于我对 OpenMP 的多线程概念缺乏了解：

for ( int i = -X/2; i < X/2; ++i )
{
    base.y = anchor + i*rho_step;
    temp = some_function( base );
    if( temp > response )
    {
        buffer.y = base.y;
        response = temp;
    }
}

这很好，然后我做了以下更改：

#pragma omp parallel for shared (buffer, response) private(base, temp)
for ( int i = -X/2; i < X/2; ++i )
{
    base.y = anchor + i*rho_step;
    temp = some_function( base );
    if( temp > response )
    {
        buffer.y = base.y;
        response = temp;
    }
}

在这段代码中，buffer.y 和 response 都没有正确的值。据我了解，每个线程都应该有自己的base.y 和temp 副本，它们只是计算的临时变量，buffer 和response 必须共享（它们将存储计算数据），但这并不像我预期的那样工作。

唯一完美的版本是以下，但显然没有性能提升：

omp_lock_t writelock;
omp_init_lock(&writelock);
omp_set_num_threads (4);

#pragma omp parallel for
for ( int i = -X/2; i < X/2; ++i )
{
    omp_set_lock(&writelock);
    base.y = anchor + i*rho_step;
    temp = some_function( base );
    if( temp > response )
    {
        buffer.y = base.y;
        response = temp;
    }
    omp_unset_lock(&writelock);
}
omp_destroy_lock(&writelock);

可能是什么问题？ （anchor 和 rho_step 是这个循环中的常量）

Answer 1

为了让您的代码处理buffer 和response 变量的跨线程，您需要为它们使用一些每个线程的局部变量，并对它们执行最终归约以更新他们的共享副本。

这是它的样子（未经测试）：

#pragma omp parallel firstprivate( base )
{
    auto localResponse = response;
    auto localBuffer = buffer;
    #pragma omp for
    for ( int i = -X/2; i < X/2; ++i )
    {
        base.y = anchor + i * rho_step;
        auto temp = some_function( base );
        if ( temp > localResponse )
        {
            localBuffer.y = base.y;
            localResponse = temp;
        }
    }
    #pragma omp critical
    {
        if ( localResponse > response )
        {
            buffer.y = localBuffer.y;
            response = localResponse;
        }
    }
}