Openmp 结果不稳健？

Question

我是openmp的新手，当我将openmp添加到我的代码中时，我发现在不同的运行中结果不一样。这是openmp的固有问题还是我的代码问题？谢谢！

#include "stdafx.h"
#include <iostream>
#include <vector>
#include <fstream>
#include <math.h>
#include<sstream>
#include <omp.h>
using namespace std;

int main()
{
    double a[1000];
    for (int i = 0; i < 1000; i++)
    {
        a[i] = 0;
    }
    for (int m = 0; m < 1000; m++)
    {
    #pragma omp parallel for shared(a)
        for (int i = 0; i < 1000000; i++)
        {
            int b = pow(i, 0.5);
            if (b < 1000)
            {
                //cout << i <<" "<<sin(i)<< endl;

                a[b] += sin(i);
            }
        }
    }

    fstream all_temp;
    all_temp.open("temperatureabcd.dat", fstream::out);
    for (int aaa = 0; aaa < 1000; aaa++)
    {
        all_temp << a[aaa] << endl;
    }
    all_temp.close();
    return 0;
}

Answer 1

您的代码正在减少数组。简单的解决方案是做

#pragma omp parallel for reduction(+:a[:1000])

但是，您使用的 MSVC（我从预编译的标头 stdafx.h 推断）不支持 OpenMP 数组缩减。您可以通过像这样更改代码来手动进行数组缩减

double a[1000] = {0};
for (int m = 0; m < 1000; m++) {
  #pragma omp parallel
  {
    double a2[1000] = {0};
    #pragma omp for nowait
    for (int i = 0; i < 1000000; i++) {
      int b = pow(i, 0.5);
      if (b < 1000) a2[b] += sin(i);
    }
    #pragma omp critical
    for(int i = 0; i < 1000; i++) a[i] += a2[i];
  }
}

另一个问题是浮点加法不是关联的，因此减少的顺序很重要。这可以通过更多的工作来解决，但这可能不是您的主要问题。