在FORTRAN90结构中使用MPI_type_contiguous和MPI_Type_CREATE_Subarray

Question

您好，我正在尝试使用fortran结构，如下所示

type some
   u                ! actual code will have 17 such scalars
end type some
TYPE(some),ALLOCATABLE,DIMENSION(:) :: metvars,newmetvars

现在，我的测试程序的目标是将10个数字从一个处理器发送到另一个处理器，但这10个数字的起点将是我的选择(例如，如果我有一个向量，假设20个数字不是必须的，我将把前10个数字带到下一个处理器，但假设我的选择是5到15)。所以首先使用mpi_type_contiguous，如下所示

 CALL MPI_TYPE_CONTIGUOUS(10,MPI_REAL,MPI_METVARS,ierr) ! declaring a derived datatype of the object to make it in to contiguous memory
 CALL MPI_TYPE_COMMIT(MPI_METVARS,ierr)

我做了发送记录，并能够将前10个数字发送到另一个处理器(我正在为2个处理器测试它)

 if(rank.EQ.0)then
     do k= 2,nz-1
     metvars(k)%u = k
     un(k)=k
     enddo
 endif

我把这个发给你

现在，对于第二部分，我使用了mpi_TYPE_CREATE_SUBARRAY

   array_size = (/20/)
   array_subsize =(/10/)
   array_start = (/5/)

   CALL MPI_TYPE_CREATE_SUBARRAY(1,array_size,array_subsize,array_start,MPI_ORDER_FORTRAN,MPI_METVARS,newtype,ierr)
   CALL MPI_TYPE_COMMIT(newtype,ierr)

   array_size = (/20/)
   array_subsize =(/10/)
   array_start = (/0/)

   CALL MPI_TYPE_CREATE_SUBARRAY(1,array_size,array_subsize,array_start,MPI_ORDER_FORTRAN,MPI_METVARS,newtype2,ierr)
   CALL MPI_TYPE_COMMIT(newtype2,ierr)

  if(rank .EQ. 0)then
     CALL MPI_SEND(metvars,1,newtype,1,19,MPI_COMM_WORLD,ierr)
  endif

  if(rank .eq. 1)then
     CALL MPI_RECV(newmetvars,1,newtype2,0,19,MPI_COMM_WORLD,MPI_STATUS_IGNORE,ierr)
  endif

我不知道该怎么做。

我收到一个错误，说

[flatm1001:14066] *** An error occurred in MPI_Recv
[flatm1001:14066] *** on communicator MPI_COMM_WORLD
[flatm1001:14066] *** MPI_ERR_TRUNCATE: message truncated
[flatm1001:14066] *** MPI_ERRORS_ARE_FATAL (your MPI job will now abort)

我在本地机器上使用openmpi。我可以在没有mpi_type_contiguous部件的情况下使用子数组命令。但是，如果我将两者结合起来，因为在本例中，我需要这样做，因为我在实际代码中有一个带有fortran的结构。我也不知道有没有更好的方法。任何形式的帮助和建议都是非常感谢的。提前感谢

Answer 1

我假设您的自定义类型包含1 real，因为没有指定它。首先构造这些变量中的10个的连续类型，即MPI_METVARS表示10个连续的实数。现在，我不知道这是否真的是问题所在，因为您发布的代码可能不完整，但现在看起来是您构造了一个包含10个MPI_METVARS类型的子数组，这意味着在newtype和newtype2中实际上有100个连续的reals。

处理该结构的“正确”方法是使用MPI_TYPE_CREATE_STRUCT为其创建一个类型，该类型应该是您的MPI_METVARS类型。

因此，请为您的自定义类型提供正确的代码，并检查newtype类型的大小。