MPI_Scatter段错误（信号11）

Question

我想并行化一些计算矩阵的代码段，但我得到了段错误错误。我已经调试了代码，但我无法找出问题所在。

这里是片段，我宣布矩阵的计算方法：

double **A, **B, **C, *tmp1,*tmp2,*tmp3; 
int N; 
int myrank, P, from, to; 

N = 100; 

tmp1 = (double *) malloc (sizeof(double) * N * N); 
tmp2 = (double *) malloc (sizeof(double) * N * N); 
tmp3 = (double *) malloc (sizeof(double) * N * N); 
A = (double **) malloc (sizeof(double *) * N); 
B = (double **) malloc (sizeof(double *) * N); 
C = (double **) malloc (sizeof(double *) * N); 

from = myrank * N/P; 
to = (myrank+1) * N/P; 

然后这里是我的inicialize MPI

MPI_Init (&argc, &argv); 
MPI_Comm_rank(MPI_COMM_WORLD, &myrank); 
MPI_Comm_size(MPI_COMM_WORLD, &P); 

但不知何故，我不能撒...

MPI_Bcast(B, N*N, MPI_INT, 0, MPI_COMM_WORLD); 

printf("Rank: %d\n", myrank); 

MPI_Scatter(A, N*N/P, MPI_INT, &A[from], N*N/P, MPI_INT, 0, MPI_COMM_WORLD); 

我不能在这里写所有的代码，因为我不被允许...但问题是在自我散布。我做错了什么？

Answer 1

您不能直接发送2-D矩阵（double**）。相反，您必须传递实际数据的地址：&A[0][0]。对于接收，您还必须将地址传递到实际存储数据的位置。在你的情况下，这将是：&A[from][0]而不是&A[from]。在其他A比那些收到仍将无效

MPI_Scatter(&A[0][0], N*N/P, MPI_INT, &A[from][0], N*N/P, MPI_INT, 0, MPI_COMM_WORLD); 

在另一方面，在接收过程中的所有值：

尝试以此为散布。您可能需要考虑创建所需大小的单独接收缓冲区：

double* local_A = (double*)malloc(sizeof(double)*N*N/P); 
MPI_Scatter(&A[0][0], N*N/P, MPI_INT, local_A, N*N/P, MPI_INT, 0, MPI_COMM_WORLD);