分配数组时出错

Question

我有一个子程序，它在程序运行期间被调用了很多次。我尝试使用尽可能多的可分配数组，并且多次调用子例程没有任何问题，但在某些时候，它终止于：

malloc.c:3790: _int_malloc: Assertion `(unsigned long)(size) >= (unsigned long)(nb)' failed.

当第一个数组被分配时，这发生在子例程的开头。

改用不可分配的数组，子例程被多次调用但再次终止，现在：

 wait: 28674: Memory fault(coredump)

我假设它在调用时终止，因为我在变量声明后立即写出一些值，没有任何计算。
呼召

do k=1, kreise
     write(*,*)k 

     call rundheit(n(k),kreis(k,1:n(k),3),kreis(k,1:n(k),2),outrnd)

end do

其中 'kreise' 的值可能高达 1500。我打印出并检查了在调用之前、子例程中和调用之后传递的参数的值。

限制 'kreise' 确实解决了问题，但限制不是一个实际的解决方案。我需要评估所有数据。不是它的断裂。

对我的环境的一些说明：
我的程序是使用英特尔 Fortran 编译器由 FEM 仿真软件编译的子例程。据我所知，我没有机会改变编译器选项，我也不能自己编译我的代码，因为它对 FEM 软件部署的子程序有很多依赖。

我在另一个更小、更简单的 Simulation 上开发并运行了这个精确的子程序，没有任何问题。只要我不使用这个特定的子程序，实际的、“更大的”模拟运行也没有任何问题。（差异主要在于节点密度，因此计算过程中考虑的数据量）其他用户子程序工作没有问题。所有子程序所做的就是在一些增量之间获取结果，进行一些分析并编写一些报告而不改变模拟。

我猜这个问题与内存处理有关，我不知道。

谢谢。

更新
我使用编译子例程-check all并发现错误发生在被指责的子例程之前。两个数组，其中一个是 n()，多次越界，但在调用时错误变得（更）严重。奇怪的部分是，当错误发生时，它是一些超出界限的迭代，例如：这里两个数组都有大小(1:72)，调用在 k=135 到 267 的某个地方中断（我在某些运行中发现的最低和最高值）。

问题是 integer Kreise，该值是在循环期间设置的：

...
allocate(n(l))
allocate(pos(l))
...
do kreise = 1,l
   pos(kreise)=minvalX+(Kreise-1)*IncX
   if(pos(kreise).gt.maxvalX) exit
end do

凡kreise八方通变为L + 1。为什么？

注意：pos(kreise).gt.maxvalX永远不应该是真的。成为真的不是问题，尽管它表明l计算错误（大）。这个出口只会在以后通过减少几个循环的迭代来节省计算时间。

Answer 1

该程序可能正在写入不应写入的内存，并损坏 Fortran allocate 使用的 malloc 内存管理结构。我建议使用 Fortran 选项进行运行时下标检查。使用 ifort，尝试-check all或-check bounds。