3

我遇到了一个分段错误，我无法用简单的代码真正理解，这只是：

  • 调用 MPI_INIT
  • 通过 MPI_COMM_DUP 复制全局通信器
  • 通过 MPI_COMM_GROUP 创建一个包含全局通信器一半进程的组
  • 最后从这个组通过 MPI_COMM_CREATE_GROUP 创建一个新的通信器

具体来说，我使用最后一个调用，而不是仅使用 MPI_COMM_CREATE，因为它仅对 group 中包含的进程组进行集体处理，而 MPI_COMM_CREATE 对 COMM 中的每个进程进行集体处理。代码如下

program mpi_comm_create_grp
  use mpi
  IMPLICIT NONE

  INTEGER :: mpi_size,  mpi_err_code
  INTEGER :: my_comm_dup, mpi_new_comm, mpi_group_world, mpi_new_group
  INTEGER :: rank_index
  INTEGER, DIMENSION(:), ALLOCATABLE :: rank_vec

  CALL mpi_init(mpi_err_code)
  CALL mpi_comm_size(mpi_comm_world, mpi_size, mpi_err_code)

  !! allocate and fill the vector for the new group
  allocate(rank_vec(mpi_size/2))
  rank_vec(:) = (/ (rank_index , rank_index=0, mpi_size/2) /)

  !! create the group directly from the comm_world: this way works
  ! CALL mpi_comm_group(mpi_comm_world, mpi_group_world, mpi_err_code)

  !! duplicating the comm_world creating the group form the dup: this ways fails
  CALL mpi_comm_dup(mpi_comm_world, my_comm_dup, mpi_err_code)
  !! creatig the group of all processes from the duplicated comm_world
  CALL mpi_comm_group(my_comm_dup, mpi_group_world, mpi_err_code)

  !! create a new group with just half of processes in comm_world
  CALL mpi_group_incl(mpi_group_world, mpi_size/2, rank_vec,mpi_new_group, mpi_err_code)

  !! create a new comm from the comm_world using the new group created
  CALL mpi_comm_create_group(mpi_comm_world, mpi_new_group, 0, mpi_new_comm, mpi_err_code)

  !! deallocate and finalize mpi
  if(ALLOCATED(rank_vec)) DEALLOCATE(rank_vec)
  CALL mpi_finalize(mpi_err_code)
end program !mpi_comm_create_grp

如果我不复制 COMM_WORLD，而是直接从全局通信器（注释行）创建组，那么一切正常。

我正在使用的并行调试器将段错误追溯到对 MPI_GROUP_TRANSLATE_RANKS 的调用，但是，据我所知，MPI_COMM_DUP 复制了复制的通信器的所有属性，包括排名编号。

我使用的是 ifort 版本 18.0.5，但我也尝试了 17.0.4 和 19.0.2，但没有更好的结果。

4

2 回答 2

3

好吧，这件事有点棘手，至少对我来说是这样，但经过一些测试和帮助，问题的根源找到了。

在代码中

CALL mpi_comm_create_group(mpi_comm_world, mpi_new_group, 0, mpi_new_comm, mpi_err_code)

mpi_new_group为之前创建的 group 创建一个新的通信器。但是mpi_comm_world，用作第一个参数的 与 不在同一上下文中mpi_new_group，如 mpich参考中所述：

MPI_COMM_DUP将在与 comm 相同的组上创建一个新的通信器，但具有新的上下文

所以正确的调用是：

CALL mpi_comm_create_group(my_comm_copy, mpi_new_group, 0, mpi_new_comm, mpi_err_code)

即，替换mpi_comm_worldfor my_comm_copy，即从中mpi_group_world创建的那个。

我仍然不确定为什么它与 OpenMPI 一起使用，但它通常对这类事情更宽容。

于 2019-04-02T07:11:09.303 回答
1

就像我写给openmpi用户列表的评论中建议的那样，他们回复了

这是完全有效的。组成该组的 MPI 进程都是通信世界的一部分。我会向英特尔 MPI 提交一个错误。

所以我尝试在英特尔论坛上发布一个问题。这是他们在 libray 的最新版本 19.3 中解决的一个错误。

于 2019-04-11T07:50:53.813 回答