/** Loading a solution file.
* Returns: an array of solution
* Impotant: This array of solution need to be desallocated
*/
_2PG_CARTESIAN_EXPORT
solution_t * loading_file_solutions(int *num_solutions_r,
int *numobj_r, const char *path_file_name){
FILE * solution_file;
char *line;
char *line_splited;
int sol;
solution_t * solutions_aux;
line = Malloc(char, MAX_LINE_SOLUTION);
//Getting initial information
initialize_for_reading_file(num_solutions_r, numobj_r, path_file_name);
//Alocating Solution
solutions_aux = allocate_solution(num_solutions_r, numobj_r);
//Reading file and set values of objective
sol = -1;
solution_file = open_file(path_file_name, fREAD);
//Removing first line that is collumn
fgets(line,MAX_LINE_SOLUTION,solution_file);
while ( fgets(line,MAX_LINE_SOLUTION,solution_file) != NULL){
sol = sol + 1;
//Removing index collumn
line_splited = strtok (line,"\t");
//Setting number of objectives
solutions_aux[sol].num_obj = *numobj_r;
for (int ob = 0; ob < *numobj_r; ob++){
line_splited = strtok (NULL,"\t");
trim(line_splited);
solutions_aux[sol].obj_values[ob] = str2float(line_splited);//str2double(line_splited);
}
}
fclose(solution_file);
free(line);
return solutions_aux;
}
int main(int argc, char *argv[]){
input_parameters_t in_param;
display_msg("Reading the configure file \n");
load_parameters_from_file(&in_param,argv[1]);
/*Represents extension of files that will be looked in directory*/
char *ext = NULL;
/*Represents how many files there is in directory*/
int num_files;
/*Represents file names that are in directory and have target extension*/
owner_file_t *file_names = NULL;
protein_t *population_p= NULL; // main population of protein
solution_t *solutions_p= NULL; // main solution
dominance_t *dominance= NULL; // main of dominance
ea_nsga2_t *nsga2_solutions_p= NULL; // main of front
/**************** START GETTING FILE NAMES *************************/
ext = Malloc(char, 4);
strcpy(ext, "pdb");
num_files = how_many_files_directory_by_extension(in_param.path_local_execute, ext);
file_names = allocate_file_t(&num_files, &in_param.number_fitness);
insert_files_directory_by_extension(file_names, in_param.path_local_execute, ext);
free(ext);
/**************** FINISHED GETTING FILE NAMES *************************/
/**************** START GETTING THE OBJECTIVES *************************/
init_gromacs_execution();
/*
// RUN PDB2GMX IN ALL FILES FOR PATTERN OF ATOM NAMES
display_msg("Run of pdb2gmx for pattern of atom names\n");
for (int ind = 0; ind < num_files; ind ++){
call_pdb2gmx_for_pattern_atom_names(file_names[ind].file_name, in_param.path_local_execute, in_param.path_gromacs_programs, in_param.force_field);
clean_gromacs_simulation(in_param.path_local_execute);
}
// FINISHED RUN PDB2GMX
*/
population_p = allocateProtein(&num_files);
solutions_p = allocate_solution(&num_files, &in_param.number_fitness);
for (int ind = 0; ind < num_files; ind ++){
printf("Computing Objectives of %s\n", file_names[ind].file_name);
char *path_PDB_file_name = path_join_file(in_param.path_local_execute, file_names[ind].file_name);
int num_atom = get_num_atom(path_PDB_file_name);
population_p[ind].p_atoms = allocate_pdbatom(&num_atom);
load_pdb_file_without_num_atom(population_p[ind].p_atoms, NULL, path_PDB_file_name);
int num_res = get_number_residues_from_atom(population_p[ind].p_atoms, &num_atom);
population_p[ind].p_topol = allocateTop_Global(&num_res, &num_atom);
renumerate_residue_number(population_p[ind].p_atoms, &num_atom);
build_topology_individual(&population_p[ind]);
solutions_p[ind].representation = &population_p[ind];
get_gromacs_objectives_of_solution(&solutions_p[ind], &in_param, &ind);
desAllocateTop_Global(population_p[ind].p_topol);
desAllocate_pdbatom(population_p[ind].p_atoms);
free(path_PDB_file_name);
}
finish_gromacs_execution();
//Deleting temporary files
system("rm PROT_IND_*.pdb");
/**************** FINISHED GETTING THE OBJECTIVES *************************/
/**************** START GETTING FRONT *************************/
in_param.size_population = num_files;
nsga2_solutions_p = allocate_nsga2_without_allocation_of_representation(&in_param);
//Setting identification
for (int i = 0; i < num_files; i++){
nsga2_solutions_p[i].sol->ID = i+1;
}
//Setting dominance
dominance = allocate_dominance(&num_files);
set_dominance(dominance, solutions_p, &num_files);
//Coping values of dominance
for (int ind = 0; ind < num_files; ind++){
// Indicates number of solutions that are dominated by me
file_names[ind].number_solutions_are_dominated = dominance[ind].max_dominated;
}
//Coping values of objective
for (int ind = 0; ind < num_files; ind++){
for (int ob = 0; ob < in_param.number_fitness; ob++)
nsga2_solutions_p[ind].sol->obj_values[ob] = solutions_p[ind].obj_values[ob];
}
//Setting front based on dominance concept
compute_fronts(nsga2_solutions_p, dominance, &num_files);
//Coping values from nsga2_solutions_p to owner_file_t
for (int ind = 0; ind < num_files; ind++){
//Coping objectives
for (int ob = 0; ob < in_param.number_fitness; ob++){
file_names[ind].obj_values[ob] = nsga2_solutions_p[ind].sol->obj_values[ob];
}
file_names[ind].front = nsga2_solutions_p[ind].front;
}
desallocate_dominance(dominance, &num_files);
desallocate_solution_nsga2(nsga2_solutions_p, &num_files);
/**************** FINISHED GETTING FRONT *************************/
/**************** START GETTING FINAL RESULTS *************************/
//Sorting solutions
sorting_solutions_by_front_dominance(file_names, &num_files, &in_param.number_fitness);
//Saving file
save_analysis_files(file_names, &num_files, &in_param.number_fitness, in_param.fitness_energies);
/**************** FINISHED FINAL RESULTS *************************/
desalocate_file_t(file_names, &num_files);
desallocate_solution(solutions_p, &num_files);
desallocateProtein(population_p, &num_files);
deAllocateload_parameters(&in_param);
display_msg("Done !!! \n");
return 0;
}
