[test-suite] Update output image dimensions in the Harris benchmark

[Meinersbur]

Summary

On AIX, I realized that the output image matrix of the harris benchmark did not match the one produced on Linux.

For instance, on Linux, the image matrix started off as:

$ more output.txt
00000000000000000000000000000000000000000000000000000000000000 . . .

But on AIX, we were getting:

$ more output.txt
25525525525500000000002552552552552552552552552552550000000000 . . .

This patch updates the dimensions of the output image in the harris benchmark in order to
resolve the issue of printing the image on AIX.

[Meinersbur]

image is meant to be a pointer to a (HEIGHT+2)*(WIDTH+2) image, where each "pixel" is a represented as a grayscale float (4 bytes) value. sizeof(float*) * (2 + HEIGHT) * (2 + WIDTH) allocates a pointer-size for each pixel (8 bytes on a 64-bit system), i.e. typically twice as much. However, the data is not supposed to be a 2-dimensional array of pointers.

I cannot claim there is no out-of-bounds array access in the program, but the size computation is not the issue. Did you try to run it under MemorySanitizer or Valgrind?

[hubert.reinterpretcast]

In D103092#2780891, @Meinersbur wrote:

image is meant to be a pointer to a (HEIGHT+2)*(WIDTH+2) image, where each "pixel" is a represented as a grayscale float (4 bytes) value. sizeof(float*) * (2 + HEIGHT) * (2 + WIDTH) allocates a pointer-size for each pixel (8 bytes on a 64-bit system), i.e. typically twice as much. However, the data is not supposed to be a 2-dimensional array of pointers.

I agree the size change is not the right fix.

I cannot claim there is no out-of-bounds array access in the program, but the size computation is not the issue. Did you try to run it under MemorySanitizer or Valgrind?

What I can say is that using calloc seems to fix the issue:

-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+      std::calloc((2 + HEIGHT) * (2 + WIDTH), sizeof(float));

That does cast suspicion on the correctness of the source program.

[hubert.reinterpretcast]

The problem is fairly obvious once found.

The kernel only initializes part of the output image (most obviously not the initial rows/columns):

