所以这很奇怪..我使用Halide(https://github.com/halide/Halide/),但我得到的错误与C ++有关。当我运行代码时,我得到了
*** glibc detected *** ./lesson_04: malloc(): memory corruption: 0x0000000000a36640 ***
它还打印"实现"但是没有打印"模糊初始"。它通过sigma等于1.0但不适用于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;
}
答案 0 :(得分:1)
有时输出卡在写缓冲区中,如果程序过早退出,则实际上永远不会写入控制台。调用fflush(stdout);
将帮助您确定是否属于这种情况。或者使用fprintf(stderr, ...)
进行调试输出,因为stderr
是无缓冲的。
然而,更有可能的是,您的程序实际上在打印“实现”和“模糊初始”...... Image<uint8_t>
复制构造函数之间发生的一项主要任务中崩溃。
为什么你通过值传递并导致副本超出我的范围,我宁愿在那里使用const
引用来避免复制大对象。尽管如此,如果你的拷贝构造函数中有一个错误,那么首先要找到并修复它是绝对值得的。