我们最近发现在Fortran中向未分配的数组进行了赋值。GNU gfortran编译器没有捕捉到错误,而且该代码可以在OSX和Linux下运行。然而,在IBM Power PC上,相同的代码会导致“段错误”。
我的问题是,下面的代码是否正确?似乎已分配给array
的数组在某些体系结构上自动分配内存,但在其他体系结构上则不会。这里有一些特定于实现的细节在起作用吗?
该代码是混合C / Fortran代码:
#include <stdlib.h>
void assign_array_(double x[], int* n);
void print_array_();
int main()
{
int n,i;
double *x;
n = 5;
x = (double*) malloc(sizeof(double)*n);
for (i = 0; i < n; i++)
x[i] = (double) i;
assign_array_(x,&n);
print_array_();
return 0;
}
而且Fortran代码:
MODULE test_mod
DOUBLE PRECISION, ALLOCATABLE, DIMENSION(:) :: array
integer :: nsize
END MODULE test_mod
SUBROUTINE assign_array(x,n)
USE test_mod
IMPLICIT NONE
INTEGER :: n
DOUBLE PRECISION :: x(n)
CALL test_allocated()
array = x
CALL test_allocated()
nsize = n
END SUBROUTINE assign_array
SUBROUTINE print_array()
USE test_mod, ONLY: nsize, array
IMPLICIT NONE
INTEGER :: i
DO i = 1,nsize
WRITE(6,'(F24.16)') array(i)
END DO
END SUBROUTINE print_array
SUBROUTINE test_allocated()
USE test_mod
IMPLICIT NONE
IF (ALLOCATED(array)) THEN
WRITE(6,*) 'Array is allocated'
WRITE(6,*) 'size is ', SIZE(array)
ELSE
WRITE(6,*) 'Array is NOT allocated'
END IF
END SUBROUTINE test_allocated
输出结果(运行时)为:
Array is NOT allocated
Array is allocated
size is 5
0.0000000000000000
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
这是在 Power PC 上的输出:
Array is NOT allocated
Segmentation fault (core dumped)
总之,该代码在GNU (GNU Fortran (MacPorts gcc5 5.4.0_0) 5.4.0) gfortran编译器下在OSX(arch:x86_64h)和Linux上运行(在OSX上托管的虚拟机中,使用GNU Fortran (Ubuntu 4.9.4-2ubuntu1~14.04.1) 4.9.4编译),但是在使用GNU Fortran (GCC) 4.4.7 20120313 (Red Hat 4.4.7-17)编译的Power PC(arch: ppc64)上无法运行。在我们的原始代码中,Power PC实现仅在稍后引用了分配数组的条目时才出现段错误,使我们的“错误”(如果它实际上是一个错误)真的很难追踪。
以上代码的正确行为是什么?
array = x
是非法的。 - stark