trying to make a test to validate the obj_function output

Author: ZocoLini

tests/CMakeLists.txt

@@ -8,6 +8,7 @@ foreach(test_name
         test_data_reader
         test_ess_solver
         test_parse_cuqdyn_conf
+        test_ess_obj_function
         test_cuqdyn_algo
 )
     add_executable(${test_name} ${test_name}.c)

tests/data/obj_function_cuqdyn_config.xml

@@ -0,0 +1,11 @@
+<?xml version="1.0" encoding="UTF-8" ?>
+
+<cuqdyn-config>
+    <tolerances>
+        <rtol>1e-8</rtol>
+        <atol>1e-12, 1e-12, 1e-12, 1e-12 1e-12</atol>
+    </tolerances>
+    <ode_expr y_count="5" p_count="5">
+        alpha-pinene
+    </ode_expr>
+</cuqdyn-config>

tests/data/obj_function_ess_config_nl2sol.dn2fb.xml

@@ -0,0 +1,66 @@
+<xml>
+    <run>
+        <typebench>customized</typebench>
+        <id>0</id>
+        <log_scale>5</log_scale>
+        <output>1</output>
+        <verbose>1</verbose>
+        <local_search>1</local_search>
+        <stopping_criteria>
+            <maxevaluation>3e3</maxevaluation>
+            <maxtime>1e10</maxtime>
+            <vtr>0</vtr>
+        </stopping_criteria>
+    </run>
+    <method name="ScatterSearch">
+        <user_options>
+            <weight>default</weight>
+            <tolc>default</tolc>
+            <prob_bound>default</prob_bound>
+            <nstuck_solution>default</nstuck_solution>
+        </user_options>
+        <global_options>
+            <dim_ref>default</dim_ref>
+            <ndiverse>default</ndiverse>
+            <combination>default</combination>
+            <n_stuck>default</n_stuck>
+        </global_options>
+        <local_options>
+            <solver>nl2sol.dn2fb</solver>
+            <tol>0</tol>
+            <evalmax>3e3</evalmax>
+            <iterprint>0</iterprint>
+            <n1>1</n1>
+            <n2>10</n2>
+            <balance>0.25</balance>
+            <bestx>default</bestx>
+        </local_options>
+    </method>
+    <parallelization name="cooperative">
+        <!-- saCeSS options -->
+        <reception_threshold>1</reception_threshold>
+        <evals_threshold>1000</evals_threshold>
+        <mult_num_sendSol>10</mult_num_sendSol>
+        <minimum_num_sendSol>20</minimum_num_sendSol>
+        <!-- CeSS option: migration time -->
+        <migration_max_time>10</migration_max_time>
+    </parallelization>
+
+    <problem>
+        <dim> 5 </dim>
+        <neq> 0 </neq>
+        <ineq> 0 </ineq>
+        <int_var> 0 </int_var>
+        <bin_var> 0 </bin_var>
+        <lb>
+            0, 0, 0, 0, 0
+        </lb>
+        <ub>
+            1, 1, 1, 1, 1
+        </ub>
+        <point>
+            1, 1, 1, 1, 1
+        </point>
+    </problem>
+
+</xml>

tests/test_ess_obj_function.c

@@ -0,0 +1,31 @@
+#include "data_reader.h"
+#include "functions.h"
+#include "matlab.h"
+#include "method_module/structure_paralleltestbed.h"
+
+#define CUQDYN_CONF "data/obj_function_cuqdyn_config.xml"
+#define SACESS_CONF "data/obj_function_ess_config_nl2sol.dn2fb.xml"
+#define DATA "data/alpha_pinene_paper_data.txt"
+
+int main()
+{
+    init_cuqdyn_conf_from_file(CUQDYN_CONF);
+
+    N_Vector texp = NULL;
+    SUNMatrix yexp = NULL;
+
+    read_data_file(DATA, &texp, &yexp);
+    N_Vector initial_values = copy_matrix_row(yexp, 0, 0, SM_COLUMNS_D(yexp));
+
+    experiment_total *exp_total = malloc(sizeof(experiment_total));
+    create_expetiment_struct(SACESS_CONF, &exp_total[0], 1, 0, "output", 1, texp, yexp, initial_values);
+
+    double x[5] = {1.0, 1.0, 1.0, 1.0, 1.0};
+
+    output_function *out = ode_model_obj_func(x, exp_total);
+
+    printf("VAlue %lf", out->value);
+
+    destroyexp(exp_total);
+    return 0;
+}