Question

我正在编写一个具有数组依赖性的OpenACC代码。内循环的每次迭代都可以更新数组的相同位置。这是一些代码：

    long unsigned int digits[d + 11];
    for (long unsigned int digit = 0; digit < d + 11; ++digit)
            digits[digit] = 0;

    for (long unsigned int i = 1; i <= n; ++i) {
            long unsigned int remainder = 1;
            for (long unsigned int digit = 0; digit < d + 11 && remainder; ++digit) {
                    long unsigned int div = remainder / i;
                    long unsigned int mod = remainder % i;
                    digits[digit] += div; // here
                    remainder = mod * 10;
            }
    }

OpenMP版本编写如下：

    #pragma omp parallel private(i)
    {
            long unsigned int digit_local[d+11];
            for(i=0;i<d+11;i++)
                    digit_local[i] = 0;

            #pragma omp for
            for (i = 1; i <= n; ++i) {
                    long unsigned int remainder = 1;
                    for (long unsigned int digit = 0; digit < d + 11 && remainder; ++digit) {
                            long unsigned int div = remainder / i;
                            long unsigned int mod = remainder % i;
                            digit_local[digit] += div;
                            remainder = mod * 10;
                    }
            }

            #pragma omp critical
            for(long unsigned int digit = 0; digit < d+11; ++digit)
                    digits[digit] += digit_local[digit];

    }

在OpenACC中，关键字private与数组一起工作，但我不知道如何将私有数组与全局数组连接。

感谢。

Answer 1

您将使用OpenACC“原子更新”指令。

    long unsigned int digit_local[d+11][n];
    #pragma acc data create(digit_local) copyout(digits)
    {

    #pragma acc parallel loop gang vector
    for (i = 1; i <= n; ++i) {
          for(j=0;j<d+11;j++) digit_local[j][i] = 0;
            long unsigned int remainder = 1;
            for (long unsigned int digit = 0; digit < d + 11 && remainder; ++digit) {
                    long unsigned int div = remainder / i;
                    long unsigned int mod = remainder % i;
                    digit_local[digit][i] += div;
                    remainder = mod * 10;
            }
    }

    #pragma acc parallel loop gang
    for(long unsigned int digit = 0; digit < d+11; ++digit) {
          long unsigned int dsum = 0;
          #pragma acc loop vector reduction(+:dsum)
          for (i = 1; i <= n; ++i) {
            dsum += digit_local[digit][i];
          }
          digits[digit] = dsum;
    }
    }

或者，您可以执行与OpenMP版本类似的操作。

{{1}}

尽管如此，我不确定这会有什么加速。

希望这有帮助，垫

在OpenACC中加入数组结果

1 个答案: