mpi - 为什么使用 `MPI_Type_vector` 来指定步幅差距？

Question

MPI_Type_vector在 MPI 文档中，我读到和之间的区别在于MPI_Type_contiguos，第一个允许指定数组元素之间的步幅间隙。

为什么人们应该对此感兴趣而不是简单地使用 a MPI_Type_contiguous？

Answer 1

您可以将其视为没有大步MPI_Type_contiguous()的特例。MPI_Type_vector()

例如，在 Fortran 中，数组是列主要的。所以如果你考虑以下数组

integer :: A(n,m)

列的派生数据类型可以通过

CALL MPI_Type_contiguous(n, MPI_INTEGER, newtype, ierr)

或者

CALL MPI_Type_vector(n, 1, 1, MPI_INTEGER, newtype, ierr)

或者

CALL MPI_Type_vector(1, n, n, MPI_INTEGER, newtype, ierr)

但描述一行的唯一方法是通过MPI_Type_vector()

CALL MPI_Type_vector(m, 1, n, MPI_INTEGER, newtype, ierr)

请注意，如果您想一次发送/接收多行，则需要调整数据类型的大小，因此完整的序列将是

CALL MPI_Type_vector(m, 1, n, MPI_INTEGER, tmptype, ierr)
CALL MPI_Type_size(MPI_INTEGER, integer_size, ierr)
CALL MPI_Type_create_resized(tmptype, 0, integer_size, newtype, ierr)
CALL MPI_Type_free(tmptype, ierr)