omp вызывает сбой R

Question

Следующая часть функции rl_compute компилируется нормально, но при запуске R выдает ошибку sh. Не знаю почему. Работает с 1 потоком.

// [[Rcpp::export]]
void rl_compute(NumericVector covparms, 
        NumericMatrix locs, 
        NumericMatrix NNarray,
        NumericVector& y, 
        NumericMatrix X)
{
    
    int n = y.length();
    int m = NNarray.ncol();
    int p = X.ncol();
    int nparms = covparms.length();
    int dim = locs.ncol();

#pragma omp parallel num_threads(2)\
    {
    int id = omp_get_thread_num();
    for(int i=id; i<n; i= i+2 ){

    int bsize = std::min(i+1,m);

    NumericMatrix locsub(bsize, dim);
    arma::mat X0( bsize, p );
    arma::vec ysub;
    NumericVector idx = na_omit(NNarray.row(i));
    ...}
    }

}