1

如何在openacc中正确地将数组设为私有？为了我的使用，我声明了一个通过递归函数并自行更新的数组。在 C++ 中，我是这样做的，

int recursion(int array[],int i)
{
   ....
   return array[i] = recursion(array,i);
}
int main()
{
  ....
   for(int run=0;run<N;++run)
   { 
    .... 
      for(int i=0;i<some N;++i)
      int f = recursion(array,i);
   }
 .....
}

现在，当我尝试使用 openacc 指令并行执行此操作时，主要问题就开始了。我想将此数组复制到并行 for 循环区域中，这样每个团伙都将拥有该数组的副本，并且能够使用递归函数更改他们自己的数组 [] 版本，而不是其他版本。我尝试这样做的方式是

#pragma acc routine seq
int recursion(int array[],int i)
{
 ....
 return array[i] = recursion(array,i);
}
int main()
{
 ....
 #pragma acc data copyin(array[0:N])
 #pragma acc parallel loop gang private(array[0:N])

 for(int run=0;run<N;++run)
 {  
  ....
    for(int i=0;i<N;++i)
    { 
     ....
     int f = recursion(array,i);
    }
  }
 }

但似乎数组没有被传递给递归函数，因为我已经检查过它没有改变。做到这一点的完美方法是什么？

ps 我也试过#pragma acc data pcreate(array[0:N])，#pragma acc parallel loop independent private(array[0:N])但结果是一样的

你可以在这里找到整个代码。它在没有指令的情况下完美运行。您唯一需要更改的是用 curand 注释该行，并在 251 268 行用 rand 取消注释该行。请帮忙！

4

1 回答 1

1

“私有”数组未初始化，但您的代码期望它们具有初始值。您要么想要在循环中初始化数组，我们使用“firstprivate”子句，它将使用主机的初始值初始化每个私有数组。

此外，您在私有子句和数据子句中都有“随机”和“ptr”。您应该从 data 指令中删除它们，因为这将使它们成为全局的。

鉴于 GPU 上的堆栈非常小（8MB），因此设备上的递归是有问题的。您可以通过将 PGI 环境变量“PGI_CUDA_STACKSIZE”设置为更大的值来增加此值，但如果递归太深，程序会崩溃。

我尝试运行您的程序，但它一直出错，可能是由于堆栈溢出。我不确定要使用什么输入值，所以它也可能是我的飞行员错误。

如果使用“firstprivate”没有帮助，请告诉我要使用的输入值，我会进一步研究。

于 2017-10-31T16:23:50.883 回答