Question

我试图学习使用OpenACC，并以简单的代码实现它来计算Mandelbrot集，但是在编译代码时遇到错误。

我发现了

OpenACC must have routine information error

，并认为在我的代码中使用起来很容易，但是似乎没有用。

这些是我想在其上使用OpenACC的例程，以及我尝试过的

#pragma acc parallel loop copyout(color)

for (int i=0;i<=n0x;i++)
{

Cx=Re_start+dx*i;

  #pragma acc loop

    for (int h=0;h<=n0y;h++)
    {

    Cy=Im_start+dy*h;
    Zx1=0.;
    Zy1=0.;

    #pragma acc loop

        for(int k=0;k<=n;k++ && mod<=4.)
        { 


        Zx=Zx1*Zx1-Zy1*Zy1+Cx;
        Zy=2*Zx1*Zy1+Cy;
        Zx1=Zx;
        Zy1=Zy; 

        mod=(Zx*Zx+Zy*Zy);

            if(mod<=4.)
            {

            if(k==n) 
            color[i][h]=-1; 


            color[i][h]=k+1;

            }

        }  

    }

}

我收到错误消息（我正在使用pgc++ -acc splot.cpp -Minfo=accel进行编译）

PGCC-S-0155-Compiler failed to translate accelerator region (see -Minfo messages): Unknown variable reference (splot.cpp: 48)
main:
     48, Generating copyout(color[:][:])
         Generating Tesla code
         52, #pragma acc loop gang /* blockIdx.x */
         59, #pragma acc loop seq
         68, #pragma acc loop vector(128) /* threadIdx.x */
     52, Accelerator restriction: scalar variable live-out from loop: color
     59, Loop is parallelizable
         Accelerator restriction: scalar variable live-out from loop: color
     68, Loop is parallelizable
         Accelerator restriction: scalar variable live-out from loop: color
PGCC-F-0704-Compilation aborted due to previous errors. (splot.cpp)
PGCC/x86-64 Linux 19.4-0: compilation aborted

我需要修改什么才能正确编译和运行我的代码？

Answer 1

问题是“ for（int k = 0; k <= n; k ++ && mod <= 4。）”中的“ && mod <= 4.0”。 for循环中的复合表达式是PGI的一个已知限制，当前在设备代码生成中不可用。

删除“ && mod <= 4.0”将解决此问题。

我可以提交一个问题报告，尽管我不清楚为什么这部分代码存在，因为这似乎是多余的。你能解释一下意图吗？

如何使用OpenACC计算Mandelbrot集？

1 个答案: