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
---
 .../benchmarks/svm/src/matlab/script_run_profile.m | 113 +++++++++++++++++++++
 1 file changed, 113 insertions(+)
 create mode 100644 SD-VBS/benchmarks/svm/src/matlab/script_run_profile.m

(limited to 'SD-VBS/benchmarks/svm/src/matlab/script_run_profile.m')

diff --git a/SD-VBS/benchmarks/svm/src/matlab/script_run_profile.m b/SD-VBS/benchmarks/svm/src/matlab/script_run_profile.m
new file mode 100644
index 0000000..b84a9d8
--- /dev/null
+++ b/SD-VBS/benchmarks/svm/src/matlab/script_run_profile.m
@@ -0,0 +1,113 @@
+function script_run_profile(dataDir, resultDir, type, commonDir,toolDir)
+
+path(path, commonDir);
+Iterations = 2;
+
+file = [dataDir, '/d16trn_1.txt'];
+trn_1 = readFile(file);
+file = [dataDir, '/d16trn_2.txt'];
+trn_2 = readFile(file);
+file = [dataDir, '/d16tst_1.txt'];
+tst_1 = readFile(file);
+file = [dataDir, '/d16tst_2.txt'];
+tst_2 = readFile(file);
+
+trn_1 = d16trn_1;
+trn_2 = d16trn_2;
+tst_1 = d16tst_1;
+tst_2 = d16tst_2;
+
+    N = 100;
+    Ntst = 100;
+    Iterations = 10;
+
+    if(strcmp(type, 'test'))
+    N = 4;
+    Ntst = 4;
+    Iterations = 2;
+    elseif(strcmp(type, 'sim_fast'))
+    N = 20;
+    Ntst = 20;
+    Iterations = 2;
+    elseif(strcmp(type, 'sim'))
+    N = 16;
+    Ntst = 16;
+    Iterations = 8;
+    elseif(strcmp(type, 'sqcif'))
+    N = 60;
+    Ntst = 60;
+    Iterations = 6;
+    elseif(strcmp(type, 'qcif'))
+    N = 72;
+    Ntst = 72;
+    Iterations = 8;
+    elseif(strcmp(type, 'vga'))
+    N = 450;
+    Ntst = 450;
+    Iterations = 15;
+    elseif(strcmp(type, 'wuxga'))
+    N = 1000;
+    Ntst = 1000;
+    Iterations = 20;
+    end
+    
+fprintf(1,'Input size\t\t- (%dx%dx%d)\n', N, Ntst,Iterations);
+
+elapsed = zeros(1,2);
+%% Timing
+start = photonStartTiming;
+
+[a,b,C,d,dim,e,eps,a_result,b_result,X,tolerance,Y,ret] = getAlphaFromTrainSet(N, trn_1, trn_2, Iterations);
+
+Yoffset=zeros(Iterations,N);
+ 
+Xtst=usps_read_partial(tst_1, tst_2 ,0,1,Ntst/Iterations, Iterations);
+Ytst=usps_read_partial(tst_1, tst_2, 0,0,Ntst/Iterations, Iterations);
+ 
+ for i=1:Iterations
+     temp=usps_read_partial(trn_1,trn_2, i,0,N/Iterations, Iterations);
+     Yoffset(i,:)=transpose(temp);
+ end
+ 
+ 
+ %error=0;
+ result=zeros(Ntst,1);
+ for n=1:Ntst
+     maxs = 0;
+     s=zeros(Iterations,1);
+     for i=1:Iterations
+         for j=1:N
+             if a_result(i,j) > 0
+                 s(i,1) = s(i,1) +a_result(i,j)*Yoffset(i,j)*polynomial(3,X(j,:),Xtst(n,:), dim);        
+             end
+         end
+         s(i,1) = s(i,1) - b_result(i,1);
+         
+         if( s(i,1) > maxs)
+             maxs = s(i,1);
+         end
+     end
+     result(n,1) = maxs;
+ end
+ 
+ %% Timing
+stop = photonEndTiming;
+
+temp = photonReportTiming(start, stop);
+elapsed(1) = elapsed(1) + temp(1);
+elapsed(2) = elapsed(2) + temp(2);
+ 
+ %% Self check
+ fWriteMatrix(result, dataDir);
+ tol = 0.001;
+ ret = fSelfCheck(result,dataDir, tol);
+ if(ret == 0)
+     disp('Error in SVM');
+ end
+
+%% Timing
+photonPrintTiming(elapsed);
+
+
+
+
-- 
cgit v1.2.2