我正在尝试使用 OpenMP 并行化以下 C++ 代码：
int np = 1000000;
double kk = 1 / pow(2 * pi, 2);
for (int kes = 1; kes <= 100; kes++) {
double E1 = 0;
#pragma omp parallel for reduction(+: E1)
for (int ies = 0; ies < np; ies++) {
for (int jes = 0; jes < np; jes++) {
if (ies != jes) {
float distanes = sqrt(pow(xp[ies] - xp[jes], 2) + pow(yp[ies] - yp[jes], 2) + pow(zp[ies] - zp[jes], 2)); …