MPI学习笔记(四)——组合通信子的管理

组和通信子的管理

组(Group)

组是一个进程的有序集合,在实现中可以看作是进程标识符的一个有序集。组内的每个进程与一个整数rank相联系,序列号从0开始并且是连续的。我们可以在通信组中使用组来描述通信空间中的参与者并对这些参与者进行分级(这样在通信空间中为它们赋予了唯一的名字)。

创建

//建立一个通信组对应的新进程组
int MPI_Comm_group(MPI_Comm comm, MPI_Group *group)

//查询调用进程在进程组的 rank
int MPI_Group_rank(MPI_Group group, int *rank)

MPI 是无法凭空构造一个组的,只能从其他以前定义的组中构造.

我们可以在初始通信子 MPI_COMM_WORLD的基础上定义组.

int MPI_Group_incl(MPI_Group old_group, int count, int *members, MPI_Group *new_group)
  
/*
MPI_Group old_group: 旧进程组
int count: members 数组中元素个数
int *members: 旧进程组中需要放入新进程组的进程的编号
MPI_Group *new_group: 新进程组
*/
创建时添加(include)

基于已经存在的进程组创建一个新的组, 并指明被包含(include)其中的成员进程.

int main(int argc, char **argv)
{
   
	int myid, numprocs, odd_rank;
	MPI_Group group_world, odd_group;
	int i;
	int members[10];
		
	MPI_Init(&argc, &argv);
	
	MPI_Comm_rank(MPI_COMM_WORLD, &myid);
    MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
	
	MPI_Comm_group(MPI_COMM_WORLD, &group_world);

	for(i=0; i<numprocs/2; i++) {
   
		members[i] = 2*i+1 ;
	}
	
	MPI_Group_incl(group_world, numprocs/2, members, &odd_group);
	
	MPI_Group_rank(odd_group, &odd_rank);
	
	printf("In process %d: odd rank is %d\n", myid, odd_rank);

	MPI_Finalize();
	return 0;
} 
/* OUTPUT */
/*
In process 2: odd rank is -32766
In process 3: odd rank is 1
In process 0: odd rank is -32766
In process 1: odd rank is 0
*/
解释
  1. 首先在一个进程中,获取了这个进程在初始通信子中的 id,以及在初始通信子中的进程数量.(myid, numprocs)

  2. 调用 MPI_Comm_group 函数, 建立新的进程组 group_world.

  3. 后面的 for 循环的意思是:
    把所有的奇数进程的 rank 存在 member 数组中.

  4. 然后调用 MPI_Group_incl 把刚刚保存的奇数进程放入新进程组(odd_group).

  5. MPI_Group_rank 获取进程在新进程组(odd_group)中的 id.

  6. 然后打印最开始的 myid 和在 odd_group 中的 rank.

  7. 在输出中为什么 2 和 0 的 rank 为-32766?
    因为这个进程不属于 odd_group 所以输出这个值.

创建时排除(exclude)
int MPI_Group_excl(MPI_group old_group, int count
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
MPI(Message Passing Interface)是一种用于并行计算的编程模型和库。归并排序是一种经典的排序算法,适合并行计算。 在MPI中,可以通过发送和接收消息来实现进程间的通信。下面是一个基于MPI的归并排序的伪代码: ```python def parallel_merge_sort(data): # 获取进程总数和当前进程编号 size = MPI.COMM_WORLD.Get_size() rank = MPI.COMM_WORLD.Get_rank() # 计算每个进程要处理的数据量 chunk_size = len(data) // size remainder = len(data) % size # 将数据分发到各个进程 if rank == 0: for i in range(size): if i < remainder: chunk = data[i * (chunk_size + 1):(i + 1) * (chunk_size + 1)] else: chunk = data[remainder + i * chunk_size:remainder + (i + 1) * chunk_size] MPI.COMM_WORLD.send(chunk, dest=i, tag=0) # 接收数据 chunk = MPI.COMM_WORLD.recv(source=0, tag=0) # 对本地数据进行排序 chunk.sort() # 归并排序 for step in range(size): # 计算要交换数据的进程编号 partner = (rank + step) % size # 发送和接收数据 sendbuf = chunk recvbuf = MPI.COMM_WORLD.recv(source=partner, tag=step) if rank < partner: sendtag = step recvtag = step + size else: sendtag = step + size recvtag = step MPI.COMM_WORLD.send(sendbuf, dest=partner, tag=sendtag) chunk = merge(chunk, recvbuf) # 将排序好的数据返回 if rank == 0: result = [] for i in range(size): chunk = MPI.COMM_WORLD.recv(source=i, tag=size) result.extend(chunk) return result else: MPI.COMM_WORLD.send(chunk, dest=0, tag=size) ``` 在这个算法中,首先将原始数据分发到各个进程,然后每个进程对本地数据进行排序,接着对每个步骤进行归并排序,并且使用MPI的send和recv函数进行交换数据。最后将排序好的数据返回到主进程。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值