如果我想定义自己的类型,并将其用作一个数据类型,MPI_Send
从一个矩阵中只取偶数行,那么这个矩阵(发送缓冲区)是否必须静态分配?
当动态分配时,我似乎遇到了问题。这是因为地址需要连续才能发送数据吗?
最佳答案
不,用MPI_Send
发送的内存不需要静态分配。
要发送数组子集,可能需要使用MPI_Type_indexed
。下面是对mpi.deino.net article on MPI_Type_indexed
示例的一个稍微修改的版本,在这里我替换了静态分配的缓冲区
int buffer[27];
到动态分配的缓冲区
int* buffer = (int*)malloc(27 * sizeof(int));
我希望它有帮助:
#include <mpi.h>
#include <stdio.h>
int main(int argc, char *argv[])
{
int rank, size, i;
MPI_Datatype type, type2;
int blocklen[3] = { 2, 3, 1 };
int displacement[3] = { 0, 3, 8 };
int* buffer = (int*)malloc(27 * sizeof(int)); //int buffer[27];
MPI_Status status;
MPI_Init(&argc, &argv);
MPI_Comm_size(MPI_COMM_WORLD, &size);
if (size < 2)
{
printf("Please run with 2 processes.\n");
MPI_Finalize();
return 1;
}
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
MPI_Type_contiguous(3, MPI_INT, &type2);
MPI_Type_commit(&type2);
MPI_Type_indexed(3, blocklen, displacement, type2, &type);
MPI_Type_commit(&type);
if (rank == 0)
{
for (i=0; i<27; i++)
buffer[i] = i;
MPI_Send(buffer, 1, type, 1, 123, MPI_COMM_WORLD);
}
if (rank == 1)
{
for (i=0; i<27; i++)
buffer[i] = -1;
MPI_Recv(buffer, 1, type, 0, 123, MPI_COMM_WORLD, &status);
for (i=0; i<27; i++)
printf("buffer[%d] = %d\n", i, buffer[i]);
fflush(stdout);
}
MPI_Finalize();
free(buffer);
return 0;
}
关于c - MPI_Send仅适用于静态分配的缓冲区,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/23552968/