ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/carte/src/main.cpp
Revision: 1053
Committed: Fri Oct 16 16:17:23 2020 UTC (4 years, 10 months ago) by francois
File size: 18470 byte(s)
Log Message:
mailleur2d_structure de quadrangle

File Contents

# Content
1 //------------------------------------------------------------
2 //------------------------------------------------------------
3 // MAGiC
4 // Jean Christophe Cuilli�re et Vincent FRANCOIS
5 // D�partement de G�nie M�canique - UQTR
6 //------------------------------------------------------------
7 // Le projet MAGIC est un projet de recherche du d�partement
8 // de g�nie m�canique de l'Universit� du Qu�bec �
9 // Trois Rivi�res
10 // Les librairies ne peuvent �tre utilis�es 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
24
25
26
27 #include "gestionversion.h"
28 #ifdef WINDOWS_VERSION
29 #include "fenetre.h"
30 #endif
31
32 #pragma hdrstop
33 #include "mg_file.h"
34 #include "mg_export.h"
35 #include "fct_generateur_frontiere.h"
36 #include "fct_generateur_constante.h"
37 #include "fct_generateur_fichier.h"
38 #include "fct_taille_fem_solution_generateur_constant.h"
39 #include "fct_taille_fem_solution_generateur_echantillon_frontiere.h"
40 #include "fct_taille_fem_solution_generateur_fichier.h"
41 #include "magic_application.h"
42 #include <string.h>
43 //---------------------------------------------------------------------------
44
45
46
47
48 int main(int argc,char **argv)
49 {
50 std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
51 MAGIC_PARAMETRE_APPLICATION p1(3,(char*)"-fichierpoint",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier d'échantillonage");
52 MAGIC_PARAMETRE_APPLICATION p2(4,(char*)"-dg",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Densité globale");
53 MAGIC_PARAMETRE_APPLICATION p3(5,(char*)"-eps",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Écart entre le maillage et la géométrie");
54 MAGIC_PARAMETRE_APPLICATION p4(6,(char*)"-taillezone",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Facteur pour la distance de variation de la carte autour de l'échantillon - Defaut 4");
55 MAGIC_PARAMETRE_APPLICATION p5(7,(char*)"-tauxtransfert",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Poid de l'échantillon virtuel place en dehors de la zone de raffinement - Defaut 0.25");
56 MAGIC_PARAMETRE_APPLICATION p6(8,(char*)"-fechantillonnage",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'échantillon par entité - Defaut 20");
57 MAGIC_PARAMETRE_APPLICATION p7(9,(char*)"-festimation",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre de cellule dans une direction pour la grille virtuelle - Defaut 10");
58 MAGIC_PARAMETRE_APPLICATION p8(10,(char*)"-constante",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille constante");
59 MAGIC_PARAMETRE_APPLICATION p9(11,(char*)"-fichierechan",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par un echantillonnage dans un fichier exterieur");
60 MAGIC_PARAMETRE_APPLICATION p10(12,(char*)"-courbure",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par approximation des courbures");
61 MAGIC_PARAMETRE_APPLICATION p11(13,(char*)"-groupe",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Groupe topologique utilise pour faire le calcul - Defaut 0");
62 MAGIC_PARAMETRE_APPLICATION p12(14,(char*)"-numgeo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Géométrie du groupe topologique - Defaut 0");
63 MAGIC_PARAMETRE_APPLICATION p13(15,(char*)"-constante2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille constante définie par une FEM_SOLUTION");
64 MAGIC_PARAMETRE_APPLICATION p14(16,(char*)"-nx",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'intervals sur x - Defaut 10");
65 MAGIC_PARAMETRE_APPLICATION p15(17,(char*)"-ny",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'intervals sur y - Defaut 10");
66 MAGIC_PARAMETRE_APPLICATION p16(18,(char*)"-nz",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'intervals sur z - Defaut 10");
67 MAGIC_PARAMETRE_APPLICATION p19(19,(char*)"-courbure2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par approximation des courbures définie par une FEM_SOLUTION");
68 MAGIC_PARAMETRE_APPLICATION p18(20,(char*)"-fichierechan2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par un echantillonnage dans un fichier exterieur définie par une FEM_SOLUTION");
69 MAGIC_PARAMETRE_APPLICATION p20(21,(char*)"-constante2mai",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille constante définie par une FEM_SOLUTION basé sur maillage");
70 MAGIC_PARAMETRE_APPLICATION p21(22,(char*)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de maillage - Defaut le premier");
71 p8.ajouter_dependance(1);p8.ajouter_dependance(2);p8.ajouter_dependance(4);p8.ajouter_dependance(8);p8.ajouter_dependance(9);p8.ajouter_dependance(13);p8.ajouter_dependance(14);
72 p9.ajouter_dependance(1);p9.ajouter_dependance(2);p9.ajouter_dependance(3);p9.ajouter_dependance(4);p9.ajouter_dependance(6);p9.ajouter_dependance(7);p9.ajouter_dependance(8);p9.ajouter_dependance(9);
73 p10.ajouter_dependance(1);p10.ajouter_dependance(2);p10.ajouter_dependance(3);p10.ajouter_dependance(4);p10.ajouter_dependance(5);p10.ajouter_dependance(6);p10.ajouter_dependance(7);p10.ajouter_dependance(8);p10.ajouter_dependance(9);p10.ajouter_dependance(13);p10.ajouter_dependance(14);
74 p13.ajouter_dependance(1);p13.ajouter_dependance(2);p13.ajouter_dependance(4);p13.ajouter_dependance(8);p13.ajouter_dependance(14);p13.ajouter_dependance(16);p13.ajouter_dependance(17);p13.ajouter_dependance(18);
75 p18.ajouter_dependance(1);p18.ajouter_dependance(2);p18.ajouter_dependance(3);p18.ajouter_dependance(4);p18.ajouter_dependance(6);p18.ajouter_dependance(7);p18.ajouter_dependance(8);p18.ajouter_dependance(14);p18.ajouter_dependance(16);p18.ajouter_dependance(17);p18.ajouter_dependance(18);
76 p19.ajouter_dependance(1);p19.ajouter_dependance(2);p19.ajouter_dependance(4);p19.ajouter_dependance(5);p19.ajouter_dependance(6);p19.ajouter_dependance(7);p19.ajouter_dependance(8);p19.ajouter_dependance(14);p19.ajouter_dependance(16);p19.ajouter_dependance(17);p19.ajouter_dependance(18);
77 p20.ajouter_dependance(1);p20.ajouter_dependance(2);p20.ajouter_dependance(4);p20.ajouter_dependance(8);p20.ajouter_dependance(22);p20.ajouter_dependance(16);p20.ajouter_dependance(17);p20.ajouter_dependance(18);
78 p19.ajouter_dependance(5);
79 lst.push_back(p1);
80 lst.push_back(p2);
81 lst.push_back(p3);
82 lst.push_back(p4);
83 lst.push_back(p5);
84 lst.push_back(p6);
85 lst.push_back(p7);
86 lst.push_back(p8);
87 lst.push_back(p9);
88 lst.push_back(p10);
89 lst.push_back(p11);
90 lst.push_back(p12);
91 lst.push_back(p13);
92 lst.push_back(p14);
93 lst.push_back(p15);
94 lst.push_back(p16);
95 lst.push_back(p19);
96 lst.push_back(p18);
97 lst.push_back(p20);
98 lst.push_back(p21);
99 MAGIC_APPLICATION app((char*)"Generation de carte de taille",argc,argv,lst,true,true);
100 if (app.get_erreur()==true) return 0;
101
102
103 if (app.get_action()==12)
104 {
105 char nomfichier[500];
106 char nomfichiersortie[500];
107 char nomfichierpoint[500];
108 double eng;
109 double eps;
110 double ind1,ind2;
111 int fechan,festim;
112 int numgroupe,numgeo;
113 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
114 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
115 app.recupere_parametre_string_avec_defaut(3,nomfichierpoint,(char*)"");
116 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
117 app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
118 app.recupere_parametre_double_avec_defaut(6,ind1,4.);
119 app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
120 app.recupere_parametre_int_avec_defaut(8,fechan,20);
121 app.recupere_parametre_int_avec_defaut(9,festim,10);
122 app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
123 app.recupere_parametre_int_avec_defaut(14,numgeo,0);
124 if (app.get_erreur()==true) return 0;
125 app.affiche((char*)"Generation de carte de taille basee sur la courbure des frontieres V1");
126 app.affiche((char*)"Lecture du fichier");
127 MG_FILE gest(nomfichier);
128 int nbptctl=0;
129 double *ptctl=NULL;
130 char mess[255];
131 if (strlen(nomfichierpoint)>0)
132 {
133 FILE* in=fopen(nomfichierpoint,"rt");
134 if (in!=NULL)
135 {
136 app.affiche((char*)"Lecture du fichier de points de controles");
137 char mess[255];
138 char* res=fgets(mess,255,in);
139 sscanf(mess,"%d",&nbptctl);
140 ptctl=new double[4*nbptctl];
141 for (int i=0;i<nbptctl;i++)
142 {
143 char* res=fgets(mess,255,in);
144 sscanf(mess,"%lf %lf %lf %lf",ptctl+4*i,ptctl+4*i+1,ptctl+4*i+2,ptctl+4*i+3);
145 }
146 fclose(in);
147 }
148 }
149 FCT_GENERATEUR_FRONTIERE carte(gest,eng,eps,ind1,ind2,ptctl,nbptctl);
150 app.affiche((char*)"Construction de la carte");
151 MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
152 if (numgroupe!=0)
153 {
154 MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
155 mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
156 }
157 carte.construit(fechan,festim,mggt);
158 app.affiche((char*)"Enregistrement de la carte");
159 carte.enregistrer(nomfichiersortie);
160 if (nbptctl!=0) delete [] ptctl;
161 app.affiche((char*)"Fin");
162 }
163 if (app.get_action()==19)
164 {
165 char nomfichier[500];
166 char nomfichiersortie[500];
167 char nomfichierpoint[500];
168 double eng;
169 double eps;
170 double ind1,ind2;
171 int fechan;
172 int numgeo;
173 int nx,ny,nz;
174 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
175 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
176 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
177 app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
178 app.recupere_parametre_double_avec_defaut(6,ind1,4.);
179 app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
180 app.recupere_parametre_int_avec_defaut(8,fechan,20);
181 app.recupere_parametre_int_avec_defaut(14,numgeo,0);
182 app.recupere_parametre_int_avec_defaut(16,nx,10);
183 app.recupere_parametre_int_avec_defaut(17,ny,10);
184 app.recupere_parametre_int_avec_defaut(18,nz,10);
185 if (app.get_erreur()==true) return 0;
186 app.affiche((char*)"Generation de carte de taille basee sur la courbure des frontieres V2");
187 app.affiche((char*)"Lecture du fichier");
188 MG_FILE gest(nomfichier);
189 char mess[255];
190 MG_GESTIONNAIRE gest1;
191 char nomsolution[500];
192 sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
193 app.affiche((char*)"Construction de la carte");
194 MG_GEOMETRIE *geo;
195 if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
196 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ECHANTILLON_FRONTIERE carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,eps,ind1,ind2);
197 carte.construit();
198 app.affiche((char*)"Enregistrement de la carte");
199 carte.enregistrer(nomfichiersortie);
200 app.affiche((char*)"Fin");
201 }
202 if (app.get_action()==10)
203 {
204 char nomfichier[500];
205 char nomfichiersortie[500];
206 double eng;
207 int fechan,festim;
208 int numgroupe,numgeo;
209 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
210 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
211 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
212 app.recupere_parametre_int_avec_defaut(8,fechan,20);
213 app.recupere_parametre_int_avec_defaut(9,festim,10);
214 app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
215 app.recupere_parametre_int_avec_defaut(14,numgeo,0);
216 if (app.get_erreur()==true) return 0;
217 app.affiche((char*)"Generation de carte de taille constante V1");
218 app.affiche((char*)"Lecture du fichier");
219 MG_FILE gest(nomfichier);
220 FCT_GENERATEUR_CONSTANTE carte(gest,eng);
221 app.affiche((char*)"Construction de la carte");
222 MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
223 if (numgroupe!=0)
224 {
225 MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
226 mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
227 }
228 carte.construit(fechan,festim,mggt);
229 app.affiche((char*)"Enregistrement de la carte");
230 carte.enregistrer(nomfichiersortie);
231 app.affiche((char*)"Fin");
232 }
233 if (app.get_action()==11)
234 {
235 char nomfichier[500];
236 char nomfichiersortie[500];
237 char nomfichierpoint[500];
238 double eng;
239 int fechan,festim;
240 double ind1,ind2;
241 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
242 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
243 app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
244 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
245 app.recupere_parametre_double_avec_defaut(6,ind1,4.);
246 app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
247 app.recupere_parametre_int_avec_defaut(8,fechan,20);
248 app.recupere_parametre_int_avec_defaut(9,festim,10);
249 if (app.get_erreur()==true) return 0;
250 app.affiche((char*)"Generation de carte de taille");
251 app.affiche((char*)" a partir du fichier d'estimation V1");
252 FCT_GENERATEUR_FICHIER carte(nomfichierpoint,eng,ind1,ind2);
253 app.affiche((char*)"Construction de la carte");
254 carte.construit(fechan,festim);
255 app.affiche((char*)"Enregistrement de la carte");
256 carte.enregistrer(nomfichiersortie);
257 app.affiche((char*)"Fin");
258 }
259 if (app.get_action()==20)
260 {
261 char nomfichier[500];
262 char nomfichiersortie[500];
263 char nomfichierpoint[500];
264 double eng;
265 int fechan,nx,ny,nz,numgeo;
266 double ind1,ind2;
267 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
268 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
269 app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
270 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
271 app.recupere_parametre_double_avec_defaut(6,ind1,4.);
272 app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
273 app.recupere_parametre_int_avec_defaut(8,fechan,20);
274 app.recupere_parametre_int_avec_defaut(14,numgeo,0);
275 app.recupere_parametre_int_avec_defaut(16,nx,10);
276 app.recupere_parametre_int_avec_defaut(17,ny,10);
277 app.recupere_parametre_int_avec_defaut(18,nz,10);
278 if (app.get_erreur()==true) return 0;
279 app.affiche((char*)"Generation de carte de taille");
280 app.affiche((char*)" a partir du fichier d'estimation V2");
281 app.affiche((char*)"Lecture du fichier");
282 MG_FILE gest(nomfichier);
283 MG_GESTIONNAIRE gest1;
284 char nomsolution[500];
285 sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
286 MG_GEOMETRIE *geo;
287 if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
288 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,nomfichierpoint,ind1,ind2);
289 app.affiche((char*)"Construction de la carte");
290 carte.construit();
291 app.affiche((char*)"Enregistrement de la carte");
292 carte.enregistrer(nomfichiersortie);
293 app.affiche((char*)"Fin");
294 }
295 if (app.get_action()==15)
296 {
297 char nomfichier[500];
298 char nomfichiersortie[500];
299 double eng;
300 int fechan,nx,ny,nz;
301 int numgeo;
302 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
303 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
304 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
305 app.recupere_parametre_int_avec_defaut(8,fechan,20);
306 app.recupere_parametre_int_avec_defaut(14,numgeo,0);
307 app.recupere_parametre_int_avec_defaut(16,nx,10);
308 app.recupere_parametre_int_avec_defaut(17,ny,10);
309 app.recupere_parametre_int_avec_defaut(18,nz,10);
310 if (app.get_erreur()==true) return 0;
311 app.affiche((char*)"Generation de carte de taille constante V2");
312 app.affiche((char*)"Lecture du fichier");
313 MG_FILE gest(nomfichier);
314 MG_GEOMETRIE* geo;
315 if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
316 MG_GESTIONNAIRE gest1;
317 char nomsolution[500];
318 sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
319 app.affiche((char*)"Construction de la carte");
320 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution);
321 carte.construit();
322 app.affiche((char*)"Enregistrement de la carte");
323 carte.enregistrer(nomfichiersortie);
324 app.affiche((char*)"Fin");
325 }
326 if (app.get_action()==21)
327 {
328 char nomfichier[500];
329 char nomfichiersortie[500];
330 double eng;
331 int fechan,nx,ny,nz;
332 int nummai;
333 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
334 app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
335 app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
336 app.recupere_parametre_int_avec_defaut(8,fechan,20);
337 app.recupere_parametre_int_avec_defaut(22,nummai,0);
338 app.recupere_parametre_int_avec_defaut(16,nx,10);
339 app.recupere_parametre_int_avec_defaut(17,ny,10);
340 app.recupere_parametre_int_avec_defaut(18,nz,10);
341 if (app.get_erreur()==true) return 0;
342 app.affiche((char*)"Generation de carte de taille constante V2");
343 app.affiche((char*)"Lecture du fichier");
344 MG_FILE gest(nomfichier);
345 MG_MAILLAGE* mai;
346 if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
347 MG_GESTIONNAIRE gest1;
348 char nomsolution[500];
349 sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
350 app.affiche((char*)"Construction de la carte");
351 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,mai,eng,fechan,nx,ny,nz,nomsolution);
352 carte.construit();
353 app.affiche((char*)"Enregistrement de la carte");
354 carte.enregistrer(nomfichiersortie);
355 app.affiche((char*)"Fin");
356 }
357
358 return 0;
359 }
360
361
362
363 #pragma package(smart_init)