在仍在使用OpenMP的功能的同时，如何获得私人向量？

Question

我试图与OpenMP并行化我的代码。

我有一个全球向量，因此我可以用自己的功能过多。 有什么方法可以使矢量的副本与每个线程相关，以便他们可以使用它？ 是描述我的问题的一些

double var = 1;
std::vector<double> vec;

void function()
{
    vec.push_back(var);
    return;
}


int main()
{
    omp_set_num_threads(2);

    #pragma omp parallel 
    {
        #pragma omp for private(vec)
        for (int i = 0; i < 4; i++)
        {        
            function(); 
        }
    }
    return 0;
}

码

• 假
• 这 在向量（更改特定值）上的工作
• （在我的原始代码中，有许多向量和功能，IVE只是试图解决问题）

IVE尝试了#pragma op threadprivate（），但这仅适用于变量而不是对矢量。 同样，对矢量进行重新启动并不有所帮助，因为我的功能始终与全局矢量一起使用，然后当不同的胎面同时称其为问题时，这会导致问题。

Answer 1

有办法可以将向量的副本分配给每个线程
他们可以用它做事吗？

是的， firstprivate 条款条款可以做到这一点：

第一个私人子句声明一个或多个列表项目是私人的
要执行任务，并以每个任务的初始化为
相应的原始物品在遇到构造时具有。

因此，它为每个线程创建一个变量的私有副本，但是该私有变量的范围是OpenMP构造之后的结构化块。在此块之外，您可以访问全局变量：

#pragma omp ... firstprivate(vec)
{                
      vec.push_back(...);  // private copy is changed here, which is threadsafe
}

void function()
{
    vec.push_back(var); // the global variable is changed here, which is not threadsafe
    return;
}

如果您希望在函数中使用变量的私有副本，则必须将其传递给您的函数：

void function(std::vector<double>& x, double y)
{
    x.push_back(y);
    return;
}
...
#pragma omp for firstprivate(vec)
   for (int i = 0; i < 4; i++)
   {        
       function(vec, 1); 

   }

但是，请注意，如@Jeromerichard所指出的那样，您应该应该不使用代码中的全局变量。