OpenMPを使用してRを囲む並列エラーを解決する

Question

RパッケージにOpenMPコードを呼び出そうとしています。 srcディレクトリの下fun.cファイル内のOpenMP地域は次のようである：

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; 

private句にR_NaNとR_NegInfを含めることはできないと思いますので、どうすればいいですか？それとも？すべての操作部分を関数にラップし、この関数をOpenMP領域で呼び出す必要がありますか？

Answer 1

この問題は、firstprivate節にR_NaNなどを追加することで解決します。その理由は、R_NaNとR_NegInfは、含まれているヘッダファイルR.hからのexternal宣言を介して呼び出されるグローバル変数です。