ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/aster/src/mgaster.cpp
Revision: 329
Committed: Thu Mar 15 19:50:10 2012 UTC (13 years, 2 months ago) by gervaislavoie
Original Path: magic/lib/optimisation/src/mgaster.cpp
File size: 2177 byte(s)
Log Message:
Sélection de la coque dans un calcul Aster

File Contents

# User Rev Content
1 francois 326 #include "gestionversion.h"
2     #include "mgaster.h"
3     #include "fem_maillage.h"
4     #include "mg_gestionnaire.h"
5     #include "mg_export.h"
6     #include "mg_import.h"
7    
8    
9    
10    
11    
12    
13    
14    
15    
16    
17     MGASTER::MGASTER()
18     {
19    
20     }
21    
22     MGASTER::MGASTER(MGASTER &mdd)
23     {
24     }
25    
26    
27     MGASTER::~MGASTER()
28     {
29     }
30    
31     void MGASTER::active_affichage(void (*fonc)(char*))
32     {
33     affiche=fonc;
34     affichageactif=1;
35     }
36    
37 francois 327 void MGASTER::calcule(FEM_MAILLAGE* fem,char *nometude,char *code)
38 francois 326 {
39     affiche((char*)"");
40     affiche((char*)"*************************");
41     affiche((char*)" Calcul par code aster");
42     affiche((char*)"*************************");
43     affiche((char*)"");
44     affiche((char*)"");
45     affiche((char*)"");
46     MG_EXPORT exp;
47 francois 327 int nb3=fem->get_nb_fem_element3();
48     if (nb3!=0)
49     {
50     FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
51     MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
52     char message[500];
53     sprintf(message," Calcul sur le volume %lu",vol->get_id());
54     affiche(message);
55     if (code!=NULL) exp.aster(vol,fem,nometude,0,code);
56     else exp.aster(vol,fem,nometude,0,(char*)"00011111");
57     }
58     else
59     {
60     FEM_ELEMENT2 *ele=fem->get_fem_element2(0);
61 gervaislavoie 329 MG_FACE* face=(MG_FACE*)ele->get_mg_element_maillage()->get_lien_topologie();
62     MG_COQUE* coq=face->get_mg_coface(0)->get_coquille()->get_mg_coque();
63 francois 327 char message[500];
64     sprintf(message," Calcul sur la coque %lu",coq->get_id());
65     affiche(message);
66     if (code!=NULL) exp.aster(coq,fem,nometude,0,code);
67     else exp.aster(coq,fem,nometude,0,(char*)"11111111");
68     }
69 francois 326 affiche((char*)" Calcul aster");
70     char nomfichiertmp[255];
71     sprintf(nomfichiertmp,"%s/as_run %s.export 1>aster.log 2>&1",getenv("PATHASTER"),nometude);
72     int codesortie=system(nomfichiertmp);
73     if (codesortie!=0)
74     {
75     sprintf(nomfichiertmp," Code de sortie aster : %d",codesortie);
76     affiche(nomfichiertmp);
77     }
78     affiche((char*)" Recuperation des resultats");
79     MG_IMPORT imp;
80     sprintf(nomfichiertmp,"%s.resu",nometude);
81     imp.aster(fem,nomfichiertmp,nometude);
82     affiche((char*)" Enregistrement MAGIC");
83     sprintf(nomfichiertmp,"%s.magic",nometude);
84     fem->get_mg_maillage()->get_gestionnaire()->enregistrer(nomfichiertmp);
85     affiche((char*)" Enregistrement GMSH");
86     exp.gmsh(fem,nometude);
87     }