并行计算实验四报告.doc

PAGE PAGE 2 广东技术师范学院实验报告 学院： 计算机科学学院 专业： 计算机科学与技术 班级： 成绩： 姓名： 学号： 组别： 组员： 实验地点： 工业中心203 实验日期： 指导教师签名： 预习情况 操作情况 考勤情况 数据处理情况 实验 （四） 项目名称： 并行程序设计(3) 一、实验目的 在一个局域网中建立能够互相通信的两台计算机，实现两台计算机并行运算。 二、实验内容： 要求：修改代码cpi.c，将任务分配改为按块分配：假设共有P个进程，将n个计算间顺序分成P块，每个进程负责一块的计算，注意当n不是P的倍数时应该尽量保持负载平衡。程序代码cpi.c参考如下: ? #include mpi.h #include stdio.h double f( double a ) { return (4.0 / (1.0 + a*a)); } ? int main( int argc, char *argv[]) { int n, myid, numprocs, i, namelen; double PI25DT = 3.141592653589793238462643; double mypi, pi, h, sum, x; double startwtime, endwtime; char processor_name[MPI_MAX_PROCESSOR_NAME]; ? MPI_Init(argc,argv); MPI_Comm_size(MPI_COMM_WORLD,numprocs); MPI_Comm_rank(MPI_COMM_WORLD,myid); MPI_Get_processor_name(processor_name,namelen); fprintf(stderr,Process %d on %s, myid, processor_name); if (myid == 0) { n=10000; startwtime = MPI_Wtime(); } MPI_Bcast(n, 1, MPI_INT, 0, MPI_COMM_WORLD); h = 1.0 / (double) n; sum = 0.0; for (i = myid; i n; i += numprocs) { x = h * ((double)i + 0.5); sum += f(x); } mypi = h * sum; MPI_Reduce(mypi, pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD); if (myid == 0) { endwtime = MPI_Wtime(); printf(pi is approximately %.16f, error is %.16f, pi, pi - PI25DT); printf(wall clock time = %f, endwtime-startwtime); } MPI_Finalize(); return 0; } 实验源程序代码： #include mpi.h #include stdio.h #include iostream.h #include fstream.h void main(int argc, char *argv[]) { int rank, size; const int MAXX = 40; const int MAXY = 40; const int MAXPROCESSOR = 64; float Data[MAXY][MAXY]; int ArraySize[2]; int i, j, k; double startwtime, endwtime; int namelen; char processor_name[MPI_MAX_PROCESSOR_NAME]; MPI_Status 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); cout线程在运行主机processor_nameendl; if(rank == 0) {