这是正确的 openMP 用法吗？ （或者：我可以相信默认设置吗？）

Question

我目前是第一次使用 openMP，并且对“数据成员不能是私有的”规则感到头疼。

我想知道下面的内容是否有效，或者它最终是否会崩溃：

class network
{
    double tau;
    void SomeFunction();
};

void network::SomeFunction()
{
    #pragma omp parallel for // <-the openMP call
    for (uint iNeu=0;iNeu<nNeurons;++iNeu)
    {
        neurons[iNeu].timeSinceSpike+=tau;  //tau is defined in some other place
        neurons[iNeu].E+=tau*tau;
    }   
}

所以，我使用最小的语法，并让 openMP 自己弄清楚一切。该版本可以编译，并且输出是正确的（到目前为止）。 我之前尝试过的是然而

void network::SomeFunction()
{
    #pragma omp parallel for default(none) shared(neurons) firstprivate(tau)  // <-the openMP call
    for (uint iNeu=0;iNeu<nNeurons;++iNeu)
    {
        neurons[iNeu].timeSinceSpike+=tau; //tau is defined in some other place
        neurons[iNeu].E+=tau*tau;
    }   
}

，正如暗示的那样，这不会编译，大概是因为 tau 和神经元是网络的数据成员。

那么问题是，如果我在第一个版本的运行中真的很幸运，以及我是否必须做类似的事情，

void network::SomeFunction()
{
    double tempTau=tau;
    vector <neuron> tempNeurons=neurons; //in realtity this copy-process would be quite involved
    #pragma omp parallel for shared(tempNeurons) firstprivate(tempTau)// <-the openMP call
    for (uint iNeu=0;iNeu<nNeurons;++iNeu)
    {
        tempNeurons[iNeu].timeSinceSpike+=tempTau;
        tempNeurons[iNeu].E+=tempTau*tempTau;
    }   
}

自然地，我更愿意坚持使用当前版本，因为它是如此简短且易于阅读，但我也想相信我的输出:) 我正在使用 gcc 4.6.1

希望有人能教育我正确的方法。

Answer 1

在此示例中，您最初所做的事情应该没问题：

• 原因是您根本没有修改 tau 成员。所以一开始就没有理由将其设为私有。如果未修改，异步共享相同的值是安全的。
• 至于神经元，您正在独立修改元素。所以这里也没有问题。

当您将变量声明为 firstprivate 时，它会被复制构造到所有线程中。所以 shared(tempNeurons) 绝对不是您想要做的。

Answer 2

http://www.openmp.org/mp-documents/OpenMP3.1.pdf

第 2.9.1 节，数据共享属性规则

• 出现在 threadprivate 指令中的变量是 threadprivate。
• 在构造内部的作用域中声明的具有自动存储期限的变量是私有的。
• 具有动态存储持续时间的对象是共享的。
• 静态数据成员是共享的。
• for 或并行 for 构造的关联 for 循环中的循环迭代变量是私有的。
• 具有 const 限定类型且没有可变成员的变量是共享的。
• 在构造内部的作用域中声明的具有静态存储持续时间的变量是共享的。

...

• for 或并行 for 构造的关联 for 循环中的循环迭代变量可以在 private 或 lastprivate 子句中列出。
• 没有可变成员的 const 限定类型的变量可以在firstprivate 子句中列出。

但是，我仍然怀念构造外部自动变量的默认共享属性。