我正在尝试使用MPI形成一个环通信,我的每个进程都将其结果发送到下一个进程,最后一个进程将结果发送到第0个进程。让我们假设我有4个进程然后我的第0个进程将结果发送到1st,1st到2nd,2nd到3rd和3rd到。
#include "mpi.h"
#include <stdio.h>
#include<stdlib.h>
#define NELEM 1000
int main (int argc, char *argv[])
{
int numtasks, rank, rc, i, dest = 1, tag = 111, source = 0, size;
double *data, result;
void *buffer;
data=(double*)malloc(sizeof(double)*NELEM);
if(data==NULL)
{
printf("Unable to allocate memory\n");
return;
}
MPI_Status status;
MPI_Init (&argc, &argv);
MPI_Comm_size (MPI_COMM_WORLD, &numtasks);
MPI_Comm_rank (MPI_COMM_WORLD, &rank);
for (i = 0; i < NELEM; i++)
data[i] = (double) random ();
if (rank == 0)
source=numtasks-1;
else
source=rank-1;
if(rank==numtasks-1)
dest=0;
else
dest=rank+1;
printf("Rank %d sending data to rank %d\n",rank,dest);
MPI_Send(data, NELEM, MPI_DOUBLE, dest, tag,MPI_COMM_WORLD);
printf("Rank %d send complete\n",rank);
printf("Rank %d receiving data from rank %d\n",rank,source);
MPI_Recv (data, NELEM, MPI_DOUBLE, source, tag, MPI_COMM_WORLD,&status);
printf("Rank %d received data from rank %d\n",rank,source);
MPI_Finalize ();
}
这里NELEM是发送或接收的元素数量。 如果我发送少于100个元素与4个线程上面代码将工作正常,但如果我增加线程数将被阻止。我不明白为什么它会被阻止。我们可以发送的数据大小是否有任何限制。
由于
的Ajay
答案 0 :(得分:3)