我正在尝试使用C中的MPI在一个“圆”中将一个进程(列)中的数组发送到另一个进程中。我可以使用单个整数使它正常工作,但是当我尝试执行以下操作时遇到分段错误通过发送数组。
{
int rank;
int numRanks;
MPI_Init(&argc, &argv);
MPI_Comm_size(MPI_COMM_WORLD, &numRanks);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
int arraysize = 100;
int *array;
if (rank == 0) {
array = malloc(sizeof(int)*arrasize);
for (int i = 0; i < arraysize; i++) {
array[i] = 1;
}
double starttime = MPI_Wtime();
MPI_Send(&array, arraysize, MPI_INT, (rank+1), 0, MPI_COMM_WORLD);
MPI_Recv(&array, arraysize, MPI_INT, numRanks-1, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
double endtime = MPI_Wtime();
printf("Rank %d got array from Rank %d in %lf", rank, numRanks-1, endtime-starttime);
} else if (rank == (numRanks-1)) {
MPI_Recv(&array, arraysize, MPI_INT, rank-1, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
printf("Rank %d got array from Rank %d", rank, rank-1);
MPI_Send(&array, arraysize, MPI_INT, 0, 0, MPI_COMM_WORLD);
} else {
MPI_Recv(&array, arraysize, MPI_INT, rank-1, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
printf("Rank %d got array from Rank %d", rank, rank-1);
MPI_Send(&array, arraysize, MPI_INT, (rank+1), 0, MPI_COMM_WORLD);
}
free(array);
MPI_Finalize();
return 0;
}
这是错误消息(为简单起见,我现在仅使用三个等级):
Rank 1 got array from Rank 0
[winston:42137] *** Process received signal ***
[winston:42137] Signal: Segmentation fault (11)
[winston:42137] Signal code: Address not mapped (1)
[winston:42137] Failing at address: 0x7add28
,然后尝试发送和接收数组的每个等级更多
任何人都能提供的任何帮助将不胜感激!
答案 0 :(得分:1)
您的array
并未指向除第0以外的每个进程(等级)中的有效内存地址。
尝试在随机(即无效)地址(即调用MPI_Recv(&array, arraysize...
)中写入内容时,这会使您的程序崩溃。
(感谢Gilles Gouaillardet在注释中指出),并且您应该将array
(而不是&array
)传递给MPI函数,因为它本身就是指向分配给您存储的地址的指针< em> some int
。
因此,解决方案将是:
array = malloc(sizeof(int)*arrasize);
。MPI_Recv(&array
(和MPI_Send(&array
)更改为MPI_Recv(array
(和MPI_Send(array
。