主页 / coding / 问题 / 79493035
Accepted
Victor Eijkhout
Asked: 2025-03-08 01:20:44 +0800 CST

OpenMP 声明减少：初始化程序引用的变量既不是“omp_priv”，也不是“omp_orig”

  • 772

省略非必要代码：

double* max_of_two(double *r,double *n);

int npoints = 1000000;
double* xy = (double*)malloc(2*npoints*sizeof(double));

double zero_coord[2] = {0.0,0.0};
double max_coord[2] = {0.0,0.0};

#pragma omp declare reduction\
  (takemax:double*:omp_out=max_of_two(omp_out,omp_in)) \
  initializer(omp_priv=zero_coord)

 #pragma omp parallel for reduction(takemax:max_coord)
  for (int i=0; i<npoints; i++) {
    max_coord = max_of_two(max_coord,xy+2*i);
  }

这会出现错误：

error: '#pragma omp declare reduction' initializer refers to variable 'zero_coord' which is not 'omp_priv' nor 'omp_orig'

这很奇怪，因为在更简单的代码中我可以初始化omp_priv=-1等等。

openmp

1 个回答

  1. Best Answer

    有几个问题：OpenMP 允许按元素减少数组。对于堆栈数组，编译器希望使用按元素减少操作并寻找类型的减少double。如果您想在双精度对上定义减少，则应使用此类类型。

    对于这种类型，您可以简单地用omp_priv = {0,0}或进行初始化omp_priv = omp_orig（适用于最大值，但不适用于一般情况）。

    omp_priv放弃不同于和 的变量的主要原因omp_orig是变量作用域的问题。一般变量不一定在初始化范围内可见。如果要从（全局）变量初始化，则需要将变量包装在函数中并从初始化子句initializer(init(&omp_priv))或调用此函数initializer(omp_priv = init())。

    以下内容使用 gcc 和 clang 进行编译和执行：

    #include <math.h>
#include <stdio.h>
#include <stdlib.h>

typedef struct {
  double x;
  double y;
} pair;

pair max_of_two(pair r, pair n) {
  r.x = r.x > n.x ? r.x : n.x;
  r.y = r.y > n.y ? r.y : n.y;
  return r;
}

int main() {
  int npoints = 100;
  pair *xy = (pair *)malloc(npoints * sizeof(pair));

  pair max_coord = {0.0, 0.0};

#pragma omp declare reduction(takemax                                          \
:pair : omp_out = max_of_two(omp_out, omp_in)) initializer(omp_priv = {0, 0})

#pragma omp parallel for reduction(takemax : max_coord)
  for (int i = 0; i < npoints; i++) {
    xy[i].x = sin(i + .5);
    xy[i].y = cos(i + .5);
    max_coord = max_of_two(max_coord, xy[i]);
  }

  printf("max_coord = %lf, %lf\n", max_coord.x, max_coord.y);
}
    • 2

相关问题