3

一个小型测试程序在 linux 64 位上使用 gfortran (4.4.5) 返回分段错误。n=2_8**22_8 时不存在故障。Gdb 指示在循环的第一次迭代期间函数 mylen 中发生分段错误。

allocate, stat=           0
size :              8388608
len, switch=false :              8388608
Segmentation fault

谢谢

function mylen(abc,n, switch) 

implicit none

logical, intent(in) :: switch
integer(kind=8), intent(in) :: n
logical, dimension(1:n), intent(in) :: abc
integer(kind=8) :: mylen
character(len=size(abc,dim=1,kind=8)) :: tmp
integer(kind=8) :: i

mylen=len(tmp,kind=8)
if (switch) then
  do i=1,len(tmp,kind=8)
    tmp(i:i)='a'
  enddo
endif

end function mylen

program test

implicit none

integer(kind=8) :: n
logical, allocatable :: abc(:)
integer(kind=8) :: mylen
integer :: ierr

n=2_8**23_8
allocate(abc(n),stat=ierr)
print *,'allocate, stat=',ierr
print *,'size :', size(abc,dim=1,kind=8)
print *,'len, switch=false :', mylen(abc,n,.false.)
print *,'len, switch=true  :', mylen(abc,n,.true.)

end program test
4

1 回答 1

4

我对此进行了测试,并意识到如果字符数组太大,您就会破坏堆栈。如果您使用可分配长度的字符串并将其分配在“mylen”的顶部,它会进入堆并且程序可以工作。

于 2013-08-21T03:42:52.297 回答