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"
#include "mailleur_fem.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);
        MAILLEUR_FEM mf;
        mf.maille(femmai,0);
        }
  app.affiche((char*)"Enregistrement");
  gest->enregistrer(fichier);
  delete gest;
  app.affiche((char*)"Fin");  
    
    
  }


return 0;
}


#pragma package(smart_init)
Revision:	883
Committed:	Thu Apr 20 13:38:18 2017 UTC (8 years, 2 months ago) by francois
File size:	7241 byte(s)
Log Message:	Creation d'un mailleur FEM pour pouvoir avoir des stratégies paramétrées de maillage. L'ancienne méthode construit disparait et est remplacée par l'utilisation d'un MAILLEUR_FEM. Stratégie de quadratisation mise en place : déplacer les noeuds pour s'assurer que les tetras quadratiques présentent une distortion au dessu d'une borne inférieure. Mais ces noeuds quittent la géométrie. Les paramètres dans ~/.magic.
#	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	francois	883	#include "mailleur_fem.h"
43	francois	432
44
45	francois	288	//---------------------------------------------------------------------------
46	francois	432	int main(int argc,char **argv)
47	francois	288
48			{
49	francois	432	std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
50			MAGIC_PARAMETRE_APPLICATION p1(2,(char)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char)"Fichier de carte de taille");
51			MAGIC_PARAMETRE_APPLICATION p2(3,(char)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Limite supérieure des mailles à optimiser - Défaut 2");
52			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");
53			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");
54			MAGIC_PARAMETRE_APPLICATION p5(6,(char)"-delaunay",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode de Delaunay");
55			MAGIC_PARAMETRE_APPLICATION p6(7,(char)"-octal",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode octal");
56			MAGIC_PARAMETRE_APPLICATION p7(8,(char)"-tailleoctree",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Taille des octants par rapport à la taille souhaitée - Défaut 2.");
57			MAGIC_PARAMETRE_APPLICATION p8(9,(char)"-fem",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Création du maillage fem");
58			MAGIC_PARAMETRE_APPLICATION p9(10,(char)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Degre du maillage fem - Défaut 1");
59	francois	532	MAGIC_PARAMETRE_APPLICATION p10(11,(char)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Valeur du pas d'integration - Defaut 5");
60			p5.ajouter_dependance(1);p5.ajouter_dependance(2);p5.ajouter_dependance(3);p5.ajouter_dependance(4),p5.ajouter_dependance(5);p5.ajouter_dependance(11);
61	francois	432	p6.ajouter_dependance(1);p6.ajouter_dependance(2);p6.ajouter_dependance(8);p6.ajouter_dependance(9);p6.ajouter_dependance(10);
62			lst.push_back(p1);
63			lst.push_back(p2);
64			lst.push_back(p3);
65			lst.push_back(p4);
66			lst.push_back(p5);
67			lst.push_back(p6);
68			lst.push_back(p7);
69			lst.push_back(p8);
70			lst.push_back(p9);
71	francois	532	lst.push_back(p10);
72	francois	288
73	francois	485	MAGIC_APPLICATION app((char*)"Mailleur d'un espace défini par une carte de taille",argc,argv,lst,false,true);
74	francois	432	if (app.get_erreur()==true) return 0;
75
76	francois	288
77
78	francois	432	if (app.get_action()==6)
79			{
80			char fichier[500];
81			char fichiercarte[500];
82			double critere;
83			bool analyse;
84			char message[500];
85			int niveau;
86	francois	532	int pas;
87	francois	432	app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
88			app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
89			app.recupere_parametre_int_avec_defaut(3,niveau,2);
90			app.recupere_parametre_double_avec_defaut(4,critere,1.5);
91			app.recupere_parametre_bool_avec_defaut(6,analyse,false);
92	francois	532	app.recupere_parametre_int_avec_defaut(11,pas,5);
93	francois	432	if (app.get_erreur()==true) return 0;
94			app.affiche((char*)"Méthode de Delaunay");
95			sprintf(message,"Fichier carte : %s",fichiercarte);
96			app.affiche(message);
97			app.affiche((char*)"");
98			app.affiche((char*)"Lecture carte de taille");
99			OT_CPU compteur;
100			FCT_GENERATEUR_3D<4> carte;
101			carte.lire(fichiercarte);
102			MG_GESTIONNAIRE* gest;
103			MG_FILE* file=new MG_FILE(fichier);
104			if (file->get_code_de_lecture()==0)
105			{
106			delete file;
107			MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
108			gest=gest2;
109			}
110			else gest=file;
111			MG_MAILLAGE* mai=new MG_MAILLAGE(NULL);
112			gest->ajouter_mg_maillage(mai);
113			MAILLEUR_DELAUNAY m(mai,&carte,critere,&compteur);
114			m.active_affichage(app.affiche);
115			m.change_niveau_optimisation(niveau);
116	francois	532	m.change_pas_integration(pas);
117	francois	432	m.maille();
118			char fich[500];
119			strcpy(fich,fichier);
120			gest->enregistrer(fich);
121			if (analyse) MG_MAILLAGE_ALGOS::calcule_caracteristiques(mai,fich,&carte);
122			strcat(fich,".magic");
123			gest->enregistrer(fich);
124			compteur.ajouter_etape("Enregistrement");
125			std::vector<std::string> tab;
126			compteur.get_tabfinal(tab);
127			for (int i=0;i<tab.size();i++)
128			app.affiche((char*)tab[i].c_str());
129			}
130			if (app.get_action()==7)
131			{
132			char fichier[500];
133			char fichiercarte[500];
134			double coef;
135			bool fem;
136			int degre;
137			app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
138			app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
139			app.recupere_parametre_double_avec_defaut(8,coef,2.);
140			app.recupere_parametre_bool_avec_defaut(9,fem,false);
141			app.recupere_parametre_int_avec_defaut(10,degre,1);
142			if (app.get_erreur()==true) return 0;
143			char message[500];
144			app.affiche((char*)"Méthode octale");
145			sprintf(message," Fichier carte : %s",fichiercarte);
146			app.affiche(message);
147			app.affiche((char*)"Lecture carte de taille");
148			FCT_GENERATEUR_3D<4> carte;
149			carte.lire(fichiercarte);
150			app.affiche((char*)"Maillage");
151			MG_GESTIONNAIRE* gest;
152			MG_FILE* file=new MG_FILE(fichier);
153			if (file->get_code_de_lecture()==0)
154			{
155			delete file;
156			MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
157			gest=gest2;
158			}
159			else gest=file;
160			MG_MAILLAGE* mgmai=new MG_MAILLAGE(NULL);
161			gest->ajouter_mg_maillage(mgmai);
162			MAILLEUR3D mailleur(&carte,mgmai);
163			mailleur.change_coef(coef);
164			mailleur.maille();
165			if (fem)
166			{
167			app.affiche((char*)"Maillage FEM");
168			FEM_MAILLAGE* femmai=new FEM_MAILLAGE(mgmai->get_mg_geometrie(),mgmai,degre);
169			gest->ajouter_fem_maillage(femmai);
170	francois	883	MAILLEUR_FEM mf;
171			mf.maille(femmai,0);
172	francois	432	}
173			app.affiche((char*)"Enregistrement");
174			gest->enregistrer(fichier);
175			delete gest;
176			app.affiche((char*)"Fin");
177
178
179
180			}
181	francois	288
182
183
184			return 0;
185			}
186
187
188
189			#pragma package(smart_init)