| 查看: 585 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 |
» 猜你喜欢
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有256人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有28人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有22人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
【答案】应助回帖
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
