I'm wondering if anyone can help. I'm trying to calculate a 2d convolution of two square arrays f(x,y) and g(x,y) using the code below, which compiles and runs fine (no error messages), but does not give the correct output for h(x,y). Instead, all values of h[n] are unchanged (e.g. if initialized before attempting the convolution) except one "column" of the output, at x=0.
I've scoured the examples included with MKL libraries, but can't spot anything different between my code and the minimal working examples of this type.