在我的代码中,我有一个循环,在其中我构build和确定的线性系统,并试图解决它: #pragma omp parallel for for (int i = 0; i < n[0]+1; i++) { for (int j = 0; j < n[1]+1; j++) { for (int k = 0; k < n[2]+1; k++) { arma::mat A(max_points, 2); arma::mat y(max_points, 1); // initialize A and y arma::vec solution = solve(A,y); } } } 有时候,程序会非常随机地挂起来,或者解决scheme向量中的结果是NaN。 如果我这样做: […]