void AlmacenaToComsol_1::aplicar(Malla *malla, string filename){ //Lineas necesarias para que scanf lea archivo en computadores seteados en otro lenguaje. setlocale(LC_NUMERIC, "POSIX"); FILE *archivo_x = fopen(filename.c_str(),"w"); char separador1[32] = "% Coordinates\n"; fprintf(archivo_x,"%s",separador1); Nodo *nodo; int* nuevoIndiceNodos=new int[malla->getMaxIndiceNodos()+1]; int nlinea=0; for(int i=0; i<=malla->getMaxIndiceNodos(); i++) { nodo = malla->getNodo(i); if(nodo != NULL) { fprintf(archivo_x,"%f %f %f\n",nodo->getPunto().getX(),nodo->getPunto().getY(),nodo->getPunto().getZ()); nuevoIndiceNodos[i]=nlinea; nlinea++; } } ostringstream myStream; myStream << "% Elements (" << malla->getTipoMalla() << ")\n"; string str_separador2 = myStream.str(); const char *separador2 = str_separador2.c_str(); fprintf(archivo_x,"%s",separador2); Cara *c; vector<int> ind_nodos; for(int i=0; i<=malla->getMaxIndiceCaras(); i++) { c = malla->getCara(i); if(c != 0) { ind_nodos = c->getNodos(); //Nota: Se invierte el orden de insercion de los puntos, para que las normales queden seteadas correctamente. //ya que en el archivo cms_1 los nodos estan al reves, dejando las normales hacia adentro. //for(int j=c->getNumElem()-1; j>=0 ;--j) { // Por bug archivos de entrada for(int j=0; j<c->getNumElem(); ++j) { // así debiera ser char *delimitador = " "; //if(j == 0) // Por bug archivos de entrada if(j == c->getNumElem()-1) // así debiera ser delimitador = "\n"; fprintf(archivo_x,"%d%s",nuevoIndiceNodos[ind_nodos[j]]+1,delimitador); } } } delete nuevoIndiceNodos; fclose(archivo_x); setlocale(LC_NUMERIC, ""); }
void AlmacenaToOff::aplicar(Malla *malla, string filename){ //Lineas necesarias para que scanf lea archivo en computadores seteados en otro lenguaje. setlocale(LC_NUMERIC, "POSIX"); FILE *archivoSalida = fopen(filename.c_str(),"w"); fprintf(archivoSalida,"OFF\n"); fprintf(archivoSalida,"%d %d %d\n",malla->getNumNodos(),malla->getNumCaras(),malla->getNumArcos()); Nodo *nodo; int* nuevoIndiceNodos=new int[malla->getMaxIndiceNodos()+1]; int nlinea=0; for(int i=0; i<=malla->getMaxIndiceNodos(); i++) { nodo = malla->getNodo(i); if(nodo != NULL) { fprintf(archivoSalida,"%f %f %f\n",nodo->getPunto().getX(),nodo->getPunto().getY(),nodo->getPunto().getZ()); nuevoIndiceNodos[i]=nlinea; nlinea++; } } Cara *c; vector<int> ind_nodos; vector<int> ind_arcos; for(int i=0; i<=malla->getMaxIndiceCaras(); i++) { c = malla->getCara(i); if(c != 0) { ind_nodos = c->getNodos(); ind_arcos = c->getArcos(); fprintf(archivoSalida,"%d ",c->getNumElem()); double conc_prom = 0; for(int j=0; j<(int)ind_nodos.size(); j++) { conc_prom = conc_prom + malla->getNodo(ind_nodos[j])->getConcentracion(); fprintf(archivoSalida,"%d ",nuevoIndiceNodos[ind_nodos[j]]); } fprintf(archivoSalida,"\n"); } } delete nuevoIndiceNodos; fclose(archivoSalida); setlocale(LC_NUMERIC, ""); }
void AlmacenaToMllDebug::aplicar(Malla *malla, string filename){ //Lineas necesarias para que scanf lea archivo en computadores seteados en otro lenguaje. setlocale(LC_NUMERIC, "POSIX"); FILE *archivoSalida = fopen(filename.c_str(),"w"); fprintf(archivoSalida,"MLL (Debug file)\n"); fprintf(archivoSalida,"%d %d %d\n",malla->getMaxIndiceNodos()+1,malla->getMaxIndiceCaras()+1,malla->getMaxIndiceArcos()+1); fprintf(archivoSalida,"%d %d %d\n",malla->getNumNodos(),malla->getNumCaras(),malla->getNumArcos()); Nodo *nodo; Vect normal; cout << "Recorriendo Nodos" << endl; for(int i=0; i<=malla->getMaxIndiceNodos(); i++) { nodo = malla->getNodo(i); if(nodo == 0) { fprintf(archivoSalida,"Nodo{%d}=Vacio \n",i); } else { fprintf(archivoSalida,"Nodo{%d}=(%f,%f,%f) ",i,nodo->getPunto().getX(),nodo->getPunto().getY(),nodo->getPunto().getZ()); vector<int> caras = nodo->getCaras(); vector<int> arcos = nodo->getArcos(); fprintf(archivoSalida,"Caras{"); for(int j=0; j<(int)caras.size(); j++) { fprintf(archivoSalida,"%d",caras[j]); if(j<(int)caras.size()-1) { fprintf(archivoSalida,","); } } fprintf(archivoSalida,"} Arcos{"); for(int j=0; j<(int)arcos.size(); j++) { fprintf(archivoSalida,"%d",arcos[j]); if(j<(int)arcos.size()-1) { fprintf(archivoSalida,","); } } fprintf(archivoSalida,"} "); normal = malla->getNodo(i)->getNormal(); fprintf(archivoSalida,"Normal=(%f,%f,%f) ",normal.getPunto().getX(),normal.getPunto().getY(),normal.getPunto().getZ()); fprintf(archivoSalida,"Concentracion={%f}",nodo->getConcentracion()); fprintf(archivoSalida,"\n"); } } cout << "Pasamos nodos" << endl; Arco *a; for(int i=0; i<=malla->getMaxIndiceArcos(); i++) { a = malla->getArco(i); //si se elimino un arco, puede que en esa posicion exista un puntero igual a 0, por eso nos aseguramos que el arco sea distinto de 0. if(a == 0) { fprintf(archivoSalida,"Arco{%d} Vacio\n",i); } else { fprintf(archivoSalida,"Arco{%d} Nodos{%d,%d} Caras{%d,%d}\n",i,a->getNodo1(),a->getNodo2(),a->getCara1(),a->getCara2()); } } Cara *c; vector<int> ind_nodos; vector<int> ind_arcos; for(int i=0; i<=malla->getMaxIndiceCaras(); i++) { c = malla->getCara(i); //si se elimino una cara, puede que en esa posicion exista un puntero igual a 0, por eso nos aseguramos que la cara sea distinta de 0. if(c == 0) { fprintf(archivoSalida,"Cara{%d} Vacia\n",i); } else { ind_nodos = c->getNodos(); ind_arcos = c->getArcos(); fprintf(archivoSalida,"Cara{%d} %d ",i,c->getNumElem()); fprintf(archivoSalida,"Nodos{"); double conc_prom = 0; for(int j=0; j<(int)ind_nodos.size(); j++) { conc_prom = conc_prom + malla->getNodo(ind_nodos[j])->getConcentracion(); fprintf(archivoSalida,"%d",ind_nodos[j]); if(j<(int)ind_nodos.size()-1) { fprintf(archivoSalida,","); } } fprintf(archivoSalida,"} "); fprintf(archivoSalida,"Arcos{"); for(int j=0; j<(int)ind_arcos.size(); j++) { if(malla->getArco(ind_arcos[j]) != 0) { fprintf(archivoSalida,"%d:%d-->%d",ind_arcos[j],malla->getArco(ind_arcos[j])->getNodo1(),malla->getArco(ind_arcos[j])->getNodo2()); if(j<(int)ind_arcos.size()-1) { fprintf(archivoSalida,","); } } } fprintf(archivoSalida,"} "); Vect normal = c->getNormal(malla); fprintf(archivoSalida,"Normal=(%f,%f,%f) ",normal.getPunto().getX(),normal.getPunto().getY(),normal.getPunto().getZ()); fprintf(archivoSalida,"Concentracion=%f",conc_prom); fprintf(archivoSalida,"\n"); } } fclose(archivoSalida); setlocale(LC_NUMERIC, ""); }