C++ 特定于循环的OpenMP并行化

我刚开始使用OpenMp，我有一个问题

我有一个如下代码：

Eigen::Matrix<double,1000,1000> mat;
for(int i=0;i<1000;i++){
   for(int j=0;j<1000;j++){
        double aux = 0;
        // code with for loops and iterators and vectors that grow and
        // finally assign value to aux....
        mat(i,j) = aux;
   }
}

Eigen:：mat矩阵；
for（int i=0；It现在显示的代码不是线程安全的？这可能是一个有趣的问题，但如果您想要一个好的答案，您需要提供更多信息。您需要更多地说明“使用for循环、迭代器和向量进行代码的增长和扩展”。提供更多代码。
Eigen::Matrix<double,1000,1000> mat;
#pragma omp parallel for
for(int i=0;i<1000;i++){
   for(int j=0;j<1000;j++){
        double aux = 0;
        // code with for loops and iterators and vectors that grow and
        // finally assign value to aux....
        mat(i,j) = aux;
   }
}