内部循环减少的openmp崩溃[重复]答案

【问题标题】：openmp collapse with inner loop reduction [duplicate]内部循环减少的openmp崩溃[重复]
【发布时间】：2021-05-06 05:06:42
【问题描述】：

我正在尝试使用#pragma omp simd 对嵌套循环进行矢量化，但内部循环对其进行了缩减。有没有办法做到这一点？

  #pragma omp simd collapse(2) reduction(+:total[i])
  for (int i=0; i<N; i++){
    for (int j=0; j<N; j++){

      float distance = sqrt(
      (x[i][0]-x[j][0]) * (x[i][0]-x[j][0]) +
      (x[i][1]-x[j][1]) * (x[i][1]-x[j][1]) +
      (x[i][2]-x[j][2]) * (x[i][2]-x[j][2])
      );

      total[i] += distance; 
    }
  }

目前，它只是在编译时出现此错误：

error: identifier "i" is undefined
#pragma omp simd collapse(2) reduction(+:total[i])
                                               ^

【问题讨论】：

标签： c++ c parallel-processing openmp simd

【解决方案1】：

您不能将变量“i”传递给+:total[i]，要么指定要减少的数组小节，要么只执行+:total。但是请记住，您需要一个支持OpenMP 4.5 数组缩减功能的编译器。

  #pragma omp simd collapse(2) reduction(+:total)
  for (int i=0; i<N; i++){
    for (int j=0; j<N; j++){

      float distance = sqrt(
      (x[i][0]-x[j][0]) * (x[i][0]-x[j][0]) +
      (x[i][1]-x[j][1]) * (x[i][1]-x[j][1]) +
      (x[i][2]-x[j][2]) * (x[i][2]-x[j][2])
      );

      total[i] += distance; 
    }
  }

当编译器不支持OpenMP 4.5 数组缩减功能时，可以选择显式实现缩减 (check this SO thread to see how)。

【讨论】：