ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/aster/src/mgaster.cpp
Revision: 327
Committed: Tue Mar 13 14:18:25 2012 UTC (13 years, 2 months ago) by francois
Original Path: magic/lib/optimisation/src/mgaster.cpp
File size: 2106 byte(s)
Log Message:
Calcul aster integre pour les volumes et pour les coques

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     MG_COQUE* coq=(MG_COQUE*)ele->get_mg_element_maillage()->get_lien_topologie();
62     char message[500];
63     sprintf(message," Calcul sur la coque %lu",coq->get_id());
64     affiche(message);
65     if (code!=NULL) exp.aster(coq,fem,nometude,0,code);
66     else exp.aster(coq,fem,nometude,0,(char*)"11111111");
67     }
68 francois 326 affiche((char*)" Calcul aster");
69     char nomfichiertmp[255];
70     sprintf(nomfichiertmp,"%s/as_run %s.export 1>aster.log 2>&1",getenv("PATHASTER"),nometude);
71     int codesortie=system(nomfichiertmp);
72     if (codesortie!=0)
73     {
74     sprintf(nomfichiertmp," Code de sortie aster : %d",codesortie);
75     affiche(nomfichiertmp);
76     }
77     affiche((char*)" Recuperation des resultats");
78     MG_IMPORT imp;
79     sprintf(nomfichiertmp,"%s.resu",nometude);
80     imp.aster(fem,nomfichiertmp,nometude);
81     affiche((char*)" Enregistrement MAGIC");
82     sprintf(nomfichiertmp,"%s.magic",nometude);
83     fem->get_mg_maillage()->get_gestionnaire()->enregistrer(nomfichiertmp);
84     affiche((char*)" Enregistrement GMSH");
85     exp.gmsh(fem,nometude);
86     }