传递时字符串数组被无效

Question

我在传递字符串数组时遇到问题。考虑以下示例代码：

! -- Module to declare variable
module my_data
   implicit none

   ! -- Declare as deferred-length allocatable array
   character(len=:), dimension(:), allocatable :: str_array
end module my_data

! -- Module to call subroutine
module my_subs
   implicit none

contains
subroutine a(str_array)
   character(len=*), dimension(:), intent(IN) :: str_array
   integer :: i, j
   character :: c

   do i=1,size(str_array)
      do j=1,len_trim(str_array(i))
         c = str_array(i)(j:j)

         ! -- Write i, j, character, and int representation
         write(*,*) 'In call: ', i, j, ' "'//c//'", ichar = ', ichar(c)
      enddo
   enddo
end subroutine a
end module my_subs

! -- Main program
program main
   use my_data, only : str_array
   use my_subs, only : a
   implicit none
   integer, parameter :: strlen = 200
   integer :: N, i, j
   character :: c

   ! -- Size of str array
   N = 2

   ! -- Allocate str_array, syntax from https://software.intel.com/en-us/forums/intel-visual-fortran-compiler-for-windows/topic/287349
   allocate(character(strlen) :: str_array(N))

   ! -- Set both to the same string
   str_array = 'abc'

   do i=1,size(str_array)
      do j=1,len_trim(str_array(i))
         c = str_array(i)(j:j)

         ! -- Write i, j, character, and int representation
         write(*,*) 'In main: ', i, j, ' "'//c//'", ichar = ', ichar(c)
      enddo
   enddo

   call a(str_array)
end program main

字符串数组被声明为假定长度元素的数组（来自the wiki）。我分配和设置字符串的值（两个元素，在这个例子中都是abc）。主程序输出有关字符串的完整详细信息，然后调用一个子程序，该子程序也输出（希望是相同的）完整详细信息。

使用 PGI、GCC 或 Intel 15.0，我得到了我期望的结果：

chaud106@ln0005 [~/Testing] % ifort --version && ifort -check all -warn all main.f90 && ./a.out
ifort (IFORT) 15.0.3 20150407
Copyright (C) 1985-2015 Intel Corporation.  All rights reserved.

 In main:            1           1  "a", ichar =           97
 In main:            1           2  "b", ichar =           98
 In main:            1           3  "c", ichar =           99
 In main:            2           1  "a", ichar =           97
 In main:            2           2  "b", ichar =           98
 In main:            2           3  "c", ichar =           99
 In call:            1           1  "a", ichar =           97
 In call:            1           2  "b", ichar =           98
 In call:            1           3  "c", ichar =           99
 In call:            2           1  "a", ichar =           97
 In call:            2           2  "b", ichar =           98
 In call:            2           3  "c", ichar =           99

但是，Intel 18.0 将字符数组的第二个元素（所有 3 个字符）设置为空字符：

chaud106@ln0005 [~/Testing] % ifort --version && ifort -check all -warn all main.f90 && ./a.out
ifort (IFORT) 18.0.0 20170811
Copyright (C) 1985-2017 Intel Corporation.  All rights reserved.

 In main:            1           1  "a", ichar =           97
 In main:            1           2  "b", ichar =           98
 In main:            1           3  "c", ichar =           99
 In main:            2           1  "a", ichar =           97
 In main:            2           2  "b", ichar =           98
 In main:            2           3  "c", ichar =           99
 In call:            1           1  "a", ichar =           97
 In call:            1           2  "b", ichar =           98
 In call:            1           3  "c", ichar =           99
 In call:            2           1  "", ichar =            0
 In call:            2           2  "", ichar =            0
 In call:            2           3  "", ichar =            0

我有几个与此行为相关的问题：

1. 为什么会发生这种情况？我在想这可能与英特尔强制执行 lhs-relocation 有关，但我不确定。添加-assume norealloc_lhs 并没有改变任何东西。

2. 传递这样的字符串数组的正确语法是什么？我可以用不同的方式声明它并避免这个问题吗？

我在这台机器上可以访问的 Intel 版本具有以下行为：

• ifort (IFORT) 15.0.2 20150121 - 不作废
• ifort (IFORT) 15.0.3 20150407 - 不作废
• ifort (IFORT) 16.0.3 20160415 - 不作废
• ifort (IFORT) 17.0.4 20170411 - 无效
• ifort (IFORT) 18.0.0 20170811 - 无效

在另一台机器上，我没有任何最新的英特尔：

• ifort (IFORT) 14.0.2 20140120 - 不作废
• ifort (IFORT) 16.0.0 20150815 - 不作废

Answer 1

您的程序在 ifort 18.0.3 上按预期运行。

我没有尝试过很多以前的版本，但我注意到 17.0.1 是 Fortran 2003 对内在赋值的自动分配成为该编译器的默认设置的点。

有问题的行似乎是

str_array = 'abc'

这里，str_array 应该首先被释放，因为右边是一个与左边不同长度参数的表达式。然后它将被分配为长度为 3 的字符（右手边的长度）和形状 [2] 和以前一样（右手边是标量）。这确实发生了，正如SIZE(str_array) 和LEN(str_array) 所见。不过，稍后将其用作实际参数时会出现一些问题。

有一些方法可以解决 18.0.1 的这个问题：

• 为虚拟参数赋予value 属性；
• str_array = ['abc','abc']（不需要之前的分配）；
• str_array(:) = 'abc'（如果您不想重新分配）。

可能还有许多其他可用，具体取决于您的需要。不过，如果可以，请将您的编译器升级到最新版本。