如何在MPI中发送具有每个处理器变体的2Darrays

我试图在根0上采用随机生成的数组，稍微随机地改变它，并将每个变体发送到另一个处理器。 这是我到目前为止的代码：

#include "stdio.h" #include "stdlib.h" #include "mpi.h" #include "math.h" int main(int argc, char **argv) { int N = 32; int dim = 3; float a = 10.0; int size, rank, i, j, k, q; float **C; float rijx, rijy, rijz, rij, Vij, E=0; float stepsize = 0.05; double Start_time, End_time, Elapse_time; MPI_Status status; MPI_Init(&argc, &argv); MPI_Comm_size(MPI_COMM_WORLD, &size); MPI_Comm_rank(MPI_COMM_WORLD, &rank); C = (float **)malloc(N * sizeof(float*)); // 32 particles for (i = 0; i < N; i++) { C[i]=(float *)malloc(dim*sizeof(float)); // x, y, z } MPI_Barrier(MPI_COMM_WORLD); if(rank == 0) { Start_time = MPI_Wtime(); } if (rank == 0) { for(i = 0; i < N; i++) { for(j = 0; j < dim; j++) { C[i][j] = (float)rand()/(float)(RAND_MAX/a); //printf("%f\n",C[i][j]); } } } // this is to generate some slight variations in the array float** randomsteps(float **matrix) { for(i = 0; i < N; i = i+(rand()%(32/size))) { for (j = 0; j  10) { C[i][j] = C[i][j] - 10; } } else { C[i][j] = C[i][j]-stepsize; if(C[i][j] < 0) { C[i][j] = C[i][j] + 10; } } } } return C; } // and here I try to send the array if(rank == 0) { for(i=0; i<size; i++) { C = randomsteps(C); MPI_Send(&C, N*3, MPI_FLOAT, i, 10+i, MPI_COMM_WORLD); } } if(rank != 0) { for(i=0; i<size; i++) { MPI_Recv(&C, N*3, MPI_FLOAT, 0, 10+i, MPI_COMM_WORLD, &status); } } MPI_Barrier(MPI_COMM_WORLD); MPI_Finalize(); return 0; } 

代码的一个明显问题是随机数的生成方式有些幼稚（每次运行程序时都会给出相同的值）。 这是我以后可以做的事情。

现在，我只是想知道 – 我发送和接收arrays的方式有什么问题？ 在使用MPI发送和接收时，如何最好地格式化数据，我遇到了很多麻烦。 我将如何修复这部分代码？

在此先感谢您的帮助！

Answer 1

这里的问题是所有MPI调用都期望内存是连续的。 你的记忆在给定的行中只是连续的，而你所指的2D数组实际上是一个指针数组。 指针不可移植，因此尝试向另一个进程发送或广播指针数组毫无意义，而且MPI本身不支持深层复制，因此这种方法不起作用。

但是，如果将arrays分配更改为以下内容：

  float** C; float* C_buff; C = (float**)malloc(N * sizeof(float*)); // 32 particles C_buff = (float*)malloc(N * dim * sizeof(float)); // buffer for particles float* p = &C_buff[0]; for (i = 0; i < N; i++) { C[i]=p; p+= dim*sizeof(float)); } 

[免责声明：用浏览器编写，完全未经测试，自担风险使用]

因此C_buff表示2D数组的连续内存， C包含C_buff连续分配内存的行指针，然后您可以使用现有代码进行初始化，但是然后执行以下操作：

 MPI_Send(&C_buff[0][0], N*DIM, MPI_FLOAT, i, 10+i, MPI_COMM_WORLD); 

即。 使用C_buff进行MPI调用，它应该可以工作。