| 查看: 763 | 回复: 5 | ||
q:1306140890木虫 (正式写手)
|
[求助]
关于并行分区中消息传递的问题？ 已有1人参与
|
|
subroutine read_b(ita) integer ny0,nz0,ii,ita if(myid==0)then allocate(du0(0:nz,1:3)) open(3,file='./inflow/du.dat',form='unformatted') read(3) read(3)((du0(k,l),k=0,nz),l=1,3) close(3) do ii=0,numprocs-1 do l=1,3 do k=0,mysizez(ii)-1 buffer(k+(l-1)*mysizez(ii))=real(du0(ibeginz(ii)+k,l)) buffer(k+(l+2)*mysizez(ii))=imag(du0(ibeginz(ii)+k,l)) enddo enddo call MPI_BSEND(buffer(0),3*mysizez(ii)*2,MPI_REAL8,ii,ii,MPI_COMM_WORLD,ierr) enddo deallocate(du0) endif call MPI_RECV(buffer(0),3*(kn+1)*2,MPI_REAL8,0,MPI_ANY_TAG,MPI_COMM_WORLD,status,ierr) do l=1,3 do k=0,kn du(k,l)=buffer(k+(l-1)*(kn+1))+ai*buffer(k+(l+2)*(kn+1)) enddo enddo endsubroutine 1）这个子程序哪儿有问题？ 2）各个进程在receive之后是如何把自己分区内的数据合成du的？接收的buffer长度和发送的不一样这不矛盾吗？（本程序是基于别的程序改的，原程序也是这么写的） 望不吝赐教，谢谢！ |
回复此楼» 猜你喜欢
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
2楼2016-08-31 10:45:14
q:1306140890
木虫 (正式写手)
- 应助: 32 (小学生)
- 金币: 3750.2
- 红花: 8
- 帖子: 648
- 在线: 290.8小时
- 虫号: 2302672
- 注册: 2013-02-27
- 专业: 流体力学
3楼2016-08-31 11:42:59
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
4楼2016-08-31 18:01:43
q:1306140890
木虫 (正式写手)
- 应助: 32 (小学生)
- 金币: 3750.2
- 红花: 8
- 帖子: 648
- 在线: 290.8小时
- 虫号: 2302672
- 注册: 2013-02-27
- 专业: 流体力学
5楼2016-08-31 19:21:47
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
6楼2016-08-31 20:13:22
10