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 faces::cargaCaras() {
// string caraNombre = "capturas/imagen_0";
// cara.setup(caraNombre);
// caras.push_back(cara);
for(int i=0; i<NCARAS; i++) {
string caraNombre = "capturas/imagen_"+ofToString(i%6);
Cara caraTmp;
caraTmp.setup(caraNombre);
caras.push_back(caraTmp);
}
swWireframe=true;
swEnLinea=true;
fuente.loadFont("Helvetica-Bold.otf", 11);
}
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 Malla::drawNormals(){
for(int i = 0; i < numCaras;i++){
Cara* c = caras[i];
int v0 = c->getIndice(0)->get_indiceVertice();
int v1 = c->getIndice(1)->get_indiceVertice();
int v2 = c->getIndice(2)->get_indiceVertice();
int v3 = c->getIndice(3)->get_indiceVertice();
int n = c->getIndice(0)->get_indiceNormal();
PV3D* normal = normales[n];
GLdouble vX = (vertices[v0]->getX()+vertices[v1]->getX()+vertices[v2]->getX()+vertices[v3]->getX())/4;
GLdouble vY = (vertices[v0]->getY()+vertices[v1]->getY()+vertices[v2]->getY()+vertices[v3]->getY())/4;
GLdouble vZ = (vertices[v0]->getZ()+vertices[v1]->getZ()+vertices[v2]->getZ()+vertices[v3]->getZ())/4;
glBegin(GL_LINES);
glVertex3d(vX,vY,vZ);
glVertex3d(normales[n]->getX()+vX,normales[n]->getY()+vY,normales[n]->getZ()+vZ);
glEnd();
}
}
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);
Nodos *nds = malla->getNodos();
Caras *crs = malla->getCaras();
Nodo *nodo;
int* nuevoIndiceNodos=new int[nds->getNumNodos()];
int nlinea=0;
for(int i=0; i<nds->getNumNodos(); i++) {
nodo = nds->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++;
}
}
char separador2[32] = "% Elements (triangular)\n";
fprintf(archivo_x,"%s",separador2);
Cara *c;
vector<int> ind_nodos;
for(int i=0; i<crs->getNumCaras(); i++) {
c = crs->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.
fprintf(archivo_x,"%d %d %d\n",nuevoIndiceNodos[ind_nodos[2]]+1,nuevoIndiceNodos[ind_nodos[1]]+1,nuevoIndiceNodos[ind_nodos[0]]+1);
}
}
delete nuevoIndiceNodos;
fclose(archivo_x);
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, "");
}
void AlmacenaToXYZU::aplicar(Malla *malla, string filename){
//Lineas necesarias para que scanf lea archivo en computadores seteados en otro lenguaje.
setlocale(LC_NUMERIC, "POSIX");
int posextension=filename.rfind(".");
string nombrefile=filename.substr(0,posextension);
string extension=filename.substr(posextension+1,filename.length());
if (posextension==0){
nombrefile=filename;
extension="txt";
}
stringstream archivox;
stringstream archivoy;
stringstream archivoz;
stringstream archivou;
archivox << nombrefile << "_nx." << extension;
archivoy << nombrefile << "_ny." << extension;
archivoz << nombrefile << "_nz." << extension;
archivou << nombrefile << "_u." << extension;
FILE *archivo_x = fopen(archivox.str().c_str(),"w");
FILE *archivo_y = fopen(archivoy.str().c_str(),"w");
FILE *archivo_z = fopen(archivoz.str().c_str(),"w");
FILE *archivo_u = fopen(archivou.str().c_str(),"w");
char separador1[32] = "% Coordinates\n";
fprintf(archivo_x,"%s",separador1);
fprintf(archivo_y,"%s",separador1);
fprintf(archivo_z,"%s",separador1);
fprintf(archivo_u,"%s",separador1);
Nodos *nds = malla->getNodos();
Caras *crs = malla->getCaras();
Nodo *nodo;
int* nuevoIndiceNodos=new int[nds->getNumNodos()];
int nlinea=0;
for(int i=0; i<nds->getNumNodos(); i++) {
nodo = nds->getNodo(i);
if(nodo != NULL) {
fprintf(archivo_x,"%f %f %f\n",nodo->getPunto()->getX(),nodo->getPunto()->getY(),nodo->getPunto()->getZ());
fprintf(archivo_y,"%f %f %f\n",nodo->getPunto()->getX(),nodo->getPunto()->getY(),nodo->getPunto()->getZ());
fprintf(archivo_z,"%f %f %f\n",nodo->getPunto()->getX(),nodo->getPunto()->getY(),nodo->getPunto()->getZ());
fprintf(archivo_u,"%f %f %f\n",nodo->getPunto()->getX(),nodo->getPunto()->getY(),nodo->getPunto()->getZ());
nuevoIndiceNodos[i]=nlinea;
nlinea++;
}
}
char separador2[32] = "% Elements (triangular)\n";
fprintf(archivo_x,"%s",separador2);
fprintf(archivo_y,"%s",separador2);
fprintf(archivo_z,"%s",separador2);
fprintf(archivo_u,"%s",separador2);
Cara *c;
vector<int> ind_nodos;
for(int i=0; i<crs->getNumCaras(); i++) {
c = crs->getCara(i);
if(c != 0) {
ind_nodos = c->getNodos();
fprintf(archivo_x,"%d %d %d\n",nuevoIndiceNodos[ind_nodos[0]]+1,nuevoIndiceNodos[ind_nodos[1]]+1,nuevoIndiceNodos[ind_nodos[2]]+1);
fprintf(archivo_y,"%d %d %d\n",nuevoIndiceNodos[ind_nodos[0]]+1,nuevoIndiceNodos[ind_nodos[1]]+1,nuevoIndiceNodos[ind_nodos[2]]+1);
fprintf(archivo_z,"%d %d %d\n",nuevoIndiceNodos[ind_nodos[0]]+1,nuevoIndiceNodos[ind_nodos[1]]+1,nuevoIndiceNodos[ind_nodos[2]]+1);
fprintf(archivo_u,"%d %d %d\n",nuevoIndiceNodos[ind_nodos[0]]+1,nuevoIndiceNodos[ind_nodos[1]]+1,nuevoIndiceNodos[ind_nodos[2]]+1);
}
}
char separador_x[32] = "% Data (nx)\n";
char separador_y[32] = "% Data (ny)\n";
char separador_z[32] = "% Data (nz)\n";
char separador_u[32] = "% Data (u)\n";
fprintf(archivo_x,"%s",separador_x);
fprintf(archivo_y,"%s",separador_y);
fprintf(archivo_z,"%s",separador_z);
fprintf(archivo_u,"%s",separador_u);
for(int i=0; i<nds->getNumNodos(); i++) {
nodo = nds->getNodo(i);
// cout << nodo->imprimir() << endl;
if(nodo != NULL) {
fprintf(archivo_x,"%f\n",nodo->getNormal()->getPunto().getX());
fprintf(archivo_y,"%f\n",nodo->getNormal()->getPunto().getY());
fprintf(archivo_z,"%f\n",nodo->getNormal()->getPunto().getZ());
fprintf(archivo_u,"%f\n",nodo->getConcentracion());
}
}
delete nuevoIndiceNodos;
fclose(archivo_x);
fclose(archivo_y);
fclose(archivo_z);
fclose(archivo_u);
setlocale(LC_NUMERIC, "");
}
