From d17b33131c14864bd1eae275f49a3f148e21cf29 Mon Sep 17 00:00:00 2001
From: Leo Chan <leochanj@live.unc.edu>
Date: Thu, 22 Oct 2020 01:53:21 -0400
Subject: Squashed commit of the sb-vbs branch.

Includes the SD-VBS benchmarks modified to:
- Use libextra to loop as realtime jobs
- Preallocate memory before starting their main computation
- Accept input via stdin instead of via argc

Does not include the SD-VBS matlab code.

Fixes libextra execution in LITMUS^RT.
---
 SD-VBS/common/c/imageBlur.c | 67 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 67 insertions(+)
 create mode 100644 SD-VBS/common/c/imageBlur.c

(limited to 'SD-VBS/common/c/imageBlur.c')

diff --git a/SD-VBS/common/c/imageBlur.c b/SD-VBS/common/c/imageBlur.c
new file mode 100644
index 0000000..8e3ad92
--- /dev/null
+++ b/SD-VBS/common/c/imageBlur.c
@@ -0,0 +1,67 @@
+/********************************
+Author: Sravanthi Kota Venkata
+********************************/
+
+#include "sdvbs_common.h"
+
+F2D* imageBlur(I2D* imageIn)
+{
+    int rows, cols;
+    F2D *imageOut, *tempOut;
+    float temp;
+    I2D *kernel;
+    int k, kernelSize, startCol, endCol, halfKernel, startRow, endRow, i, j, kernelSum;
+
+    rows = imageIn->height;
+    cols = imageIn->width;
+
+    imageOut = fSetArray(rows, cols, 0);
+    tempOut = fSetArray(rows, cols, 0);
+    kernel = iMallocHandle(1, 5);
+
+    asubsref(kernel,0) = 1;
+    asubsref(kernel,1) = 4;
+    asubsref(kernel,2) = 6;
+    asubsref(kernel,3) = 4;
+    asubsref(kernel,4) = 1;
+    kernelSize = 5;
+    kernelSum = 16;
+
+    startCol = 2;  
+    endCol = cols - 2;  
+    halfKernel = 2;   
+
+    startRow = 2;    
+    endRow = rows - 2;  
+
+    for(i=startRow; i<endRow; i++){
+        for(j=startCol; j<endCol; j++)
+        {
+            temp = 0;
+            for(k=-halfKernel; k<=halfKernel; k++)
+            {
+                temp += subsref(imageIn,i,j+k) * asubsref(kernel,k+halfKernel);
+            }
+            subsref(tempOut,i,j) = temp/kernelSum;
+        }
+    }
+    
+    for(i=startRow; i<endRow; i++)
+    {
+        for(j=startCol; j<endCol; j++)
+        {
+            temp = 0;
+            for(k=-halfKernel; k<=halfKernel; k++)
+            {
+                temp += subsref(tempOut,(i+k),j) * asubsref(kernel,k+halfKernel);
+            }
+            subsref(imageOut,i,j) = temp/kernelSum;
+        }
+    }
+
+    fFreeHandle(tempOut);
+    iFreeHandle(kernel);
+    return imageOut;
+}
+             
+
-- 
cgit v1.2.2