3

これは、「MPI_Send/Recv を使用して Fortran 90 で多次元配列のチャンクを処理する」に関する以前のスレッドに関連しています。私の配列は実数で、5 x 5 x 5、つまり x(1:5,1:5,1:5) と言う 3 次元です。配列の次の部分を送信したい場合は、

x(2:3,2:5,4:5)

0からProc 1まで、次のテストプログラムを使用しています

program mpi
implicit none
include "mpif.h"

integer :: ierr,myid,nprocs,status(mpi_status_size),i,j,k,&
&          starts(3),newsize(3),oldsize(3)
real    :: x(1:5,1:5,1:5),y(1:5,1:5,1:5),z(2:3,2:5,4:5)
integer :: arr

call mpi_init(ierr)
call mpi_comm_rank(mpi_comm_world,myid,ierr)
call mpi_comm_size(mpi_comm_world,nprocs,ierr)

if(myid == 0) then
  x = 0.0
  call random_number(x)
  starts = (/2,2,4/)
  newsize = (/2,4,2/)
  oldsize = (/5,5,5/)
  call mpi_type_create_subarray(3,oldsize,newsize,starts,mpi_order_fortran, &
  & mpi_real,arr,ierr)
  call mpi_type_commit(arr,ierr)
  call mpi_send(x,1,arr,1,1,mpi_comm_world,ierr)
  do i = 2,3
    do j = 2,5
      do k = 4,5
        print*,'#1',x(i,j,k)
      enddo
    enddo
  enddo
  print*,' '
else
  y = 0.0
  call mpi_recv(z,16,mpi_real,0,1,mpi_comm_world,status,ierr)
  do i = 2,3
    do j = 2,5
      do k = 4,5
        print*,'#2',z(i,j,k)
      enddo
    enddo
  enddo
endif

call mpi_finalize(ierr)

stop
end

「開始」配列に関連するランタイム エラーが発生します。要素は 0 または 1 でなければなりません。ここで何が欠けていますか? これの正しい形は何ですか?これを使用する FORTRAN の例が見つかりませんでした。

4

1 に答える 1