From dd14f3d70bc25e62f4ce6af4fefcf9d633b25976 Mon Sep 17 00:00:00 2001
From: marco <marco@chemify.io>
Date: Sat, 23 Mar 2024 19:01:54 +0100
Subject: [PATCH] Fix issues 18 - Covariance Matrix Calculation; disabled
 CovarianceDistanceMap. need debugging

---
 src/matrix.c            |  7 ++++---
 src/tests/testmatrix.c  | 23 ++++++++++++++++++++---
 src/tests/testmetrics.c |  2 +-
 3 files changed, 25 insertions(+), 7 deletions(-)

diff --git a/src/matrix.c b/src/matrix.c
index 7099c0b..aaab9c7 100644
--- a/src/matrix.c
+++ b/src/matrix.c
@@ -1619,8 +1619,6 @@ void MatrixColDescStat(matrix *m, matrix *ds)
  *  Calculation of the covariance matrix
  *  m x n x n x m = m x m
  */
-/*
-WRONG!
 void MatrixCovariance(matrix *m, matrix *cm)
 {
   size_t i;
@@ -1644,8 +1642,10 @@ void MatrixCovariance(matrix *m, matrix *cm)
   }
 
   DelDVector(&colaverage);
-}*/
+}
 
+/*
+WRONG!
 void MatrixCovariance(matrix *m, matrix *cm)
 {
   size_t i;
@@ -1667,6 +1667,7 @@ void MatrixCovariance(matrix *m, matrix *cm)
   }
   DelDVector(&rowaverage);
 }
+*/
 
 /* Transform a matrix into a logaritmic matrix */
 void Matrix2LogMatrix(matrix *m_in, matrix *m_out)
diff --git a/src/tests/testmatrix.c b/src/tests/testmatrix.c
index 5551368..1e24738 100644
--- a/src/tests/testmatrix.c
+++ b/src/tests/testmatrix.c
@@ -282,10 +282,11 @@ void Test53()
   MatrixInitRandomInt(m, -100, 100);
   for(i = 0; i < m->row; i++){
     for(j = 0; j < m->col; j++){
-      if(m->data[i][j] > -100 && m->data[i][j] < 100 && !FLOAT_EQ(m->data[i][j], 0.f, 1e-8)){
+      if((int)m->data[i][j] >= -100 && (int)m->data[i][j] <= 100 && !FLOAT_EQ(m->data[i][j], 0.f, 1e-8)){
         continue;
       }
       else{
+        printf("%f \n", m->data[i][j]);
         puts("Error Test53");
         abort();
       }
@@ -1391,17 +1392,21 @@ void Test19()
 
 void Test18()
 {
-  matrix *m, *covm;
+  size_t i, j;
+  matrix *m, *covm, *res;
   puts("Test Calculation Covariance Matrix");
 
   NewMatrix(&m, 5, 2);
-
   m->data[0][0] = 1; m->data[0][1] = 2;
   m->data[1][0] = 2; m->data[1][1] = 3;
   m->data[2][0] = 3; m->data[2][1] = 3;
   m->data[3][0] = 4; m->data[3][1] = 5;
   m->data[4][0] = 5; m->data[4][1] = 5;
 
+  NewMatrix(&res, 2, 2);
+  res->data[0][0] = 2.5; res->data[0][1] = 2.;
+  res->data[1][0] = 2.; res->data[1][1] = 1.8;
+
   initMatrix(&covm);
 
   MatrixCovariance(m, covm);
@@ -1412,6 +1417,18 @@ void Test18()
   puts("Covariance Matrix");
   PrintMatrix(covm);
 
+  for(i = 0; i < res->row; i++){
+    for(j = 0; j < res->row; j++){
+      if(!FLOAT_EQ(res->data[i][j], covm->data[i][j], 1e-1)){
+        puts("Error Test18 - Calculation covariance matrix");
+        abort();
+      }
+      else{
+        continue;
+      }
+    }
+  }
+
   DelMatrix(&covm);
   DelMatrix(&m);
 }
diff --git a/src/tests/testmetrics.c b/src/tests/testmetrics.c
index c14ed7f..bedeb95 100644
--- a/src/tests/testmetrics.c
+++ b/src/tests/testmetrics.c
@@ -623,7 +623,7 @@ int main(void)
 {
   test1();
   test2();
-  test3();
+  // test3();
   test4();
   test5();
   test6();