c/c++语言开发共享使用OpenMP计算PI的值

我正在尝试通过并行化monte carlo代码来学习如何使用OpenMP，该代码使用给定的迭代次数来计算PI的值。 代码的内容是这样的：

int chunk = CHUNKSIZE; count=0; #pragma omp parallel shared(chunk,count) private(i) { #pragma omp for schedule(dynamic,chunk) for ( i=0; i<niter; i++) { x = (double)rand()/RAND_MAX; y = (double)rand()/RAND_MAX; z = x*x+y*y; if (z<=1) count++; } } pi=(double)count/niter*4; printf("# of trials= %d , estimate of pi is %g n",niter,pi); 

虽然在10,000次迭代时，这并没有产生适当的pi值。 如果取出所有OpenMP东西，它就可以正常工作。 我应该提一下，我使用了monte carlo代码： http ： //www.dartmouth.edu/~rc/classes/soft_dev/C_simple_ex.html

我只是用它来尝试学习OpenMP。 有什么想法为什么会聚在1.4ish？ 我不能用多个线程增加变量吗？ 我猜这个问题是变量count 。

谢谢！

 #pragma omp parallel shared(chunk,count) private(i) 

 #pragma omp parallel shared(chunk) private(i,x,y,z) reduction(+:count) 

以上就是c/c++开发分享使用OpenMP计算PI的值相关内容,想了解更多C/C++开发(异常处理)及C/C++游戏开发关注计算机技术网(www.ctvol.com)!)。