ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/carte/src/main.cpp
Revision: 523
Committed: Tue Jun 10 19:02:38 2014 UTC (11 years, 2 months ago) by francois
File size: 16042 byte(s)
Log Message:
integration de la nouvelle carte de taille dans les outils magic V4. Reste un cas d'echantillon eventuellement superposé a faire.

File Contents

# User Rev Content
1 francois 283 //------------------------------------------------------------
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 francois 469 #include "fct_taille_fem_solution_generateur_constant.h"
39 francois 523 #include "fct_taille_fem_solution_generateur_echantillon_frontiere.h"
40     #include "fct_taille_fem_solution_generateur_fichier.h"
41 francois 432 #include "magic_application.h"
42 francois 283 #include <string.h>
43     //---------------------------------------------------------------------------
44    
45    
46    
47    
48     int main(int argc,char **argv)
49     {
50 francois 432 std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
51 francois 523 MAGIC_PARAMETRE_APPLICATION p1(3,(char*)"-fichierpoint",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier d'échantillonage");
52 francois 432 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 francois 469 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 francois 523 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 francois 432 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);
70     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);
71 francois 523 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);
72     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);
73     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);
74     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);
75    
76    
77 francois 432 lst.push_back(p1);
78     lst.push_back(p2);
79     lst.push_back(p3);
80     lst.push_back(p4);
81     lst.push_back(p5);
82     lst.push_back(p6);
83     lst.push_back(p7);
84     lst.push_back(p8);
85     lst.push_back(p9);
86     lst.push_back(p10);
87     lst.push_back(p11);
88     lst.push_back(p12);
89 francois 469 lst.push_back(p13);
90     lst.push_back(p14);
91     lst.push_back(p15);
92     lst.push_back(p16);
93 francois 523 lst.push_back(p19);
94     lst.push_back(p18);
95 francois 432 MAGIC_APPLICATION app((char*)"Generation de carte de taille",argc,argv,lst,true,true);
96     if (app.get_erreur()==true) return 0;
97    
98    
99     if (app.get_action()==12)
100 francois 283 {
101 francois 432 char nomfichier[500];
102     char nomfichiersortie[500];
103     char nomfichierpoint[500];
104     double eng;
105     double eps;
106     double ind1,ind2;
107     int fechan,festim;
108     int numgroupe,numgeo;
109     app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
110     app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
111     app.recupere_parametre_string_avec_defaut(3,nomfichierpoint,(char*)"");
112     app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
113     app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
114     app.recupere_parametre_double_avec_defaut(6,ind1,4.);
115     app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
116     app.recupere_parametre_int_avec_defaut(8,fechan,20);
117     app.recupere_parametre_int_avec_defaut(9,festim,10);
118     app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
119     app.recupere_parametre_int_avec_defaut(14,numgeo,0);
120     if (app.get_erreur()==true) return 0;
121     app.affiche((char*)"Generation de carte de taille");
122     app.affiche((char*)"Lecture du fichier");
123 francois 283 MG_FILE gest(nomfichier);
124     int nbptctl=0;
125     double *ptctl=NULL;
126     char mess[255];
127 francois 432 if (strlen(nomfichierpoint)>0)
128     {
129     FILE* in=fopen(nomfichierpoint,"rt");
130     if (in!=NULL)
131     {
132     app.affiche((char*)"Lecture du fichier de points de controles");
133     char mess[255];
134     char* res=fgets(mess,255,in);
135     sscanf(mess,"%d",&nbptctl);
136     ptctl=new double[4*nbptctl];
137     for (int i=0;i<nbptctl;i++)
138     {
139 francois 283 char* res=fgets(mess,255,in);
140     sscanf(mess,"%lf %lf %lf %lf",ptctl+4*i,ptctl+4*i+1,ptctl+4*i+2,ptctl+4*i+3);
141 francois 432 }
142     fclose(in);
143     }
144     }
145 francois 283 FCT_GENERATEUR_FRONTIERE carte(gest,eng,eps,ind1,ind2,ptctl,nbptctl);
146 francois 432 app.affiche((char*)"Construction de la carte");
147 francois 283 MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
148     if (numgroupe!=0)
149     {
150     MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
151     mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
152     }
153     carte.construit(fechan,festim,mggt);
154 francois 432 app.affiche((char*)"Enregistrement de la carte");
155 francois 283 carte.enregistrer(nomfichiersortie);
156     if (nbptctl!=0) delete [] ptctl;
157 francois 432 app.affiche((char*)"Fin");
158 francois 283 }
159 francois 523 if (app.get_action()==19)
160     {
161     char nomfichier[500];
162     char nomfichiersortie[500];
163     char nomfichierpoint[500];
164     double eng;
165     double eps;
166     double ind1,ind2;
167     int fechan;
168     int numgeo;
169     int nx,ny,nz;
170     app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
171     app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
172     app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
173     app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
174     app.recupere_parametre_double_avec_defaut(6,ind1,4.);
175     app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
176     app.recupere_parametre_int_avec_defaut(8,fechan,20);
177     app.recupere_parametre_int_avec_defaut(14,numgeo,0);
178     app.recupere_parametre_int_avec_defaut(16,nx,10);
179     app.recupere_parametre_int_avec_defaut(17,ny,10);
180     app.recupere_parametre_int_avec_defaut(18,nz,10);
181     if (app.get_erreur()==true) return 0;
182     app.affiche((char*)"Generation de carte de taille");
183     app.affiche((char*)"Lecture du fichier");
184     MG_FILE gest(nomfichier);
185     char mess[255];
186     MG_GESTIONNAIRE gest1;
187     char nomsolution[500];
188     sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
189     app.affiche((char*)"Construction de la carte");
190     MG_GEOMETRIE *geo;
191     if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
192     FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ECHANTILLON_FRONTIERE carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,eps,ind1,ind2);
193     carte.construit();
194     app.affiche((char*)"Enregistrement de la carte");
195     carte.enregistrer(nomfichiersortie);
196     app.affiche((char*)"Fin");
197     }
198 francois 469 if (app.get_action()==10)
199 francois 283 {
200 francois 432 char nomfichier[500];
201     char nomfichiersortie[500];
202     double eng;
203     int fechan,festim;
204     int numgroupe,numgeo;
205     app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
206     app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
207     app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
208     app.recupere_parametre_int_avec_defaut(8,fechan,20);
209     app.recupere_parametre_int_avec_defaut(9,festim,10);
210     app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
211     app.recupere_parametre_int_avec_defaut(14,numgeo,0);
212     if (app.get_erreur()==true) return 0;
213     app.affiche((char*)"Generation de carte de taille");
214     app.affiche((char*)"Lecture du fichier");
215 francois 283 MG_FILE gest(nomfichier);
216     FCT_GENERATEUR_CONSTANTE carte(gest,eng);
217 francois 432 app.affiche((char*)"Construction de la carte");
218 francois 283 MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
219     if (numgroupe!=0)
220     {
221     MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
222     mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
223     }
224     carte.construit(fechan,festim,mggt);
225 francois 432 app.affiche((char*)"Enregistrement de la carte");
226 francois 283 carte.enregistrer(nomfichiersortie);
227 francois 432 app.affiche((char*)"Fin");
228 francois 283 }
229 francois 432 if (app.get_action()==11)
230 francois 283 {
231 francois 432 char nomfichier[500];
232     char nomfichiersortie[500];
233     char nomfichierpoint[500];
234     double eng;
235     int fechan,festim;
236     double ind1,ind2;
237     app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
238     app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
239     app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant");
240     app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
241     app.recupere_parametre_double_avec_defaut(6,ind1,4.);
242     app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
243     app.recupere_parametre_int_avec_defaut(8,fechan,20);
244     app.recupere_parametre_int_avec_defaut(9,festim,10);
245     if (app.get_erreur()==true) return 0;
246     app.affiche((char*)"Generation de carte de taille");
247     app.affiche((char*)" a partir du fichier d'estimation");
248 francois 394 FCT_GENERATEUR_FICHIER carte(nomfichierpoint,eng,ind1,ind2);
249 francois 432 app.affiche((char*)"Construction de la carte");
250 francois 283 carte.construit(fechan,festim);
251 francois 432 app.affiche((char*)"Enregistrement de la carte");
252 francois 283 carte.enregistrer(nomfichiersortie);
253 francois 432 app.affiche((char*)"Fin");
254 francois 283 }
255 francois 523 if (app.get_action()==20)
256     {
257     char nomfichier[500];
258     char nomfichiersortie[500];
259     char nomfichierpoint[500];
260     double eng;
261     int fechan,nx,ny,nz,numgeo;
262     double ind1,ind2;
263     app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
264     app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
265     app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant");
266     app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
267     app.recupere_parametre_double_avec_defaut(6,ind1,4.);
268     app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
269     app.recupere_parametre_int_avec_defaut(8,fechan,20);
270     app.recupere_parametre_int_avec_defaut(14,numgeo,0);
271     app.recupere_parametre_int_avec_defaut(16,nx,10);
272     app.recupere_parametre_int_avec_defaut(17,ny,10);
273     app.recupere_parametre_int_avec_defaut(18,nz,10);
274     if (app.get_erreur()==true) return 0;
275     app.affiche((char*)"Generation de carte de taille");
276     app.affiche((char*)" a partir du fichier d'estimation");
277     app.affiche((char*)"Lecture du fichier");
278     MG_FILE gest(nomfichier);
279     MG_GESTIONNAIRE gest1;
280     char nomsolution[500];
281     sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
282     MG_GEOMETRIE *geo;
283     if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
284     FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,nomfichierpoint,ind1,ind2);
285     app.affiche((char*)"Construction de la carte");
286     carte.construit();
287     app.affiche((char*)"Enregistrement de la carte");
288     carte.enregistrer(nomfichiersortie);
289     app.affiche((char*)"Fin");
290     }
291 francois 469 if (app.get_action()==15)
292     {
293     char nomfichier[500];
294     char nomfichiersortie[500];
295     double eng;
296     int fechan,nx,ny,nz;
297     int numgeo;
298     app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
299     app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
300     app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
301     app.recupere_parametre_int_avec_defaut(8,fechan,20);
302     app.recupere_parametre_int_avec_defaut(14,numgeo,0);
303     app.recupere_parametre_int_avec_defaut(16,nx,10);
304     app.recupere_parametre_int_avec_defaut(17,ny,10);
305     app.recupere_parametre_int_avec_defaut(18,nz,10);
306     if (app.get_erreur()==true) return 0;
307     app.affiche((char*)"Generation de carte de taille");
308     app.affiche((char*)"Lecture du fichier");
309     MG_FILE gest(nomfichier);
310     MG_GEOMETRIE* geo;
311     if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
312     MG_GESTIONNAIRE gest1;
313 francois 523 char nomsolution[500];
314     sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
315 francois 469 app.affiche((char*)"Construction de la carte");
316 francois 523 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution);
317     carte.construit();
318 francois 469 app.affiche((char*)"Enregistrement de la carte");
319     carte.enregistrer(nomfichiersortie);
320     app.affiche((char*)"Fin");
321     }
322 francois 283
323     return 0;
324     }
325    
326    
327    
328     #pragma package(smart_init)