9 files changed, 168 insertions, 0 deletions
diff --git a/src/c/statisticsFunctions/moment/dmomentcola.c b/src/c/statisticsFunctions/moment/dmomentcola.c
new file mode 100644
index 00000000..538328e3
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/dmomentcola.c
@@ -0,0 +1,18 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+
+void dmomentcola (double* inp, int row, int col, double ord, double* out)
+{
+    double vect[col];
+    for(int i = 0; i < row ; i++)
+    {
+	for(int j = 0; j < col ; j++)
+	{
+	    vect[j] = inp[i + j*row];
+	}
+	out[i] = dmoments(vect, col, ord);
+    }
+
+}
diff --git a/src/c/statisticsFunctions/moment/dmomentrowa.c b/src/c/statisticsFunctions/moment/dmomentrowa.c
new file mode 100644
index 00000000..811aac39
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/dmomentrowa.c
@@ -0,0 +1,18 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+
+void dmomentrowa (double* inp, int row, int col, double ord, double* out)
+{
+    double vect[row];
+    for(int i = 0; i < col ; i++)
+    {
+	for(int j = 0; j < row ; j++)
+	{
+	    vect[j] = inp[j + i*row];
+	}
+	out[i] = dmoments(vect, row, ord);
+    }
+
+}
diff --git a/src/c/statisticsFunctions/moment/dmoments.c b/src/c/statisticsFunctions/moment/dmoments.c
new file mode 100644
index 00000000..779ffa0f
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/dmoments.c
@@ -0,0 +1,17 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+
+
+double dmoments (double* inp, int size, double ord)
+{
+    double sum = 0;
+
+    for(int i = 0; i < size; i++)		// Elements are raised to the order and then their mean is calculated to give moment
+    {
+        sum = sum + pow(inp[i], ord);
+    }
+
+    return sum/size;
+}
diff --git a/src/c/statisticsFunctions/moment/smomentcola.c b/src/c/statisticsFunctions/moment/smomentcola.c
new file mode 100644
index 00000000..4255782c
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/smomentcola.c
@@ -0,0 +1,18 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+
+void smomentcola (float* inp, int row, int col, double ord, float* out)
+{
+    float vect[col];
+    for(int i = 0; i < row ; i++)
+    {
+	for(int j = 0; j < col ; j++)
+	{
+	    vect[j] = inp[i + j*row];
+	}
+	out[i] = smoments(vect, col, ord);
+    }
+
+}
diff --git a/src/c/statisticsFunctions/moment/smomentrowa.c b/src/c/statisticsFunctions/moment/smomentrowa.c
new file mode 100644
index 00000000..574ebe76
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/smomentrowa.c
@@ -0,0 +1,18 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+
+void smomentrowa (float* inp, int row, int col, double ord, float* out)
+{
+    float vect[row];
+    for(int i = 0; i < col ; i++)
+    {
+	for(int j = 0; j < row ; j++)
+	{
+	    vect[j] = inp[j + i*row];
+	}
+	out[i] = smoments(vect, row, ord);
+    }
+
+}
diff --git a/src/c/statisticsFunctions/moment/smoments.c b/src/c/statisticsFunctions/moment/smoments.c
new file mode 100644
index 00000000..c900d815
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/smoments.c
@@ -0,0 +1,17 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+
+
+float smoments (float* inp, int size, double ord)
+{
+    float sum = 0;
+
+    for(int i = 0; i < size; i++)			// Elements are raised to the order and then their mean is calculated to give moment
+    {
+        sum = sum + pow(inp[i], ord);
+    }
+
+    return sum/size;
+}
diff --git a/src/c/statisticsFunctions/moment/zmomentcola.c b/src/c/statisticsFunctions/moment/zmomentcola.c
new file mode 100644
index 00000000..2b3cebb1
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/zmomentcola.c
@@ -0,0 +1,20 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+#include "doubleComplex.h"
+#include "floatComplex.h"
+
+void zmomentcola (doubleComplex* inp, int row, int col, double ord, doubleComplex* out)
+{
+    double vect[col];
+    for(int i = 0; i < row ; i++)
+    {
+	for(int j = 0; j < col ; j++)
+	{
+	    vect[j] = inp[i + j*row];
+	}
+	out[i] = zmoments(vect, col, ord);
+    }
+
+}
diff --git a/src/c/statisticsFunctions/moment/zmomentrowa.c b/src/c/statisticsFunctions/moment/zmomentrowa.c
new file mode 100644
index 00000000..e1e1b47d
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/zmomentrowa.c
@@ -0,0 +1,20 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+#include "doubleComplex.h"
+#include "floatComplex.h"
+
+void zmomentrowa (doubleComplex* inp, int row, int col, double ord, doubleComplex* out)
+{
+    doubleComplex vect[row];
+    for(int i = 0; i < col ; i++)
+    {
+	for(int j = 0; j < row ; j++)
+	{
+	    vect[j] = inp[j + i*row];
+	}
+	out[i] = zmoments(vect, row, ord);
+    }
+
+}
diff --git a/src/c/statisticsFunctions/moment/zmoments.c b/src/c/statisticsFunctions/moment/zmoments.c
new file mode 100644
index 00000000..b70df707
--- /dev/null
+++ b/src/c/statisticsFunctions/moment/zmoments.c
@@ -0,0 +1,22 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "moment.h"
+#include "doubleComplex.h"
+#include "floatComplex.h"
+#include "pow.h"
+#include "addition.h"
+#include "division.h"
+
+
+doubleComplex zmoments (doubleComplex* inp, int size, double ord)
+{
+    doubleComplex sum = DoubleComplex(0,0);
+
+    for(int i = 0; i < size; i++)			// Elements are raised to the order and then their mean is calculated to give moment
+    {
+        sum = zadds(sum,zpows(inp[i], ord));
+    }
+
+    return zrdivs(sum,size);
+}