Question

最近我一直在研究C ++。素数查找器是一个练习。在我编写了串行版本代码之后，我尝试使用openmp来并行化代码。但是，串行版本工作正常但是一旦我尝试了openmp版本，输出（1到100之间的素数）是非常错误的。有什么建议吗？

这是我的代码（相关功能）：

    void prime2_bits(const int& n)
    {
       omp_set_num_threads(4);
       int i,j,tp,sq=sqrt(n+1);
       int p[n/32+1];
       initint(p,n/32+1);
    #pragma omp parallel for default(shared) private(i,j,tp) schedule(dynamic)
       for(i=2;i<sq;i++){
         if(!((p[i/32]>>(i%32))&1)){ // see whether the ith bit is 0
           tp=i*i;
           for(j=tp;j<n+1;j+=i){
             p[j/32]|=(1<<(j%32)); // set the jth bit as 1
           }
         }
       }
    // Do the printing job
       for(i=2;i<n+1;i++)
         if(!((p[i/32]>>(i%32))&1)) cout<<i<<' ';
    }

，输出如下：

using method 2 (bits to save space):
2 3 5 7 11 13 15 17 19 21 23 27 29 31 35 37 41 43 47 49 53 55 59 61 67 69 71 73 77 79 81 83 87 89 91 93 97 99 0 msec used!

另外，为了提高效率，我再次修改了串行代码。然后尝试将它与openmp并行而没有成功。有什么建议吗？

修改后的代码：

void prime3_norev(const int& n)
{
   int i,j,tp,m=n+1;
   int pi=0;
   bool pFlag[m];
   int p[n/3];
   omp_set_num_threads(4);
   initbool(pFlag,m);
   initint(p,n/3);
#pragma omp parallel for default(shared) private(i,j) schedule(dynamic) 
   for(i=2;i<m;i++){
     if(!pFlag[i]) p[pi++]=i;
     for(j=0;(j<pi)&&(i*p[j]<m);j++){
       pFlag[i*p[j]]=true;
       if(i%p[j]==0) break;
     }
   }
// Do the printing job
   for(i=0;i<pi;i++)
     cout<<p[i]<<' ';
}

和相应的输出：

using method 3 (no revisits, space for time):
2 3 6 7 11 13 15 19 23 29 31 35 37 41 43 47 53 55 59 61 65 67 71 73 79 83 85 89 95 97 0 msec used!

代码编译为：

g++ -O2 -fopenmp fib.cc

非常感谢！

Answer 1

要修复代码，您可以使用有序的子句

#pragma omp parallel for schedule(dynamic) ordered
for(int i=2;i<m;i++) {
    #pragma omp ordered
    if(!pFlag[i]) p[pi++]=i;
    //rest of code

有关工作示例，请参阅http://coliru.stacked-crooked.com/a/fa1eebf126940fb0。我使用schedule（static）只是为了在没有ordered子句的情况下更容易显示错误。

我不确定订购后会对性能产生什么影响。在我看来，你正试图将Eratosthenes的筛子并行化。在这种情况下，我建议您查看http://create.stephan-brumme.com/eratosthenes/。它使用OpenMP在大约1秒钟内找到所有素数达到10 ^ 9。

基于并行位操作的寻星器与openMP？

1 个答案: