2D Convolution return wrong result

Emad_B_ — Mon, 05 Oct 2015 23:43:23 GMT

I am new to MKL, so this might be a stupid mistake. I am trying to convolve 2D 1024x1024 array with its elements set to 5 with 3x3 kernel with its element set to 1. The expected result is 45 everywhere except the boundary, however, I am getting positive and negative very small fraction close to zero, and all returned status code are ok.

Here a simple repro:

int status = VSL_STATUS_OK;
int mode = VSL_CONV_MODE_FFT; //VSL_CONV_MODE_AUTO; // VSL_CONV_MODE_DIRECT;
int xstride[] = {1, 1024};
int ystride[] = {1, 3};
int zstride[] = {1, 1026};
int xshape[] = {1024, 1024};
int yshape[] = {3, 3};
int zshape[] = {1026, 1026};
int offset[] = {0, 0};

double* x = new double[1024*1024];
for (int i = 0; i < (1024*1024); i++) {
x = 5.0;
}

double* y = new double[3*3];
for (int i = 0; i < (3*3); i++) {
y = 1.0;
}

double* z = new double[1026*1026];
memset((void*)z, 0, 1026*1026*sizeof(double));

VSLConvTaskPtr task;
status = vslsConvNewTask(&task, mode, 2, xshape, yshape, zshape);
assert(status == VSL_STATUS_OK);

status = vslConvSetStart(task, offset);
assert(status == VSL_STATUS_OK);

status = vsldConvExec(task, x, xstride, y, ystride, z, zstride);
assert(status == VSL_STATUS_OK);

vslConvDeleteTask(&task);

delete[] x;

delete[] y;

delete[] z;

Hi Emad B.,

Evgueni_P_Intel — Tue, 06 Oct 2015 04:19:12 GMT

Hi Emad B.,

Since d means double, s means single, please use vsldConvNewTask.

Evgueni

Thanks Evgueni, that's did it

Emad_B_ — Tue, 06 Oct 2015 16:16:25 GMT

Thanks Evgueni, that's did it.

topic 2D Convolution return wrong result in Intel® oneAPI Math Kernel Library

2D Convolution return wrong result

Hi Emad B.,

Thanks Evgueni, that's did it