我强迫自己用C语言MPI中的个人结构来做广播基本上我已经创建了struct Vector3d。
typedef struct {
double x, y, z;
} Vector3d;
然后我一直在阅读并为每个MPI进程编写代码。
int rank, size;
MPI_Init(&argc, &argv);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
MPI_Comm_size(MPI_COMM_WORLD, &size);
// MPI Struct para Vector3d
int nroItems = 3;
int blockLengths[3] = { 1, 1, 1 };
MPI_Datatype types[3] = { MPI_DOUBLE, MPI_DOUBLE, MPI_DOUBLE };
MPI_Datatype MPI_Vector3d;
MPI_Aint offsets[3];
offsets[0] = offsetof(Vector3d, x);
offsets[1] = offsetof(Vector3d, y);
offsets[3] = offsetof(Vector3d, z);
MPI_Type_create_struct(nroItems, blockLengths, offsets, types, &MPI_Vector3d);
MPI_Type_commit(&MPI_Vector3d);
然后我用这个做一个矢量3D阵列的广播。
Vector3d * num = (Vector3d *) malloc(sizeof(Vector3d) * 10);
if(rank == 0) {
...
...
MPI_Bcast(num, 10, MPI_Vector3d, 0, MPI_COMM_WORLD);
} else {
MPI_Bcast(num, 10, MPI_Vector3d, 0, MPI_COMM_WORLD);
}
当我运行它的时候我得到了这个。
[mario-elementary:24020] *** Process received signal ***
[mario-elementary:24020] Signal: Segmentation fault (11)
[mario-elementary:24020] Signal code: Address not mapped (1)
[mario-elementary:24020] Failing at address: 0x56fae13e2cc8
[mario-elementary:24020] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0xfcb0) [0x2b7d6fce7cb0]
[mario-elementary:24020] [ 1] /lib/x86_64-linux-gnu/libc.so.6(+0x14ae90) [0x2b7d70040e90]
[mario-elementary:24020] [ 2] /usr/lib/libmpi.so.0(+0x3f812) [0x2b7d6fa67812]
...
...etc
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 24020 on node mario-elementary exited on signal 11 (Segmentation fault).
--------------------------------------------------------------------------
请帮忙!
最佳答案
你的偏移有误
offsets[3] = offsetof(Vector3d, z);
应改为
offset[2]
我想这就是解决BCAST问题的方法。另外,您的offsetof函数的结果是什么它是如何工作的?
关于c - MPI Bcast与MPI_Type_create_struct段错误,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/20776797/