程序并行化使用OpenMP C++

Question

我试图计算积分

#include <iostream> 
#include <omp.h> 

using namespace std; 

double my_exp(double x) { 
    double res = 1., term = 1.; 
    for(int n=1; n<=1000; n++) { 
    term *= x/n; 
    res += term; 
    } 
    return res; 
} 

double f(double x) { 
    return x*my_exp(x); 
} 


int main() { 
    double a=0., b=1., result = 0.; 
    int nsteps = 1000000; 
    double h = (b - a)/nsteps; 


    for(int i=1; i<nsteps; i++) result += f(a + i*h); 
    result = (result + .5*(f(a) + f(b)))*h; 

    cout.precision(15); 
    cout << "Result: " << result << endl; 
    return 0; 
} 

这项计划数积分和返回结果Result: 1.00000000000035 来算积分。但执行时间很多。 我应该平行我的计划，我想我应该补充#pragma omp parallel for，但它不工作

Answer 1

改变你的主要功能

#pragma omp parallel 
    { 
    double localresult = 0.0; 
#pragma omp for 
    for(int i=1; i<nsteps; i++) 
     localresult += f(a + i*h); 
#pragma omp critical 
    { 
     result += localresult; 
    } 
    } 

    result = (result + .5*(f(a) + f(b)))*h; 

编辑：沿muXXmit2X的线条更简单的解决办法是

#pragma omp parallel for reduction(+:result) 
for(int i=1; i<nsteps; i++) result += f(a + i*h); 

result = (result + .5*(f(a) + f(b)))*h;