From f618466c25d43f3bae9e40920273bf77de1e1149 Mon Sep 17 00:00:00 2001
From: leochanj105 <leochanj@live.unc.edu>
Date: Mon, 19 Oct 2020 23:09:30 -0400
Subject: initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs
---
 SD-VBS/benchmarks/svm/src/c/polynomial.c | 42 ++++++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100644 SD-VBS/benchmarks/svm/src/c/polynomial.c

(limited to 'SD-VBS/benchmarks/svm/src/c/polynomial.c')

diff --git a/SD-VBS/benchmarks/svm/src/c/polynomial.c b/SD-VBS/benchmarks/svm/src/c/polynomial.c
new file mode 100644
index 0000000..cce9190
--- /dev/null
+++ b/SD-VBS/benchmarks/svm/src/c/polynomial.c
@@ -0,0 +1,42 @@
+/********************************
+Author: Sravanthi Kota Venkata
+********************************/
+
+#include "svm.h"
+
+float polynomial(int d, F2D* a, F2D* b, int dim)
+{
+    float ret;
+    F2D *bt, *bt1;
+    int i,j,r,c;
+
+    r = b->height;
+    c = b->width;
+
+    bt = fMallocHandle(c, r);
+
+    for(i=0; i<r; i++)
+    {
+        for(j=0; j<c; j++)
+        {
+            subsref(bt,j,i) = subsref(b,i,j);
+        }
+    }
+
+    bt1 = fMtimes(a, bt);
+    fFreeHandle(bt);
+    
+    if(bt1->height == 1 && bt1->width ==1)
+        ret = pow(asubsref(bt1,0),d)/dim;
+    else
+    {
+        fFreeHandle(bt1);
+        return -1;
+    }
+
+    fFreeHandle(bt1);
+    return ret;
+}
+
+
+
-- 
cgit v1.2.2