mailleur/src/main2.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       main2.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2025
//####//  Derniere modification par francois
//####//  mer 14 mai 2025 17:54:47 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------


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