MPI错误 - 发生了什么？

Question

这可能是一个非常愚蠢的问题，但我看不到我的错误在此代码...输出是错误的，1个打印：

3.000000 - 3.000000 - 3.000000 - 3.000000 - 3.000000 -
3.000000 - 3.000000 - 3.000000 - 3.000000 - 3.000000 -
3.000000 - 3.000000 - 3.000000 - 3.000000 - 3.000000 -
3.000000 - 3.000000 - 3.000000 - 3.000000 - 3.000000 -
3.000000 - 0.000000 - 0.000000 - 0.000000 - 0.000000 -

我需要用MPI_Recv来代替Bcast ... B发生了什么事？ = /这是我的malloc？或者我不能使用MPI_Recv发送整个矩阵？为什么整个数组不会进入另一个进程？

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

int main(int argc, char **argv){ 
    int rank, size; 
    int lines, cols; 
    int i, j; 

    MPI_Init(&argc, &argv); 
    MPI_Comm_rank(MPI_COMM_WORLD, &rank); 
    MPI_Comm_size(MPI_COMM_WORLD, &size); 

    MPI_Status status; 

    lines = 5; 
    cols = 5; 

    if(rank == 0){ 
     double** array = (double**) malloc(lines*sizeof(double*)); 
     for(i=0; i<lines; i++) 
      array[i] = (double*) malloc(cols*sizeof(double)); 

     for(i=0; i<lines; i++) 
      for(j=0; j<cols; j++) 
        array[i][j] = 3; 

     for(i=0; i<lines; i++){ 
      for(j=0; j<cols; j++) 
        printf("%f - ", array[i][j]); 
      printf("\n"); 
     } 

     MPI_Send(&array[0][0], lines*cols, MPI_DOUBLE, 1, 0, MPI_COMM_WORLD); 
    } 
    else{ 
     double** arrayNew = (double**)malloc(lines*sizeof(double*)); 
     for (i=0; i<lines; i++) 
       arrayNew[i] = (double*) malloc(cols*sizeof(double)); 

     MPI_Recv(&arrayNew[0][0], lines*cols, MPI_DOUBLE, 0, 0, MPI_COMM_WORLD, &status); 

     for(i=0; i<lines; i++){ 
       for(j=0; j<cols; j++) 
         printf("%f - ", arrayNew[i][j]); 
       printf("\n"); 
     } 
    } 
    MPI_Finalize(); 
} 

Answer 1

这一直出现，问题是C和多维数组。

的问题是，这条线

MPI_Send(&array[0][0], lines*cols, MPI_DOUBLE, 1, 0, MPI_COMM_WORLD); 

告诉MPI发送lines*cols双打从&(array[0][0])位置开始，同样这条线

MPI_Recv(&arrayNew[0][0], lines*cols, MPI_DOUBLE, 0, 0, MPI_COMM_WORLD, &status); 

告诉MPI接收lines*cols双打到内存中的起始从位置&(array[0][0])。但是这种配置：

double** arrayNew = (double**)malloc(lines*sizeof(double*)); 
    for (i=0; i<lines; i++) 
      arrayNew[i] = (double*) malloc(cols*sizeof(double)); 

不会使lines*cols的连续数组变成双数;它使lines阵列的cols双打，并且这些行可能分散在整个内存。你需要做的是这样的：

double** array = (double**) malloc(lines*sizeof(double*)); 
    array[0] = (double *)malloc(lines*cols*sizeof(double)); 
    for(i=1; i<lines; i++) 
     array[i] = &(array[0][i*cols]); 

    /* ... */ 

    free(array[0]); 
    free(array); 

分配和空闲内存连续lines*cols块，你可以发送和接收到。