欢迎来到天天文库
浏览记录
ID:8838450
大小:696.34 KB
页数:43页
时间:2018-04-09
《多核并行编程实验报告》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、多核并行编程实验报告42一、实验目的使用三种编程方法实现积分法或蒙特卡洛法计算pi值(选择实现其中一种,建议蒙特卡洛法)·串行方法·Windows环境多线程编程·Windows环境OpenMP库并行编程二、实验原理·积分法计算pi值的基本原理因此可以用上述公式来计算π的值。·蒙特卡罗法计算pi值的基本原理·蒙特卡洛算法可理解为通过大量实验,模拟实际行为,来收集统计数据。本例中,算法随机产生一系列点,模拟这些点落在如下图所示的正方形区域内的情况。其几何解释如下·如图所示,正方形边长为1,左下顶点与原点重合,两边分别与x轴重合。曲线为1/4圆弧,圆心位于原点,与正方形
2、左下定点重合,半径为1。正方形面积S1=1,圆弧内面积S2=p/4。算法模拟大量点随机落在此正方形区域内,落在圆弧内的点的数量(n2)与点的总数(n1)的比例与面积成正比关系。即42·三、实验环境·系统硬件环境:通过CPU-Z获取本机CPU型号,二级cache容量,处理器内核数和逻辑线程数等信息·操作系统:winxp32位(学校机房电脑)·开发工具:Visualstudio2010四、实验任务和步骤系统硬件环境4242积分法1.传统单线程代码#include#include#include#include
3、staticlongnum_steps=1000000000;//块数可以修改doublestep;usingnamespacestd;intmain(intargc,char*argv[]){inti;doublex,pi,starttime,endtime,sum=0.0;step=1.0/(double)num_steps;time_tt;starttime=clock();srand((unsigned)time(&t));for(i=0;i4、x);}pi=step*sum;endtime=clock();printf("块数=%d",num_steps);printf("Pi=%f",pi);printf("time=%f",endtime-starttime);return0;}42块数1000000块数1000000042块数1000000042块数为1000000002.WinAPI多线程编程代码#include#include#include#include#includeusingn5、amespacestd;staticlongnum_step;constintThreadnum=10;//修改此处线程数为1,4,8doublestep,pi=0,sum=0.0;CRITICAL_SECTIONg;DWORDWINAPIthreadFunc(LPVOIDpParam){42intnum=*((int*)pParam);doublex,sum1=0;step=1.0/(double)num_step;for(inti=num;i6、EnterCriticalSection(&g);sum+=sum1;LeaveCriticalSection(&g);return0;}voidmain(){cout<<"输入分块的数目:";cin>>num_step;LARGE_INTEGERtemp;doubledFreq;//系统时钟QueryPerformanceFrequency(&temp);dFreq=(double)temp.QuadPart;//获得计数器的时刻LONGLONGStart,End;QueryPerformanceCounter(&temp);Start=temp.QuadPar7、t;//获得计数器的初始值HANDLEhthread[Threadnum];InitializeCriticalSection(&g);for(inti=0;i8、uadPa
4、x);}pi=step*sum;endtime=clock();printf("块数=%d",num_steps);printf("Pi=%f",pi);printf("time=%f",endtime-starttime);return0;}42块数1000000块数1000000042块数1000000042块数为1000000002.WinAPI多线程编程代码#include#include#include#include#includeusingn
5、amespacestd;staticlongnum_step;constintThreadnum=10;//修改此处线程数为1,4,8doublestep,pi=0,sum=0.0;CRITICAL_SECTIONg;DWORDWINAPIthreadFunc(LPVOIDpParam){42intnum=*((int*)pParam);doublex,sum1=0;step=1.0/(double)num_step;for(inti=num;i6、EnterCriticalSection(&g);sum+=sum1;LeaveCriticalSection(&g);return0;}voidmain(){cout<<"输入分块的数目:";cin>>num_step;LARGE_INTEGERtemp;doubledFreq;//系统时钟QueryPerformanceFrequency(&temp);dFreq=(double)temp.QuadPart;//获得计数器的时刻LONGLONGStart,End;QueryPerformanceCounter(&temp);Start=temp.QuadPar7、t;//获得计数器的初始值HANDLEhthread[Threadnum];InitializeCriticalSection(&g);for(inti=0;i8、uadPa
6、EnterCriticalSection(&g);sum+=sum1;LeaveCriticalSection(&g);return0;}voidmain(){cout<<"输入分块的数目:";cin>>num_step;LARGE_INTEGERtemp;doubledFreq;//系统时钟QueryPerformanceFrequency(&temp);dFreq=(double)temp.QuadPart;//获得计数器的时刻LONGLONGStart,End;QueryPerformanceCounter(&temp);Start=temp.QuadPar
7、t;//获得计数器的初始值HANDLEhthread[Threadnum];InitializeCriticalSection(&g);for(inti=0;i8、uadPa
8、uadPa
此文档下载收益归作者所有