| 查看: 657 | 回复: 2 | |||
[求助]
初学MPI，以下Fortran程序请教
|
|
以下程序编译后运行： mpiexec -n 4 test.exe 出现如下结果： myid= 2 results( 3 )= 10 myid= 1 results( 2 )= 10 myid= 3 results( 4 )= 10 myid= 0 results( 1 )= 10 results( 1 )= 10 results( 2 )= -842150451 results( 3 )= -842150451 results( 4 )= -842150451 期望所有输出结果均为10，请问哪里错了？ 谢谢！ 附：Fortran程序： program test implicit none include 'mpif.h' integer::i integer rank, size, myid, ierror,commandsize integer,allocatable::results(  commandsize=4 allocate(results(commandsize)) call MPI_INIT(ierror) call MPI_COMM_SIZE(MPI_COMM_WORLD, size, ierror) call MPI_COMM_RANK(MPI_COMM_WORLD, myid, ierror) write(*,*)'myid=',myid results(myid+1)=10 write(*,*)'results(',myid+1,')=',results(myid+1) call MPI_ALLGATHER(results,commandsize,MPI_INTEGER,results,commandsize,MPI_INTEGER,MPI_COMM_WORLD,ierror) if(myid==0)then do i=1,commandsize write(*,*)'results(',i,')=',results(i) enddo endif call MPI_FINALIZE(ierror) deallocate(results) end |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 |
» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有185人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran和c混合编译问题？
已经有9人回复
- ABAQUS中的用户子程序fortran编写问题
已经有4人回复
- 求教，我没安装FORTRAN，能不能用其它简单的方法打开后缀为.f的FORTRAN源程序呢？
已经有13人回复
- 急求fortran运行错误原因，在线等
已经有7人回复
- 请教 fortran 运行错误的原因
已经有13人回复
- 配置一台1万左右的计算机跑Fortran程序以及高斯的计算
已经有8人回复
- 关于Fortran程序运行的求助
已经有2人回复
- fortran程序怎样调用lapack库进行矩阵对角化？
已经有7人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 关于MPI并行编程的基本问题
已经有8人回复
- 请教关于in cell形式向original形式坐标转化的fortran算法
已经有3人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】初学fortran，
已经有15人回复
- 请教fortran导txt数据进origin问题
已经有1人回复
baghnac
金虫 (小有名气)
- 应助: 14 (小学生)
- 金币: 1706.4
- 红花: 6
- 帖子: 177
- 在线: 44小时
- 虫号: 1365578
- 注册: 2011-08-11
- 性别: GG
- 专业: 宇宙学
【答案】应助回帖
tang76(金币+1): 2011-10-26 00:08:22
|
⑴发送与接收缓冲区如果一样，会冲突 ⑵接收缓冲区应大于各发送缓冲区的总大小 虽然书上没有找到这么说，但是我试出来的是这个结果 结果如下 myid=0 results( 1 )= 10 results( 2 )= 0 results( 3 )= 0 results( 4 )= 0 recvs( 1 )= 10 recvs( 2 )= 0 recvs( 3 )= 0 recvs( 4 )= 0 recvs( 5 )= 0 recvs( 6 )= 10 recvs( 7 )= 0 recvs( 8 )= 0 recvs( 9 )= 0 recvs( 10 )= 0 recvs( 11 )= 10 recvs( 12 )= 0 recvs( 13 )= 0 recvs( 14 )= 0 recvs( 15 )= 0 recvs( 16 )= 10 myid=1 results( 1 )= 0 results( 2 )= 10 results( 3 )= 0 results( 4 )= 0 myid=2 results( 1 )= 0 results( 2 )= 0 results( 3 )= 10 results( 4 )= 0 myid=3 results( 1 )= 0 results( 2 )= 0 results( 3 )= 0 results( 4 )= 10 附程序： program test implicit none include 'mpif.h' integer::i integer rank, size, myid, ierror,commandsize,sendcnt,recvcnt,sendssize integer,allocatable::results( integer,allocatable::recvs( commandsize=4 sendssize=16 allocate(results(commandsize)) allocate(recvs(sendssize)) call MPI_INIT(ierror) call MPI_COMM_SIZE(MPI_COMM_WORLD, size, ierror) call MPI_COMM_RANK(MPI_COMM_WORLD, myid, ierror) ! write(*,*)'myid=',myid results(myid+1)=10 ! write(*,*)'results(',myid+1,')=',results(myid+1) call MPI_ALLGATHER(results,commandsize,MPI_INTEGER,recvs,commandsize,MPI_INTEGER,MPI_COMM_WORLD,ierror) if(myid==0)then do i=1,commandsize write(*,*)'results(',i,')=',results(i) enddo do i=1,sendssize write(*,*)'recvs(',i,')=',recvs(i) enddo endif call MPI_FINALIZE(ierror) deallocate(results) end  运行结果 |
2楼2011-10-25 14:29:53
3楼2011-10-26 00:09:47
5
投票:
