并行计算_实验二_矩阵乘法的OpenMP实现及性能分析.docVIP

深 圳 大 学 实 验 报 告 课程名称：并行计算 实验名称：矩阵乘法的OpenMP实现及性能分析 姓 名： 学 号： 班 级： 实验日期：2011年10月21日、11月4日  一. 实验目的 1) 用OpenMP实现最基本的数值算法“矩阵乘法” 2) 掌握for编译制导语句 3) 对并行程序进行简单的性能 二. 实验环境 1) 硬件环境：32核CPU、32G内存计算机； 2) 软件环境：Linux、Win2003、GCC、MPICH、VS2008； 4) Windows登录方式：通过远程桌面连接192.168.150.197，用户名和初始密码都是自己的学号。 三. 实验内容 1. 用OpenMP编写两个n阶的方阵a和b的相乘程序，结果存放在方阵c中，其中乘法用for编译制导语句实现，写出代码。 初始值： 输入： 方阵的阶并行域的线程数 输出：c中所有元素之和程序的执行时间 提示： 用中的clock_t clock( void )函数得到当前程序执行的时间#include sys/time.h timeval start,end; gettimeofday(start,NULL); gettimeofday(end,NULL); coutexecution time: (end.tv_sec-start.tv_sec)+(double)(end.tv_usec-start.tv_usec)/ 1000000seconds endl; 答： 在windows下使用Microsofe Visual Studio编程，源代码如下： #include omp.h #include stdio.h #include time.h #define NN 2000 int a[NN][NN], b[NN][NN]; long long c[NN][NN]; void solve(int n, int num_thread) { int i, j, t, k, time; clock_t startTime, endTime; long long sum; omp_set_num_threads(num_thread); for(i=0;in;i++)//对矩阵a和矩阵b进行初始化 { t=i+1; for(j=0;jn;j++) { a[i][j]=t++; b[i][j]=1; } } startTime=clock(); sum=0; #pragma omp parallel shared(a,b,c) private(i,j,k) { #pragma omp for schedule(dynamic) for(i=0;in;i++) { for(j=0;jn;j++) { c[i][j]=0; for(k=0;kn;k++) { c[i][j]+=a[i][k]*b[k][j]; } } } } for(i=0;in;i++)for(j=0;jn;j++) sum+=c[i][j]; endTime=clock(); time=endTime-startTime; printf(sum=%lld time=%dms\n,sum,time); } int main() { int n, num_thread; while(scanf(%d%d,n,num_thread)!=EOF) { solve(n,num_thread); } return 0; } 2. 分析矩阵相乘程序的执行时间、加速比和效率：方阵阶固定为1000，节点数分别取1、2、4、8、16和32时，为减少误差，每项实验进行5次，取平均值作为实验结果。 答： 串行执行时程序的执行时间为：T = 15.062s 加速比=顺序执行时间/并行执行时间 效率=加速比/节点数 表1 不同节点数下程序的执行时间（秒） 节点数 实验结果 1 2 4 8 16 32 第1次 16.640 8.172 4.078 2.125 1.093 0.594 第2次 16.422 8.156 4.172 2.141 1.078 0.578 第3次 16.406 8.266 4.078 2.125 1.094 0.563 第4次 16.781 8.172 4.079 2.109 1.094 0.563 第5次 16.422 8.171 4.078 2.125 1.093 0.578 平均值 16.5342 8.1874 4.0970 2.1250 1.0904 0.5752 图1 不同节点数下程序的执