Cross Added

10 files changed, 348 insertions, 0 deletions

diff --git a/src/c/matrixOperations/cross/ccrossa.c b/src/c/matrixOperations/cross/ccrossa.c
new file mode 100644
index 0000000..5f74bb3
--- /dev/null
+++ b/src/c/matrixOperations/cross/ccrossa.c
@@ -0,0 +1,37 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+#include "floatComplex.h"
+#include "stdlib.h"
+#include "string.h"
+#include "cat.h"
+
+void ccrossa (floatComplex* inp1, int size1_r, int size1_c, floatComplex* inp2, int size2_r, int size2_c ,floatComplex* 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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/cross/dcrossa.c b/src/c/matrixOperations/cross/dcrossa.c
new file mode 100644
index 0000000..203ec24
--- /dev/null
+++ b/src/c/matrixOperations/cross/dcrossa.c
@@ -0,0 +1,33 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+
+void dcrossa (double* inp1, int size1_r, int size1_c, double* inp2, int size2_r, int size2_c ,double* 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)];
+
+            }
+        }
+
+    }
+}
+
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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/cross/i8crossa.c b/src/c/matrixOperations/cross/i8crossa.c
new file mode 100644
index 0000000..adfbd53
--- /dev/null
+++ b/src/c/matrixOperations/cross/i8crossa.c
@@ -0,0 +1,34 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+#include "int8.h"
+
+void i8crossa (int8* inp1, int size1_r, int size1_c, int8* inp2, int size2_r, int size2_c ,int8* 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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/cross/scrossa.c b/src/c/matrixOperations/cross/scrossa.c
new file mode 100644
index 0000000..414b85a
--- /dev/null
+++ b/src/c/matrixOperations/cross/scrossa.c
@@ -0,0 +1,33 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+
+void scrossa (float* inp1, int size1_r, int size1_c, float* inp2, int size2_r, int size2_c ,float* 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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/cross/u16crossa.c b/src/c/matrixOperations/cross/u16crossa.c
new file mode 100644
index 0000000..e0a32bc
--- /dev/null
+++ b/src/c/matrixOperations/cross/u16crossa.c
@@ -0,0 +1,34 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+#include "uint16.h"
+
+void u16crossa (uint16* inp1, int size1_r, int size1_c, uint16* inp2, int size2_r, int size2_c ,uint16* 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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/cross/u8crossa.c b/src/c/matrixOperations/cross/u8crossa.c
new file mode 100644
index 0000000..7a8b857
--- /dev/null
+++ b/src/c/matrixOperations/cross/u8crossa.c
@@ -0,0 +1,34 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+#include "uint8.h"
+
+void u8crossa (uint8* inp1, int size1_r, int size1_c, uint8* inp2, int size2_r, int size2_c ,uint8* 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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/cross/zcrossa.c b/src/c/matrixOperations/cross/zcrossa.c
new file mode 100644
index 0000000..85a1960
--- /dev/null
+++ b/src/c/matrixOperations/cross/zcrossa.c
@@ -0,0 +1,37 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cross.h>
+#include "doubleComplex.h"
+#include "stdlib.h"
+#include "string.h"
+#include "cat.h"
+
+void zcrossa (doubleComplex* inp1, int size1_r, int size1_c, doubleComplex* inp2, int size2_r, int size2_c ,doubleComplex* 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)];
+
+            }
+        }
+
+    }
+}
+
diff --git a/src/c/matrixOperations/includes/cross.h b/src/c/matrixOperations/includes/cross.h
new file mode 100644
index 0000000..0abd001
--- /dev/null
+++ b/src/c/matrixOperations/includes/cross.h
@@ -0,0 +1,37 @@
+ /* Copyright (C) 2017 - 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: Brijesh Gupta C R
+ Organization: FOSSEE, IIT Bombay
+ Email: toolbox@scilab.in
+ */
+
+#ifndef __CROSS_H__
+#define __CROSS_H__
+#include "types.h"
+#include "doubleComplex.h"
+#include "floatComplex.h"
+
+#ifdef  __cplusplus
+extern "C" {
+#endif
+
+void dcrossa (double* inp1, int size1_r, int size1_c, double* inp2, int size2_r, int size2_c ,double* oup);
+void zcrossa (doubleComplex* inp1, int size1_r, int size1_c, doubleComplex* inp2, int size2_r, int size2_c ,doubleComplex* oup);
+void ccrossa (floatComplex* inp1, int size1_r, int size1_c, floatComplex* inp2, int size2_r, int size2_c ,floatComplex* oup);
+void scrossa (float* inp1, int size1_r, int size1_c, float* inp2, int size2_r, int size2_c ,float* oup);
+void u8crossa (uint8* inp1, int size1_r, int size1_c, uint8* inp2, int size2_r, int size2_c ,uint8* oup);
+void u16crossa (uint16* inp1, int size1_r, int size1_c, uint16* inp2, int size2_r, int size2_c ,uint16* oup);
+void i8crossa (int8* inp1, int size1_r, int size1_c, int8* inp2, int size2_r, int size2_c ,int8* oup);
+void i16crossa (int16* inp1, int size1_r, int size1_c, int16* inp2, int size2_r, int size2_c ,int16* oup);
+
+
+#ifdef  __cplusplus 
+} /* extern "C" */
+#endif
+
+#endif /*__CROSS_H__*/
diff --git a/src/c/matrixOperations/interfaces/int_cross.h b/src/c/matrixOperations/interfaces/int_cross.h
new file mode 100644
index 0000000..8f4e782
--- /dev/null
+++ b/src/c/matrixOperations/interfaces/int_cross.h
@@ -0,0 +1,35 @@
+ /* Copyright (C) 2017 - 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: Brijesh Gupta C R
+ Organization: FOSSEE, IIT Bombay
+ Email: toolbox@scilab.in
+ */
+#ifndef __INT_CROSS_H__
+#define __INT_CROSS_H__
+
+#ifdef  __cplusplus
+extern "C" {
+#endif
+
+
+
+#define d2d2crossd2(in1,size1,in2,size2,out) dcrossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define z2z2crossz2(in1,size1,in2,size2,out) zcrossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define c2c2crossc2(in1,size1,in2,size2,out) ccrossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define s2s2crosss2(in1,size1,in2,size2,out) scrossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define u82u82crossu82(in1,size1,in2,size2,out) u8crossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define i82i82crossi82(in1,size1,in2,size2,out) i8crossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define u162u162crossu162(in1,size1,in2,size2,out) u16crossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+#define i162i162crossi162(in1,size1,in2,size2,out) i16crossa(in1,size1[0],size1[1],in2,size2[0],size2[1],out)
+
+
+#ifdef  __cplusplus
+} /* extern "C" */
+#endif
+
+#endif /*__INT_CROSS_H__*/