如何在 OpenACC 中循环内的 std::vectors？

Question

我有以下代码

int main(int argc, char** argv )
{
    std::vector< std::vector< std::vector<double> > > vec { 
        {{1,2},{3,4}, {5,6},{7,8}}, 
        {{9,10}, {11,12}}, 
        {{13,14}, {15,16}, {17,18}} };

    #pragma acc parallel loop
    for (int k = 0; k <3; k++) {

        std::vector<std::vector<double>>& vec2d = vec[k];
        int L = vec2d.size();

        //std::vector<int>dVec{67,51,1,0,50};
        std::vector<double>dVec(L, 0.0);

        for (int i = 0; i < L; i++)
        {
            dVec[i] = vec2d[i][1] - vec2d[i][0];
        }

        for (int j=0; j<2; j++) {
            printf("k: %d j: %d vec0: %f, vec1: %f\n", k, j, vec2d[j][0], vec2d[j][1]);
        }
    }
    std::cout<<"finished\n";

    return 0;
}

我用pgc++ -fast -ta=tesla:cuda9.2,managed -o runEx runEx.cpp -std=c++17 && ./runEx编译

如果我注释掉#pragma acc parallel loop，那么它就可以工作。但是如果我把它留在那里，那么我会得到错误

PGCC-S-0155-Procedures called in a compute region must have acc routine information: operator delete (void *) (runEx.cpp: 425)
PGCC-S-0155-Accelerator region ignored; see -Minfo messages  (runEx.cpp: 6)
PGCC/x86-64 Linux 19.10-0: compilation completed with severe errors

另外，如果我注释掉 std::vector<int>dVec 和包含它的 for 循环，那么即使使用 #pragma acc parallel loop，代码也可以工作

但是，如果我改变循环，它就变成了：

#pragma acc parallel loop
for (int k = 0; k <3; k++) {
    std::vector<int>dVec{67,51,1,0,50};
}

然后我得到同样的错误

这是为什么？

Answer 1

来自网站：https://docs.computecanada.ca/wiki/OpenACC_Tutorial_-_Optimizing_loops

我也能看到东西：

1. 您的声明：

#pragma acc parallel loop

缺少参数。

1. 站点调用函数时使用：#pragma acc 并行循环存在（row_offsets,cols,Acoefs,xcoefs,ycoefs） 这意味着参数为 0 的函数可能是非法的。

Answer 2

这里的问题是 std::vector 有一些在设备上不可用的成员函数。编译错误是专门调出delete，但是我怀疑size这个函数也会有问题，构造函数也一样。由于您不控制std::vector 的来源，因此您无法向它们添加acc routine。我过去所做的解决方法是在 OpenACC 区域之外严格使用向量，并将指向其数据的原始指针传递到区域中。这肯定很麻烦，尤其是对于大型代码，但它确实有效。否则，您也可以尝试实现自己的最小向量类，该类使用acc routine 装饰成员函数。我也看到这个成功了。

对于一个非常大的 L，您可以通过将 acc parallel loop 单独放在您的 i 循环中来获得，但除非您将数组提升到 @ 之外，否则您将大量来回复制数据987654328@ 循环以启用重用。