我是 MPI 的新手,我想为 Residence
struct
创建一个新的数据类型。我只是想看看我是否可以正确地创建新类型。
结构住宅
{
双 x;
双y;
};
我的新 MPI 类型
MPI_Datatype createRecType()
{
// Set-up the arguments for the type constructor
MPI_Datatype new_type;
int count = 2;
int blocklens[] = { 1,1 };
MPI_Aint indices[2];
//indices[0]=0;
MPI_Type_extent( MPI_DOUBLE, &indices[0] );
MPI_Type_extent( MPI_DOUBLE, &indices[1] );
MPI_Datatype old_types[] = {MPI_DOUBLE,MPI_DOUBLE};
MPI_Type_struct(count,blocklens,indices,old_types,&new_type);
MPI_Type_commit(&new_type);
}
最佳答案
除了 indices
应该从结构的 开始 以字节为单位给出每个结构字段的偏移量之外,您几乎已经猜对了。构造这种类型的正确方法是使用 offsetof
运算符,在 stddef.h
中定义:
#include <stddef.h> // or <cstddef> for C++
struct Residence
{
double x;
double y;
};
MPI_Datatype createRecType()
{
// Set-up the arguments for the type constructor
MPI_Datatype new_type;
int count = 2;
int blocklens[] = { 1,1 };
MPI_Aint indices[2];
indices[0] = (MPI_Aint)offsetof(struct Residence, x);
indices[1] = (MPI_Aint)offsetof(struct Residence, y);
MPI_Datatype old_types[] = {MPI_DOUBLE,MPI_DOUBLE};
MPI_Type_struct(count,blocklens,indices,old_types,&new_type);
MPI_Type_commit(&new_type);
return new_type;
}
虽然这对于该特定结构就足够了,但通常必须调整结构化类型的长度,以考虑编译器可能在结构末尾插入的任何尾随填充。仅当一个人想要发送该结构类型的多个项目时才需要,即结构元素数组。旧的方法是将第三个成员添加到
MPI_UB
类型的结构中(UB 来自 Upper Bound)并将该成员的偏移量设置为等于 sizeof(struct Residence)
(填充计入 sizeof
返回的结构大小) .现代方法是使用 MPI_Type_create_resized
,它创建一个新的 MPI 类型,其类型签名与原始类型相同,但范围不同: MPI_Type_struct(count,blocklens,indices,old_types,&new_type);
// Create a resized type
MPI_Type resized_new_type;
MPI_Type_create_resized(new_type,
// lower bound == min(indices) == indices[0]
indices[0],
(MPI_Aint)sizeof(struct Residence),
&resized_new_type);
MPI_Type_commit(&resized_new_type);
// Free new_type as it is no longer needed
MPI_Type_free(&new_type);
return resized_new_type;
仅显示相关代码行。上面的代码假设
indices[0]
给出了第一个结构元素的偏移量。可以改用 MPI_Type_get_extent
来获得真正的下限,这适用于具有负偏移量的结构类型。没有必要提交 new_type
因为它仅用于构造调整大小的类型。也没有必要保留它,这就是为什么在创建 resized_new_type
后将其释放。关于c++ - 如何在 MPI 中创建新类型,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/20040663/