#include #include #include #include #include "read_ini.h" #include "../malleability/ProcessDist.h" #include "ini.h" void malloc_config_arrays(configuration *user_config, int resizes); void def_struct_config_file(configuration *config_file, MPI_Datatype *config_type); void def_struct_config_file_array(configuration *config_file, MPI_Datatype *config_type); /* * Funcion utilizada para leer el fichero de configuracion * y guardarlo en una estructura para utilizarlo en el futuro. * * Primero lee la seccion "general" y a continuacion cada una * de las secciones "resize%d". */ static int handler(void* user, const char* section, const char* name, const char* value) { configuration* pconfig = (configuration*)user; char *resize_name = malloc(10 * sizeof(char)); int act_resize = pconfig->actual_resize; snprintf(resize_name, 10, "resize%d", act_resize); #define MATCH(s, n) strcmp(section, s) == 0 && strcmp(name, n) == 0 if (MATCH("general", "resizes")) { pconfig->resizes = atoi(value) + 1; malloc_config_arrays(pconfig, pconfig->resizes); } else if (MATCH("general", "matrix_tam")) { pconfig->matrix_tam = atoi(value); } else if (MATCH("general", "comm_tam")) { pconfig->comm_tam = atoi(value); } else if (MATCH("general", "SDR")) { pconfig->sdr = atoi(value); } else if (MATCH("general", "ADR")) { pconfig->adr = atoi(value); } else if (MATCH("general", "AIB")) { //TODO Refactor cambiar nombre pconfig->aib = atoi(value); } else if (MATCH("general", "CST")) { pconfig->cst = atoi(value); } else if (MATCH("general", "CSS")) { pconfig->css = atoi(value); } else if (MATCH("general", "time")) { pconfig->general_time = atof(value); // Resize } else if (MATCH(resize_name, "iters")) { pconfig->iters[act_resize] = atoi(value); } else if (MATCH(resize_name, "procs")) { pconfig->procs[act_resize] = atoi(value); } else if (MATCH(resize_name, "factor")) { pconfig->factors[act_resize] = atof(value); } else if (MATCH(resize_name, "physical_dist")) { char *aux = strdup(value); if (strcmp(aux, "node") == 0) { pconfig->phy_dist[act_resize] = COMM_PHY_NODES; } else { pconfig->phy_dist[act_resize] = COMM_PHY_CPU; } free(aux); pconfig->actual_resize = pconfig->actual_resize+1; // Ultimo elemento del grupo } else { return 0; /* unknown section or name, error */ } free(resize_name); return 1; } /* * Crea y devuelve una estructura de configuracion a traves * de un nombre de fichero dado. * * La memoria de la estructura se reserva en la funcion y es conveniente * liberarla con la funcion "free_config()" */ configuration *read_ini_file(char *file_name) { configuration *config = NULL; config = malloc(sizeof(configuration) * 1); if(config == NULL) { printf("Error when reserving configuration structure\n"); return NULL; } config->actual_resize=0; if(ini_parse(file_name, handler, config) < 0) { // Obtener configuracion printf("Can't load '%s'\n", file_name); return NULL; } return config; } /* * Reserva de memoria para los vectores de la estructura de configuracion * * Si se llama desde fuera de este fichero, la memoria de la estructura * tiene que reservarse con la siguiente linea: * "configuration *config = malloc(sizeof(configuration));" * * Sin embargo se puede obtener a traves de las funciones * - read_ini_file * - recv_config_file */ void malloc_config_arrays(configuration *user_config, int resizes) { if(user_config != NULL) { user_config->iters = malloc(sizeof(int) * resizes); user_config->procs = malloc(sizeof(int) * resizes); user_config->factors = malloc(sizeof(float) * resizes); user_config->phy_dist = malloc(sizeof(int) * resizes); } } /* * Libera toda la memoria de una estructura de configuracion */ void free_config(configuration *user_config) { if(user_config != NULL) { free(user_config->iters); free(user_config->procs); free(user_config->factors); free(user_config->phy_dist); free(user_config); } } /* * Imprime por salida estandar toda la informacion que contiene * la configuracion pasada como argumento */ void print_config(configuration *user_config, int grp) { if(user_config != NULL) { int i; printf("Config loaded: resizes=%d, matrix=%d, comm_tam=%d, sdr=%d, adr=%d, aib=%d, css=%d, cst=%d, time=%f || grp=%d\n", user_config->resizes, user_config->matrix_tam, user_config->comm_tam, user_config->sdr, user_config->adr, user_config->aib, user_config->css, user_config->cst, user_config->general_time, grp); for(i=0; iresizes; i++) { printf("Resize %d: Iters=%d, Procs=%d, Factors=%f, Phy=%d\n", i, user_config->iters[i], user_config->procs[i], user_config->factors[i], user_config->phy_dist[i]); } } } /* * Imprime por salida estandar la informacion relacionada con un * solo grupo de procesos en su configuracion. */ void print_config_group(configuration *user_config, int grp) { if(user_config != NULL) { int parents, sons; parents = sons = 0; if(grp > 0) { parents = user_config->procs[grp-1]; } if(grp < user_config->resizes - 1) { sons = user_config->procs[grp+1]; } printf("Config: matrix=%d, comm_tam=%d, sdr=%d, adr=%d, aib=%d, css=%d, cst=%d, time=%f\n", user_config->matrix_tam, user_config->comm_tam, user_config->sdr, user_config->adr, user_config->aib, user_config->css, user_config->cst, user_config->general_time); printf("Config Group: iters=%d, factor=%f, phy=%d, procs=%d, parents=%d, sons=%d\n", user_config->iters[grp], user_config->factors[grp], user_config->phy_dist[grp], user_config->procs[grp], parents, sons); } } //||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| || //||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| || //| FUNCIONES DE INTERCOMUNICACION DE ESTRUCTURA DE CONFIGURACION || //||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| || //||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| |/ /* * Envia una estructura de configuracion al grupo de procesos al que se * enlaza este grupo a traves del intercomunicador pasado como argumento. * * Esta funcion tiene que ser llamada por todos los procesos del mismo grupo * e indicar cual es el proceso raiz que se encargara de enviar la * configuracion al otro grupo. */ void send_config_file(configuration *config_file, int root, MPI_Comm intercomm) { MPI_Datatype config_type, config_type_array; // Obtener un tipo derivado para enviar todos los // datos escalares con una sola comunicacion def_struct_config_file(config_file, &config_type); // Obtener un tipo derivado para enviar los tres vectores // de enteros con una sola comunicacion def_struct_config_file_array(config_file, &config_type_array); MPI_Bcast(config_file, 1, config_type, root, intercomm); MPI_Bcast(config_file, 1, config_type_array, root, intercomm); MPI_Bcast(config_file->factors, config_file->resizes, MPI_FLOAT, root, intercomm); //Liberar tipos derivados MPI_Type_free(&config_type); MPI_Type_free(&config_type_array); } /* * Recibe una estructura de configuracion desde otro grupo de procesos * y la devuelve. La memoria de la estructura se reserva en esta funcion. * * Esta funcion tiene que ser llamada por todos los procesos del mismo grupo * e indicar cual es el proceso raiz del otro grupo que se encarga de enviar * la configuracion a este grupo. * * La memoria de la configuracion devuelta tiene que ser liberada con * la funcion "free_config". */ configuration *recv_config_file(int root, MPI_Comm intercomm) { MPI_Datatype config_type, config_type_array; configuration *config_file = malloc(sizeof(configuration) * 1); // Obtener un tipo derivado para recibir todos los // datos escalares con una sola comunicacion def_struct_config_file(config_file, &config_type); MPI_Bcast(config_file, 1, config_type, root, intercomm); // Obtener un tipo derivado para enviar los tres vectores // de enteros con una sola comunicacion malloc_config_arrays(config_file, config_file->resizes); // Reserva de memoria de los vectores def_struct_config_file_array(config_file, &config_type_array); MPI_Bcast(config_file, 1, config_type_array, root, intercomm); MPI_Bcast(config_file->factors, config_file->resizes, MPI_FLOAT, root, intercomm); //Liberar tipos derivados MPI_Type_free(&config_type); MPI_Type_free(&config_type_array); return config_file; } /* * Tipo derivado para enviar 6 elementos especificos * de la estructura de configuracion con una sola comunicacion. */ void def_struct_config_file(configuration *config_file, MPI_Datatype *config_type) { int i, counts = 11; int blocklengths[11] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1}; MPI_Aint displs[counts], dir; MPI_Datatype types[counts]; // Rellenar vector types types[0] = types[1] = types[2] = types[3] = types[4] = types[5] = types[6] = types[7] = types[8] = MPI_INT; types[9] = MPI_FLOAT; types[10] = MPI_DOUBLE; // Rellenar vector displs MPI_Get_address(config_file, &dir); MPI_Get_address(&(config_file->resizes), &displs[0]); MPI_Get_address(&(config_file->actual_resize), &displs[1]); MPI_Get_address(&(config_file->matrix_tam), &displs[2]); MPI_Get_address(&(config_file->comm_tam), &displs[3]); MPI_Get_address(&(config_file->sdr), &displs[4]); MPI_Get_address(&(config_file->adr), &displs[5]); MPI_Get_address(&(config_file->aib), &displs[6]); MPI_Get_address(&(config_file->css), &displs[7]); MPI_Get_address(&(config_file->cst), &displs[8]); MPI_Get_address(&(config_file->general_time), &displs[9]); MPI_Get_address(&(config_file->Top), &displs[10]); for(i=0;iresizes; //Rellenar vector displs MPI_Get_address(config_file, &dir); MPI_Get_address(config_file->iters, &displs[0]); MPI_Get_address(config_file->procs, &displs[1]); MPI_Get_address(config_file->phy_dist, &displs[2]); for(i=0;i