[test-suite] Add Image Processing Kernels Using Benchmark Library: Anisotropic diffusion

[proton]

Summary

Runtime:
With Polly

Benchmark	Time	CPU	Iterations
BENCHMARK_ANISTROPIC_DIFFUSION/32	242 us	242 us	2832
BENCHMARK_ANISTROPIC_DIFFUSION/64	1120 us	1120 us	614
BENCHMARK_ANISTROPIC_DIFFUSION/128	4680 us	4680 us	131
BENCHMARK_ANISTROPIC_DIFFUSION/256	19263 us	19261 us	36

Without Polly

Benchmark	Time	CPU	Iterations
BENCHMARK_ANISTROPIC_DIFFUSION/32	2390 us	2383 us	292
BENCHMARK_ANISTROPIC_DIFFUSION/64	10961 us	10949 us	64
BENCHMARK_ANISTROPIC_DIFFUSION/128	46711 us	46707 us	15
BENCHMARK_ANISTROPIC_DIFFUSION/256	194217 us	194172 us	4

[Meinersbur]

Wow, that's 10x fast with Polly!

MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/anisotropicDiffusionKernel.c

10 ↗

(On Diff #157799)

[suggestion] Parenthesize to avoid precedence issues: #define TIME_STEP (1.0 / 8)

MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/main.cpp

41–42 ↗

(On Diff #157799)

[nit] Can be on one line: int *outputImage = (int*) malloc(...

[Meinersbur]

LGTM

test-suite/trunk/MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/CMakeLists.txt

+set(IMAGEPROC_UTILS MicroBenchmarks/ImageProcessing/utils)
+
+list(APPEND CPPFLAGS -I ${CMAKE_SOURCE_DIR}/${IMAGEPROC_UTILS} -std=c++11)
+
+llvm_test_verify("${CMAKE_SOURCE_DIR}/HashProgramOutput.sh ${CMAKE_CURRENT_BINARY_DIR}/anisotropicDiffusionOutput.txt")
+llvm_test_verify("${FPCMP} ${CMAKE_CURRENT_BINARY_DIR}/anisotropicDiffusionOutput.txt ${CMAKE_CURRENT_SOURCE_DIR}/anisotropicDiffusion.reference_output")
+
+llvm_test_run(WORKDIR ${CMAKE_CURRENT_BINARY_DIR})
+
+llvm_test_executable(AnisotropicDiffusion ../utils/ImageHelper.cpp ../utils/glibc_compat_rand.c main.cpp anisotropicDiffusionKernel.c)
+
+target_link_libraries(AnisotropicDiffusion benchmark)

test-suite/trunk/MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/anisotropicDiffusion.reference_output

+df287cfd3bd0d85a2bffd28b5bd37efd

test-suite/trunk/MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/anisotropicDiffusionKernel.c

+/*
+  Source:
+  https://cboard.cprogramming.com/c-programming/164080-anisotropic-diffusion-implementation-c.html
+  Modified by Pankaj Kukreja (github.com/proton0001)
+  Indian Institute of Technology Hyderabad
+*/
+#include "diffusion.h"
+#include <math.h>
+
+#define TIME_STEP (1.0 / 8)
+#define DIFFUSION_CONST 0.5
+void anisotropicDiffusionKernel(int height, int width,
+                                int inputImage[HEIGHT][WIDTH],
+                                int outputImage[height][width],
+                                int no_iteration) {
+  int i, j, k;
+  for (k = 0; k < no_iteration; k++) {
+    for (i = 2; i < height - 2; i++) {
+      for (j = 2; j < width - 2; j++) {
+        int gradient_x, gradient_y, grad_final;
+        float result;
+        float ne, nw, nn, ns;
+        gradient_x = inputImage[i][j + 2] - inputImage[i][j];
+        gradient_y = inputImage[i + 1][j + 1] - inputImage[i - 1][j + 1];
+        grad_final = gradient_x + gradient_y;
+        ne = exp(-(grad_final * grad_final) / DIFFUSION_CONST *
+                 DIFFUSION_CONST) *
+             (inputImage[i][j + 1] - inputImage[i][j]);
+
+        gradient_x = inputImage[i][j] - inputImage[i][j - 2];
+        gradient_y = inputImage[i + 1][j - 1] - inputImage[i - 1][j - 1];
+        grad_final = gradient_x + gradient_y;
+        nw = exp(-(grad_final * grad_final) / DIFFUSION_CONST *
+                 DIFFUSION_CONST) *
+             (inputImage[i][j - 1] - inputImage[i][j]);
+
+        gradient_x = inputImage[i - 1][j + 1] - inputImage[i - 1][j - 1];
+        gradient_y = inputImage[i][j] - inputImage[i - 2][j];
+        grad_final = gradient_x + gradient_y;
+        nn = exp(-(grad_final * grad_final) / DIFFUSION_CONST *
+                 DIFFUSION_CONST) *
+             (inputImage[i - 1][j] - inputImage[i][j]);
+
+        gradient_x = inputImage[i + 1][j + 1] - inputImage[i + 1][j - 1];
+        gradient_y = inputImage[i + 2][j] - inputImage[i][j];
+        grad_final = gradient_x + gradient_y;
+        ns = exp(-(grad_final * grad_final) / DIFFUSION_CONST *
+                 DIFFUSION_CONST) *
+             (inputImage[i + 1][j] - inputImage[i][j]);
+
+        outputImage[i][j] = inputImage[i][j] + TIME_STEP * (ne + nw + nn + ns);
+      }
+    }
+  }
+}

test-suite/trunk/MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/diffusion.h

+/**
+  Pankaj Kukreja
+  github.com/proton0001
+  Indian Institute of Technology Hyderabad
+*/
+#ifndef _DIFFUSSION_H
+#define _DIFFUSSION_H
+
+#define HEIGHT 512
+#define WIDTH 512
+
+#define ITERATION 10
+
+#endif /* _DIFFUSSION_H */

test-suite/trunk/MicroBenchmarks/ImageProcessing/AnisotropicDiffusion/main.cpp

+/**
+  Pankaj Kukreja
+  github.com/proton0001
+  Indian Institute of Technology Hyderabad
+*/
+
+#include "ImageHelper.h"
+#include "diffusion.h"
+#include <iostream> // std::cerr
+
+#define BENCHMARK_LIB
+#ifdef BENCHMARK_LIB
+#include "benchmark/benchmark.h"
+#endif
+
+int *inputImage;
+extern "C" void anisotropicDiffusionKernel(int height, int width, int *source,
+                                           int *destination, int num_iteration);
+int main(int argc, char *argv[]) {
+
+#ifdef BENCHMARK_LIB
+  ::benchmark::Initialize(&argc, argv);
+#endif
+
+  const char *anisotropicDiffusionOutputFileName =
+      (const char *)"./anisotropicDiffusionOutput.txt";
+
+  inputImage = (int *)malloc(sizeof(int) * (HEIGHT) * (WIDTH));
+
+  if (inputImage == NULL) {
+    std::cerr << "Insufficient memory\n";
+    exit(EXIT_FAILURE);
+  }
+
+  initializeRandomImage(inputImage, HEIGHT, WIDTH);
+
+#ifdef BENCHMARK_LIB
+  ::benchmark::RunSpecifiedBenchmarks();
+#endif
+
+  int *outputImage = (int *)malloc(sizeof(int) * (HEIGHT) * (WIDTH));
+  if (outputImage == NULL) {
+    std::cerr << "Insufficient memory\n";
+    exit(EXIT_FAILURE);
+  }
+
+  anisotropicDiffusionKernel(HEIGHT, WIDTH, inputImage, outputImage, ITERATION);
+
+  for (int j = 0; j < WIDTH; j++) {
+    outputImage[0 + j] = 0;
+    outputImage[1 * WIDTH + j] = 0;
+    outputImage[(HEIGHT - 1) * WIDTH + j] = 0;
+    outputImage[(HEIGHT - 2) * WIDTH + j] = 0;
+  }
+
+  for (int i = 0; i < HEIGHT; i++) {
+    outputImage[i * WIDTH + 0] = 0;
+    outputImage[i * WIDTH + 1] = 0;
+    outputImage[i * WIDTH + (WIDTH - 1)] = 0;
+    outputImage[i * WIDTH + (WIDTH - 2)] = 0;
+  }
+
+  saveImage(outputImage, anisotropicDiffusionOutputFileName, HEIGHT, WIDTH);
+
+  free(outputImage);
+  free(inputImage);
+
+  return EXIT_SUCCESS;
+}
+
+#ifdef BENCHMARK_LIB
+void BENCHMARK_ANISTROPIC_DIFFUSION(benchmark::State &state) {
+
+  int height = state.range(0);
+  int width = state.range(0);
+
+  int *outputImage = (int *)malloc(sizeof(int) * (height) * (width));
+  if (outputImage == NULL) {
+    std::cerr << "Insufficient memory\n";
+    exit(EXIT_FAILURE);
+  }
+  /* first call made to warm up cache */
+  anisotropicDiffusionKernel(height, width, inputImage, outputImage, ITERATION);
+
+  for (auto _ : state) {
+    anisotropicDiffusionKernel(height, width, inputImage, outputImage,
+                               ITERATION);
+  }
+  /* Since we are not passing state.range as 20 this if case will always be
+   * false.
+   * This call is to make compiler think that outputImage may be used
+   * later so that above kernel calls will not get optimize out */
+  if (state.range(0) == 20) {
+    saveImage(outputImage, (const char *)"testFailed.txt", height, width);
+  }
+
+  free(outputImage);
+}
+
+BENCHMARK(BENCHMARK_ANISTROPIC_DIFFUSION)
+    ->Arg(32)
+    ->Arg(64)
+    ->Arg(128)
+    ->Arg(256)
+    ->Unit(benchmark::kMicrosecond);
+#endif

test-suite/trunk/MicroBenchmarks/ImageProcessing/CMakeLists.txt

 add_subdirectory(Dither)
+add_subdirectory(AnisotropicDiffusion)