openmp咋就这么不好调试呢？

hjnln85

  
double x,y,A=0;
#pragma omp parallel for num_threads(2) private(x,y) reduction(+:A)
      for(m=1;m<=num;m++)
        {
          x=m*1.0;         
          for(n=1;n<=num;n++)
                {
               y=n*1.0;
              A=A+(x+y);        
              }
        }
      cout<<"A: "<<A<<endl;

各位大侠，小弟刚学openmp，我就想求个加和的值怎么就不对呢？

井蛙夏虫

我测试了一下，结果没问题。

井蛙夏虫

我测试了3个数:
2为12；10为1100； 100 为1.01e+06；

1. #include <stdlib.h>
   
2. #include <iostream>
   
3. 
   
4. int main(int argc, char*argv[])
   
5. {
   
6.     if (argc != 2)
   
7.     {
   
8.         std::cout << "use: " << argv[0] << " num" << std::endl;
   
9.         return 1;
   
10.     }
    
11. 
    
12.     int num = atoi(argv[1]);
    
13.     int m, n;
    
14.     double x,y,A=0;
    
15. #pragma omp parallel for num_threads(2) private(x,y) reduction(+:A)
    
16.     for(m=1;m<=num;m++)
    
17.     {
    
18.         x=m*1.0;
    
19.         for(n=1;n<=num;n++)
    
20.         {
    
21.             y=n*1.0;
    
22.             A=A+(x+y);
    
23.         }
    
24.     }
    
25.     std::cout<< "A: " << A << std::endl;
    
26.     return 0;
    
27. }
    

复制代码

hjnln85

回复 3# 井蛙夏虫
您是怎么编译的呢？要是单核运行，结果没有错误，但是多核心的话，每次结果都不一样！

   

井蛙夏虫

回复 4# hjnln85
我编译的时候忘了带-fopenmp，结果被忽略了。
你的程序是有问题，你的private(x,y)应当换为private(x,y,n)。
n也应当为私有，否则所有线程共用一个n。


   

hjnln85

回复 5# 井蛙夏虫
你是对的！得把n设置为私有的