MPI_Send在大数据大小的环通信中阻塞

时间:2013-05-02 11:52:37

标签: parallel-processing mpi

我正在尝试使用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

1 个答案:

答案 0 :(得分:3)

您的所有进程都在尝试发送。但他们不能因为没有人愿意倾听。

对于较小的元素大小,我希望消息适合缓冲区。

正如Jonathan建议的那样,答案是使用MPI_Sendrecv()或非阻塞通信。