modifying x0 values inside the experiment_total struct works for the mpi version

Author: ZocoLini

example-files/lotka_volterra_ess_mpi_config.xml

@@ -7,7 +7,7 @@
         <verbose>0</verbose>
         <local_search>1</local_search>
         <stopping_criteria>
-            <maxevaluation>2e3</maxevaluation>
+            <maxevaluation>3e3</maxevaluation>
             <maxtime>1e10</maxtime>
             <vtr>0</vtr>
         </stopping_criteria>

modules/cuqdyn-c/src/cuqdyn.c

@@ -97,16 +97,11 @@ CuqdynResult *cuqdyn_algo(const char *data_file, const char *sacess_conf_file, c
             execute_ess_solver(sacess_conf_file, output_file, texp, yexp, tmp_initial_condition, NULL);
 
         memcpy(NV_DATA_S(initial_params), NV_DATA_S(predicted_params), NV_LENGTH_S(predicted_params) * sizeof(sunrealtype));
-
         N_VDestroy(predicted_params);
-
-#ifdef MPI
-        MPI_Bcast(NV_DATA_S(initial_params), NV_LENGTH_S(initial_params), MPI_DOUBLE, 0, MPI_COMM_WORLD);
-#endif
     }
 
 #ifdef MPI
-    MPI_Barrier(MPI_COMM_WORLD);
+    MPI_Bcast(NV_DATA_S(initial_params), NV_LENGTH_S(initial_params), MPI_DOUBLE, 0, MPI_COMM_WORLD);
 
     long iterations;
     long start_index;

modules/cuqdyn-c/src/ess_solver.c

@@ -82,9 +82,7 @@ N_Vector execute_ess_solver(const char *file, const char *path, N_Vector texp, S
     {
         for (int i = 0; i < NV_LENGTH_S(initial_params); i++)
         {
-            printf("QASDFASDFASDFASDF\n");
             exptotal[0].test.bench.X0[0][i] = NV_Ith_S(initial_params, i);
-            printf("%lf", exptotal[0].test.bench.X0[0][i]);
         }
     }
 

modules/sacess/src/input_module/input_module.c

@@ -1043,7 +1043,6 @@ int extract_element_problem(xmlDocPtr doc, xmlNodePtr *root, experiment_total *e
     		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++;
 	    	}

mpi_toolchain.cmake

@@ -1,8 +1,8 @@
 set(CMAKE_C_COMPILER mpicc)
 set(CMAKE_CXX_COMPILER mpic++)
 set(CMAKE_Fortran_COMPILER mpifort)
-set(CMAKE_C_FLAGS "-O3 -cpp -DGNU -fPIC -no-pie -DOPENMP -DMPI -fopenmp -w")
-set(CMAKE_Fortran_FLAGS "-O3 -cpp -DGNU -w -fallow-argument-mismatch -std=gnu -DOPENMP -DMPI -DGNU -fopenmp")
+set(CMAKE_C_FLAGS "-O3 -cpp -DGNU -fPIC -no-pie -DMPI -w")
+set(CMAKE_Fortran_FLAGS "-O3 -cpp -DGNU -w -fallow-argument-mismatch -std=gnu -DGNU")
 
 set(MISQP_LIBRARY_DIR "${CMAKE_CURRENT_SOURCE_DIR}/deps/misqp/gnu")
 set(LIBRARIES mpi)