2

我正在用 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<float> 的 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);

谢谢!

崩溃

4

2 回答 2

5

也许之前分配的内存块有一个缓冲区溢出,它破坏了堆?

于 2009-08-05T06:13:02.317 回答
0

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

查看kiss_fftndr.h 头文件:

/* 输入时间数据有 dims[0] X dims[1] X ... X dims[ndims-1] 标量点

输出频率数据有 dims[0] X dims[1] X ... X dims[ndims-1]/2+1 复点 *

于 2009-08-13T00:26:47.480 回答