mailleur/src/main2.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�e et Vincent FRANCOIS
//  D�artement de G�ie M�anique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�artement
//  de g�ie m�anique de l'Universit�du Qu�ec �
//  Trois Rivi�es
//  Les librairies ne peuvent �re utilis�s sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       main.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 �11H25
//------------------------------------------------------------
//------------------------------------------------------------
#include "gestionversion.h"
#ifdef WINDOWS_VERSION
#include "fenetre.h"
#endif

#pragma hdrstop

#include "mg_file.h"
#include "vct_multi_modele.h"
#include "tpl_octree.h"
#include "tpl_grille.h"
#include "ot_decalage_parametre.h"
#include "ot_cpu.h"
#include "fct_generateur_3d.h"
#include "mailleur_delaunay.h"
#include <fstream>
#include <string.h>
#include "mg_maillage_algos.h"
#include "magic_application.h"


//---------------------------------------------------------------------------
int main(int argc,char **argv)

{
std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
MAGIC_PARAMETRE_APPLICATION p1(2,(char*)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier de carte de taille");
MAGIC_PARAMETRE_APPLICATION p2(3,(char*)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Limite supérieure des mailles à optimiser - Défaut 2");
MAGIC_PARAMETRE_APPLICATION p3(4,(char*)"-criterearret",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Valeur maximale de la longueur des arêtes par rapport à la longueur initiale - Défaut 1.5");
MAGIC_PARAMETRE_APPLICATION p4(5,(char*)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ajout d'une analyse de la qualité du maillage en fin de programme");
MAGIC_PARAMETRE_APPLICATION p5(6,(char*)"-delaunay",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Maillage d'une boite par la méthode de Delaunay");
MAGIC_PARAMETRE_APPLICATION p6(7,(char*)"-octal",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Maillage d'une boite par la méthode octal");
MAGIC_PARAMETRE_APPLICATION p7(8,(char*)"-tailleoctree",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Taille des octants par rapport à la taille souhaitée - Défaut 2.");
MAGIC_PARAMETRE_APPLICATION p8(9,(char*)"-fem",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Création du maillage fem");
MAGIC_PARAMETRE_APPLICATION p9(10,(char*)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Degre du maillage fem - Défaut 1");
MAGIC_PARAMETRE_APPLICATION p10(11,(char*)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Valeur du pas d'integration - Defaut 5");
p5.ajouter_dependance(1);p5.ajouter_dependance(2);p5.ajouter_dependance(3);p5.ajouter_dependance(4),p5.ajouter_dependance(5);p5.ajouter_dependance(11);
p6.ajouter_dependance(1);p6.ajouter_dependance(2);p6.ajouter_dependance(8);p6.ajouter_dependance(9);p6.ajouter_dependance(10);
lst.push_back(p1);
lst.push_back(p2);
lst.push_back(p3);
lst.push_back(p4);
lst.push_back(p5);
lst.push_back(p6);
lst.push_back(p7);
lst.push_back(p8);
lst.push_back(p9);
lst.push_back(p10);

MAGIC_APPLICATION app((char*)"Mailleur d'un espace défini par une carte de taille",argc,argv,lst,false,true);
if (app.get_erreur()==true) return 0;  
 

if (app.get_action()==6)
  {
  char fichier[500];
  char fichiercarte[500];
  double critere;
  bool analyse;
  char message[500];
  int niveau;
  int pas;
  app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
  app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
  app.recupere_parametre_int_avec_defaut(3,niveau,2);
  app.recupere_parametre_double_avec_defaut(4,critere,1.5);
  app.recupere_parametre_bool_avec_defaut(6,analyse,false);
  app.recupere_parametre_int_avec_defaut(11,pas,5);
  if (app.get_erreur()==true) return 0;
  app.affiche((char*)"Méthode de Delaunay");
  sprintf(message,"Fichier carte : %s",fichiercarte);
  app.affiche(message);
  app.affiche((char*)"");
  app.affiche((char*)"Lecture carte de taille");
  OT_CPU compteur;  
  FCT_GENERATEUR_3D<4> carte;
  carte.lire(fichiercarte);
  MG_GESTIONNAIRE* gest;
  MG_FILE* file=new MG_FILE(fichier);
  if (file->get_code_de_lecture()==0)
        {
        delete file;
        MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
        gest=gest2;
        }
  else gest=file;
  MG_MAILLAGE* mai=new MG_MAILLAGE(NULL);
  gest->ajouter_mg_maillage(mai);
  MAILLEUR_DELAUNAY m(mai,&carte,critere,&compteur);  
  m.active_affichage(app.affiche);
  m.change_niveau_optimisation(niveau);
  m.change_pas_integration(pas);
  m.maille();
  char fich[500];
  strcpy(fich,fichier);
  gest->enregistrer(fich);
  if (analyse) MG_MAILLAGE_ALGOS::calcule_caracteristiques(mai,fich,&carte);
  strcat(fich,".magic");
  gest->enregistrer(fich);
  compteur.ajouter_etape("Enregistrement");
  std::vector<std::string> tab;
  compteur.get_tabfinal(tab);
  for (int i=0;i<tab.size();i++)
    app.affiche((char*)tab[i].c_str());
  }
if (app.get_action()==7)
  {
  char fichier[500];
  char fichiercarte[500];
  double coef;
  bool fem;
  int degre;
  app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
  app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
  app.recupere_parametre_double_avec_defaut(8,coef,2.);
  app.recupere_parametre_bool_avec_defaut(9,fem,false);
  app.recupere_parametre_int_avec_defaut(10,degre,1);
  if (app.get_erreur()==true) return 0;
  char message[500];
  app.affiche((char*)"Méthode octale");
  sprintf(message,"   Fichier carte : %s",fichiercarte);
  app.affiche(message);
  app.affiche((char*)"Lecture carte de taille");
  FCT_GENERATEUR_3D<4> carte;
  carte.lire(fichiercarte);
  app.affiche((char*)"Maillage");
  MG_GESTIONNAIRE* gest;
  MG_FILE* file=new MG_FILE(fichier);
  if (file->get_code_de_lecture()==0)
        {
        delete file;
        MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
        gest=gest2;
        }
  else gest=file;
  MG_MAILLAGE* mgmai=new MG_MAILLAGE(NULL);
  gest->ajouter_mg_maillage(mgmai);
  MAILLEUR3D mailleur(&carte,mgmai);
  mailleur.change_coef(coef);
  mailleur.maille();
  if (fem)
        {
        app.affiche((char*)"Maillage FEM");
        FEM_MAILLAGE* femmai=new FEM_MAILLAGE(mgmai->get_mg_geometrie(),mgmai,degre);
        gest->ajouter_fem_maillage(femmai);
        femmai->construire(0);
        }
  app.affiche((char*)"Enregistrement");
  gest->enregistrer(fichier);
  delete gest;
  app.affiche((char*)"Fin");  
    
    
  }


return 0;
}


#pragma package(smart_init)
Revision:	747
Committed:	Tue Oct 6 19:53:22 2015 UTC (9 years, 8 months ago) by francois
File size:	7197 byte(s)
Log Message:	gestion du multivolume a l'importation des fichiers OCC avec triangulation associée
#	User	Rev	Content
1	francois	288	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�e et Vincent FRANCOIS
5			// D�artement de G�ie M�anique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�artement
8			// de g�ie m�anique de l'Universit�du Qu�ec �
9			// Trois Rivi�es
10			// Les librairies ne peuvent �re utilis�s sans l'accord
11			// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// main.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 �11H25
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23			#include "gestionversion.h"
24			#ifdef WINDOWS_VERSION
25			#include "fenetre.h"
26			#endif
27
28			#pragma hdrstop
29
30			#include "mg_file.h"
31			#include "vct_multi_modele.h"
32			#include "tpl_octree.h"
33			#include "tpl_grille.h"
34			#include "ot_decalage_parametre.h"
35			#include "ot_cpu.h"
36	francois	481	#include "fct_generateur_3d.h"
37	francois	288	#include "mailleur_delaunay.h"
38			#include <fstream>
39			#include <string.h>
40			#include "mg_maillage_algos.h"
41	francois	432	#include "magic_application.h"
42
43
44	francois	288	//---------------------------------------------------------------------------
45	francois	432	int main(int argc,char **argv)
46	francois	288
47			{
48	francois	432	std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
49			MAGIC_PARAMETRE_APPLICATION p1(2,(char)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char)"Fichier de carte de taille");
50			MAGIC_PARAMETRE_APPLICATION p2(3,(char)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Limite supérieure des mailles à optimiser - Défaut 2");
51			MAGIC_PARAMETRE_APPLICATION p3(4,(char)"-criterearret",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Valeur maximale de la longueur des arêtes par rapport à la longueur initiale - Défaut 1.5");
52			MAGIC_PARAMETRE_APPLICATION p4(5,(char)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Ajout d'une analyse de la qualité du maillage en fin de programme");
53			MAGIC_PARAMETRE_APPLICATION p5(6,(char)"-delaunay",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode de Delaunay");
54			MAGIC_PARAMETRE_APPLICATION p6(7,(char)"-octal",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode octal");
55			MAGIC_PARAMETRE_APPLICATION p7(8,(char)"-tailleoctree",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Taille des octants par rapport à la taille souhaitée - Défaut 2.");
56			MAGIC_PARAMETRE_APPLICATION p8(9,(char)"-fem",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Création du maillage fem");
57			MAGIC_PARAMETRE_APPLICATION p9(10,(char)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Degre du maillage fem - Défaut 1");
58	francois	532	MAGIC_PARAMETRE_APPLICATION p10(11,(char)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Valeur du pas d'integration - Defaut 5");
59			p5.ajouter_dependance(1);p5.ajouter_dependance(2);p5.ajouter_dependance(3);p5.ajouter_dependance(4),p5.ajouter_dependance(5);p5.ajouter_dependance(11);
60	francois	432	p6.ajouter_dependance(1);p6.ajouter_dependance(2);p6.ajouter_dependance(8);p6.ajouter_dependance(9);p6.ajouter_dependance(10);
61			lst.push_back(p1);
62			lst.push_back(p2);
63			lst.push_back(p3);
64			lst.push_back(p4);
65			lst.push_back(p5);
66			lst.push_back(p6);
67			lst.push_back(p7);
68			lst.push_back(p8);
69			lst.push_back(p9);
70	francois	532	lst.push_back(p10);
71	francois	288
72	francois	485	MAGIC_APPLICATION app((char*)"Mailleur d'un espace défini par une carte de taille",argc,argv,lst,false,true);
73	francois	432	if (app.get_erreur()==true) return 0;
74
75	francois	288
76
77	francois	432	if (app.get_action()==6)
78			{
79			char fichier[500];
80			char fichiercarte[500];
81			double critere;
82			bool analyse;
83			char message[500];
84			int niveau;
85	francois	532	int pas;
86	francois	432	app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
87			app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
88			app.recupere_parametre_int_avec_defaut(3,niveau,2);
89			app.recupere_parametre_double_avec_defaut(4,critere,1.5);
90			app.recupere_parametre_bool_avec_defaut(6,analyse,false);
91	francois	532	app.recupere_parametre_int_avec_defaut(11,pas,5);
92	francois	432	if (app.get_erreur()==true) return 0;
93			app.affiche((char*)"Méthode de Delaunay");
94			sprintf(message,"Fichier carte : %s",fichiercarte);
95			app.affiche(message);
96			app.affiche((char*)"");
97			app.affiche((char*)"Lecture carte de taille");
98			OT_CPU compteur;
99			FCT_GENERATEUR_3D<4> carte;
100			carte.lire(fichiercarte);
101			MG_GESTIONNAIRE* gest;
102			MG_FILE* file=new MG_FILE(fichier);
103			if (file->get_code_de_lecture()==0)
104			{
105			delete file;
106			MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
107			gest=gest2;
108			}
109			else gest=file;
110			MG_MAILLAGE* mai=new MG_MAILLAGE(NULL);
111			gest->ajouter_mg_maillage(mai);
112			MAILLEUR_DELAUNAY m(mai,&carte,critere,&compteur);
113			m.active_affichage(app.affiche);
114			m.change_niveau_optimisation(niveau);
115	francois	532	m.change_pas_integration(pas);
116	francois	432	m.maille();
117			char fich[500];
118			strcpy(fich,fichier);
119			gest->enregistrer(fich);
120			if (analyse) MG_MAILLAGE_ALGOS::calcule_caracteristiques(mai,fich,&carte);
121			strcat(fich,".magic");
122			gest->enregistrer(fich);
123			compteur.ajouter_etape("Enregistrement");
124			std::vector<std::string> tab;
125			compteur.get_tabfinal(tab);
126			for (int i=0;i<tab.size();i++)
127			app.affiche((char*)tab[i].c_str());
128			}
129			if (app.get_action()==7)
130			{
131			char fichier[500];
132			char fichiercarte[500];
133			double coef;
134			bool fem;
135			int degre;
136			app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
137			app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
138			app.recupere_parametre_double_avec_defaut(8,coef,2.);
139			app.recupere_parametre_bool_avec_defaut(9,fem,false);
140			app.recupere_parametre_int_avec_defaut(10,degre,1);
141			if (app.get_erreur()==true) return 0;
142			char message[500];
143			app.affiche((char*)"Méthode octale");
144			sprintf(message," Fichier carte : %s",fichiercarte);
145			app.affiche(message);
146			app.affiche((char*)"Lecture carte de taille");
147			FCT_GENERATEUR_3D<4> carte;
148			carte.lire(fichiercarte);
149			app.affiche((char*)"Maillage");
150			MG_GESTIONNAIRE* gest;
151			MG_FILE* file=new MG_FILE(fichier);
152			if (file->get_code_de_lecture()==0)
153			{
154			delete file;
155			MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
156			gest=gest2;
157			}
158			else gest=file;
159			MG_MAILLAGE* mgmai=new MG_MAILLAGE(NULL);
160			gest->ajouter_mg_maillage(mgmai);
161			MAILLEUR3D mailleur(&carte,mgmai);
162			mailleur.change_coef(coef);
163			mailleur.maille();
164			if (fem)
165			{
166			app.affiche((char*)"Maillage FEM");
167			FEM_MAILLAGE* femmai=new FEM_MAILLAGE(mgmai->get_mg_geometrie(),mgmai,degre);
168			gest->ajouter_fem_maillage(femmai);
169			femmai->construire(0);
170			}
171			app.affiche((char*)"Enregistrement");
172			gest->enregistrer(fichier);
173			delete gest;
174			app.affiche((char*)"Fin");
175
176
177
178			}
179	francois	288
180
181
182			return 0;
183			}
184
185
186
187			#pragma package(smart_init)