并行处理实验报告：用MPI实现的矩阵乘法的加速比分析(推荐文档).doc

《并行处理实验报告：用MPI实现的矩阵乘法的加速比分析(推荐文档).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、华中科技大学课程名称并行处理实验名称矩阵乘法的实现及加速比分析考生姓名李佩佩考生学号M系、年级计算机软件与理论2013级类别硕士研究生考试日期2014年1月3日一.实验目的1)学会如何使用集群2)掌握怎么用并行或分布式的方式编程3)掌握如何以并行的角度分析一个特定的问题二.实验环境1)硬件环境：4核CPU、2GB内存计算机；2)软件环境：WindowsXP、MPICH2、VS2010、XmanagerEnterprise3；3)集群登录方式：通过远程桌面连接211.69.198.2，用户名：pppusr，密码：AE2Q3P0。三.实验内容1.实验代码编写四个.c文件，分

2、别为DenseMulMatrixMPI.c、DenseMulMatrixSerial.c、SparseMulMatrixMPI.c和SparseMulMatrixSerial.c，用于比较并行和串行矩阵乘法的加速比，以及稀疏矩阵和稠密矩阵的加速比。这里需要说明一下，一开始的时候我是把串、并行放在一个程序中，那么就只有两个.c文件DenseMulMatrix.c和SparseMulMatrix.c，把串行计算矩阵乘的部分放到了主进程中，即procsID=0的进程，但是结果发现执行完串行后，再执行并行就特别的慢。另外，对于稀疏矩阵的处理方面可能不太好，在生成稀疏矩阵的过程中

3、非0元素位置的生成做到了随机化，但是在进行稀疏矩阵乘法时没有对矩阵压缩，所以跟稠密矩阵乘法在计算时间上没多大区别。方阵A和B的初始值是利用rand()和srand()函数随机生成的。根据稀疏矩阵和稠密矩阵的定义，对于稀疏矩阵和稠密矩阵的初始化方法InitMatrix(int*M,int*N,intlen)会有所不同。这里需要说明一下，一开始对于矩阵A和B的初始化是两次调用InitMatrix(int*M,intlen)，生成A和B矩阵，但是随后我发现，由于两次调用方法InitMatrix的时间间隔非常短，又由于srand()函数的特点，导致生成的矩阵A和B完全一样；然后

4、，我就在两次调用之间加入了语句“Sleep(1000)；”，加入头文件“#include”，这样生成的A、B矩阵就不一样了，但很快问题又出现了，在Xshell中不能识别头文件“#include”。所以，最后决定用下面的方法生成矩阵A和B，B是A的转置。//稠密矩阵的生成方法voidInitMatrix(int*M,int*N,intlen){srand((unsigned)time(NULL));for(i=0;i

5、

6、进程数procsNum，对于串行计算，只需要np=1；输出：程序的执行时间。在WindowsXP下使用MicrosoftVisualStudio2010编程，由于稀疏矩阵和稠密矩阵的代码只是初始化部分不同，所以以稠密矩阵乘法为例，列出并行和串行的源代码。并行计算的矩阵乘法源代码：DenseMulMatrixMPI.c#include#include#include#include#defineLength1000int*A,*B,*C,*buffer,*ans;inttemp,i,j,k;intproc

7、sID,procsNum,line;doublestartTime,endTime,totalTime;voidInitMatrix(int*M,int*N,intlen);//实现部分见上面voiddel(){free(A);free(B);free(C);free(buffer);free(ans);}intmain(intargc,char*argv[]){MPI_Statusstatus;MPI_Init(&argc,&argv);MPI_Comm_rank(MPI_COMM_WORLD,&procsID);//获取当前进程号MP