1

考虑这个简单的程序

program main
implicit none

integer :: array(2,3) = transpose(reshape((/1,2,3,4,5,6/),(/ size(array, 2), size(array, 1) /)))
integer i,j
do i=1,2
   write(*,*) array(i,:)
end do
end program main

当我使用 intel fortran 编译器(版本 13.0.0)编译时出现错误:

main.f90(4):错误 #6263:此内在函数在常量表达式中无效。[移调]

integer :: array(2,3) = transpose(reshape((/1,2,3,4,5,6/),(/ size(array, 2), size(array, 1) /)))

似乎转置不能用于常量表达式(?)。那么有没有办法在定义时初始化数组?在我的问题中,数组非常大,因此手动转置不是一种选择。

4

2 回答 2

3

在初始化表达式(或 F2008 术语中的常量表达式)中使用 TRANSPOSE 是该编译器尚不支持的 Fortran 2003 功能。

使用带有适当 ORDER 参数的 RESHAPE 是 Fortran 95 的一部分,由该编译器支持,并且可以提供您想要的等价物:

integer :: array(2,3) = reshape([1,2,3,4,5,6], shape(array), ORDER=[2,1])
于 2013-10-12T23:48:57.370 回答
0

这意味着你不能在那里定义 array。如果你这样做

integer :: array(2,3)
integer :: i,j
array = transpose....

它会编译。

于 2013-10-12T23:38:52.493 回答