使用FORTRAN 90中的MPI_Type_Create_Subarray 90

这与“使用MPI_Send/Recv处理Fortran 90中的多灰数组的块”上的较早线程相关。我的阵列是真实的，3调光说5 5 5，即X（1：5,1：5,1：5）如果我想给我的数组的下一部分说使用FORTRAN 90中的MPI_Type_Create_Subarray 90

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示例。

来源

2013-10-21 Madhurjya

该标准规定array_of_starts（您的starts）开始索引为零，而不是一个！因此，如果Fortran编号从1开始，则必须减去1。因此，您的starts阵列应该是(/ 1, 1, 3 /)。

从MPI 3.0 Standard, Ch. 4.1.3：

给用户的建议。在Fortran程序中，数组索引从1开始，如果子数组的特定维度的起始坐标为N，则该维度的array_of_starts中的条目为 n-1。（向用户建议结束）

来源

2013-10-21 12:09:49

感谢您的信息。如果数组从一个负数索引开始，对于父数组x（-n：n，-m：m），我想发送数据块x（-2：0，:)？现在应该是什么开始数组？ – Madhurjya

@Madhurjya，'-2'与'-n'的偏移量是'（-2） - （ - n）= n-2'。对于第二维，偏移量为“0”。 –

谢谢Alexander和Hristo！讨论对我有很大的帮助。 – Madhurjya

使用FORTRAN 90中的MPI_Type_Create_Subarray 90

回答

相关问题