并行计算实验三报告

PAGEPAGE2广东技术师范学院实验报告学院：计算机科学学院专业：计算机科学与技术班级：成绩：姓名：学号：组别：组员：实验地点：工业中心203实验日期：指导教师签名：预习情况操作情况考勤情况数据处理情况实验（三）项目名称：并行程序设计(2)一、实验目的在一个局域网中建立能够互相通信的两台计算机，实现两台计算机并行运算。二、实验内容：要求：用MPI_Isend和MPI_Irecv代替代码deadlock.c中的MPI_send和MPI_recv来避免出现程序死锁。程序代码deadlock.c参考如下#include<stdio.h>#include<string.h>#include"mpi.h"#defineSIZE16intmain(intargc,char**argv){staticintbuf1[SIZE],buf2[SIZE];intnprocs,rank,tag,src,dst;MPI_Statusstatus;MPI_Init(&argc,&argv);MPI_Comm_size(MPI_COMM_WORLD,&nprocs);/*获取总进程数*/MPI_Comm_rank(MPI_COMM_WORLD,&rank);/*获取本进程的进程号*//*初始化buf1*/memset(buf1,1,SIZE);tag=123;dst=(rank>=nprocs-1)?0:rank+1;src=(rank==0)?nprocs-1:rank-1;MPI_Send(buf1,SIZE,MPI_INT,dst,tag,MPI_COMM_WORLD);MPI_Recv(buf2,SIZE,MPI_INT,src,tag,MPI_COMM_WORLD,&status);MPI_Finalize();return0;}实验源程序代码：#include"mpi.h"#include<stdio.h>#include<iostream.h>#include<fstream.h>voidmain(intargc,char*argv[]){ intrank,size; constintMAXX=8; constintMAXY=3; constintMAXPROCESSOR=64; intmyid,namelen; floatData[8][8]; intArraySize[2]; inti,j,k; doublestartwtime,endwtime; charprocessor_name[MPI_MAX_PROCESSOR_NAME]; MPI_Request*Status; MPI_Init(&argc,&argv); MPI_Comm_rank(MPI_COMM_WORLD,&rank); MPI_Comm_size(MPI_COMM_WORLD,&size); MPI_Get_processor_name(processor_name,&namelen); if(rank==0) { ifstreamin("input.txt"); in>>ArraySize[0]>>ArraySize[1]; for(i=0;i<ArraySize[0];i++) for(j=0;j<ArraySize[1];j++) in>>Data[i][j]; cout<<"二维数组如下:"<<endl; for(i=0;i<ArraySize[0];i++) { for(j=0;j<ArraySize[1];j++) cout<<Data[i][j]<<""; cout<<endl; } } MPI_Bcast(ArraySize,2,MPI_INT,0,MPI_COMM_WORLD); intAverageLineNumber,HeavyProcessorNumber,MyLineNumber; intCurrentLine,StartLine,SendSize; floatSendArray[8*3]; MPI_Request*IERROR; AverageLineNumber=ArraySize[0]/size; HeavyProcessorNumber=ArraySize[0]%size; if(rank<HeavyProcessorNumber) MyLineNumber=AverageLineNumber+1; else MyLineNumber=AverageLineNumber; if(rank==0) { CurrentLine=ArraySize[0]; for(i=size-1;i>=0;i--) { SendSize=0; if(i==0) startwtime=MPI_Wtime(); if(i<HeavyProcessorNumber) StartLine=CurrentLine-AverageLineNumber; else StartLine=CurrentLine-AverageLineNumber+1; for(j=StartLine;j<CurrentLine;j++) for(k=0;k<ArraySize[1];k++) SendArray[SendSize++]=Data[j][k]; if(i!=0) MPI_Send(SendArray,SendSize,MPI_FLOAT,i,10,MPI_COMM_WORLD); CurrentLine=StartLine-1; } } else MPI_Irecv(SendArray,MyLineNumber*ArraySize[1],MPI_FLOAT,0,10,MPI_COMM_WORLD,Status); floatSum; Sum=0; for(i=0;i<MyLineNumber*ArraySize[1];i++) Sum+=SendArray[i]; floatALLSum[64]; MPI_Gather(&Sum,1,MPI_FLOAT,ALLSum,1,MPI_FLOAT,0,MPI_COMM_WORLD); Sum=0; if(rank==0) { for(i=0;i<size;i++) Sum+=ALLSum[i]; endwtime