added test to validate the state transformation but has an unknown but. despite that, the state transformation works as expected

Author: ZocoLini

modules/cuqdyn-c/src/states_transformer.c

@@ -4,7 +4,6 @@
 #include <nvector/nvector_serial.h>
 #include <stdio.h>
 #include <sunmatrix/sunmatrix_dense.h>
-#include <unistd.h>
 
 #include "config.h"
 #include "matlab.h"
@@ -40,7 +39,7 @@ SUNMatrix transform_states(SUNMatrix states)
 
         N_VDestroy(input);
     }
-
+    
     free(output);
     SUNMatDestroy(states);
     return transformed_result;

tests/CMakeLists.txt

@@ -10,6 +10,7 @@ foreach(test_name
         test_ess_obj_function
         test_cuqdyn_algo
         test_context_and_config
+        test_transform_states
 )
     add_executable(${test_name} ${test_name}.c)
     target_link_libraries(${test_name} cuqdyn-c)

tests/test_transform_states.c

@@ -0,0 +1,45 @@
+#include <assert.h>
+#include <sunmatrix/sunmatrix_dense.h>
+#include "config.h"
+#include "cuqdyn.h"
+
+#define CUQDYN_CONF "data/nfkb_cuqdyn_config.xml"
+
+int main()
+{
+    CuqDynContext context = init_cuqdyn_context_from_file(CUQDYN_CONF);
+    CuqdynConf *conf = get_cuqdyn_conf(context);
+
+    sunrealtype states_values[] = {
+        0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
+        1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2
+    };
+    sunrealtype expected_transformation_values[] = {
+        0, 1, 2, 1, 3, 1, 1,
+        1, 2, 4, 2, 6, 2, 2
+    };
+
+    assert(6 == conf->states_transformer.count);
+
+    SUNMatrix states = NewDenseMatrix(2, 16);
+
+    for (int i = 0; i < 2; ++i) {
+        for (int j = 0; j < 16; ++j) {
+            SM_ELEMENT_D(states, i, j) = states_values[16 * i + j];
+        }
+    }
+
+    SUNMatrix transformation = transform_states(states);
+
+    for (int i = 0; i < 2; ++i) {
+        for (int j = 0; j < 7; ++j) {
+            // TODO: This transformation works as expected but SM_ELEMENT_D thrown seg fault
+            assert(expected_transformation_values[7 * i + j] == SM_ELEMENT_D(transformation, i, j));
+        }
+    }
+
+    destroy_cuqdyn_context(context);
+    SUNMatDestroy(transformation);
+
+    return 0;
+}