如何使用 OpenMP 并行化这个内循环？

Question

我有多个嵌套循环，我在并行第二个内部循环时遇到问题。部分原因是我正在使用 Eigen 库，并且在第二个内部循环中，我正在计算复杂矩阵的总和，该矩阵将在第一个内部循环中使用。这是一个代码。

for (int i = 0; i < 4; ++i) {
#pragma omp parallel for
    for (int n = 0; n < 1000000; ++n) {
        //initialize sum to zero 
        for (int j = 0; j < 4; ++j) {
            if (j != i) {
                //sum matrices
            }
        }
        //use it here
    }
}

并行化此循环的最佳方法是什么，我可以在那里执行复杂矩阵的并行求和吗？