MPI + FORTRAN 发送接收数组操作

在FORTRAN 中，以二维数组为例，arr(i, j)既可以表示第i行j列的元素值，也可以表示其地址，这样，我们就可以发送接收连续的序列。

下面的代码，是将0号进程中数组的第一列发送到1号进程中数组的第二列，并且输出数组的值。采用标准的阻塞式通信。

发送接收前的大部分输出操作，只是为了验证数组的存取细节，并无必要。

 

program maininclude "mpif.h"integer, parameter:: n = 3integer myid, numprocs, i, j, ierr, tmpinteger arr(n, n),status(mpi_status_size)call MPI_INIT(ierr)call MPI_COMM_RANK(MPI_COMM_WORLD, myid, ierr)call MPI_COMM_SIZE(MPI_COMM_WORLD, numprocs, ierr)!set values for arraytmp = 0do 10, i = 1, ndo 10, j = 1, narr(i, j) = tmptmp = tmp + 1 10continue!print the array, unnecessaryif(myid.eq.0) thendo 20, i = 1, n do 20, j = 1, n print *, arr(i, j)20 continueprint *, '#:'do i = 1, nprint *,arr(1,i)enddoprint *, '##:' do i = 1, n print *,arr(2,i)enddoprint *, '###:' do i = 1, n print *,arr(3,i)enddo!print *, arrendifif(myid.eq.0) then! send and viewcall mpi_send(arr(1,1), 3, mpi_integer, 1, 100, &mpi_comm_world, ierr)print *, 'sending...'print *, arrelse!viewcall mpi_recv(arr(1,2), 3, mpi_integer, 0, 100, &mpi_comm_world, status, ierr)print *, 'receiving...'print *, arrendifcall MPI_FINALIZE(ierr)stopend

 

 

 

结果：

[root@c001 zlt]# mpiexec -n 2 ./send 0 1 2 3 4 5 6 7 8 #: 0 1 2 ##: 3 4 5 ###: 6 7 8 sending... 0 3 6 1 4 7 2 5 8 receiving... 0 3 6 0 3 6 2 5 8[root@c001 zlt]#

 

 

 

 

三维数组，发送一个页面：

program main include "mpif.h" integer, parameter:: n = 3 integer myid, numprocs, i, j, k, ierr, tmp integer arr(n,n,n),status(mpi_status_size) call MPI_INIT(ierr) call MPI_COMM_RANK(MPI_COMM_WORLD, myid, ierr) call MPI_COMM_SIZE(MPI_COMM_WORLD, numprocs, ierr) tmp = 1 do 10, k = 1, n do 10, j = 1, n do 10, i = 1, n arr(i, j, k) = tmp tmp = tmp + 110 continue if(myid.eq.0) then do i = 1, n do j = 1, n print *, arr(i,j,2) enddo enddo endif if(myid.eq.0) then call mpi_send(arr(1,1,1), n*n, mpi_integer, 1, 100, & mpi_comm_world, ierr) print *, 'sending...' print *, arr else call mpi_recv(arr(1,1,2), n*n, mpi_integer, 0, 100, & mpi_comm_world, status, ierr) print *, 'receiving...' print *, arr endif call MPI_FINALIZE(ierr) stop end

 result:

[root@c0109 zlt]# mpiexec -n 2 ./send 10 13 16 11 14 17 12 15 18 sending... 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 receiving... 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 9 19 20 21 22 23 24 25 26 27[root@c0109 zlt]#