Question

我面临以下问题，无法弄清楚发生了什么。

我有一个例程，在我的代码开头分配一些工作数组。这些工作数组是数据结构的一部分。

结构定义为：

 Type bond_stat
  Real*8,Allocatable,Dimension( : , : ) ::data_space
  Real*8,Allocatable,Dimension( : , : ) ::data_time
  Real*8,Allocatable,Dimension( : )     ::bin_width         
  Real*8,Allocatable,Dimension( : )     ::time
  Integer,Allocatable,Dimension( : , : )  ::connection_table
 End Type bond_stat
 Type( bond_stat ),allocatable,dimension( : ) ::Hbonds

这是在模块内部完成的。我正在做的下一步是将此结构分配为

Allocate( Hbonds( 1:2 ) )

然后我调用一个子程序，其中分配了结构中包含的数组。现在奇怪的是，如果我尝试用

分配数组

 Allocate( Hbonds( i )%data_space( 1 : 3 , 1 : Nbins ) )
 Allocate( Hbonds( i )%data_time( 1 : 1 , 1 : Nbins ) )
 Allocate( Hbonds( i )%time( 1 : N ) )
 Allocate( Hbonds( i )%bin_width( 1 : 2 ) )
 Allocate( Hbonds( i )%connection_table( 1 : 2 , 1 : N ) ) )

如果我执行了代码状态，则执行时已经分配了时间数组。我100％确定我之前没有分配这个数组。我也确定之前没有使用该数组，也许编译器会分配它，因为它已经在使用中。所以我想我会检查数组是否已经分配，如果已分配，则解除分配并重新分配。如果我这样做，我会在释放行中收到分段错误：

    If ( Allocated( Hbonds( i )%time ) ) then
       Deallocate( Hbonds( i )%time )
    End If
    Allocate( Hbonds( i )%time( 1 : N ) )

我完全不知道这里发生了什么。请任何人帮忙。输出显示为：

Program received signal SIGSEGV: Segmentation fault - invalid memory 
reference.
Backtrace for this error:
#0  0x2b54685d3f0d in ???
#1  0x2b54685d314d in ???
#2  0x2b5468e4acaf in ???
#3  0x2b5468e9715c in ???
#4  0x4d27cf in __bond_statistics_MOD_bond_stat_init
at/fs/home/jona/programs/gfortran/eval_new01/playground/H_bond_analysis.f90:35
#5  0x4c8f5f in __bond_statistics_MOD_prep_bond_stat
at/fs/home/jona/programs/gfortran/eval_new01/playground/H_bond_analysis.f90:
350
#6  0x44477a in ???
#7  0x446724 in ???
#8  0x2b5468e35f44 in ???
#9  0x4016c8 in ???
#10  0xffffffffffffffff in ???
Segmentation fault (core dumped)

当没有尝试在它看起来之前解除分配数组时：

At line 38 of file H_bond_analysis.f90
Fortran runtime error: Attempting to allocate already allocated variable 'hbonds'
Error termination. Backtrace:
#0  0x2ae8991bef0d in ???
#1  0x2ae8991bfa45 in ???
#2  0x2ae8991bfdfa in ???
#3  0x4d28e5 in __bond_statistics_MOD_bond_stat_init
at/fs/home/jona/programs/gfortran/eval_new01/playground/H_bond_analysis.f90:38
#4  0x4c8f5f in __bond_statistics_MOD_prep_bond_stat
at/fs/home/jona/programs/gfortran/eval_new01/playground/H_bond_analysis.f90:350
#5  0x44477a in ???
#6  0x446724 in ???
#7  0x2ae899a20f44 in ???
#8  0x4016c8 in ???
#9  0xffffffffffffffff in ???

第38行是分配时间的行。我正在使用带有标志的gfortran编译器 FLAGS = -fopenmp -g -Wall -fcheck = all -fbounds-check

Answer 1

我现在想出了问题所在。我正在使用makefile来编译代码。今天我想在bond_stat模块中添加更多数组，并在执行makefile后再次遇到同样的错误。我现在做的是手动重新编译所有模块而不改变代码本身的任何内容。这解决了这个问题。

Fortran认为已经分配了未分配的数组

1 个答案: