| 查看: 602 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 |
» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有249人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有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人回复
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
