arrays - 通过指针引用重叠数组

Question

我需要使用调整到重叠数组的指针执行 memmove() 风格的操作。更准确地说，我需要将 b(:) 复制到 a(:) 中。

在 C 编程语言中，我会调用 memmove(a, b, N)。但是，据我所知，Fortran 标准在指针别名方面非常严格。

那么，以下三个选项中哪一个可以安全使用(根据标准)，哪个会导致未定义的行为:

• 带有别名指针的向量语法(选项 1)，
• 使用别名指针的显式 do 循环(选项 2)，
• 调用将重叠数组作为参数的例程(选项 3)

?

program ptr
  implicit none
  integer, parameter :: ARRAY_SIZE = 10, N = 6
  integer, target :: array(ARRAY_SIZE)
  integer, dimension(:), pointer :: a, b
  integer :: i

  a => array(1: N)
  b => array(3: N+2)

  ! Option 1
  a(1: N) = b(1: N)

  ! Option 2
  do i = 1, N
    a(i) = b(i)
  end do

  ! Option 3
  call foobar(a, b, N)

  contains

    subroutine foobar(a, b, length)
      integer, dimension(:), intent(out) :: a
      integer, dimension(:), intent(in) :: b
      integer, intent(in) :: length

      a(1: length) = b(1: length)
    end subroutine foobar
end program ptr

Answer 1

选项 1 和 2 都可以，编译器知道指针可以别名。

子例程本身没问题，但您不能向其传递两个别名参数，这确实违反了标准(未定义行为是 C 术语)。您可以使参数成为指针，然后可以使用别名并且程序将符合标准。

  integer, dimension(:), pointer :: a
  integer, dimension(:), pointer :: b

我还删除了意图，因为它们与指针关联状态相关，而不是目标值，正如 IanH 评论的那样。此外，他们还需要 Fortran2003。