summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/matrixOperations/kron/skrona.c
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/c/matrixOperations/kron/skrona.c')
-rw-r--r--2.3-1/src/c/matrixOperations/kron/skrona.c38
1 files changed, 38 insertions, 0 deletions
diff --git a/2.3-1/src/c/matrixOperations/kron/skrona.c b/2.3-1/src/c/matrixOperations/kron/skrona.c
new file mode 100644
index 00000000..368d04b5
--- /dev/null
+++ b/2.3-1/src/c/matrixOperations/kron/skrona.c
@@ -0,0 +1,38 @@
+// Copyright (C) 2016 - IIT Bombay - FOSSEE
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution. The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+// Author: Siddhesh Wani
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+
+//Function for kroneker product of two matrices
+
+#include "kron.h"
+
+void skrona (float *in1, int row1, int col1, float *in2, int row2, \
+ int col2, float *out)
+{
+ int row1_count, col1_count,row2_count, col2_count;
+ int row = row1*row2, col = col1*col2;
+ int temp = 0;
+
+ for(col1_count = 0;col1_count < col1; col1_count++)
+ {
+ for(row1_count = 0;row1_count < row1; row1_count++)
+ {
+ for(col2_count = 0;col2_count < col2; col2_count++)
+ {
+ for(row2_count = 0;row2_count < row2; row2_count++)
+ {
+ temp = (col1_count*col2+col2_count)*row+(row1_count*row2+row2_count);
+ out[temp] = in1[col1_count*row1+row1_count]*in2[col2_count*row2+row2_count];
+ }
+ }
+ }
+ }
+
+} \ No newline at end of file