for (int _i0 = 2; (_i0 < height); _i0++) {
  for (int _i1 = 2; (_i1 < width); _i1++) {
    outputImg[_i0][_i1] =

The check output does print from the initial rows/columns:

void printImage(int height, int width, Image &arr,
                int dummy) {
  std::ofstream myfile;
  myfile.open("output.txt");
  for (int i = 0; i < height - 2; i++) {
    for (int j = 0; j < width - 2; j++) {
      if (arr[i][j] < 0) {

if (argc == 2) {
  printImage(HEIGHT + 2, WIDTH + 2, *imageOutput, sum);
} else {
  printImage(HEIGHT + 2, WIDTH + 2, *imageOutput, -1);
}

Likewise, the "sum" is calculated using additional rows and columns not initialized by the kernel:

for (int i = 0; i < height + 2; i++) {
  for (int j = 0; j < width + 2; j++) {
    sum = (sum + 1) & (int)(*imageOutput)[i][j];
  }
}

[Meinersbur]

Thank you @hubert.reinterpretcast for the analysis.

@amyk are you going to update that patch? Either bey use of calloc, or (preferably) fixing the array bounds of what is being checked.

In D103092#2786318, @hubert.reinterpretcast wrote:

The kernel only initializes part of the output image (most obviously not the initial rows/columns):

for (int _i0 = 2; (_i0 < height); _i0++) {
  for (int _i1 = 2; (_i1 < width); _i1++) {
    outputImg[_i0][_i1] =

Should this be

for (int _i0 = 1; (_i0 < height + 1); _i0++) {
  for (int _i1 = 1; (_i1 < width + 1); _i1++) {
    outputImg[_i0][_i1] =

(i.e. a 1-pixel margin left/right/top/bottom instead a 2-pixel margin left/top)? Alternatively, the output image doesn't even need the margin.

I remember having reviewed the original patch, seem I didn't pay enough attention.

[amyk]

Thanks @hubert.reinterpretcast and @Meinersbur! I was also further investigating this and I can confirm that changing the malloc calls to calloc resolves the issue for the harris benchmark.

I noticed that this change (a 1-pixel margin left/right/top/bottom), still did not provide the correct result.

for (int _i0 = 1; (_i0 < height + 1); _i0++) {
  for (int _i1 = 1; (_i1 < width + 1); _i1++) {
    outputImg[_i0][_i1] =

However, you did mention that the output image does not require a margin. Did you mean the following?

diff --git a/MicroBenchmarks/harris/harrisKernel.cpp b/MicroBenchmarks/harris/harrisKernel.cpp
index 0fb9c894..09d8aacc 100644
--- a/MicroBenchmarks/harris/harrisKernel.cpp
+++ b/MicroBenchmarks/harris/harrisKernel.cpp
@@ -103,8 +103,8 @@ void harrisKernel(
     }
   }

-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; (_i0 < height); _i0++) {
+    for (int _i1 = 0; (_i1 < width); _i1++) {
       outputImg[_i0][_i1] =
           (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
     }

If so, then this change works to produce the correct output on AIX, as well.

And yes, I will be updating the patch - either by calloc or the above change (if that is preferred).

[hubert.reinterpretcast]

In D103092#2786379, @amyk wrote:

And yes, I will be updating the patch - either by calloc or the above change (if that is preferred).

Instead of calloc, we can use new:

diff --git a/MicroBenchmarks/harris/main.cpp b/MicroBenchmarks/harris/main.cpp
index 6d43c930..f9833588 100644
--- a/MicroBenchmarks/harris/main.cpp
+++ b/MicroBenchmarks/harris/main.cpp
@@ -66,4 +66,3 @@ void BENCHMARK_HARRIS(benchmark::State &state) {
   float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
-  imageOutput = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+  imageOutput = new float[1][2 + HEIGHT][2 + WIDTH]();

@@ -129,3 +128,3 @@ void BENCHMARK_HARRIS(benchmark::State &state) {

-  free((void *)imageOutput);
+  delete [] imageOutput;
   free((void *)image);
@@ -157,4 +156,3 @@ int main(int argc, char *argv[]) {
   float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
-  imageOutput = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+  imageOutput = new float[1][2 + HEIGHT][2 + WIDTH]();

@@ -213,3 +211,3 @@ int main(int argc, char *argv[]) {
   free((void *)image);
-  free((void *)imageOutput);
+  delete [] imageOutput;
   return 0;

[Meinersbur]

In D103092#2786379, @amyk wrote:

However, you did mention that the output image does not require a margin. Did you mean the following?

diff --git a/MicroBenchmarks/harris/harrisKernel.cpp b/MicroBenchmarks/harris/harrisKernel.cpp
index 0fb9c894..09d8aacc 100644
--- a/MicroBenchmarks/harris/harrisKernel.cpp
+++ b/MicroBenchmarks/harris/harrisKernel.cpp
@@ -103,8 +103,8 @@ void harrisKernel(
     }
   }

-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; (_i0 < height); _i0++) {
+    for (int _i1 = 0; (_i1 < width); _i1++) {
       outputImg[_i0][_i1] =
           (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
     }

I don't have an overview on what Harris edge detection does with trace whether trace[0][0], trace[0][1] etc contains useful information. The output should also be a WIDTH*HEIGHT-sized image, i.e.

- float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
+ float(*imageOutput)[HEIGHT][WIDTH];

It seems that outputImg[WIDTH], outputImg[WIDTH+1][ etc is ignored anyway. Same with trace[WIDTH], trace[WIDTH+1]`, etc.

[amyk]

In D103092#2787051, @hubert.reinterpretcast wrote:

Instead of calloc, we can use new:

Thanks Hubert. I can con confirm that new also works for this test and we can possibly move forward with this, too.

In D103092#2787107, @Meinersbur wrote:`

I don't have an overview on what Harris edge detection does with trace whether trace[0][0], trace[0][1] etc contains useful information. The output should also be a WIDTH*HEIGHT-sized image, i.e.

- float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
+ float(*imageOutput)[HEIGHT][WIDTH];

It seems that outputImg[WIDTH], outputImg[WIDTH+1][ etc is ignored anyway. Same with trace[WIDTH], trace[WIDTH+1]`, etc.

@Meinersbur Sorry - just a clarification question. From your comment, did you mean that imageOutput/trace/det should be [HEIGHT][WIDTH] since the WIDTH and WIDTH+1, etc. is ignored anyway?

[amyk]

In D103092#2787107, @Meinersbur wrote:

In D103092#2786379, @amyk wrote:

However, you did mention that the output image does not require a margin. Did you mean the following?

diff --git a/MicroBenchmarks/harris/harrisKernel.cpp b/MicroBenchmarks/harris/harrisKernel.cpp
index 0fb9c894..09d8aacc 100644
--- a/MicroBenchmarks/harris/harrisKernel.cpp
+++ b/MicroBenchmarks/harris/harrisKernel.cpp
@@ -103,8 +103,8 @@ void harrisKernel(
     }
   }

-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; (_i0 < height); _i0++) {
+    for (int _i1 = 0; (_i1 < width); _i1++) {
       outputImg[_i0][_i1] =
           (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
     }

I don't have an overview on what Harris edge detection does with trace whether trace[0][0], trace[0][1] etc contains useful information. The output should also be a WIDTH*HEIGHT-sized image, i.e.

- float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
+ float(*imageOutput)[HEIGHT][WIDTH];

It seems that outputImg[WIDTH], outputImg[WIDTH+1][ etc is ignored anyway. Same with trace[WIDTH], trace[WIDTH+1]`, etc.

Additionally, do you mean that we should be doing this change, along with this change too:

diff --git a/MicroBenchmarks/harris/harrisKernel.cpp b/MicroBenchmarks/harris/harrisKernel.cpp
index 0fb9c894..09d8aacc 100644
--- a/MicroBenchmarks/harris/harrisKernel.cpp
+++ b/MicroBenchmarks/harris/harrisKernel.cpp
@@ -103,8 +103,8 @@ void harrisKernel(
     }
   }

-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; (_i0 < height); _i0++) {
+    for (int _i1 = 0; (_i1 < width); _i1++) {
       outputImg[_i0][_i1] =
           (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
     }

If we decide not to go with the new/delete or calloc approach?

[Meinersbur]

In D103092#2789477, @amyk wrote:

@Meinersbur Sorry - just a clarification question. From your comment, did you mean that imageOutput/trace/det should be [HEIGHT][WIDTH] since the WIDTH and WIDTH+1, etc. is ignored anyway?

Would be nice. However, for a minimal fix, just not printing the irrelevant border would be sufficient.

[amyk]

• Update the dimensions of outputImg and the loop indices for outputImg in the harris kernel

@Meinersbur It's possible that I may have misunderstood your comment, but could you please take another look at this patch?

[Meinersbur]

LGTM, thank you

MicroBenchmarks/harris/harrisKernel.cpp

3–4

[nit] remove parens around HEIGHT

There wasn't really a reason for parens before, but now its even more obvious

MicroBenchmarks/harris/main.cpp

32

[nit] remove parens around HEIGHT

[amyk]

In D103092#2809760, @Meinersbur wrote:

LGTM, thank you

Thanks for reviewing. My apologies in the delay on getting to this patch. I wanted to address the review comment this morning - upon realizing that I was not testing this patch with the latest revision of test-suite, it looks like this patch no longer addresses the issue on AIX.
It seems since the D101475 revision, the harris benchmark will fail on AIX even with this change.

The output image from this patch is slightly different with this patch (as seen in D103092-output.txt, although it still looks a more similar to the reference output when comparing with the output without this change.

$ diff reference-output.out D103092-output.out
10,12c10,13
< 0000000033330000003333000000 . . . 
< 0000000033330000003333000000 . . .
< 0000000033330000003333000000 . . . 
---
> 33000000333300000033330000003333 . . . 
> 33000000333300000033330000003333 . . . 
> 33000000333300000033330000003333 . . . 
> 33000000000000000000000000000000 . . . 

So, I think if I update the array bounds, this patch needs more work. As a side note, doing the new/delete or calloc approach will still produce the correct output image.

[Meinersbur]

If the sizes of the output changes a mismatch with the old reference output is expected (the removed border contains undefined values anyway). In that case, just replace the reference-output.out with D103092-output.txt.

[amyk]

In D103092#2830621, @Meinersbur wrote:

If the sizes of the output changes a mismatch with the old reference output is expected (the removed border contains undefined values anyway). In that case, just replace the reference-output.out with D103092-output.txt.

If it's okay, I just wanted to double check if my understanding is clear: your suggestion is to add a new reference output file for AIX to reflect these changes, right? (Since the Linux result still matches the original reference output)

Also, ever since the D101475 revision (updating the benchmark dependency of Microbenchmarks), the different output image is produced even just by printing the margin like so:

diff --git a/MicroBenchmarks/harris/harrisKernel.cpp b/MicroBenchmarks/harris/harrisKernel.cpp
index 0fb9c894..09d8aacc 100644
--- a/MicroBenchmarks/harris/harrisKernel.cpp
+++ b/MicroBenchmarks/harris/harrisKernel.cpp
@@ -103,8 +103,8 @@ void harrisKernel(
     }
   }

-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; (_i0 < height); _i0++) {
+    for (int _i1 = 0; (_i1 < width); _i1++) {
       outputImg[_i0][_i1] =
           (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
     }

Without changing the dimensions of outputImg to [HEIGHT][WIDTH]. Is adding a new reference output still okay in this situation?
I just wanted to make sure since your previous comment stated "If the sizes of the output changes a mismatch with the old reference output is expected", but the output is different even without changing the dimensions of outputImg.

[Meinersbur]

In D103092#2832411, @amyk wrote:

If it's okay, I just wanted to double check if my understanding is clear: your suggestion is to add a new reference output file for AIX to reflect these changes, right? (Since the Linux result still matches the original reference output)

How can the Linux output still be the same? harrisKernel assigns:

for (int _i0 = 2; (_i0 < height); _i0++) {
  for (int _i1 = 2; (_i1 < width); _i1++) {
    outputImg[_i0][_i1] = ...

printImage accesses:

for (int i = 0; i < height - 2; i++) {
  for (int j = 0; j < width - 2; j++) {
     ... arr[i][j] ...

There is a gap of the indices 0 and 1 that is never assigned, and the last two indices never printed. Output of Linux and AIX must be the same. It seems that the platform difference between Linux and AIX is that on Linux the memory is somehow initialized to zero. Please don't rely on that on either platform. There seems to be an off-by-two error here. Print all the values that have been assigned in harrisKernel, and only those.

[amyk]

In D103092#2841458, @Meinersbur wrote:

How can the Linux output still be the same? harrisKernel assigns:

for (int _i0 = 2; (_i0 < height); _i0++) {
  for (int _i1 = 2; (_i1 < width); _i1++) {
    outputImg[_i0][_i1] = ...

printImage accesses:

for (int i = 0; i < height - 2; i++) {
  for (int j = 0; j < width - 2; j++) {
     ... arr[i][j] ...

There is a gap of the indices 0 and 1 that is never assigned, and the last two indices never printed. Output of Linux and AIX must be the same. It seems that the platform difference between Linux and AIX is that on Linux the memory is somehow initialized to zero. Please don't rely on that on either platform. There seems to be an off-by-two error here. Print all the values that have been assigned in harrisKernel, and only those.

I apologize, I misunderstood your previous comment - I agree that the output of Linux and AIX should be the same.

In the current patch, harrisKernel doesn't assign outputImg starting at 2, since I changed it to start at 0:

for (int _i0 = 0; (_i0 < height); _i0++) {
  for (int _i1 = 0; (_i1 < width); _i1++) {
    outputImg[_i0][_i1] =
        (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
  }
}

Moreover, unless I am missing something, just by looking at the code, it looks like printImage does actually print the last two indices, despite going up to height - 2 and width - 2 because the arguments of height and width to printImage are HEIGHT + 2 and WIDTH + 2:

if (argc == 2) {
  printImage(HEIGHT + 2, WIDTH + 2, *imageOutput, sum);
} else {
  printImage(HEIGHT + 2, WIDTH + 2, *imageOutput, -1);
}

[amyk]

Since the revision of updating the benchmark dependency of Microbenchmarks affected this change I made to harris, I tried to redo this patch and noticed something.
If I start off with the change of making imageOutput to have the size of [HEIGHT][WIDTH] instead of [HEIGHT + 2][WIDTH + 2], the image output I get on Linux and AIX is almost the same, except for the first line that is printed:

$ diff AIX_output.txt Linux_output.txt
1c1
< 0000002552552552552552552552552552550000000000255 . . . 
> 0000000000000000000000000000000000000000000000000 . . . 

So I think this relates to the height of the image. Updating harrisKernel to assign outputImg from 0 until height resolves the issue of the first line being different.

However, I found that the I can no longer change the width to also start from 0.

-  for (int _i0 = 2; (_i0 < height); _i0++) {
+  for (int _i0 = 0; (_i0 < height); _i0++) {
     for (int _i1 = 2; (_i1 < width); _i1++) {
       outputImg[_i0][_i1] =
           (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
     }
   }
 }

It looks like currently, width needs to stay the same, or there is more variation in the lines that are produced.

$ diff AIX_output.txt Linux_output.txt
10,13c10,12
< 3300000033330000003333000000333300000033330000003333000000 . . . 
< 3300000033330000003333000000333300000033330000003333000000 . . . 
< 3300000033330000003333000000333300000033330000003333000000 . . .
< 3300000000000000000000000000000000000000000000000000000000 . . . 
---
> 0000000033330000003333000000333300000033330000003333000000 . . .
> 0000000033330000003333000000333300000033330000003333000000 . . . 
> 0000000033330000003333000000333300000033330000003333000000 . . .
. . . 

I can update the patch based on this change that I've observed.

[amyk]

Update patch to:

• change imageOutput to have a size of [HEIGHT][WIDTH].
• Remove unnecessary parentheses around HEIGHT.
• Update harrisKernel to assign height from 0.

This updated revision reflects the change that is discussed in the comment above.

[Meinersbur]

• Fix use of unitialized memory in MicroBenchmarks/harris

[Meinersbur]

Since we did not seem to advance here, I took the initiative and fixed the test according to what I had in mind.

Changes:

• Increase border size of input image to 2 in each direction
• Remove unused border of Sxx,Syy,Sxz,det,trace
• Normalize loop bounds to start with 0
• Remove superfluous parentheses

Valgind memcheck, AddressSanitizer and MemorySanitizer(*)-clean

(*) msan still complains about use of initialized memory in Google Benchmark

[Meinersbur]

• Remove some more unnecessary parens

[amyk]

In D103092#2869208, @Meinersbur wrote:

Since we did not seem to advance here, I took the initiative and fixed the test according to what I had in mind.

Changes:

• Increase border size of input image to 2 in each direction
• Remove unused border of Sxx,Syy,Sxz,det,trace
• Normalize loop bounds to start with 0
• Remove superfluous parentheses

Valgind memcheck, AddressSanitizer and MemorySanitizer(*)-clean

(*) msan still complains about use of initialized memory in Google Benchmark

Hi, I'd like to apologize for previously misunderstanding your review comments earlier, and sorry for a lot of back and forth on this patch.
Nonetheless, thank you very much for updating the patch, I greatly appreciate it.

I've applied this patch and this also runs successfully on AIX.

Sorry if this is a silly question, but may I ask behind the justification of increasing the border size for the input image?

[Meinersbur]

In D103092#2870379, @amyk wrote:

Sorry if this is a silly question, but may I ask behind the justification of increasing the border size for the input image?

Computing Ix and Iy each uses 1 neighbor in east/west respectively north/south direction from the original image, which requires the input to be +2 larger than the output. In application where input/output must be the same size, it is common to mirror the pixels along the image borders. Here, we just ask initCheckboardImage to create a larger image.

Similarly computing Iyy, Sxy, Sxx, Syy requires a larger border from Ix and Iy. Both accumulate to a +4 for both stencils.

[amyk]

In D103092#2871419, @Meinersbur wrote:

In D103092#2870379, @amyk wrote:

Sorry if this is a silly question, but may I ask behind the justification of increasing the border size for the input image?

Computing Ix and Iy each uses 1 neighbor in east/west respectively north/south direction from the original image, which requires the input to be +2 larger than the output. In application where input/output must be the same size, it is common to mirror the pixels along the image borders. Here, we just ask initCheckboardImage to create a larger image.

Similarly computing Iyy, Sxy, Sxx, Syy requires a larger border from Ix and Iy. Both accumulate to a +4 for both stencils.

Ah, I see. Thanks a lot for the explanation.
I would also like to ask would it make sense for you to commandeer this revision since you've updated the patch?

[Meinersbur]

OK. I took over the patch. Please try it on AIX and LGTM it if you think it is good.

[amyk]

In D103092#2872718, @Meinersbur wrote:

OK. I took over the patch. Please try it on AIX and LGTM it if you think it is good.

I've tried this on AIX and the output I'm getting matches the reference output. Thank you again for your assistance on this patch.
LGTM.

MicroBenchmarks/harris/harris.h

 #ifndef __HARRIS_H__
 #define __HARRIS_H__
 
 #include <cstdlib>
 #include <cstring>
 #include <fstream>
 #include <iomanip>
 #include <iostream>
 #include <string>
 
 // ============================================================================
 
 // This parameter is used to initialize the input image
 #define BOX_SIZE 10
 
 /* Comment this to not use google benchmark library */
-#define BENCHMARK_LIB 
+#define BENCHMARK_LIB
 
 // Image Size
-#define HEIGHT 2050
-#define WIDTH  2050
+#define HEIGHT 2048
+#define WIDTH 2048
 
 // Parameters For harris kernel
-#define THRESHOLD 0.1 
+#define THRESHOLD 0.1
 
 // ============================================================================
 
 // Initialize a checkboard image
-void initCheckboardImage(int height, int width); 
+void initCheckboardImage(int height, int width);
 // print image to output.txt
-void printImage(int height, int width, float img[(2 + HEIGHT)][2 + WIDTH]);
+void printImage(int height, int width, float img[HEIGHT][WIDTH]);
 
 // harris kernel from polymage_naive.cpp
-void harrisKernel(int height 
-                  , int width 
-                  , float inputImg[2 + HEIGHT][2 + WIDTH] 
-                  , float outputImg[(2 + HEIGHT)][2 + WIDTH] 
-                  , float Ix [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Iy [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Ixx [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Ixy [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Iyy [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Sxx [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Sxy [(2 + HEIGHT)][2 + WIDTH] 
-                  , float Syy [(2 + HEIGHT)][2 + WIDTH] 
-                  , float det [(2 + HEIGHT)][2 + WIDTH] 
-                  , float trace [(2 + HEIGHT)][2 + WIDTH]);
+void harrisKernel(int height, int width, float inputImg[4 + HEIGHT][4 + WIDTH],
+                  float outputImg[HEIGHT][WIDTH],
+                  float Ix[2 + HEIGHT][2 + WIDTH],
+                  float Iy[2 + HEIGHT][2 + WIDTH],
+                  float Ixx[2 + HEIGHT][2 + WIDTH],
+                  float Ixy[2 + HEIGHT][2 + WIDTH],
+                  float Iyy[2 + HEIGHT][2 + WIDTH], float Sxx[HEIGHT][WIDTH],
+                  float Sxy[HEIGHT][WIDTH], float Syy[HEIGHT][WIDTH],
+                  float det[HEIGHT][WIDTH], float trace[HEIGHT][WIDTH]);
 
 // ============================================================================
 #endif

MicroBenchmarks/harris/harris.reference_output

-8e390d93bb5942d09e958f5c791c94ad
+36d1e6a1490ed3da2e3c323cb57852ba

MicroBenchmarks/harris/harrisKernel.cpp

 #include "harris.h"
 
 // harris kernel from polymage_naive.cpp
-void harrisKernel(
-    int height, int width, float inputImg[2 + HEIGHT][2 + WIDTH],
-    float outputImg[(2 + HEIGHT)][2 + WIDTH], float Ix[(2 + HEIGHT)][2 + WIDTH],
-    float Iy[(2 + HEIGHT)][2 + WIDTH], float Ixx[(2 + HEIGHT)][2 + WIDTH],
-    float Ixy[(2 + HEIGHT)][2 + WIDTH], float Iyy[(2 + HEIGHT)][2 + WIDTH],
-    float Sxx[(2 + HEIGHT)][2 + WIDTH], float Sxy[(2 + HEIGHT)][2 + WIDTH],
-    float Syy[(2 + HEIGHT)][2 + WIDTH], float det[(2 + HEIGHT)][2 + WIDTH],
-    float trace[(2 + HEIGHT)][2 + WIDTH]) {
-  for (int _i0 = 0; (_i0 - height < 0); _i0++) {
-    for (int _i1 = 0; (_i1 - width < 0); _i1++) {
-      Iy[_i0 + 1][_i1 + 1] =
-          (((((((inputImg[_i0][_i1]) * -0.0833333333333f) +
-               ((inputImg[_i0][_i1 + 2]) * 0.0833333333333f)) +
-              ((inputImg[_i0 + 1][_i1]) * -0.166666666667f)) +
-             ((inputImg[_i0 + 1][_i1 + 2]) * 0.166666666667f)) +
-            ((inputImg[_i0 + 2][_i1]) * -0.0833333333333f)) +
-           ((inputImg[_i0 + 2][_i1 + 2]) * 0.0833333333333f));
+void harrisKernel(int height, int width, float inputImg[4 + HEIGHT][4 + WIDTH],
+                  float outputImg[HEIGHT][WIDTH],
+                  float Ix[2 + HEIGHT][2 + WIDTH],
+                  float Iy[2 + HEIGHT][2 + WIDTH],
+                  float Ixx[2 + HEIGHT][2 + WIDTH],
+                  float Ixy[2 + HEIGHT][2 + WIDTH],
+                  float Iyy[2 + HEIGHT][2 + WIDTH], float Sxx[HEIGHT][WIDTH],
+                  float Sxy[HEIGHT][WIDTH], float Syy[HEIGHT][WIDTH],
+                  float det[HEIGHT][WIDTH], float trace[HEIGHT][WIDTH]) {
+  for (int _i0 = 0; _i0 < height + 2; _i0++) {
+    for (int _i1 = 0; _i1 < width + 2; _i1++) {
+      Iy[_i0][_i1] = inputImg[_i0][_i1] * -0.0833333333333f +
+                     inputImg[_i0][_i1 + 2] * 0.0833333333333f +
+                     inputImg[_i0 + 1][_i1] * -0.166666666667f +
+                     inputImg[_i0 + 1][_i1 + 2] * 0.166666666667f +
+                     inputImg[_i0 + 2][_i1] * -0.0833333333333f +
+                     inputImg[_i0 + 2][_i1 + 2] * 0.0833333333333f;
     }
   }
 
-  for (int _i0 = 0; (_i0 - height < 0); _i0++) {
-    for (int _i1 = 0; (_i1 - width < 0); _i1++) {
-      Ix[_i0 + 1][_i1 + 1] =
-          ((((((inputImg[_i0][_i1] * -0.0833333333333f) +
-               (inputImg[_i0 + 2][_i1] * 0.0833333333333f)) +
-              (inputImg[_i0][_i1 + 1] * -0.166666666667f)) +
-             (inputImg[_i0 + 2][_i1 + 1] * 0.166666666667f)) +
-            (inputImg[_i0][_i1 + 2] * -0.0833333333333f)) +
-           (inputImg[_i0 + 2][_i1 + 2] * 0.0833333333333f));
+  for (int _i0 = 0; _i0 < height + 2; _i0++) {
+    for (int _i1 = 0; _i1 < width + 2; _i1++) {
+      Ix[_i0][_i1] = inputImg[_i0][_i1] * -0.0833333333333f +
+                     inputImg[_i0 + 2][_i1] * 0.0833333333333f +
+                     inputImg[_i0][_i1 + 1] * -0.166666666667f +
+                     inputImg[_i0 + 2][_i1 + 1] * 0.166666666667f +
+                     inputImg[_i0][_i1 + 2] * -0.0833333333333f +
+                     inputImg[_i0 + 2][_i1 + 2] * 0.0833333333333f;
     }
   }
 
-  for (int _i0 = 1; (_i0 - height - 1 < 0); _i0++) {
-    for (int _i1 = 1; (_i1 - width - 1 < 0); _i1++) {
+  for (int _i0 = 0; _i0 < height + 2; _i0++) {
+    for (int _i1 = 0; _i1 < width + 2; _i1++) {
       Iyy[_i0][_i1] = Iy[_i0][_i1] * Iy[_i0][_i1];
     }
   }
 
-  for (int _i0 = 1; (_i0 - height - 1 < 0); _i0++) {
-    for (int _i1 = 1; (_i1 - width - 1 < 0); _i1++) {
+  for (int _i0 = 0; _i0 < height + 2; _i0++) {
+    for (int _i1 = 0; _i1 < width + 2; _i1++) {
       Ixy[_i0][_i1] = Ix[_i0][_i1] * Iy[_i0][_i1];
     }
   }
 
-  for (int _i0 = 1; (_i0 - height - 1 < 0); _i0++) {
-    for (int _i1 = 1; (_i1 - width - 1 < 0); _i1++) {
+  for (int _i0 = 0; _i0 < height + 2; _i0++) {
+    for (int _i1 = 0; _i1 < width + 2; _i1++) {
       Ixx[_i0][_i1] = Ix[_i0][_i1] * Ix[_i0][_i1];
     }
   }
 
-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
-      Syy[_i0][_i1] = ((((((((Iyy[-1 + _i0][-1 + _i1] + Iyy[-1 + _i0][_i1]) +
-                             Iyy[-1 + _i0][1 + _i1]) +
-                            Iyy[_i0][-1 + _i1]) +
-                           Iyy[_i0][_i1]) +
-                          Iyy[_i0][1 + _i1]) +
-                         Iyy[1 + _i0][-1 + _i1]) +
-                        Iyy[1 + _i0][_i1]) +
-                       Iyy[1 + _i0][1 + _i1]);
+  for (int _i0 = 0; _i0 < height; _i0++) {
+    for (int _i1 = 0; _i1 < width; _i1++) {
+      Syy[_i0][_i1] =
+          Iyy[_i0 + 1][_i1] + Iyy[_i0][_i1 + 1] + Iyy[_i0][_i1 + 2] +
+          Iyy[_i0 + 1][_i1] + Iyy[_i0 + 1][_i1 + 1] + Iyy[_i0 + 1][_i1 + 2] +
+          Iyy[_i0 + 2][_i1] + Iyy[_i0 + 2][_i1 + 1] + Iyy[_i0 + 2][_i1 + 2];
     }
   }
 
-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
-      Sxy[_i0][_i1] = ((((((((Ixy[-1 + _i0][-1 + _i1] + Ixy[-1 + _i0][_i1]) +
-                             Ixy[-1 + _i0][1 + _i1]) +
-                            Ixy[_i0][-1 + _i1]) +
-                           Ixy[_i0][_i1]) +
-                          Ixy[_i0][1 + _i1]) +
-                         Ixy[1 + _i0][-1 + _i1]) +
-                        Ixy[1 + _i0][_i1]) +
-                       Ixy[1 + _i0][1 + _i1]);
+  for (int _i0 = 0; _i0 < height; _i0++) {
+    for (int _i1 = 0; _i1 < width; _i1++) {
+      Sxy[_i0][_i1] = Ixy[_i0][_i1] + Ixy[_i0][_i1 + 1] + Ixy[_i0][_i1 + 2] +
+                      Ixy[_i0 + 1][_i1] + Ixy[_i0 + 1][_i1 + 1] +
+                      Ixy[_i0 + 1][_i1 + 2] + Ixy[_i0 + 2][_i1] +
+                      Ixy[_i0 + 2][_i1 + 1] + Ixy[_i0 + 2][_i1 + 2];
     }
   }
 
-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
-      Sxx[_i0][_i1] = ((((((((Ixx[-1 + _i0][-1 + _i1] + Ixx[-1 + _i0][_i1]) +
-                             Ixx[-1 + _i0][1 + _i1]) +
-                            Ixx[_i0][-1 + _i1]) +
-                           Ixx[_i0][_i1]) +
-                          Ixx[_i0][1 + _i1]) +
-                         Ixx[1 + _i0][-1 + _i1]) +
-                        Ixx[1 + _i0][_i1]) +
-                       Ixx[1 + _i0][1 + _i1]);
+  for (int _i0 = 0; _i0 < height; _i0++) {
+    for (int _i1 = 0; _i1 < width; _i1++) {
+      Sxx[_i0][_i1] = Ixx[_i0][_i1] + Ixx[_i0][_i1 + 1] + Ixx[_i0][_i1 + 2] +
+                      Ixx[_i0 + 1][_i1] + Ixx[_i0 + 1][_i1 + 1] +
+                      Ixx[_i0 + 1][_i1 + 2] + Ixx[_i0 + 2][_i1] +
+                      Ixx[_i0 + 2][_i1 + 1] + Ixx[_i0 + 2][_i1 + 2];
     }
   }
 
-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
-      trace[_i0][_i1] = (Sxx[_i0][_i1] + Syy[_i0][_i1]);
+  for (int _i0 = 0; _i0 < height; _i0++) {
+    for (int _i1 = 0; _i1 < width; _i1++) {
+      trace[_i0][_i1] = Sxx[_i0][_i1] + Syy[_i0][_i1];
     }
   }
 
-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; _i0 < height; _i0++) {
+    for (int _i1 = 0; _i1 < width; _i1++) {
       det[_i0][_i1] =
-          ((Sxx[_i0][_i1] * Syy[_i0][_i1]) - (Sxy[_i0][_i1] * Sxy[_i0][_i1]));
+          Sxx[_i0][_i1] * Syy[_i0][_i1] - Sxy[_i0][_i1] * Sxy[_i0][_i1];
     }
   }
 
-  for (int _i0 = 2; (_i0 < height); _i0++) {
-    for (int _i1 = 2; (_i1 < width); _i1++) {
+  for (int _i0 = 0; _i0 < height; _i0++) {
+    for (int _i1 = 0; _i1 < width; _i1++) {
       outputImg[_i0][_i1] =
-          (det[_i0][_i1] - ((0.04f * trace[_i0][_i1]) * trace[_i0][_i1]));
+          det[_i0][_i1] - 0.04f * trace[_i0][_i1] * trace[_i0][_i1];
     }
   }
 }

MicroBenchmarks/harris/main.cpp

 #include "harris.h"
 int sum = 0;
 
 #ifdef BENCHMARK_LIB
 #include "benchmark/benchmark.h"
 #endif
 
 // This function initializes the input image to checkbox image
 void initCheckboardImage(int height, int width,
-                         float image[(2 + HEIGHT)][2 + WIDTH]) {
+                         float image[4 + HEIGHT][4 + WIDTH]) {
   int last_pixel_x = 0;
   int last_pixel_y = 0;
   for (int i = 0; i < height; i++) {
     if (i % BOX_SIZE == 0) {
       last_pixel_y = (last_pixel_y + 1) % 2;
     }
     last_pixel_x = last_pixel_y;
     for (int j = 0; j < width; j++) {
       if (j % BOX_SIZE == 0) {
         last_pixel_x = (last_pixel_x + 1) % 2;
       }
       if (last_pixel_x == 0) {
         image[i][j] = 255;
       } else {
         image[i][j] = 0;
       }
     }
   }
 }
 
 // Writes image matrix to a file.
-void printImage(int height, int width, float arr[(2 + HEIGHT)][2 + WIDTH],
+void printImage(int height, int width, float arr[HEIGHT][WIDTH],
                 int dummy) {
   std::ofstream myfile;
   myfile.open("output.txt");
-  for (int i = 0; i < height - 2; i++) {
-    for (int j = 0; j < width - 2; j++) {
+  for (int i = 0; i < height; i++) {
+    for (int j = 0; j < width; j++) {
       if (arr[i][j] < 0) {
         myfile << 0;
       } else if (arr[i][j] > 255) {
         myfile << 3;
       } else {
         myfile << (int)(arr[i][j]);
       }
     }
     myfile << "\n";
   }
 
   // Dummy code to make sure the allocated ImageOutput Array is not optimized
   // out
   if (dummy > 0) {
     myfile << sum;
   }
 }
 
 #ifdef BENCHMARK_LIB
 void BENCHMARK_HARRIS(benchmark::State &state) {
   int height = state.range(0);
   int width = state.range(1);
 
-  float(*image)[HEIGHT + 2][WIDTH + 2];
-  image = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  initCheckboardImage((HEIGHT + 2), (WIDTH + 2), *image);
-
-  float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
-  imageOutput = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+  float(*image)[HEIGHT + 4][WIDTH + 4];
+  image = (float(*)[HEIGHT + 4][WIDTH + 4])
+      malloc(sizeof(float) * (4 + HEIGHT) * (4 + WIDTH));
+  initCheckboardImage(HEIGHT + 4, WIDTH + 4, *image);
+
+  float(*imageOutput)[HEIGHT][WIDTH];
+  imageOutput = (float(*)[HEIGHT][WIDTH])
+      malloc(sizeof(float) * HEIGHT * WIDTH);
 
   float(*Ix)[2 + HEIGHT][2 + WIDTH];
   float(*Iy)[2 + HEIGHT][2 + WIDTH];
   float(*Ixx)[2 + HEIGHT][2 + WIDTH];
   float(*Ixy)[2 + HEIGHT][2 + WIDTH];
   float(*Iyy)[2 + HEIGHT][2 + WIDTH];
-  float(*Sxx)[2 + HEIGHT][2 + WIDTH];
-  float(*Sxy)[2 + HEIGHT][2 + WIDTH];
-  float(*Syy)[2 + HEIGHT][2 + WIDTH];
-  float(*det)[2 + HEIGHT][2 + WIDTH];
-  float(*trace)[2 + HEIGHT][2 + WIDTH];
+  float(*Sxx)[HEIGHT][WIDTH];
+  float(*Sxy)[HEIGHT][WIDTH];
+  float(*Syy)[HEIGHT][WIDTH];
+  float(*det)[HEIGHT][WIDTH];
+  float(*trace)[HEIGHT][WIDTH];
 
   Ix = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Iy = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Ixx = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Ixy = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Iyy = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  Sxx = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  Sxy = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  Syy = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  det = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  trace = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+  Sxx = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  Sxy = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  Syy = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  det = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  trace = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
 
   harrisKernel(height, width, *image, *imageOutput, *Ix, *Iy, *Ixx, *Ixy, *Iyy,
                *Sxx, *Sxy, *Syy, *det, *trace);
 
   for (auto _ : state) {
     harrisKernel(height, width, *image, *imageOutput, *Ix, *Iy, *Ixx, *Ixy,
                  *Iyy, *Sxx, *Sxy, *Syy, *det, *trace);
   }
 
   free((void *)Ix);
   free((void *)Iy);
   free((void *)Ixx);
   free((void *)Ixy);
   free((void *)Iyy);
   free((void *)Sxx);
   free((void *)Sxy);
   free((void *)Syy);
   free((void *)det);
   free((void *)trace);
 
-  for (int i = 0; i < height + 2; i++) {
-    for (int j = 0; j < width + 2; j++) {
+  for (int i = 0; i < height; i++) {
+    for (int j = 0; j < width; j++) {
       sum = (sum + 1) & (int)(*imageOutput)[i][j];
     }
   }
 
-  state.SetBytesProcessed(sizeof(float) * (height + 2) * (width + 2) *
-                          state.iterations());
+  state.SetBytesProcessed(sizeof(float) * height * width * state.iterations());
 
   free((void *)imageOutput);
   free((void *)image);
 }
 BENCHMARK(BENCHMARK_HARRIS)
     ->Args({256, 256})
     ->Args({512, 512})
     ->Args({1024, 1024})
     ->Args({2048, 2048})
     ->Unit(benchmark::kMicrosecond);
 
 #endif
 
 int main(int argc, char *argv[]) {
   sum = 1;
 #ifdef BENCHMARK_LIB
   ::benchmark::Initialize(&argc, argv);
   if (::benchmark::ReportUnrecognizedArguments(argc, argv))
     return 1;
   ::benchmark::RunSpecifiedBenchmarks();
 #endif
 
   // Extra Call to verify output of kernel
-  float(*image)[HEIGHT + 2][WIDTH + 2];
-  image = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  initCheckboardImage((HEIGHT + 2), (WIDTH + 2), *image);
-
-  float(*imageOutput)[2 + HEIGHT][2 + WIDTH];
-  imageOutput = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+  float(*image)[HEIGHT + 4][WIDTH + 4];
+  image = (float(*)[HEIGHT + 4][WIDTH + 4])
+      malloc(sizeof(float) * (HEIGHT + 4) * (WIDTH + 4));
+  initCheckboardImage(HEIGHT + 4, WIDTH + 4, *image);
+
+  float(*imageOutput)[HEIGHT][WIDTH];
+  imageOutput = (float(*)[HEIGHT][WIDTH])
+      malloc(sizeof(float) * HEIGHT * WIDTH);
 
   float(*Ix)[2 + HEIGHT][2 + WIDTH];
   float(*Iy)[2 + HEIGHT][2 + WIDTH];
   float(*Ixx)[2 + HEIGHT][2 + WIDTH];
   float(*Ixy)[2 + HEIGHT][2 + WIDTH];
   float(*Iyy)[2 + HEIGHT][2 + WIDTH];
-  float(*Sxx)[2 + HEIGHT][2 + WIDTH];
-  float(*Sxy)[2 + HEIGHT][2 + WIDTH];
-  float(*Syy)[2 + HEIGHT][2 + WIDTH];
-  float(*det)[2 + HEIGHT][2 + WIDTH];
-  float(*trace)[2 + HEIGHT][2 + WIDTH];
+  float(*Sxx)[HEIGHT][WIDTH];
+  float(*Sxy)[HEIGHT][WIDTH];
+  float(*Syy)[HEIGHT][WIDTH];
+  float(*det)[HEIGHT][WIDTH];
+  float(*trace)[HEIGHT][WIDTH];
 
   Ix = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Iy = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Ixx = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Ixy = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
   Iyy = (float(*)[2 + HEIGHT][2 + WIDTH])
       malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  Sxx = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  Sxy = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  Syy = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  det = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
-  trace = (float(*)[2 + HEIGHT][2 + WIDTH])
-      malloc(sizeof(float) * (2 + HEIGHT) * (2 + WIDTH));
+  Sxx = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  Sxy = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  Syy = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  det = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
+  trace = (float(*)[HEIGHT][WIDTH])malloc(sizeof(float) * HEIGHT * WIDTH);
 
   harrisKernel(HEIGHT, WIDTH, *image, *imageOutput, *Ix, *Iy, *Ixx, *Ixy, *Iyy,
                *Sxx, *Sxy, *Syy, *det, *trace);
 
   free((void *)Ix);
   free((void *)Iy);
   free((void *)Ixx);
   free((void *)Ixy);
   free((void *)Iyy);
   free((void *)Sxx);
   free((void *)Sxy);
   free((void *)Syy);
   free((void *)det);
   free((void *)trace);
 
   if (argc == 2) {
-    printImage(HEIGHT + 2, WIDTH + 2, *imageOutput, sum);
+    printImage(HEIGHT, WIDTH, *imageOutput, sum);
   } else {
-    printImage(HEIGHT + 2, WIDTH + 2, *imageOutput, -1);
+    printImage(HEIGHT, WIDTH, *imageOutput, -1);
   }
 
   free((void *)image);
   free((void *)imageOutput);
   return 0;
 }