Openmp结果不稳健

Openmp results are not robust?

本文关键字：结果 Openmp      更新时间：2023-10-16

我是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;
}

您的代码正在进行数组缩减。简单的解决方案是进行

#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];
}
}

另一个问题是浮点加法是不关联的，所以进行归约的顺序很重要。这可以通过更多的工作来解决，但这可能不是你的主要问题。

相关文章：