-ta=tesla:deepcopy 标志和#pragma acc 形状

Question

我刚刚发现了 deepcopy 标志。 直到此刻，我一直使用-ta=tesla:managed来处理深拷贝，我想探索替代方案。 我读了这篇文章： https://www.pgroup.com/blogs/posts/deep-copy-beta.htm 写得很好，但我认为它不包括我的情况。 我有这种类型的结构：

typedef struct Data_{
  double ****Vc;   
  double ****Uc;
} Data

这些到数组的形状不是由结构本身的元素定义的，而是由另一个结构的元素定义的，并且它们本身仅在程序执行期间定义。

在这种情况下如何使用#pragma acc shape(Vc, Uc) ？

如果没有此编译指示并复制结构如下：

int main(){
  Data     data;

  Initialize(&data);

}

int Initialize(Data *data){

  data->Uc = ARRAY_4D(ntot[KDIR], ntot[JDIR], ntot[IDIR], NVAR, double);
  data->Vc = ARRAY_4D(NVAR, ntot[KDIR], ntot[JDIR], ntot[IDIR], double);
  #pragma acc enter data copyin(data)

  PrimToCons3D(data->Vc, data->Uc, grid, NULL);

}

void PrimToCons3D(double ****V, double ****U, Grid *grid, RBox *box){
  #pragma acc parallel loop collapse(3) present(V[:NVAR][:nx3_tot][:nx2_tot][:nx1_tot])
  for (k = kbeg; k <= kend; k++){
  for (j = jbeg; j <= jend; j++){
  for (i = ibeg; i <= iend; i++){
    double v[NVAR];

    #pragma acc loop
    for (nv = 0; nv < NVAR; nv++) v[nv] = V[nv][k][j][i];
}

我收到致命错误：在设备 1 上找不到 PRESENT 子句中的数据：名称=V 主机：0x1fd2b80 文件：/home/Prova/Src/mappers3D.c PrimToCons3D 行：140

顺便说一句，同样的代码适用于-ta=tesla:managed 。

Answer 1

由于您没有提供完整的复制示例，因此我无法对此进行测试，但它看起来像：

typedef struct Data_{
  int i,j,k,l;
  double ****Vc;
  double ****Uc;
  #pragma acc shape(Vc[0:k][0:j][0:i][0:l])
  #pragma acc shape(Uc[0:k][0:j][0:i][0:l])
} Data;

int Initialize(Data *data){
  data->Vc.i = ntot[IDIR];
  data->Vc.j = ntot[JDIR];
  data->Vc.k = ntot[KDIR];
  data->Vc.l = NVAR;
  data->Uc.i = ntot[IDIR];
  data->Uc.j = ntot[JDIR];
  data->Uc.k = ntot[KDIR];
  data->Uc.l = NVAR;
  data->Uc = ARRAY_4D(ntot[KDIR], ntot[JDIR], ntot[IDIR], NVAR, double);
  data->Vc = ARRAY_4D(NVAR, ntot[KDIR], ntot[JDIR], ntot[IDIR], double);
  #pragma acc enter data copyin(data)
  PrimToCons3D(data->Vc, data->Uc, grid, NULL);
}

void PrimToCons3D(double ****V, double ****U, Grid *grid, RBox *box){
  int kbeg, jbeg, ibeg, kend, jend, iend;


  #pragma acc parallel loop collapse(3) present(V, U)
  for (int k = kbeg; k <= kend; k++){
  for (int j = jbeg; j <= jend; j++){
  for (int i = ibeg; i <= iend; i++){

尽管请记住，OpenACC 标准没有采用我们的“形状”和“策略”指令，我们（NVHPC 编译器团队）只做了一个 Beta 版本，我们没有维护。

可能最好进行手动深度复制，这将符合标准，如果您可以提供包含您如何进行数组分配的复制器，我可以提供帮助，即“ARRAY_4D”。