| 查看: 481 | 回复: 1 | |||
| 【悬赏金币】回答本帖问题，作者fo3mt2将赠送您 10 个金币 | |||
fo3mt2金虫 (小有名气)
|
[求助]
MPI并行收集数据
|
||
|
大家好， 我想请教一段用fortran写的mpi并行程序。代码如下： do i = 1, nomi, 1 if (mod(i-1,world_size) /= world_rank) cycle kkim = kkimp(i,1:3) call proce(i,ndegen,nbc,npo,kkim,hil,thil,tecv,teve,bcb,eni,er,pi,nne,hb,bcp1) if (world_rank == 0) then bcp2 = bcp1 ！the 6th line do j = 1, world_size-1, 1 if (i-1+j == nomi) exit k = i + j + 1000 call mpi_recv(bcp3,nomi*nne*6,mpi_double_complex,j,k,mpi_comm_world,mpi_status_ignore,ierr) bcp2 = bcp2 + bcp3 end do else k = i + 1000 call mpi_send(bcp1,nomi*nne*6,mpi_double_complex,0,k,mpi_comm_world,ierr) end if end do 假设根进程（root processor）在外循环里一共被使用了两次，每次都会调用子程序proce，并从子程序得到传回的数组bcp1。 当用根进程（root processor）使用数组bcp2，从各个进程（包括其自身）收集各自得到的数组bcp1时，根进程（root processor）需要先收集自己得到的数组bcp1。 因为假设根进程（root processor）在外循环里一共被使用了两次，根进程（root processor）一共要从自身收集两次bcp1数组。那程序里的第6行（标记为！the 6th line），写成bcp2 = bcp1是否合理呢？ 我觉得这样写，会使得根进程（root processor）在循环里第二次被使用得到的数组bcp1，覆盖掉第一次被使用得到的数组bcp1，进而导致丢失了第一次使用时得到的数组bcp1信息。 如果不丢失信息，应该进行叠加，也就是把第6行写成bcp2 = bcp2 + bcp1；但程序里却不是这样写的。 我想不通是什么原因，能否请大家给些建议呢？或者是否是程序写错了呢？ 谢谢啦。 盼复。 |
回复此楼» 猜你喜欢
海南大学田新龙教授课题组招24级博士生或科研助理 (氢燃料电池与器件，重点支持方向)
已经有6人回复
-
谁要做个人网站，可以交流啊！
已经有28人回复
-
申请基金代表性成果
已经有11人回复
-
推荐转投（ transfer pending）是否有用？
已经有6人回复
-
博三一直没文章怎么办
已经有23人回复
-
蹲一个男朋友
已经有26人回复
-
基金委也是打工的
已经有22人回复
-
河北省基金
已经有18人回复
-
换导师
已经有10人回复
-
中文投稿
已经有7人回复
taichiman
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 660.6
- 沙发: 5
- 帖子: 346
- 在线: 10.4小时
- 虫号: 32291817
- 注册: 2022-11-12
- 专业: 石油、天然气地质学
2楼2022-12-05 12:05:02
30