用valgrind检查fftw3
在我的程序的一个步骤中，我需要卷积图像。 为此，我使用fftw3提供的function。 当我在我的程序上运行valgrind ，我得到了这个堆栈跟踪。 我的函数叫做fftw3 ，它运行fftw3的fftw_plan_dft_r2c_2d两次（一次在图像上，一次在卷积内核上。为了使它更具可读性，我删除了所有地址和进程ID。
HEAP SUMMARY: in use at exit: 62,280 bytes in 683 blocks total heap usage: 178,271 allocs, 177,588 frees, 36,617,058 bytes allocated 3,304 (24 direct, 3,280 indirect) bytes in 1 blocks are definitely lost in loss record 129 of 131 at : malloc (vg_replace_malloc.c:291) by : fftw_malloc_plain (in ./prog) by : fftw_mkapiplan (in ./prog) by : fftw_plan_many_dft_r2c (in ./prog) by : fftw_plan_dft_r2c (in ./prog) by : fftw_plan_dft_r2c_2d (in ./prog) by : convolve (freqdomain.c:199) by : convolve (conv.c:290) by : main (main.c:332) 3,304 (24 direct, 3,280 indirect) bytes in 1 blocks are definitely lost in loss record 130 of 131 at : malloc (vg_replace_malloc.c:291) by : fftw_malloc_plain (in ./prog) by : fftw_mkapiplan (in ./prog) by : fftw_plan_many_dft_r2c (in ./prog) by : fftw_plan_dft_r2c (in ./prog) by : fftw_plan_dft_r2c_2d (in ./prog) by : convolve (freqdomain.c:203) by : convolve (conv.c:290) by : main (main.c:332) LEAK SUMMARY: definitely lost: 48 bytes in 2 blocks indirectly lost: 6,560 bytes in 60 blocks possibly lost: 0 bytes in 0 blocks still reachable: 55,672 bytes in 621 blocks suppressed: 0 bytes in 0 blocks Reachable blocks (those to which a pointer was found) are not shown. To see them, rerun with: --leak-check=full --show-leak-kinds=all For counts of detected and suppressed errors, rerun with: -v ERROR SUMMARY: 3 errors from 3 contexts (suppressed: 4 from 4)
根据手册的建议，完成后我使用了fftw_free和fftw_destroy_plan 。 我想看看有什么我可以用这个做的，或者这是fftw3的内部问题？ 由于未释放的malloc位于FFTW源代码的深处。
编辑：包括fftw_cleanup（）之后
在下面你可以看到我添加fftw_cleanup()后的diff 。
[me@mycomputer]$ diff NoCleanup WithCleanup 2,3c2,3 < in use at exit: 62,280 bytes in 683 blocks in use at exit: 9,008 bytes in 66 blocks > total heap usage: 178,271 allocs, 178,205 frees, 36,617,058 bytes allocated 5c5 3,304 (24 direct, 3,280 indirect) bytes in 1 blocks are definitely lost in loss record 39 of 40 16c16 3,304 (24 direct, 3,280 indirect) bytes in 1 blocks are definitely lost in loss record 40 of 40 31c31 still reachable: 2,400 bytes in 4 blocks
退出时使用的still reachable字节数明显减少，释放的malloc数量也增加了。 但主要的错误（ 3,304 (24 direct, 3,280 indirect) bytes in 1 blocks are definitely lost ）仍然存在。
我相信您的代码可能有错误，但我建议您进行测试以缩小问题范围。
我（随机）发现了一些调用fftw_plan_dft_r2c_2d() 示例代码 ，以及FFTW库函数的其他部分。
在我的环境中：基于amd64的Debian 7.4“wheezy”，libfftw 3.3.2 + Debian补丁（3.3.2-3.1）和gcc 4.7.2（Debian 4.7.2-5）在valgrind下运行此测试代码似乎没有泄露内存一旦你明确地将fftw_cleanup()添加到每个test0#函数的结尾或者在main（）的结尾处，在返回之前。
因此，当你多次调用fftw_plan_dft_r2c_2d()时，除非FFTW库有内存泄漏，我怀疑错误是在你自己的代码中。
如果您认为库中存在内存泄漏而不是使用它，请随意调整John Burkardt编写并在LGPL下许可的示例代码重复调用fftw_plan_dft_r2c_2d() 。
从FFTW文档 ：
FFTW的计划器保存了一些其他持久性数据，例如累积的智慧和当前配置中可用的算法列表。 如果要释放所有这些并将FFTW重置为原始状态，那么当您启动程序时，可以调用：
void fftw_cleanup(void);
我发现只需调用* fftw_destroy_plan *也会从valgrind输出中删除丢失的块消息。
当然，调用* fftw_cleanup *会将所有东西都拖掉。