error related to the initial point of the eSS found

Author: ZocoLini

example-files/lotka_volterra_ess_serial_config.xml

@@ -58,9 +58,9 @@
         <ub>
             1.000, 1.000, 1.000, 1.000
         </ub>
-        <initial_point>
+        <point>
             0.3, 0.3, 0.3, 0.3
-        </initial_point>
+        </point>
     </problem>
 
 </xml>

modules/cuqdyn-c/src/cuqdyn.c

@@ -77,6 +77,18 @@ CuqdynResult *cuqdyn_algo(const char *data_file, const char *sacess_conf_file,
     {
         resid_loo = NewDenseMatrix(m, n);
         predicted_params_matrix = NewDenseMatrix(m, config->ode_expr.p_count);
+
+        N_Vector init_vals = New_Serial(NV_LENGTH_S(initial_values));
+        memcpy(NV_DATA_S(init_vals), NV_DATA_S(initial_values), NV_LENGTH_S(initial_values) * sizeof(sunrealtype));
+
+        N_Vector texp = copy_vector_remove_indices(times, create_array((long[]) {}, 0));
+
+        SUNMatrix yexp = copy_matrix_remove_rows(observed_data, create_array((long[]) {}, 0));
+
+        N_Vector predicted_params =
+            execute_ess_solver(sacess_conf_file, output_file, texp, yexp, init_vals, rank, nproc);
+
+
     }
 
 #ifdef MPI

modules/sacess/src/input_module/input_module.c

@@ -1025,6 +1025,7 @@ int extract_element_problem(xmlDocPtr doc, xmlNodePtr *root, experiment_total *e
     }
 
 // INIT POINT
+    printf("asdfasdf");
     if (extract_element_uniq(doc, cur, name_sol) != NULL){
 	test->bench.number_init_sol = count_element_multi(doc, cur, name_sol);
 
@@ -1036,12 +1037,14 @@ int extract_element_problem(xmlDocPtr doc, xmlNodePtr *root, experiment_total *e
 	     vector_char[i] = (char *) malloc(100000* sizeof(char));
 	     test->bench.X0[i] = (double *) malloc(test->bench.dim*sizeof(double));
 	}
+
 	extract_element_multi(doc,cur,name_sol,vector_char);
         for (i=0;i<test->bench.number_init_sol;i++) {
     		counter=0;
     		pt = strtok (vector_char[i],",");
     		while (pt != NULL) {
        	     		test->bench.X0[i][counter] = atof(pt);
+    		    printf("%lf", test->bench.X0[i][counter]);
 	        	pt = strtok (NULL, ",");
 	        	counter++;
 	    	}