Pragma omp parallel + ntl

Question

我尝试使用openmp运行以下代码，但是我遇到了分段错误

void modKeyGenPrs(mat_GF2E *&Prs, mat_GF2E Lst[], mat_GF2E L1, mat_GF2E L2) {
    Prs = new mat_GF2E[m];
    mat_GF2E L1_trans = transpose(L1);
        #pragma omp parallel shared(L1_trans,L2,Lst,Prs,L1) 
        {
               #pragma omp for
               for (int i = 0; i < m; i++) {
                   (Prs[i]).SetDims(n, n);
                   for (int j = 0; j < m; j++) {
                       Prs[i] = Prs[i] + (L2[i][j] * (L1_trans * (Lst[i]) * L1));
                   }
               }
       }
}