ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/mgoperation/src/main.cpp
Revision: 326
Committed: Mon Mar 12 15:34:00 2012 UTC (13 years, 3 months ago) by francois
File size: 28148 byte(s)
Log Message:
integration du calcul par code aster dans magic. Seul les volumes sont faits

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