summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/statisticsFunctions/center/scentercola.c
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/c/statisticsFunctions/center/scentercola.c')
-rw-r--r--2.3-1/src/c/statisticsFunctions/center/scentercola.c37
1 files changed, 37 insertions, 0 deletions
diff --git a/2.3-1/src/c/statisticsFunctions/center/scentercola.c b/2.3-1/src/c/statisticsFunctions/center/scentercola.c
new file mode 100644
index 00000000..d5ef5eeb
--- /dev/null
+++ b/2.3-1/src/c/statisticsFunctions/center/scentercola.c
@@ -0,0 +1,37 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include <ones.h>
+#include "subtraction.h"
+#include "matrixMultiplication.h"
+#include "center.h"
+
+void scentercola (float* inp, int row, int col, float* out)
+{
+
+ float sum[row], xbar[row];
+ for(int i = 0; i < row; i++)
+ sum[i] = 0;
+
+
+ for(int i = 0; i < row; i++)
+ {
+ for(int j = 0; j < col; j++)
+ {
+ sum[i] += inp[i + j*row];
+ }
+ }
+
+ for(int i = 0; i < row; i++)
+ xbar[i] = sum[i]/col;
+
+ float one[col]; // Creating a matrix of ones
+ sonesa(one,1,col);
+
+
+ float prod[row*col];
+ smulma(xbar, row, 1, one, 1, col, prod);
+
+ for(int i = 0; i< row*col; i++)
+ out[i] = inp[i] - prod[i];
+}