我正在尝试用于mtrace
检测 fortran 程序中的内存泄漏。我正在使用 gfortran 编译器。有关 mtrace 的(工作)C 示例,请参阅维基百科条目:http ://en.wikipedia.org/wiki/Mtrace
我尝试了两种方法,即包装 mtrace() 和 muntrace() 并从 fortran 程序中调用它们,以及创建一个直接调用 mtrace() 和 muntrace() 的 C 程序,除了中间的泄漏 fortran 代码。这两种方法都无法检测到内存泄漏,但这里我只介绍后者。
例子.c
#include <stdlib.h>
#include <mcheck.h>
extern void leaky_(); // this might be different on your system
// if it doesn't work, try to run:
// 1) gfortran leaky.f90 -c
// 2) nm leaky.o
// and then change this declaration and its use below
void main() {
mtrace();
leaky_();
muntrace();
}
泄漏的.f90
subroutine leaky()
real, allocatable, dimension(:) :: tmp
integer :: error
allocate (tmp(10), stat=error)
if (error /= 0) then
print*, "subroutine leaky could not allocate space for array tmp"
endif
tmp = 1
!of course the actual code makes more...
print*, ' subroutine leaky run '
return
end subroutine leaky
我编译:
gfortran -g example.c leaky.f90
然后我运行:
export MALLOC_TRACE=`pwd`/raw.txt; ./a.out
然后我解析raw.txt
mtrace
输出:
mtrace a.out raw.txt
并得到:
没有内存泄漏。
有什么我做错了,或者我可以做些什么来mtrace
找到泄漏的fortran内存分配?我猜 gfortran 正在使用不同的malloc
调用,它mtrace
不会跟踪......事实上,正如我在上面所写的,如果我编写一个调用 (wrapped)mtrace()
和muntrace()
.
已编辑:我考虑了其他选项(包括此处尚未提及的一些选项),但正在调试的实际代码在 P6/AIX 上运行,因此 Valgrind 将“只是”不方便(它需要在不同的机器上运行),而 Forcheck 将是不方便(它需要在不同的机器上运行)和昂贵的(〜3k $)。目前 mtrace 将是最好的解决方案,如果它有效的话。
再次编辑:我的猜测
我猜 gfortran 正在使用不同的
malloc
调用,它mtrace
不会跟踪......
是正确的。查看可执行文件(使用nm
或readelf
)没有任何malloc()
调用,而是调用_gfortran_allocate_array
- 可能会调用 malloc)。还有其他想法吗?
再次编辑:我发布了答案,但我不能接受(访问http://stackoverflow.uservoice.com/pages/general/suggestions/39426并请求该功能,它确实需要 - 不需要获得声誉)