FORTRAN中的释放和内存分配问题[关闭]

Question

我的部分 FORTRAN 代码的解除分配和分配方面存在问题。特别是，我认为这个问题与在网上搜索我的错误消息时的内存分配有关。错误消息涉及无效指针，但是，我没有在我的程序中使用任何指针

在完成我的 f 循环的第 2 次迭代后（见下文），程序崩溃了，或者说大部分时间它崩溃了，有时它只是冻结了。我相信这是错误所在。随着程序运行到这一点。

我有未显示的子程序，但由于它们适用于其他模拟组合，我有理由相信它们不是问题。我正在程序中的其他地方使用 deallocate 和 allocate（成功），所以我很惊讶它在这里不起作用。

为了便于阅读，我只展示了程序的一部分。特别是，我删除了对我编写的子例程的调用。我希望我已经为你们程序员提供了足够的信息来帮助我找出问题所在。如果不是，请指定您想要的其他信息，我很乐意遵守。我已经使用各种编译器选项编译了该程序，并修复了一些错误并删除了任何警告。然而，此时编译器选项并没有给我更多信息。

allocate(poffvect(1:6))
allocate(phi1out(1:1))
allocate(phi2out(1:1))
allocate(phi1outs1(1:1))
allocate(phi2outs1(1:1))

     dummy allocation
   allocate(phi1outind(1:1))
   allocate(phi2outind(1:1))
   allocate(phi1outinds1(1:1))
   allocate(phi2outinds1(1:1))

   do e = 1, 6
     print *,"e", e
     do f = 1, 3

       print *,"f", f, iteratst1(f), trim(filenumcharimp)

       deallocate(phi1outinds1, STAT = AllocateStatus)
     if (AllocateStatus /= 0) stop "Error during  deallocation of phi1outinds1"
     print *, "Allocatestatus of phi1outinds1 is",  AllocateStatus
     deallocate(phi2outinds1, STAT = AllocateStatus)
    print *, "DeAllocatestatus of phi1outinds2 is", AllocateStatus

     if (AllocateStatus /= 0) stop "Error during deallocation of phi2outinds1"
    print *, "we deallocate f loop ok", iteratst1(f)

      allocate(phi1outinds1(1:iteratst1(f)), STAT = AllocateStatus)
     if (AllocateStatus /= 0) stop "Error during allocation of phi1outinds1"
     allocate(phi2outinds1(1:iteratst1(f)), STAT = AllocateStatus)
    if (AllocateStatus /= 0) stop "Error during deallocation of phi1outinds1"

 end do
 end do

编译器选项

ifort -free -check  -traceback -o adatptmultistage1new.out adatptmultistage1new.f90

输出

       e           1
       f           1        5000 43
     DeAllocatestatus of phi1outinds1 is           0
    DeAllocatestatus of phi1outinds2 is           0
    we deallocate f loop ok        5000
   f loop done           1
  f           2       10000 43
 Allocatestatus of phi1outinds1 is           0
 DeAllocatestatus of phi1outinds2 is           0
 we deallocate f loop ok       10000
   f loop done           2
   f           3       15000 43
  Allocatestatus of phi1outinds1 is           0

错误信息

*** glibc detected *** ./adatptmultistage1new.out: munmap_chunk(): invalid pointer:   0x0000000000d3ddd0 ***
    ======= Backtrace: =========
    /lib/libc.so.6(+0x77806)[0x7f5863b7b806]
   .    /adatptmultistage1new.out[0x43247c]
   .    /adatptmultistage1new.out[0x404368] 
   ./adatptmultistage1new.out[0x4031ec]
   /lib/libc.so.6(__libc_start_main+0xfd)[0x7f5863b22c4d]
  .    /adatptmultistage1new.out[0x4030e9]
   ======= Memory map: ========
00400000-004d4000 r-xp 00000000 08:03    9642201 
/home/jgold/smwcv/error_infect/test/surfaces/multistage/adaptonly/adatptmultistage1new.out
006d4000-006dc000 rw-p 000d4000 08:03 9642201

[为简洁起见，未显示其余错误消息]

  7fffb004d000-7fffb00bc000 rw-p 00000000 00:00 0                          [stack]
  7fffb01d7000-7fffb01d8000 r-xp 00000000 00:00 0                          [vdso]
  ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0                  [vsyscall]
  Aborted

Answer 1

这是我们要尝试找出的大量代码。您是否使用尽可能多的编译器调试选项对其进行编译？特别是，您是否使用数组边界检查？你用的是什么编译器？我没有看到“使用”语句......最好将您的子例程放入一个模块并“使用”该模块，以便编译器可以检查实际参数和虚拟参数之间的参数一致性。

编辑：“双重释放或损坏”表明内存已损坏。由于您似乎没有任何指针，因此有三种可能的方式来破坏内存：

1. 使用尚未分配的可分配变量。如果 allocate 语句失败，程序可能会在此时抛出错误。您可能正在使用忘记分配的变量。
2. 对过程的调用中的参数与过程期望的参数之间存在分歧，即实际参数和伪参数之间存在分歧。使用模块将允许编译器对此进行更好的检查。
3. 使用非法下标值写入数组大小以外的内容——这将覆盖“随机”内存，例如描述下一个数组的内部结构。打开运行时下标或数组绑定检查将对此进行测试。使用 ifort：-check bounds 或 -check all。对于非常通过检查尝试：-O2 -stand f03 -check all -traceback -warn all -fstack-protector -assume protect_parens -implicitnone