奇怪的回溯 - 错误在哪里？

Question

我正在用 C++ 开发图像处理应用程序。 我见过很多编译器错误和回溯，但这对我来说是新的。

#0  0xb80c5430 in __kernel_vsyscall ()
#1  0xb7d1b6d0 in raise () from /lib/tls/i686/cmov/libc.so.6
#2  0xb7d1d098 in abort () from /lib/tls/i686/cmov/libc.so.6
#3  0xb7d5924d in ?? () from /lib/tls/i686/cmov/libc.so.6
#4  0xb7d62276 in ?? () from /lib/tls/i686/cmov/libc.so.6
#5  0xb7d639c5 in malloc () from /lib/tls/i686/cmov/libc.so.6
#6  0xb7f42f47 in operator new () from /usr/lib/libstdc++.so.6
#7  0x0805bd20 in Image<Color>::fft (this=0xb467640) at ../image_processing/image.cpp:545

这里发生了什么事？ new 操作符崩溃了，好吧。 但为什么？ 这并不是内存不足（它尝试分配大约 128Kb，一个 128x64 像素，每个像素有两个浮点数）。 另外，它也没有接缝，因为它是我自己的代码中的错误（构造函数没有被触及！）。

上述行 (#7) 中的代码是：

Image<Complex> *result = new Image<Complex>(this->resX, resY); 
// this->resX = 128, resY = 64 (both int), Complex is a typedef for std::complex<float>

几乎相同的实例化适用于我的代码中的其他位置。 如果我注释掉这部分代码，稍后它会在类似的部分崩溃。 我不明白，我也没有任何想法，如何调试它。 有什么帮助吗？

编译器是 gcc 4.3.3，libc 是 2.9（均来自 Ubuntu Jaunty）

更新：

我在相同的方法和 main() 中的错误行上方包含了以下几行

    Image<Complex> *test = new Image<Complex>(128, 64);
    delete test;

奇怪的事情：在同样的方法中它会崩溃，在 main() 中它不会。 正如我所提到的，Complex 是 std::complex的 typedef。 构造函数不会被调用，我在该行之前和构造函数本身中插入了一个 cout 。

更新 2：

感谢 KPexEA 提供此提示！ 我试过这个：

Image<Complex> *test = new Image<Complex>(128, 64);
delete test;

kiss_fft_cpx *output = (kiss_fft_cpx*) malloc( this->resX * this->resY/2 * sizeof(kiss_fft_cpx) );
kiss_fftndr( cfg, input, output );

Image<Complex> *test2 = new Image<Complex>(128, 64);
delete test2;

它崩溃了 - 你猜？ - 测试2！ 所以我的 Kissfft 的 malloc 似乎是有问题的。 我会看一下。

最终更新：

好的，完成了！ 感谢大家！

其实我早该注意到的。 上周，我注意到 Kissfft（一个快速傅里叶变换库）从 128x128 像素源图像生成了 130x64 像素 FFT 图像。 是的，130 像素宽，而不是 128。不要问我为什么，我不知道！ 因此，必须分配 130x64x2xsizeof(float) 字节，而不是我之前想到的 128x64x...。 奇怪的是，它并不是在我修复该错误后就崩溃了，而是几天后。

作为记录，我的最终代码是：

int resY = (int) ceil(this->resY/2);

kiss_fft_cpx *output = (kiss_fft_cpx*) malloc( (this->resX+2) * resY * sizeof(kiss_fft_cpx) );
kiss_fftndr( cfg, input, output );

Image<Complex> *result = new Image<Complex>(this->resX, resY);

谢谢！

克雷什

Answer 1

也许先前分配的内存块存在缓冲区溢出，从而损坏了堆？

Answer 2

您没有分配足够的内存。 Kissfft 的半谱格式（以及 FFTW 和 IMKL）包含 X*(Y/2+1) 复杂元素。

参见kiss_fftndr.h头文件：

/*
输入 timedata 具有 dims[0] X dims[1] X ... X dims[ndims-1] 标量点

输出 freqdata 具有 dims[0] X dims[1] X ... X dims[ndims-1]/2 +1 复杂点
*