ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/mgoperation/src/main.cpp
Revision: 318
Committed: Fri Feb 24 17:52:15 2012 UTC (13 years, 3 months ago) by gervaislavoie
File size: 27220 byte(s)
Log Message:
Correction d'un bug sur la saisie des ccf

File Contents

# User Rev Content
1 francois 283
2     //------------------------------------------------------------
3     //------------------------------------------------------------
4     // MAGiC
5     // Jean Christophe Cuilli�e et Vincent FRANCOIS
6     // D�artement de G�ie M�anique - UQTR
7     //------------------------------------------------------------
8     // Le projet MAGIC est un projet de recherche du d�artement
9     // de g�ie m�anique de l'Universit�du Qu�ec �
10     // Trois Rivi�es
11     // Les librairies ne peuvent �re utilis�s sans l'accord
12     // des auteurs (contact : francois@uqtr.ca)
13     //------------------------------------------------------------
14     //------------------------------------------------------------
15     //
16     // main.cpp
17     //
18     //------------------------------------------------------------
19     //------------------------------------------------------------
20     // COPYRIGHT 2000
21     // Version du 02/03/2006 �11H25
22     //------------------------------------------------------------
23     //------------------------------------------------------------
24     #include "gestionversion.h"
25     #ifdef WINDOWS_VERSION
26     #include "fenetre.h"
27     #endif
28    
29    
30     #include <string.h>
31     #include "main.h"
32     #include "mg_file.h"
33     #include "mg_export.h"
34     #include "mg_import.h"
35     #include "mg_lissage.h"
36     #include "step_import.h"
37     #include "acis_import.h"
38     #include "step_import.h"
39     #include "occ_import.h"
40     #include "aide.h"
41 francois 296 #include "mc_gestionnaire.h"
42 francois 283 #include "vct_multi_modele.h"
43     #include <ctype.h>
44     #include <math.h>
45    
46    
47     //---------------------------------------------------------------------------
48    
49     void affiche(char* message)
50     {
51     #ifdef WINDOWS_VERSION
52     Form1->Memo1->Lines->Add(message);
53     #else
54     std::cout << message << std::endl;
55     #endif
56     }
57    
58     void minuscule(char * mess)
59     {
60     int nb=strlen(mess);
61     for (int i=0;i<nb;i++)
62     mess[i]=tolower(mess[i]);
63     }
64    
65     #ifdef WINDOWS_VERSION
66     int amain(int argc,char **argv)
67     #else
68     int main(int argc,char **argv)
69     #endif
70     {
71     affiche((char*)"Operation dans l'environnement MAGiC\n\n");
72     char nomfichier[3000];
73     char fichierout[3000];
74     char fichieraster[3000];
75     char fichieroptis[3000];
76     char fichierflottant[3000];
77     char fichierdensite[3000];
78     char coderesu[10];
79     int nummai;
80     int degre;
81     int action=0;
82     int numgeo=0;
83     int numgt=0;
84     int numentite=0;
85     int expand=0;
86     int numsol1=0;
87     int numsol2=0;
88     int numsol3=0;
89     int numchamp1=0;
90     int numchamp2=0;
91     int numchamp3=0;
92     int flottant=0;
93 gervaislavoie 304 int desactivedeforme=0;
94 francois 283 int consimpose=0;
95     int consmailleurauto=0;
96     int consopti=0;
97     double unite=1;
98     double seuil=0.8;
99     double eps=1e-6;
100     double eps2=1.;
101     double niveau=150.;
102     int iter_max = 10;
103     double sigma;
104     double gamma=1.; //par d�faut 1.0
105     double sigmaf=1.; //par d�faut 1.0
106     double sigmag=1.; //par d�faut 1.0
107     int reactivation = 0;
108     int bruitage = 0;
109     int lissage = 0;
110     int importtri=0;
111     sigma = 0.038/sqrt(10.); //valeur par d�faut
112 francois 292 char ccf[500];
113 francois 283 INIAIDE;
114     if (argc==2) strcpy(nomfichier,argv[1]);
115     else if (argc>2)
116     {
117     for (int i=0;i<argc;i++)
118     {
119     if (strcmp(argv[i],"-file")==0) strcpy(nomfichier,argv[i+1]);
120     if (strcmp(argv[i],"-fileout")==0) strcpy(fichierout,argv[i+1]);
121     if (strcmp(argv[i],"-mesh")==0) nummai=atol(argv[i+1]);
122     if (strcmp(argv[i],"-numsol1")==0) numsol1=atol(argv[i+1]);
123     if (strcmp(argv[i],"-numsol2")==0) numsol2=atol(argv[i+1]);
124     if (strcmp(argv[i],"-numsol3")==0) numsol3=atol(argv[i+1]);
125     if (strcmp(argv[i],"-numchamp1")==0) numchamp1=atol(argv[i+1]);
126     if (strcmp(argv[i],"-numchamp2")==0) numchamp2=atol(argv[i+1]);
127     if (strcmp(argv[i],"-numchamp3")==0) numchamp3=atol(argv[i+1]);
128     if (strcmp(argv[i],"-geo")==0) numgeo=atol(argv[i+1]);
129     if (strcmp(argv[i],"-groupe")==0) numgt=atol(argv[i+1]);
130     if (strcmp(argv[i],"-entite")==0) numentite=atol(argv[i+1]);
131     if (strcmp(argv[i],"-degre")==0) degre=atoi(argv[i+1]);
132     if (strcmp(argv[i],"-unite")==0) unite=atof(argv[i+1]);
133     if (strcmp(argv[i],"-seuil")==0) seuil=atof(argv[i+1]);
134     if (strcmp(argv[i],"-niveau")==0) niveau=atof(argv[i+1]);
135     if (strcmp(argv[i],"-aster")==0) strcpy(fichieraster,argv[i+1]);
136     if (strcmp(argv[i],"-densite")==0) strcpy(fichierdensite,argv[i+1]);
137     if (strcmp(argv[i],"-optis")==0) strcpy(fichieroptis,argv[i+1]);
138     if (strcmp(argv[i],"-ficflottant")==0) strcpy(fichierflottant,argv[i+1]);
139     if (strcmp(argv[i],"-coderesu")==0) strcpy(coderesu,argv[i+1]);
140     if (strcmp(argv[i],"-flottant")==0) flottant = 1;
141     if (strcmp(argv[i],"-eps")==0) eps=atof(argv[i+1]);
142     if (strcmp(argv[i],"-iter")==0) iter_max=atol(argv[i+1]);
143     if (strcmp(argv[i],"-sigma")==0) sigma=atof(argv[i+1]);
144     if (strcmp(argv[i],"-gamma")==0) gamma=atof(argv[i+1]);
145     if (strcmp(argv[i],"-sigmaf")==0) sigmaf=atof(argv[i+1]);
146     if (strcmp(argv[i],"-sigmag")==0) sigmag=atof(argv[i+1]);
147     if (strcmp(argv[i],"-fem")==0) action=1;
148     if (strcmp(argv[i],"-femdef")==0) action=15;
149     if (strcmp(argv[i],"-expand")==0) expand=1;
150     if (strcmp(argv[i],"-stepocc")==0) action=2;
151     if (strcmp(argv[i],"-magicgmsh")==0) action=3;
152     if (strcmp(argv[i],"-magicaster")==0) action=4;
153     if (strcmp(argv[i],"-stepmagic")==0) action=5;
154     if (strcmp(argv[i],"-satmagic")==0) action=6;
155     if (strcmp(argv[i],"-cttgmsh")==0) action=7;
156     if (strcmp(argv[i],"-occmagic")==0) action=8;
157     if (strcmp(argv[i],"-astermagic")==0) action=9;
158     if (strcmp(argv[i],"-magicoptis")==0) action=10;
159     if (strcmp(argv[i],"-optismagic")==0) action=11;
160     if (strcmp(argv[i],"-selectionne")==0) action=12;
161     if (strcmp(argv[i],"-deselectionne")==0) action=13;
162     if (strcmp(argv[i],"-examine")==0) action=14;
163     if (strcmp(argv[i],"-fusiongeo")==0) action=17;
164     if (strcmp(argv[i],"-importdensite")==0) action=18;
165     if (strcmp(argv[i],"-lissage")==0) action=16;
166     if (strcmp(argv[i],"-liss1")==0) lissage = 1;
167     if (strcmp(argv[i],"-liss2")==0) lissage = 2;
168     if (strcmp(argv[i],"-liss3")==0) lissage = 3;
169     if (strcmp(argv[i],"-bruitage")==0) bruitage = 1;
170     if (strcmp(argv[i],"-reactivation")==0) reactivation = 1;
171     if (strcmp(argv[i],"-cons_optimise")==0) consopti = 1;
172     if (strcmp(argv[i],"-cons_impose")==0) consimpose = 1;
173     if (strcmp(argv[i],"-cons_mailleurauto")==0) consmailleurauto = 1;
174     if (strcmp(argv[i],"-importtriangulation")==0) importtri=1;
175     if (strcmp(argv[i],"-importstl")==0) importtri=2;
176     if (strcmp(argv[i],"-deflexion")==0) eps2=atof(argv[i+1]);
177 francois 293 if (strcmp(argv[i],"-ccf_id")==0) {action=19;strcpy(ccf,argv[i+1]);};
178     if (strcmp(argv[i],"-ccft_id")==0) {action=20;strcpy(ccf,argv[i+1]);};
179     if (strcmp(argv[i],"-ccf_num")==0) {action=21;strcpy(ccf,argv[i+1]);};
180     if (strcmp(argv[i],"-ccft_num")==0) {action=22;strcpy(ccf,argv[i+1]);};
181 francois 296 if (strcmp(argv[i],"-ccfhelp")==0) action=23;
182 gervaislavoie 304 if (strcmp(argv[i],"-nodeforme")==0) desactivedeforme=1;
183 francois 292
184 francois 283 }
185     }
186     else {
187     AFFAIDE;
188     return 0;
189     }
190 francois 296 char mess[3000];
191     char *p;
192     char extension[50];
193     char chaine[3000];
194    
195     if (strcmp(nomfichier,"-ccfhelp")==0) action=23;
196     else
197     {
198 francois 283
199     sprintf(mess," Fichier : %s",nomfichier);
200     affiche(mess);
201 francois 296 p=strrchr(nomfichier,'.');
202 francois 283 strcpy(extension,p);
203     minuscule(extension);
204     for (int i=0;i<3000;i++) chaine[i]=0;
205 francois 296 }
206 francois 283 if ( ((strcmp(extension,".magic")==0) && (action==0)) || (action==3))
207     {
208     sprintf(mess," Conversion MAGiC vers GMSH");
209     affiche(mess);
210     MG_FILE gest(nomfichier);
211     int nb=gest.get_nb_mg_maillage();
212     char chaine[3000];
213     for (int i=0;i<nb;i++)
214     {
215     MG_MAILLAGE* mai=gest.get_mg_maillage(i);
216     unsigned long id=mai->get_id();
217     strncpy(chaine,nomfichier,p-nomfichier);
218     char chaine2[3000];
219     sprintf(chaine2,"%s%lu",chaine,id);
220     std::string namefic=chaine2;
221     MG_EXPORT exp;
222     exp.gmsh(mai,namefic);
223     }
224     nb=gest.get_nb_fem_maillage();
225     for (int i=0;i<nb;i++)
226     {
227     FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
228 gervaislavoie 304 if (desactivedeforme==1) mai->desactive_deforme();
229     unsigned long id=mai->get_id();
230 francois 283 strncpy(chaine,nomfichier,p-nomfichier);
231     char chaine2[3000];
232     sprintf(chaine2,"%s%lu",chaine,id);
233     std::string namefic=chaine2;
234     MG_EXPORT exp;
235     double coef=eps;
236     if (coef<0.)
237     {
238     coef=mai->calcul_coef();
239     sprintf(mess," Coefficient de deforme pour le maillage %lu est de %le",mai->get_id(),coef);
240     affiche(mess);
241     }
242     exp.gmsh(mai,namefic,coef);
243     }
244     }
245     if (action==4)
246     {
247     sprintf(mess," Conversion MAGiC vers code ASTER");
248     affiche(mess);
249     MG_FILE gest(nomfichier);
250     int nb=gest.get_nb_fem_maillage();
251     char chaine[1000];
252     for (int i=0;i<nb;i++)
253     {
254     FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
255     unsigned long id=mai->get_id();
256     strncpy(chaine,nomfichier,p-nomfichier);
257     char chaine2[3000];
258     sprintf(chaine2,"%s%lu",chaine,id);
259     std::string namefic=chaine2;
260     MG_EXPORT exp;
261 gervaislavoie 302 exp.aster2(mai,namefic,0,coderesu);
262 francois 283 }
263     }
264     if (((strcmp(extension,".ctt")==0) && (action==0)) || (action==7))
265     {
266     sprintf(mess," Conversion carte de taille vers GMSH");
267     affiche(mess);
268     FCT_GENERATEUR_3D<4> carte;
269     carte.lire(nomfichier);
270     strncpy(chaine,nomfichier,p-nomfichier);
271     std::string namefic=chaine;
272     MG_EXPORT exp;
273     exp.gmsh(&carte,namefic);
274     }
275 francois 295 #ifdef BREP_STEP
276 francois 283 if (((strcmp(extension,".step")==0) &&(action==0)) || (action==5))
277     {
278     sprintf(mess," Conversion STEP vers MAGiC");
279     affiche(mess);
280     STEP_IMPORT stepimport;
281     MG_GESTIONNAIRE gest;
282     stepimport.importer(gest,nomfichier) ;
283     strncpy(chaine,nomfichier,p-nomfichier);
284     std::string namefic=chaine;
285     namefic=namefic + ".magic";
286     int nb=gest.get_nb_mg_geometrie();
287     for (int i=0;i<nb;i++)
288     gest.get_mg_geometrie(i)->change_valeur_unite(unite);
289     gest.enregistrer(namefic.c_str());
290     }
291 francois 295 #endif
292     #ifdef BREP_SAT
293 francois 283 if (((strcmp(extension,".sat")==0) && (action==0)) || (action==6))
294     {
295     sprintf(mess," Conversion SAT vers MAGiC");
296     affiche(mess);
297     MG_GESTIONNAIRE gest;
298     ACIS_IMPORT acisimport;
299     acisimport.importer(gest,nomfichier) ;
300     strncpy(chaine,nomfichier,p-nomfichier);
301     std::string namefic=chaine;
302     namefic=namefic + ".magic";
303     int nb=gest.get_nb_mg_geometrie();
304     for (int i=0;i<nb;i++)
305     gest.get_mg_geometrie(i)->change_valeur_unite(unite);
306     gest.enregistrer(namefic.c_str());
307     }
308 francois 295 #endif
309 francois 283 if (((strcmp(extension,".brep")==0) && (action==0)) ||(action==8))
310     {
311     #ifdef BREP_OCC
312     sprintf(mess," Conversion OpenCascade vers MAGiC");
313     affiche(mess);
314     MG_GESTIONNAIRE gest;
315     OCC_IMPORT occimport;
316     MG_GEOMETRIE* mggeo=occimport.importer(gest,nomfichier,FICHIEROCC,unite,eps) ;
317     if (importtri>0) occimport.importer(gest,mggeo,eps2,importtri);
318     strncpy(chaine,nomfichier,p-nomfichier);
319     std::string namefic=chaine;
320     namefic=namefic + ".magic";
321     int nb=gest.get_nb_mg_geometrie();
322     for (int i=0;i<nb;i++)
323     gest.get_mg_geometrie(i)->change_valeur_unite(unite);
324     gest.enregistrer(namefic.c_str());
325     #endif
326     }
327     if (action==2)
328     {
329     #ifdef BREP_OCC
330     sprintf(mess," Conversion STEP via OpenCascade vers MAGiC avec epsilon=%lf",eps);
331     affiche(mess);
332     MG_GESTIONNAIRE gest;
333     OCC_IMPORT occimport;
334     MG_GEOMETRIE* mggeo=occimport.importer(gest,nomfichier,FICHIERSTEP,unite,eps) ;
335     if (importtri>0) occimport.importer(gest,mggeo,eps2,importtri);
336     strncpy(chaine,nomfichier,p-nomfichier);
337     std::string namefic=chaine;
338     namefic=namefic + ".magic";
339     int nb=gest.get_nb_mg_geometrie();
340     for (int i=0;i<nb;i++)
341     gest.get_mg_geometrie(i)->change_valeur_unite(unite);
342     gest.enregistrer(namefic.c_str());
343     #endif
344     }
345     if (action==1)
346     {
347     sprintf(mess," Creation d'un maillage FEM");
348     affiche(mess);
349     MG_FILE gest(nomfichier);
350     MG_MAILLAGE* mai=gest.get_mg_maillageid(nummai);
351     FEM_MAILLAGE* fem=new FEM_MAILLAGE(mai->get_mg_geometrie(),mai,degre);
352     gest.ajouter_fem_maillage(fem);
353     fem->construire(0);
354     gest.enregistrer(nomfichier);
355     }
356     if (action==9)
357     {
358     sprintf(mess," Importation des résultats ASTER vers MAGIC");
359     affiche(mess);
360     MG_FILE gest(nomfichier);
361     FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
362     MG_IMPORT imp;
363 gervaislavoie 302 imp.aster2(mai,fichieraster,nomfichier);
364 francois 283 gest.enregistrer(nomfichier);
365     }
366     if (action==10)
367     {
368     sprintf(mess," Conversion MAGiC vers OPTIS");
369     affiche(mess);
370     MG_FILE gest(nomfichier);
371     int nb=gest.get_nb_fem_maillage();
372     char chaine[1000];
373     for (int i=0;i<nb;i++)
374     {
375     FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
376     unsigned long id=mai->get_id();
377     strncpy(chaine,nomfichier,p-nomfichier);
378     char chaine2[3000];
379     sprintf(chaine2,"%s%lu",chaine,id);
380     std::string namefic=chaine2;
381     MG_EXPORT exp;
382 gervaislavoie 318 exp.aster2(mai,namefic,1,coderesu);
383 francois 283 }
384     }
385     if (action==11)
386     {
387     sprintf(mess," Importation des résultats d'optimisation vers MAGIC");
388     affiche(mess);
389     MG_FILE gest(nomfichier);
390     FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
391     MG_IMPORT imp;
392     imp.optis(mai,fichieroptis,fichierflottant,seuil,niveau,flottant);
393     gest.enregistrer(nomfichier);
394     }
395     if (action==12)
396     {
397     sprintf(mess," Selection dans MAGIC");
398     affiche(mess);
399     MG_FILE gest(nomfichier);
400     MG_GEOMETRIE* mggeo=gest.get_mg_geometrieid(numgeo);
401     MG_GROUPE_TOPOLOGIQUE* mggt;
402     if (numgt==0)
403     {
404     mggt=new MG_GROUPE_TOPOLOGIQUE;
405     mggeo->ajouter_mg_groupe_topologique(mggt);
406     }
407     else mggt=mggeo->get_mg_groupe_topologiqueid(numgt);
408     MG_ELEMENT_TOPOLOGIQUE* ele;
409     ele=mggeo->get_mg_sommetid(numentite);
410     if (ele==NULL) ele=mggeo->get_mg_areteid(numentite);
411     if (ele==NULL) ele=mggeo->get_mg_faceid(numentite);
412     if (ele==NULL) ele=mggeo->get_mg_volumeid(numentite);
413     if (ele!=NULL)
414     {
415     mggt->ajouter(ele);
416     if (expand)
417     {
418     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
419     ele->get_topologie_sousjacente(&lst);
420     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
421     for (MG_ELEMENT_TOPOLOGIQUE *ele2=lst.get_premier(it);ele2!=NULL;ele2=lst.get_suivant(it))
422     mggt->ajouter(ele2);
423     }
424     }
425     gest.enregistrer(nomfichier);
426     }
427     if (action==13)
428     {
429     sprintf(mess," Deselection dans MAGIC");
430     affiche(mess);
431     MG_FILE gest(nomfichier);
432     MG_GEOMETRIE* mggeo=gest.get_mg_geometrieid(numgeo);
433     MG_GROUPE_TOPOLOGIQUE* mggt=mggeo->get_mg_groupe_topologiqueid(numgt);
434     MG_ELEMENT_TOPOLOGIQUE* ele;
435     ele=mggeo->get_mg_sommetid(numentite);
436     if (ele==NULL) ele=mggeo->get_mg_areteid(numentite);
437     if (ele==NULL) ele=mggeo->get_mg_faceid(numentite);
438     if (ele==NULL) ele=mggeo->get_mg_volumeid(numentite);
439     if (ele!=NULL)
440     {
441     mggt->supprimer(ele);
442     if (expand)
443     {
444     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
445     ele->get_topologie_sousjacente(&lst);
446     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
447     for (MG_ELEMENT_TOPOLOGIQUE *ele2=lst.get_premier(it);ele2!=NULL;ele2=lst.get_suivant(it))
448     mggt->supprimer(ele2);
449     }
450     }
451     gest.enregistrer(nomfichier);
452     }
453     if (action==16)
454     {
455     sprintf(mess," Lissage de l'enveloppe externe du maillage optimisee");
456     affiche(mess);
457     MG_FILE gest(nomfichier);
458     FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
459     MG_GESTIONNAIRE gest2;
460     MG_LISSAGE liss;
461     liss.active_affichage(affiche);
462     liss.lisse(mai,gest2,eps,sigma,sigmaf,sigmag,gamma,iter_max,reactivation,bruitage,lissage,consopti,consimpose,consmailleurauto);
463     gest2.enregistrer(fichierout);
464     }
465     sprintf(mess,"Conversion reussie");
466    
467     if (action==14)
468     {
469     sprintf(mess," Analyse MAGIC");
470     affiche(mess);
471     MG_FILE gest(nomfichier);
472     char message [5000];
473     int nbgeo=gest.get_nb_mg_geometrie();
474     sprintf(message," %d geometrie :",nbgeo);
475     for (int i=0;i<nbgeo;i++)
476     {
477     MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
478     sprintf(message,"%s %lu ",message,geo->get_id());
479     }
480     affiche(message);
481     for (int i=0;i<nbgeo;i++)
482     {
483     MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
484     sprintf(message," geometrie %lu",geo->get_id());
485     affiche(message);
486     int nbmg=geo->get_nb_mg_groupe_topologique();
487     sprintf(message," %d groupe topologique",nbmg);
488     for (int j=0;j<nbmg;j++)
489     {
490     MG_GROUPE_TOPOLOGIQUE* mggp=geo->get_mg_groupe_topologique(j);
491     sprintf(message,"%s %lu ",message,mggp->get_id());
492     }
493     affiche(message);
494 francois 299 int nbvolume=geo->get_nb_mg_volume();
495     if (nbvolume==0) sprintf(message," %d volumes",geo->get_nb_mg_volume());
496     else sprintf(message," %d volumes : ",geo->get_nb_mg_volume());
497     for (int i=0;i<nbvolume;i++) sprintf(message,"%s %lu",message,geo->get_mg_volume(i)->get_id());
498 francois 290 affiche(message);
499 francois 299 int nbcoque=geo->get_nb_mg_coque();
500     if (nbcoque==0) sprintf(message," %d coques",geo->get_nb_mg_coque());
501     else sprintf(message," %d coques : ",geo->get_nb_mg_coque());
502     for (int i=0;i<nbcoque;i++) sprintf(message,"%s %lu",message,geo->get_mg_coque(i)->get_id());
503 francois 290 affiche(message);
504 francois 299 int nbpoutre=geo->get_nb_mg_poutre();
505     if (nbpoutre==0) sprintf(message," %d poutres",geo->get_nb_mg_poutre());
506     else sprintf(message," %d poutres : ",geo->get_nb_mg_poutre());
507     for (int i=0;i<nbpoutre;i++) sprintf(message,"%s %lu",message,geo->get_mg_poutre(i)->get_id());
508     affiche(message);
509 francois 290 sprintf(message," %d coquilles",geo->get_nb_mg_coquille());
510     affiche(message);
511 francois 283 sprintf(message," %d faces",geo->get_nb_mg_face());
512     affiche(message);
513     sprintf(message," %d aretes",geo->get_nb_mg_arete());
514     affiche(message);
515     sprintf(message," %d sommets",geo->get_nb_mg_sommet());
516     affiche(message);
517    
518     }
519     int nbmai=gest.get_nb_mg_maillage();
520     sprintf(message," %d maillage geometrique :",nbmai);
521     affiche(message);
522     for (int i=0;i<nbmai;i++)
523     {
524     MG_MAILLAGE* mai=gest.get_mg_maillage(i);
525     sprintf(message," maillage geometrique %lu ",mai->get_id());
526     affiche(message);
527     sprintf(message," %d noeuds",mai->get_nb_mg_noeud());
528     affiche(message);
529     sprintf(message," %d segments",mai->get_nb_mg_segment());
530     affiche(message);
531     sprintf(message," %d triangles",mai->get_nb_mg_triangle());
532     affiche(message);
533     sprintf(message," %d tetras",mai->get_nb_mg_tetra());
534     affiche(message);
535     }
536     int nbmaif=gest.get_nb_fem_maillage();
537     sprintf(message," %d maillage FEM :",nbmaif);
538     affiche(message);
539     for (int i=0;i<nbmaif;i++)
540     {
541     FEM_MAILLAGE* maif=gest.get_fem_maillage(i);
542     sprintf(message," maillage FEM %lu ",maif->get_id());
543     affiche(message);
544     sprintf(message," %d noeuds",maif->get_nb_fem_noeud());
545     affiche(message);
546 francois 309 sprintf(message," %d elements 1D",maif->get_nb_fem_element1());
547 francois 283 affiche(message);
548 francois 309 sprintf(message," %d elements 2D",maif->get_nb_fem_element2());
549 francois 283 affiche(message);
550 francois 309 sprintf(message," %d elements 3D",maif->get_nb_fem_element3());
551 francois 283 affiche(message);
552     }
553     int nbsol=gest.get_nb_mg_solution();
554     sprintf(message," %d solution maillage geometrique :",nbsol);
555     for (int i=0;i<nbsol;i++)
556     {
557     MG_SOLUTION* sol=gest.get_mg_solution(i);
558     sprintf(message,"%s %lu ",message,sol->get_id());
559     }
560     affiche(message);
561     int nbsolf=gest.get_nb_fem_solution();
562     sprintf(message," %d solution maillage FEM :",nbsolf);
563     for (int i=0;i<nbsolf;i++)
564     {
565     FEM_SOLUTION* sol=gest.get_fem_solution(i);
566     sprintf(message,"%s %lu ",message,sol->get_id());
567     }
568     affiche(message);
569     }
570     if (action==15)
571     {
572     sprintf(mess," Calcul de deforme");
573     affiche(mess);
574     MG_FILE gest(nomfichier);
575     FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
576     FEM_SOLUTION* sol1=gest.get_fem_solutionid(numsol1);
577     FEM_SOLUTION* sol2=gest.get_fem_solutionid(numsol2);
578     FEM_SOLUTION* sol3=gest.get_fem_solutionid(numsol3);
579     affiche(mess);
580     mai->calcul_deforme(sol1,numchamp1,sol2,numchamp2,sol3,numchamp3);
581     gest.enregistrer(nomfichier);
582     }
583     if (action==17)
584     {
585     sprintf(mess," Fusion de géométrie");
586     affiche(mess);
587     MG_FILE gest(nomfichier);
588     MG_GEOMETRIE* mggeo=gest.get_mg_geometrieid(numgeo);
589     VCT_MULTI_MODELE fusion(mggeo);
590     fusion.recherche_identite();
591     gest.enregistrer(nomfichier);
592     }
593     if (action==18)
594     {
595 francois 292 sprintf(mess," Création d'une pseudo-solution de la répartition de densite");
596 francois 283 affiche(mess);
597     MG_FILE gest(nomfichier);
598     FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
599     MG_IMPORT imp;
600 gervaislavoie 302 imp.import_densite(mai,fichierdensite,nomfichier,3);
601 francois 283 gest.enregistrer(nomfichier);
602     }
603 francois 293 if ((action==19) || (action==20))
604 francois 292 {
605 francois 296 unsigned long identity;
606     char typeccf[3];
607 gervaislavoie 318 typeccf[2]=0;
608 gervaislavoie 302 double valccf=0.5; // Valeur par défaut
609 francois 296 sscanf(ccf,"%lu:%c%c:%le\n",&identity,&typeccf[0],&typeccf[1],&valccf);
610     MC_GESTIONNAIRE mcgest;
611     int res=mcgest.existe(typeccf);
612     if (res==0)
613     {
614     sprintf(mess," ERREUR!!!! Code conditions aux limites inconnu");
615     affiche(mess);
616     return 0;
617     }
618     sprintf(mess," Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
619 francois 292 affiche(mess);
620     MG_FILE gest(nomfichier);
621     MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
622     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
623 francois 296 MG_VOLUME* ele;
624     ele=geo->get_mg_volumeid(identity);
625 francois 292 if (ele!=NULL) liste.ajouter(ele);
626 francois 296 else
627     {
628     MG_COQUE* ele;
629     ele=geo->get_mg_coqueid(identity);
630     if (ele!=NULL) liste.ajouter(ele);
631     else
632 francois 292 {
633 francois 296 MG_POUTRE* ele;
634     ele=geo->get_mg_poutreid(identity);
635 francois 292 if (ele!=NULL) liste.ajouter(ele);
636     else
637 francois 296 {
638     MG_FACE* ele;
639     ele=geo->get_mg_faceid(identity);
640     if (ele!=NULL) liste.ajouter(ele);
641     else
642     {
643     MG_ARETE* ele;
644     ele=geo->get_mg_areteid(identity);
645     if (ele!=NULL) liste.ajouter(ele);
646     else
647     {
648     MG_SOMMET* ele;
649     ele=geo->get_mg_sommetid(identity);
650     if (ele!=NULL) liste.ajouter(ele);
651     }
652     }
653     }
654 francois 292 }
655 francois 296 }
656 francois 292 if (action==20)
657     {
658     int nb=liste.get_nb();
659     for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
660     }
661     int nb=liste.get_nb();
662     for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf);;
663     gest.enregistrer(nomfichier);
664     }
665 francois 293 if ((action==21)|| (action==22))
666     {
667 francois 296 MC_GESTIONNAIRE mcgest;
668     unsigned long identity;
669 francois 293 char typeccf[3];
670     double valccf=0.5;
671     char c;
672     sscanf(ccf,"%c:%lu:%c%c:%le\n",&c,&identity,&typeccf[0],&typeccf[1],&valccf);
673 francois 296 int res=mcgest.existe(typeccf);
674     if (res==0)
675     {
676     sprintf(mess," ERREUR!!!! Code conditions aux limites inconnu");
677     affiche(mess);
678     return 0;
679     }
680     sprintf(mess," Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
681     affiche(mess);
682     MG_FILE gest(nomfichier);
683     MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
684     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
685 francois 293 MG_FACE* ele=NULL;
686     if (c=='F') ele=geo->get_mg_face(identity);
687     if (ele!=NULL) liste.ajouter(ele);
688     else
689     {
690     MG_ARETE* ele=NULL;
691     if (c=='A') ele=geo->get_mg_arete(identity);
692     if (ele!=NULL) liste.ajouter(ele);
693     else
694     {
695     MG_SOMMET* ele=NULL;
696     if (c=='S') ele=geo->get_mg_sommet(identity);
697     if (ele!=NULL) liste.ajouter(ele);
698 francois 296 else
699     {
700     MG_VOLUME* ele=NULL;
701     if (c=='V') ele=geo->get_mg_volume(identity);
702     if (ele!=NULL) liste.ajouter(ele);
703     else
704     {
705     MG_COQUE* ele=NULL;
706     if (c=='C') ele=geo->get_mg_coque(identity);
707     if (ele!=NULL) liste.ajouter(ele);
708     else
709     {
710     MG_POUTRE* ele=NULL;
711     if (c=='P') ele=geo->get_mg_poutre(identity);
712     if (ele!=NULL) liste.ajouter(ele);
713    
714     }
715     }
716     }
717 francois 293 }
718     }
719     if (action==22)
720     {
721     int nb=liste.get_nb();
722     for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
723     }
724     int nb=liste.get_nb();
725     for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf);;
726     gest.enregistrer(nomfichier);
727     }
728 francois 296 if (action==23)
729     {
730     sprintf(mess," Liste des conditions aux limites disponibles");
731     affiche(mess);
732     MC_GESTIONNAIRE mcgest;
733     int ok=0;
734     do
735     {
736     std::string code;
737     std::string description;
738     int res=mcgest.get_description(code,description);
739     if (res)
740     {
741     sprintf(mess," %s : %s",code.c_str(),description.c_str());
742     affiche(mess);
743     }
744     else ok=1;
745     }
746     while (ok==0);
747     }
748 francois 293
749 francois 296 sprintf(mess,"Fin");
750 francois 283 affiche(mess);
751    
752     return 0;
753     }
754    
755    
756    
757     #pragma package(smart_init)