mailleur/src/main3.cpp

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       main3.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:57:20 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "magic_application.h"
#include <string.h>
#include "mg_file.h"
#include "fct_generateur_3d.h"
#include "fct_taille_fem_solution.h"
#include "mailleur3d_peau.h"
#include "mailleur3d_maillage_partiel.h"

int main(int argc,char **argv)
{
std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
MAGIC_PARAMETRE_APPLICATION p1(3,(char*)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier de carte de taille");lst.push_back(p1);
MAGIC_PARAMETRE_APPLICATION p2(4,(char*)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Limite supérieure des mailles à optimiser - Défaut 2");lst.push_back(p2);
MAGIC_PARAMETRE_APPLICATION p3(5,(char*)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de maillage à mailler - Défaut le premier");lst.push_back(p3);
MAGIC_PARAMETRE_APPLICATION p4(6,(char*)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ajout d'une analyse de la qualité du maillage en fin de programme");lst.push_back(p4);
MAGIC_PARAMETRE_APPLICATION p5(7,(char*)"-enregistrement",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Enregistrement de la dernière étape (pour debbugage) - Defaut false");lst.push_back(p5);
MAGIC_PARAMETRE_APPLICATION p6(8,(char*)"-numvol",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de volume à mailler - Défaut le premier ou NULL s il n y a en pas");lst.push_back(p6);
MAGIC_PARAMETRE_APPLICATION p7(9,(char*)"-carte2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Utilisation de la nouvelle version de carte de taille");lst.push_back(p7);
MAGIC_PARAMETRE_APPLICATION p8(10,(char*)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Valeur du pas d'integration - Defaut 32");lst.push_back(p8);
MAGIC_PARAMETRE_APPLICATION p9(11,(char*)"-continue",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Continuation d'un maillage peau arrete -Default non");lst.push_back(p9);
MAGIC_PARAMETRE_APPLICATION p10(12,(char*)"-prioritemetrique",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Valeur de la prise en compte de la métrique versus la qualité - Défaut 0.65");lst.push_back(p10);

MAGIC_APPLICATION app((char*)"Mailleur 3D d'une enveloppe triangulee",argc,argv,lst,true,true);
if (app.get_erreur()==true) return 0;  
 
char fichier[500];
char fichiersortie[500];
char fichiersortie2[500];
char fichiercarte[500];
int niveau;
int nummai;
int numvol;
bool analyse;
bool save;
bool version;
int pas;
bool recover;
double prioritemetrique;
app.recupere_parametre_string(1,fichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
app.recupere_parametre_string_avec_defaut(2,fichiersortie,fichier);
app.recupere_parametre_string(3,fichiercarte,(char*)"Fichier carte manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
app.recupere_parametre_int_avec_defaut(4,niveau,2);
app.recupere_parametre_int_avec_defaut(5,nummai,0);
app.recupere_parametre_bool_avec_defaut(6,analyse,false);
app.recupere_parametre_bool_avec_defaut(7,save,false);
app.recupere_parametre_int_avec_defaut(8,numvol,0);
app.recupere_parametre_bool_avec_defaut(9,version,false);
app.recupere_parametre_int_avec_defaut(10,pas,32);
app.recupere_parametre_bool_avec_defaut(11,recover,false);
app.recupere_parametre_double_avec_defaut(12,prioritemetrique,0.65);
 
  strcpy(fichiersortie2,fichiersortie);
  strcat(fichiersortie2,"dbg");
if (app.get_erreur()==true) return 0; 
if (recover==true)
      {
        int versioncarte;
        if (version==false) versioncarte=MAGIC::VERSION_CARTE_TAILLE::ANCIENNE; else versioncarte=MAGIC::VERSION_CARTE_TAILLE::NOUVELLE;
        MAILLEUR3D_MAILLAGE_PARTIEL m3d(fichier,fichiersortie,fichiercarte,versioncarte,nummai,numvol,save);
        m3d.active_affichage(app.affiche);
        m3d.change_niveau_optimisation(niveau);
        m3d.change_pas_integration(pas);
        m3d.change_priorite_metrique(prioritemetrique);
        if (save==true) m3d.change_nom_fichierdbg(fichiersortie2);
        m3d.maille();
      }
else 
{
  app.affiche((char*)"Lecture fichier magic");
MG_FILE gest(fichier);
app.affiche((char*)"Lecture carte de taille");
FCT_TAILLE *metrique;
if (!version)
  {
  FCT_GENERATEUR_3D<4> *carte=new FCT_GENERATEUR_3D<4>;
  carte->lire(fichiercarte);
  metrique=carte;
  }
else
  {
  FCT_TAILLE_FEM_SOLUTION *carte=new FCT_TAILLE_FEM_SOLUTION(fichiercarte);
  metrique=carte;
  }
MG_MAILLAGE* mai;
if (nummai==0) mai=gest.get_mg_maillage(0);
else mai=gest.get_mg_maillageid(nummai);
MG_GEOMETRIE* mggeo=mai->get_mg_geometrie();
MG_VOLUME *mgvol=NULL;
if (mggeo!=NULL)
  {
  if (numvol==0) mgvol=mggeo->get_mg_volume(0);
  else mgvol=mggeo->get_mg_volumeid(numvol);
  }
MAILLEUR3D_PEAU m3d(mai,NULL,NULL,metrique,save);
if (save) 
          {
          m3d.change_nom_fichierdbg(fichiersortie2);
          }
m3d.active_affichage(app.affiche);
m3d.change_niveau_optimisation(niveau);
m3d.change_pas_integration(pas);
m3d.change_priorite_metrique(prioritemetrique);
app.affiche((char*)"Maillage 3D");
m3d.maille(mgvol,NULL,NULL);
if (analyse==true)
                {
                app.affiche((char*)"Analyse maillage");
                double vol;
                int tab[4];
                 m3d.analyse_maillage_obtenu(vol,tab);
                char mess[100];
                sprintf(mess,"    1) Respect de la carte de taille");
                app.affiche(mess);
                sprintf(mess,"    nombre de tetra theorique %f, nombre de tetra pratique %d ",vol,mai->get_nb_mg_tetra());
                app.affiche(mess);
                sprintf(mess,"    erreur : %.2f%%",(mai->get_nb_mg_tetra()-vol)*100./vol);
                app.affiche(mess);
                sprintf(mess,"    2) Qualite du maillage");
                app.affiche(mess);
                sprintf(mess,"    nombre de tetra bon (qualite>0.5)        : %d (%.2f%%)",tab[3],tab[3]*100./mai->get_nb_mg_tetra());
                app.affiche(mess);
                sprintf(mess,"    nombre de tetra acceptable (qualite>0.2) : %d (%.2f%%) ",tab[2],tab[2]*100./mai->get_nb_mg_tetra());
                app.affiche(mess);
                sprintf(mess,"    nombre de tetra moyen (qualite>0.1)      : %d (%.2f%%) ",tab[1],tab[1]*100./mai->get_nb_mg_tetra());
                app.affiche(mess);
                sprintf(mess,"    nombre de tetra mauvais (qualite>0.)     : %d (%.2f%%) ",tab[0],tab[0]*100./mai->get_nb_mg_tetra());
                app.affiche(mess);
                }
app.affiche((char*)"Enregistrement");
MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO=3;
gest.enregistrer(fichiersortie);
delete metrique;
}
app.affiche((char*)"Fin");
return 0;
}


Revision:	1157
Committed:	Thu Jun 13 22:18:27 2024 UTC (15 months, 3 weeks ago) by francois
File size:	7411 byte(s)
Log Message:	compatibilité Ubuntu 22.04 Suppression des refeences à Windows Ajout d'une banière
#	User	Rev	Content
1	francois	1157	//####//------------------------------------------------------------
2			//####//------------------------------------------------------------
3			//####// MAGiC
4			//####// Jean Christophe Cuilliere et Vincent FRANCOIS
5			//####// Departement de Genie Mecanique - UQTR
6			//####//------------------------------------------------------------
7			//####// MAGIC est un projet de recherche de l equipe ERICCA
8			//####// du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
9			//####// http://www.uqtr.ca/ericca
10			//####// http://www.uqtr.ca/
11			//####//------------------------------------------------------------
12			//####//------------------------------------------------------------
13			//####//
14			//####// main3.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:57:20 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22	francois	422	#include "gestionversion.h"
23	francois	432	#include "magic_application.h"
24			#include <string.h>
25	francois	422	#include "mg_file.h"
26	francois	481	#include "fct_generateur_3d.h"
27	francois	531	#include "fct_taille_fem_solution.h"
28	francois	422	#include "mailleur3d_peau.h"
29	francois	558	#include "mailleur3d_maillage_partiel.h"
30	francois	422
31	francois	432	int main(int argc,char **argv)
32	francois	422	{
33	francois	432	std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
34			MAGIC_PARAMETRE_APPLICATION p1(3,(char)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char)"Fichier de carte de taille");lst.push_back(p1);
35			MAGIC_PARAMETRE_APPLICATION p2(4,(char)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Limite supérieure des mailles à optimiser - Défaut 2");lst.push_back(p2);
36	francois	531	MAGIC_PARAMETRE_APPLICATION p3(5,(char)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Numero de maillage à mailler - Défaut le premier");lst.push_back(p3);
37	francois	432	MAGIC_PARAMETRE_APPLICATION p4(6,(char)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Ajout d'une analyse de la qualité du maillage en fin de programme");lst.push_back(p4);
38	francois	494	MAGIC_PARAMETRE_APPLICATION p5(7,(char)"-enregistrement",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Enregistrement de la dernière étape (pour debbugage) - Defaut false");lst.push_back(p5);
39	francois	531	MAGIC_PARAMETRE_APPLICATION p6(8,(char)"-numvol",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Numero de volume à mailler - Défaut le premier ou NULL s il n y a en pas");lst.push_back(p6);
40			MAGIC_PARAMETRE_APPLICATION p7(9,(char)"-carte2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Utilisation de la nouvelle version de carte de taille");lst.push_back(p7);
41	francois	532	MAGIC_PARAMETRE_APPLICATION p8(10,(char)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Valeur du pas d'integration - Defaut 32");lst.push_back(p8);
42	francois	558	MAGIC_PARAMETRE_APPLICATION p9(11,(char)"-continue",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Continuation d'un maillage peau arrete -Default non");lst.push_back(p9);
43			MAGIC_PARAMETRE_APPLICATION p10(12,(char)"-prioritemetrique",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Valeur de la prise en compte de la métrique versus la qualité - Défaut 0.65");lst.push_back(p10);
44	francois	422
45	francois	432	MAGIC_APPLICATION app((char*)"Mailleur 3D d'une enveloppe triangulee",argc,argv,lst,true,true);
46			if (app.get_erreur()==true) return 0;
47
48	francois	422	char fichier[500];
49			char fichiersortie[500];
50	francois	494	char fichiersortie2[500];
51	francois	422	char fichiercarte[500];
52	francois	432	int niveau;
53			int nummai;
54	francois	531	int numvol;
55	francois	432	bool analyse;
56	francois	494	bool save;
57	francois	531	bool version;
58	francois	532	int pas;
59	francois	558	bool recover;
60			double prioritemetrique;
61	francois	1051	app.recupere_parametre_string(1,fichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
62	francois	462	app.recupere_parametre_string_avec_defaut(2,fichiersortie,fichier);
63	francois	1051	app.recupere_parametre_string(3,fichiercarte,(char*)"Fichier carte manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
64	francois	462	app.recupere_parametre_int_avec_defaut(4,niveau,2);
65			app.recupere_parametre_int_avec_defaut(5,nummai,0);
66			app.recupere_parametre_bool_avec_defaut(6,analyse,false);
67	francois	494	app.recupere_parametre_bool_avec_defaut(7,save,false);
68	francois	531	app.recupere_parametre_int_avec_defaut(8,numvol,0);
69			app.recupere_parametre_bool_avec_defaut(9,version,false);
70	francois	532	app.recupere_parametre_int_avec_defaut(10,pas,32);
71	francois	558	app.recupere_parametre_bool_avec_defaut(11,recover,false);
72			app.recupere_parametre_double_avec_defaut(12,prioritemetrique,0.65);
73
74			strcpy(fichiersortie2,fichiersortie);
75			strcat(fichiersortie2,"dbg");
76	francois	462	if (app.get_erreur()==true) return 0;
77	francois	558	if (recover==true)
78			{
79			int versioncarte;
80			if (version==false) versioncarte=MAGIC::VERSION_CARTE_TAILLE::ANCIENNE; else versioncarte=MAGIC::VERSION_CARTE_TAILLE::NOUVELLE;
81			MAILLEUR3D_MAILLAGE_PARTIEL m3d(fichier,fichiersortie,fichiercarte,versioncarte,nummai,numvol,save);
82			m3d.active_affichage(app.affiche);
83			m3d.change_niveau_optimisation(niveau);
84			m3d.change_pas_integration(pas);
85			m3d.change_priorite_metrique(prioritemetrique);
86			if (save==true) m3d.change_nom_fichierdbg(fichiersortie2);
87			m3d.maille();
88			}
89			else
90			{
91			app.affiche((char*)"Lecture fichier magic");
92	francois	422	MG_FILE gest(fichier);
93	francois	432	app.affiche((char*)"Lecture carte de taille");
94	francois	531	FCT_TAILLE *metrique;
95			if (!version)
96			{
97			FCT_GENERATEUR_3D<4> *carte=new FCT_GENERATEUR_3D<4>;
98			carte->lire(fichiercarte);
99			metrique=carte;
100			}
101			else
102			{
103			FCT_TAILLE_FEM_SOLUTION *carte=new FCT_TAILLE_FEM_SOLUTION(fichiercarte);
104			metrique=carte;
105			}
106	francois	422	MG_MAILLAGE* mai;
107			if (nummai==0) mai=gest.get_mg_maillage(0);
108			else mai=gest.get_mg_maillageid(nummai);
109	francois	531	MG_GEOMETRIE* mggeo=mai->get_mg_geometrie();
110			MG_VOLUME *mgvol=NULL;
111			if (mggeo!=NULL)
112			{
113			if (numvol==0) mgvol=mggeo->get_mg_volume(0);
114			else mgvol=mggeo->get_mg_volumeid(numvol);
115			}
116			MAILLEUR3D_PEAU m3d(mai,NULL,NULL,metrique,save);
117	francois	494	if (save)
118			{
119			m3d.change_nom_fichierdbg(fichiersortie2);
120			}
121	francois	432	m3d.active_affichage(app.affiche);
122	francois	422	m3d.change_niveau_optimisation(niveau);
123	francois	532	m3d.change_pas_integration(pas);
124	francois	558	m3d.change_priorite_metrique(prioritemetrique);
125	francois	432	app.affiche((char*)"Maillage 3D");
126	francois	531	m3d.maille(mgvol,NULL,NULL);
127	francois	432	if (analyse==true)
128	francois	424	{
129	francois	432	app.affiche((char*)"Analyse maillage");
130	francois	424	double vol;
131			int tab[4];
132			m3d.analyse_maillage_obtenu(vol,tab);
133			char mess[100];
134			sprintf(mess," 1) Respect de la carte de taille");
135	francois	432	app.affiche(mess);
136	francois	424	sprintf(mess," nombre de tetra theorique %f, nombre de tetra pratique %d ",vol,mai->get_nb_mg_tetra());
137	francois	432	app.affiche(mess);
138	francois	424	sprintf(mess," erreur : %.2f%%",(mai->get_nb_mg_tetra()-vol)*100./vol);
139	francois	432	app.affiche(mess);
140	francois	424	sprintf(mess," 2) Qualite du maillage");
141	francois	432	app.affiche(mess);
142	francois	424	sprintf(mess," nombre de tetra bon (qualite>0.5) : %d (%.2f%%)",tab[3],tab[3]*100./mai->get_nb_mg_tetra());
143	francois	432	app.affiche(mess);
144	francois	424	sprintf(mess," nombre de tetra acceptable (qualite>0.2) : %d (%.2f%%) ",tab[2],tab[2]*100./mai->get_nb_mg_tetra());
145	francois	432	app.affiche(mess);
146	francois	424	sprintf(mess," nombre de tetra moyen (qualite>0.1) : %d (%.2f%%) ",tab[1],tab[1]*100./mai->get_nb_mg_tetra());
147	francois	432	app.affiche(mess);
148	francois	424	sprintf(mess," nombre de tetra mauvais (qualite>0.) : %d (%.2f%%) ",tab[0],tab[0]*100./mai->get_nb_mg_tetra());
149	francois	432	app.affiche(mess);
150	francois	424	}
151	francois	432	app.affiche((char*)"Enregistrement");
152	francois	422	MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO=3;
153			gest.enregistrer(fichiersortie);
154	francois	531	delete metrique;
155	francois	558	}
156	francois	432	app.affiche((char*)"Fin");
157	francois	422	return 0;
158			}
159
160
161
162	francois	432