1

所以这很奇怪..我正在使用 Halide ( https://github.com/halide/Halide/ )，但我得到的错误是 C++ 相关的。当我运行我得到的代码时

*** glibc detected *** ./lesson_04: malloc(): memory corruption: 0x0000000000a36640 ***

它还打印“realizing”，但不打印“blur initial”。当我通过等于 1.0 的 sigma 但不适用于 2.0 或 3.0 时，它可以工作。这是一段较长的代码。更长的代码即使对于 1 也不起作用。我不明白一个值如何影响任何事情以及它在哪里卡住。有什么帮助吗？谢谢

#include <Halide.h>
#include <stdio.h>
#include <math.h>

#include "armadillo"
using namespace arma;

using namespace Halide;
#include "../apps/support/image_io.h"


Image<float> calculateGaussian(double sigmaD){
    int sigma=(int)sigmaD;
    Image<float> convolution(sigma*6+1,sigma*6+1); float sum=0;

    for (int i=-sigma*3; i<sigma*3+2; i++){
        for (int j=-sigma*3; j<sigma*3+2; j++){
            float r=std::sqrt(i*i+j*j);
            convolution(i+sigma*3,j+sigma*3)=std::exp(-r*r/(2*sigma*sigma));
            sum+=convolution(i+sigma*3,j+sigma*3);
        }
    }
    for (int i=0; i<sigma*6+1; i++){
        for (int j=0; j<sigma*6+1; j++){
            convolution(i,j)/=sum;
        }
    }
    return convolution;
}

Image<uint8_t> convolve(Image<uint8_t> im, Image<float> convolution){
    RDom r(convolution);
    int w=convolution.width(), h=convolution.height();
    int W=im.width(), H=im.height();
    Var x,y,c;
    Func blur, input;
    input(x, y,c) = im(clamp(x, 0, W-1), clamp(y, 0, H-1),c);   
    blur(x, y,c) = sum(cast<uint8_t>(convolution(r.x, r.y) * input(x + r.x - w/2, y + r.y - h/2,c)));
    Image<uint8_t> out=blur.realize(im.width(),im.height(),im.channels());
    return out;
}

Image<uint8_t> blurGaussian(Image<uint8_t> im, float sigma){
    printf("blur initial");
    Image<float> gaussian=calculateGaussian(sigma);

    return convolve(im,gaussian);
}

int main(int argc, char **argv) {

    Image<uint8_t> im1 = load<uint8_t>("../apps/images/stata-1.png");
    Image<uint8_t> output;
    printf("realizing \n");
    output=blurGaussian(im1,3.0);
    save(output, "blur.png");

    printf("Success!\n");
    return 0;
}
4

1 回答 1

1

有时，如果程序过早退出，输出会卡在写缓冲区中，并且永远不会真正写入控制台。致电fflush(stdout);将帮助您确定是否属于这种情况。或者fprintf(stderr, ...)用于调试输出，因为stderr没有缓冲。

然而，更有可能的是，您的程序实际上确实在打印“实现”和“模糊初始”之间发生的一项主要任务中崩溃......Image<uint8_t>复制构造函数。

为什么你按值传递并导致复制超出了我的范围，我更喜欢const那里的引用以避免复制大对象。尽管如此，如果您的复制构造函数中存在错误，则绝对值得首先找到并修复它。

于 2013-06-27T03:35:58.850 回答