【openmp】利用openmp实现蒙特卡洛方法计算pi

邱奇迹

于 2022-10-03 09:32:42 发布

阅读量1.4k

点赞数 2

CC 4.0 BY-SA版权

文章标签： c++ 算法开发语言

本文链接：https://blog.csdn.net/qq_44495605/article/details/127149325

本文介绍了一种使用OpenMP并行框架与蒙特卡洛方法相结合来高效计算圆周率π的方法。通过随机生成大量点并判断这些点是否位于单位圆内，可以估算出π的大致数值。该程序利用多线程技术显著提高了计算效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

利用openmp实现蒙特卡洛方法计算pi

#include "stdafx.h"
#include<stdio.h>
#include<omp.h>
#include<stdlib.h>
#include<time.h>

#define MaxSize 204800
int _tmain(int argc, _TCHAR* argv[])
{
	double x,y;
	clock_t start,finish;
	double count=0;
	double total_time;
	srand((unsigned)time(NULL));
	start = clock();
#pragma omp parallel for reduction(+:count)
	for(int i=0;i<MaxSize;i++){
		//printf("%d\n",omp_get_thread_num());
		x=rand()/(double)RAND_MAX;
		y=rand()/(double)RAND_MAX;
		if((x*x+y*y)<1){
			count++;
		}
	}
	finish = clock();
	total_time = (double)(finish-start);
	double pi=(double)(count/MaxSize)*4.0;
	printf("%lf\n",pi);
	printf("%lf",total_time);
	return 0;
}