MPI：散射后的值已损坏

Question

我只想用进程 0 创建一个具有随机值的矩阵。 创建它后，我将值分散到不同的过程中。 我有问题，我收到损坏的值，例如

1980714952

32716

-672521120

-672521856

21858

-310078720

21858

0

0

如何分散进程0创建的矩阵？

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

#define N 3

int main(int argc, char** argv)
{
    
    int rank;
    int size;
    
    int aa[N];
    int a[N][N];
    
    MPI_Init (&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    
    if (rank = 0)
    {
        srand(time(0));
        for (int i = 0; i < N; i++)
        {
            for (int j = 0; j < N; j++)
            {

                    a[i][j] = (rand() % 5);
            }
        }
        

    }
    
    MPI_Scatter(a, N*N/size, MPI_INT, aa, N*N/size, MPI_INT,0,MPI_COMM_WORLD);

    for (int x = 0; x < N; x++)
    {

            printf("%d\n", aa[x]);  
    }
        
        
    

    MPI_Finalize();


}


float maxValue(float* array, int _size)
{
    float v = 0.0;
    
    for (int i = 0; i < _size; i++)
    {
        if (array[i] >= v)
        {
            v = array[i];
        }   
    }
    
    return v;
} 

Answer 1

这是因为N*N/size 定义。您将N 定义为3，如果进程数超过3，假设进程数为4，那么N*N/Size 将是9/4，因此值将是2。所以2 元素只会被rank 0 分散到所有进程中。结果，最后，所有进程都将打印N（此处为3）值，但2 值仅从等级0 接收，第三个将是垃圾值。同样，此行为发生在其他数量大于N 的进程上。

您可以通过 a[i][j] = (rand() % 5); 将其更改为 a[i][j] = i; 并从进程数为 3 开始进行测试。它将按预期工作。如果您增加进程的大小，您会看到垃圾值将与输出一起打印。

if (rank = 0) 也应该是(rank == 0)。

回答评论中的问题：

要将N (x,y) 值发送到P 进程，更改将是。

    #define N 9 // this value must be divisible by P number of processes, because scatter is used for distributing uniform size. for 9 processes, the N can be any multiple of 9
    int aa[2*N/size]; // the receiving process only needs N/size * 2 values
    int a[N][2]; // since 2 values for each N
    MPI_Scatter(a, N*2/size, MPI_INT, aa, N*2/size, MPI_INT,0,MPI_COMM_WORLD);
    for (int x = 0; x < N*2/size; x++) // each processes print its local values
      {
        printf("rank %d, %d\n", rank,aa[x]);
      }

简而言之：要使用scatter 在P 进程之间平均分配一对N (x,y) 值，N 应该可以被P 整除。否则你必须使用scatterv