Matrix 矩阵乘以MPI（分段故障）
我最近开始学习MPI，你可能已经猜到了，我已经遇到了一个我自己无法解决的错误 我想写一个将两个矩阵相乘的程序。然而，我还没有走到那一步，事实上，我一开始就坚持播放矩阵Matrix 矩阵乘以MPI（分段故障）,matrix,segmentation-fault,mpi,Matrix,Segmentation Fault,Mpi,我最近开始学习MPI，你可能已经猜到了，我已经遇到了一个我自己无法解决的错误 我想写一个将两个矩阵相乘的程序。然而，我还没有走到那一步，事实上，我一开始就坚持播放矩阵 #define MASTER 0 if (rank == MASTER) { A = (double *) malloc(N * N * sizeof(double)); B = (double *) malloc(N * N * sizeof(double)); matFillRa
#define MASTER 0
if (rank == MASTER) {
A = (double *) malloc(N * N * sizeof(double));
B = (double *) malloc(N * N * sizeof(double));
matFillRand(N, A);
matFillRand(N, B);
}
if (rank == MASTER) {
P = (double *) malloc(N * N * sizeof(double));
}
matMulMPI(N, A, B, P);
if (rank == MASTER) {
printMatrix(N, P);
}
void matMulMPI(long N, double *a, double *b, double *c) {
long i, j, k;
int rank, size;
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
MPI_Comm_size(MPI_COMM_WORLD, &size);
MPI_Bcast(&N, 1, MPI_LONG, MASTER, MPI_COMM_WORLD);
MPI_Bcast(b, N*N, MPI_DOUBLE, MASTER, MPI_COMM_WORLD);
printMatrix(N, b);
//TO-DO: Broadcast A
//TO-DO: Do Math
}
我已经弄明白了。所有进程（不仅仅是主进程）都需要首先分配内存 因此，丢失的行是
void matMulMPI(long N, double *a, double *b, double *c) {
...
MPI_Bcast(&N, 1, MPI_LONG, MASTER, MPI_COMM_WORLD);
b = (double *) malloc(N * N * sizeof(double));
MPI_Bcast(b, N*N, MPI_DOUBLE, MASTER, MPI_COMM_WORLD);
...
}
我已经弄明白了。所有进程（不仅仅是主进程）都需要首先分配内存 因此，丢失的行是
void matMulMPI(long N, double *a, double *b, double *c) {
...
MPI_Bcast(&N, 1, MPI_LONG, MASTER, MPI_COMM_WORLD);
b = (double *) malloc(N * N * sizeof(double));
MPI_Bcast(b, N*N, MPI_DOUBLE, MASTER, MPI_COMM_WORLD);
...
}
为了澄清为什么会出现这种情况：广播意味着广播后的每个进程都将在
MPI\u BcastMPI\u BcastMPI\u BcastMPI\u Bcast