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 "fct_taille_fem_solution.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 avec ancienne carte de taille");
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");
MAGIC_PARAMETRE_APPLICATION p11(12,(char*)"-delaunay2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Maillage d'une boite par la méthode de Delaunay avec nouvelle carte de taille");
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);
p11.ajouter_dependance(1);p11.ajouter_dependance(2);p11.ajouter_dependance(3);p11.ajouter_dependance(4),p11.ajouter_dependance(5);p11.ajouter_dependance(11);
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);
lst.push_back(p11);

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) || (app.get_action()==12))
  {
  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(5,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_TAILLE* metrique;
  if (app.get_action()==6)
          {
          FCT_GENERATEUR_3D<4> *carte=new FCT_GENERATEUR_3D<4>;
          carte->lire(fichiercarte);
          metrique=carte;
          }
   //MG_FILE *gestcarte;
   if (app.get_action()==12)
          {
          FCT_TAILLE_FEM_SOLUTION *carte=new FCT_TAILLE_FEM_SOLUTION(fichiercarte);
          metrique=carte;
          }
  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,metrique,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,metrique);
    strcat(fich,".magic");
    gest->enregistrer(fich);
    }
  delete metrique;
  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:	925
Committed:	Wed Apr 25 21:39:05 2018 UTC (7 years, 1 month ago) by francois
File size:	8027 byte(s)
Log Message:	bug dans le mailleur de delaunay resultant du changement de carte taille (bug deja ancien)
#	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	925	#include "fct_taille_fem_solution.h"
38	francois	288	#include "mailleur_delaunay.h"
39			#include <fstream>
40			#include <string.h>
41			#include "mg_maillage_algos.h"
42	francois	432	#include "magic_application.h"
43	francois	883	#include "mailleur_fem.h"
44	francois	432
45
46	francois	288	//---------------------------------------------------------------------------
47	francois	432	int main(int argc,char **argv)
48	francois	288
49			{
50	francois	432	std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
51			MAGIC_PARAMETRE_APPLICATION p1(2,(char)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char)"Fichier de carte de taille");
52			MAGIC_PARAMETRE_APPLICATION p2(3,(char)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Limite supérieure des mailles à optimiser - Défaut 2");
53			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");
54			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");
55	francois	925	MAGIC_PARAMETRE_APPLICATION p5(6,(char)"-delaunay",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode de Delaunay avec ancienne carte de taille");
56	francois	432	MAGIC_PARAMETRE_APPLICATION p6(7,(char)"-octal",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode octal");
57			MAGIC_PARAMETRE_APPLICATION p7(8,(char)"-tailleoctree",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Taille des octants par rapport à la taille souhaitée - Défaut 2.");
58			MAGIC_PARAMETRE_APPLICATION p8(9,(char)"-fem",MAGIC_PARAMETRE_APPLICATION::BOOL,(char)"Création du maillage fem");
59			MAGIC_PARAMETRE_APPLICATION p9(10,(char)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Degre du maillage fem - Défaut 1");
60	francois	532	MAGIC_PARAMETRE_APPLICATION p10(11,(char)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Valeur du pas d'integration - Defaut 5");
61	francois	925	MAGIC_PARAMETRE_APPLICATION p11(12,(char)"-delaunay2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Maillage d'une boite par la méthode de Delaunay avec nouvelle carte de taille");
62	francois	532	p5.ajouter_dependance(1);p5.ajouter_dependance(2);p5.ajouter_dependance(3);p5.ajouter_dependance(4),p5.ajouter_dependance(5);p5.ajouter_dependance(11);
63	francois	432	p6.ajouter_dependance(1);p6.ajouter_dependance(2);p6.ajouter_dependance(8);p6.ajouter_dependance(9);p6.ajouter_dependance(10);
64	francois	925	p11.ajouter_dependance(1);p11.ajouter_dependance(2);p11.ajouter_dependance(3);p11.ajouter_dependance(4),p11.ajouter_dependance(5);p11.ajouter_dependance(11);
65	francois	432	lst.push_back(p1);
66			lst.push_back(p2);
67			lst.push_back(p3);
68			lst.push_back(p4);
69			lst.push_back(p5);
70			lst.push_back(p6);
71			lst.push_back(p7);
72			lst.push_back(p8);
73			lst.push_back(p9);
74	francois	532	lst.push_back(p10);
75	francois	925	lst.push_back(p11);
76	francois	288
77	francois	485	MAGIC_APPLICATION app((char*)"Mailleur d'un espace défini par une carte de taille",argc,argv,lst,false,true);
78	francois	432	if (app.get_erreur()==true) return 0;
79
80	francois	288
81
82	francois	925	if ((app.get_action()==6) \|\| (app.get_action()==12))
83	francois	432	{
84			char fichier[500];
85			char fichiercarte[500];
86			double critere;
87			bool analyse;
88			char message[500];
89			int niveau;
90	francois	532	int pas;
91	francois	432	app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
92			app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
93			app.recupere_parametre_int_avec_defaut(3,niveau,2);
94			app.recupere_parametre_double_avec_defaut(4,critere,1.5);
95	francois	925	app.recupere_parametre_bool_avec_defaut(5,analyse,false);
96	francois	532	app.recupere_parametre_int_avec_defaut(11,pas,5);
97	francois	432	if (app.get_erreur()==true) return 0;
98			app.affiche((char*)"Méthode de Delaunay");
99			sprintf(message,"Fichier carte : %s",fichiercarte);
100			app.affiche(message);
101			app.affiche((char*)"");
102			app.affiche((char*)"Lecture carte de taille");
103			OT_CPU compteur;
104	francois	925	FCT_TAILLE* metrique;
105			if (app.get_action()==6)
106			{
107			FCT_GENERATEUR_3D<4> *carte=new FCT_GENERATEUR_3D<4>;
108			carte->lire(fichiercarte);
109			metrique=carte;
110			}
111			//MG_FILE *gestcarte;
112			if (app.get_action()==12)
113			{
114			FCT_TAILLE_FEM_SOLUTION *carte=new FCT_TAILLE_FEM_SOLUTION(fichiercarte);
115			metrique=carte;
116			}
117	francois	432	MG_GESTIONNAIRE* gest;
118			MG_FILE* file=new MG_FILE(fichier);
119			if (file->get_code_de_lecture()==0)
120			{
121			delete file;
122			MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
123			gest=gest2;
124			}
125			else gest=file;
126			MG_MAILLAGE* mai=new MG_MAILLAGE(NULL);
127			gest->ajouter_mg_maillage(mai);
128	francois	925	MAILLEUR_DELAUNAY m(mai,metrique,critere,&compteur);
129	francois	432	m.active_affichage(app.affiche);
130			m.change_niveau_optimisation(niveau);
131	francois	532	m.change_pas_integration(pas);
132	francois	432	m.maille();
133			char fich[500];
134			strcpy(fich,fichier);
135			gest->enregistrer(fich);
136	francois	925	if (analyse)
137			{
138			MG_MAILLAGE_ALGOS::calcule_caracteristiques(mai,fich,metrique);
139			strcat(fich,".magic");
140			gest->enregistrer(fich);
141			}
142			delete metrique;
143	francois	432	compteur.ajouter_etape("Enregistrement");
144			std::vector<std::string> tab;
145			compteur.get_tabfinal(tab);
146			for (int i=0;i<tab.size();i++)
147			app.affiche((char*)tab[i].c_str());
148			}
149			if (app.get_action()==7)
150			{
151			char fichier[500];
152			char fichiercarte[500];
153			double coef;
154			bool fem;
155			int degre;
156			app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
157			app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant");
158			app.recupere_parametre_double_avec_defaut(8,coef,2.);
159			app.recupere_parametre_bool_avec_defaut(9,fem,false);
160			app.recupere_parametre_int_avec_defaut(10,degre,1);
161			if (app.get_erreur()==true) return 0;
162			char message[500];
163			app.affiche((char*)"Méthode octale");
164			sprintf(message," Fichier carte : %s",fichiercarte);
165			app.affiche(message);
166			app.affiche((char*)"Lecture carte de taille");
167			FCT_GENERATEUR_3D<4> carte;
168			carte.lire(fichiercarte);
169			app.affiche((char*)"Maillage");
170			MG_GESTIONNAIRE* gest;
171			MG_FILE* file=new MG_FILE(fichier);
172			if (file->get_code_de_lecture()==0)
173			{
174			delete file;
175			MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
176			gest=gest2;
177			}
178			else gest=file;
179			MG_MAILLAGE* mgmai=new MG_MAILLAGE(NULL);
180			gest->ajouter_mg_maillage(mgmai);
181			MAILLEUR3D mailleur(&carte,mgmai);
182			mailleur.change_coef(coef);
183			mailleur.maille();
184			if (fem)
185			{
186			app.affiche((char*)"Maillage FEM");
187			FEM_MAILLAGE* femmai=new FEM_MAILLAGE(mgmai->get_mg_geometrie(),mgmai,degre);
188			gest->ajouter_fem_maillage(femmai);
189	francois	883	MAILLEUR_FEM mf;
190			mf.maille(femmai,0);
191	francois	432	}
192			app.affiche((char*)"Enregistrement");
193			gest->enregistrer(fichier);
194			delete gest;
195			app.affiche((char*)"Fin");
196
197
198
199			}
200	francois	288
201
202
203			return 0;
204			}
205
206
207
208			#pragma package(smart_init)