MPI将大量数据发送给自我

Question

我尝试使用MPI发送并接收一个大型阵列。该代码适用于小阵列（阵列尺寸＆lt; 10000）时，当我将阵列大小进一步增加到100000时，它将死去。 Here is my code:

#include <stdio.h>
#include <stdlib.h>
#include <mpi.h>

int main(int argc, char *argv[])
{

MPI_Init(&argc, &argv);

int size;
MPI_Comm_size(MPI_COMM_WORLD, &size);

int rank;
MPI_Comm_rank(MPI_COMM_WORLD, &rank);

MPI_Status stat;


    int N = 100000;

    // Allocate memory for A on CPU
    double *A = (double*)malloc(N*sizeof(double));
    double *B = (double*)malloc(N*sizeof(double));

    // Initialize all elements of A to 0.0
    for(int i=0; i<N; i++){
        A[i] = 0.0;
    }

    int tag1 = 10;
    int tag2 = 20;

//printf("rank=%d",rank);


    MPI_Send(A, N, MPI_DOUBLE, 0, tag1, MPI_COMM_WORLD);
    
    MPI_Recv(B, N, MPI_DOUBLE, 0, tag1, MPI_COMM_WORLD, &stat);

MPI_Finalize();

return 0;
}

I compile the code by: mpicc
运行的命令：Mpirun -NP 1 ./

Answer 1

这就是发送和接收的定义：它们正在阻止。在发送发送后，发送后的声明将不会执行。

1. 解决此问题的安全方法是使用mpi_isend和mpi_irecv。
2. 一个不阻止的小消息的情况是系统依赖的优化，您无法指望它。