0

嗨，我正在尝试使用这样的 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 部分的情况下使用 subarray 命令。但是，如果我将两者结合起来，因为在这种情况下我需要这样做，因为我在实际代码中有一个带有 fortran 的结构。我不知道是否有更好的方法来做到这一点。任何形式的帮助和建议表示赞赏。提前致谢

4

1 回答 1

0

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

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

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

于 2011-03-20T19:19:37.453 回答