1 files changed, 34 insertions, 0 deletions
diff --git a/src/c/matrixOperations/cross/i16crossa.c b/src/c/matrixOperations/cross/i16crossa.c
new file mode 100644
index 0000000..91cdb7f
--- /dev/null
+++ b/src/c/matrixOperations/cross/i16crossa.c
@@ -0,0 +1,34 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+#include "int16.h"
+
+void i16crossa (int16* inp1, int size1_r, int size1_c, int16* inp2, int size2_r, int size2_c ,int16* oup)
+{
+    int i,j;
+
+    if (size1_c == 3 && size1_r == 1)               // Calculating the cross product for a 1X3 Matrix
+    {
+        oup[0] = inp1[1]*inp2[2] - inp1[2]*inp2[1];
+        oup[1] = inp1[2]*inp2[0] - inp1[0]*inp2[2];
+        oup[2] = inp1[0]*inp2[1] - inp1[1]*inp2[0];
+        printf("%lf ", oup[0]);
+        printf("%lf ", oup[1]);
+        printf("%lf \n", oup[2]);
+    }
+    else					// Calculating the cross product for a 3XN Matrix
+    {
+        for(i = 0; i < size1_r; i++)
+        {
+            for(j = 0; j < size1_c; j = j+3)
+            {
+               oup[j+(i*size1_r)] = inp1[(j+1)+(i*size1_r)] * inp2[(j+2)+(i*size1_r)] - inp1[(j+2)+(i*size1_r)] * inp2[(j+1)+(i*size1_r)];
+               oup[(j+1)+(i*size1_r)] = inp1[(j+2)+(i*size1_r)] * inp2[(j+0)+(i*size1_r)] - inp1[(j+0)+(i*size1_r)] * inp2[(j+2)+(i*size1_r)];
+               oup[(j+2)+(i*size1_r)] = inp1[(j+0)+(i*size1_r)] * inp2[(j+1)+(i*size1_r)] - inp1[(j+1)+(i*size1_r)] * inp2[(j+0)+(i*size1_r)];
+
+            }
+        }
+
+    }
+}
+