使用 ifort 进行数组复制时程序崩溃

这个程序崩溃了Illegal instruction: 4在 MacOSX Lion 和 ifort (IFORT) 12.1.0 20111011 上

program foo
      real, pointer :: a(:,:), b(:,:)
      allocate(a(5400, 5400))
      allocate(b(5400, 3600))
      a=1.0
      b(:, 1:3600) = a(:, 1:3600)

      print *, a
      print *, b

      deallocate(a)
      deallocate(b)

end program 

相同的程序适用于 gfortran。我没有看到任何问题。有任何想法吗 ？展开副本并对列执行显式循环在两个编译器中都有效。

请注意，使用可分配而不是指针我没有问题。

无论语句是否位于模块内部，行为都是相同的。 我在 ifort (IFORT) 12.1.3 20120130 上确认了相同的行为。

显然，Linux和ifort 12.1.5没有出现问题

我尝试使用以下链接选项增加堆栈大小

ifort -Wl,-stack_size,0x40000000,-stack_addr,0xf0000000 test.f90

但我仍然遇到同样的错误。增加 ulimit -s 来解决同样的问题。

编辑2：我做了一些更多的调试，显然问题发生在数组拼接操作时

      b(:, 1:3600) = a(:, 1:3600)

涉及可疑接近 16 M 数据的值。

我正在比较生成的操作码，但如果有一种方法可以查看更具沟通性的中间代码形式，我将很乐意欣赏。

您的程序是正确的（尽管如果您不需要能够重新指向它，我更喜欢可分配到指针）。问题是 ifort 默认将所有临时数组放置在堆栈上，无论它们有多大。它似乎需要一个临时数组来完成您在这里执行的复制操作。为了解决 ifort 愚蠢的默认行为，always编译时使用 -heap-arrays 标志。 IE。

ifort -o test test.f90 -heap-arrays 1600

-heap-arrays 后面的数字是应该开始使用堆的阈值。对于低于此的大小，使用堆栈。我在这里选择了一个相当低的数字 - 您可能可以安全地使用更高的数字。理论上，堆栈数组更快，但差异通常完全可以忽略不计。我希望英特尔能解决这个问题。每个其他编译器对此设置都有合理的默认值。