我试着做一个简单的排序,其中每个非零进程都将一个文件(filename=proc)读入一个缓冲区。zero proc收集所有这些buf,对其进行排序,然后打印出来。但是,在下面的代码中,proc 0收集proc1的缓冲区,而不是proc 2的缓冲区。有什么建议吗?
我用mpirun-np 3 a执行它
我的输入文件是
文件名:“1”
四十
十
一百
文件名:“2”
九十
二十
二十五
代码是:
#include <stdio.h>
#include <stdlib.h>
#include <mpi.h>
#define DEBUG 1
//#undef DEBUG
int compare (const void * a, const void * b)
{
return ( *(int*)a - *(int*)b );
}
int main (int argc, char *argv[])
{
int values[3];
int recv[3];
int n, i=0, temp;
FILE *in1, *in2;
int rank, size;
MPI_Init(&argc, &argv);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
char filename[20];
if(rank!=0){
// read files with filename"<proc#>" into the buffer
sprintf(filename, "%d", rank);
in1=fopen(filename,"r");
while(fscanf(in1,"%d",&values[i]) != EOF){
printf("rank %d Read data %d\n", rank,values[i]);
i++;
}
}
// gather values from all procs.
MPI_Gather(values,i,MPI_INT,recv,i,MPI_INT,0,MPI_COMM_WORLD);
printf("Gather done!");
if(rank==0){
// sort
qsort (recv, 6, sizeof(int), compare);
// print results
for (n=0; n<6; n++)
printf ("%d ",recv[n]);
printf("\n");
}
if(rank!=0)
fclose(in1);
MPI_Finalize();
return 0;
}
最佳答案
int recv[3];
接收缓冲区的大小应为要在收集节点上收集的来自所有处理器的元素总数。
所以recv[9]应该在这里工作。
查看此示例以进行收集
example
关于c - MPI收集和qsort,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/7857574/