快速排序代码.docVIP

#include stdio.h #include stdlib.h #include mpi.h #define TRUE 1 /* * 函数名: main * 功能：实现快速排序的主程序 * 输入：argc为命令行参数个数； * argv为每个命令行参数组成的字符串数组。 * 输出：返回0代表程序正常结束 */ main(int argc,char *argv[]) { int DataSize; int *data; /*MyID表示进程标志符；SumID表示组内进程数*/ int MyID, SumID; int i, j; int m, r; MPI_Status status; /*启动MPI计算*/ MPI_Init(argc,argv); /*MPI_COMM_WORLD是通信子*/ /*确定自己的进程标志符MyID*/ MPI_Comm_rank(MPI_COMM_WORLD,MyID); /*组内进程数是SumID*/ MPI_Comm_size(MPI_COMM_WORLD,SumID); /*根处理机(MyID=0)获取必要信息，并分配各处理机进行工作*/ if(MyID==0) { /*获取待排序数组的长度*/ DataSize=GetDataSize(); data=(int *)malloc(DataSize*sizeof(int)); /*内存分配错误*/ if(data==0) ErrMsg(Malloc memory error!); /*动态生成待排序序列*/ srand(396); for(i=0;iDataSize;i++) { data[i]=(int)rand(); printf(%10d,data[i]); } printf(\n); } m=log2(SumID); //调用函数:求以2为底的SumID的对数 /* 从根处理器将数据序列广播到其他处理器*/ /*{1表示传送的输入缓冲中的元素的个数, */ /* MPI_INT表示输入元素的类型, */ /* 0表示root processor的ID } */ MPI_Bcast(DataSize,1,MPI_INT,0,MPI_COMM_WORLD); /*ID号为0的处理器调度执行排序*/ para_QuickSort(data,0,DataSize-1,m,0,MyID); /*ID号为0的处理器打印排序完的有序序列*/ if(MyID==0) { printf(实际使用到的处理器个数：%d个\n,exp2(m)); printf(快速排序后的序列：\n); for(i=0;iDataSize;i++) { printf(%10d,data[i]); } printf(\n); } MPI_Finalize(); //结束计算 } /* * 函数名: para_QuickSort * 功能：并行快速排序，对起止位置为start和end的序列，使用2的m次幂个处理器进行排序 * 输入：无序数组data[1,n]，使用的处理器个数2^m * 输出：有序数组data[1,n] */ para_QuickSort(int *data,int start,int end,int m,int id,int MyID) { int i, j; int r; int MyLength; int *tmp; MPI_Status status; MyLength=-1; /*如果可供选择的处理器只有一个，那么由处理器id调用串行排序，对应于算法步骤(1.1)*/ /*(1.1) Pid call quicksort(data,i,j) */ if(m==0) { if(MyID==id) QuickSort(data,start,end); return; } /*由第id号处理器划分数据，并将后一部分数据发送到处理器id+exp2(m-1)，对应于算法步骤(1.2,1.3)*/ /*(1.2) Pid: r=patrition(data,i,j)*/ if(MyID==id) { /*将当前的无序区R[1，n]划分成左右两个无序的子区R[1，i-1]和R[i，n](1≤i≤n)*/ r=partition(data,start,end); MyLength=end-r;