c - 解决在R中使用OpenMP时的封闭并行错误

Question

我尝试调用OpenMP我的R中的代码包裹。 fun.c 中的 OpenMP 区域src下的文件目录如下:

    static void mp_func(double *x, double *a, double *b,
                           double *den, int *P) {
        #pragma omp parallel for num_threads(P[0]) default(none)          \
            firstprivate(length, q, a, b, isLog, lowTail, cdf, R_NaN) private(i) \
                reduction(||:naflag)
                for (i = 0; i < length; i++){
                    ...
                    if(some condition) den[i] = R_NaN;
                    ...
                    if(some condition) den[i] = R_NegInf;
                }
    }

但是，当我运行 R CMD CHECK 时，`它给出的错误消息是:

 enclosing parallel;

和

R_NaN not specified in enclosing parallel;

我不认为包括R_NaN和R_NegInf在 private条款是可能的，那么我该如何解决呢？我是否应该将所有操作部分包装到一个函数中并在 OpenMP 区域中调用该函数？

Answer 1

好吧，这个问题可以通过在firstprivate子句中添加R_NaN等来解决。这样做的原因是 R_NaN 和 R_NegInf 是通过包含的头文件 R.h 中的 external 声明调用的全局变量>.