为什么在英特尔编译器上使用 OpenMP 的 firstprivate 和 std::vector 时会出现未定义的行为？

Question

在英特尔 c++ 编译器上将 OpenMP 与 firstprivate 和 std::vector 结合使用时出现问题。取以下三个函数：

#include <omp.h>

void pass_vector_by_value(std::vector<double> p) {
#pragma omp parallel
{
    //do sth
}
}

void pass_vector_by_value_and_use_firstprivate(std::vector<double> p) {
#pragma omp parallel firstprivate(p)
{
    //do sth
}
}

void create_vector_locally_and_use_firstprivate() {
std::vector<double> p(3, 7);
#pragma omp parallel firstprivate(p)
{
    //do sth
}
}

代码编译时没有警告：

icc filename.cpp -openmp -Wall -pedantic

（icc 14.0.1版（gcc 4.7.0版兼容性））

或：

g++ filename.cpp -fopenmp -Wall -pedantic

(gcc 版本 4.7.2 20130108 [gcc-4_7-branch 修订版 195012] (SUSE Linux))

但在使用 icc 编译后，我遇到了运行时错误，例如：

*** Error in `./a.out': munmap_chunk(): invalid pointer: 0x00007fff31bcc980 ***

调用第二个函数时（pass_vector_by_value_and_use_firstprivate）

因此，只有在使用 firstprivate 子句（应该调用复制构造函数）并且向量按值传递给函数（也应该调用复制构造函数）时才会发生错误。当不传递向量但在函数中本地创建它或不使用 firstprivate 时，没有错误！在 gcc 上我没有收到任何错误。

我想知道代码是否会以某种方式产生未定义的行为，或者这是否是 icc 中的错误？

Answer 1

我在 ICC 上遇到了同样的问题，但在 GCC 上没有。看起来像一个错误。这是一个解决方法

void pass_vector_by_value2(std::vector<double> p) {
    #pragma omp parallel
    {
        std::vector<double> p_private = p;
        //do sth with p_private      
    }
}

另一方面，总的来说，我不会将非 POD 按值传递给函数。我会使用参考，但如果你这样做，你会得到错误

error: ‘p’ has reference type for ‘firstprivate’

解决方案是我在上面发布的代码。通过值或引用传递它，然后在并行区域内定义一个私有副本，就像我在上面的代码中所做的那样。