| 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 |
francois |
470 |
MGASTER::MGASTER():affichageactif(0) |
| 18 |
francois |
326 |
{ |
| 19 |
|
|
|
| 20 |
|
|
} |
| 21 |
|
|
|
| 22 |
francois |
470 |
MGASTER::MGASTER(MGASTER &mdd):affichageactif(mdd.affichageactif) |
| 23 |
francois |
326 |
{ |
| 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 |
706 |
int MGASTER::calcule(char* nomfichierparam,FEM_MAILLAGE* fem,char *nometude,int typeetude,char *code,bool avecenreg) |
| 38 |
francois |
326 |
{ |
| 39 |
gervaislavoie |
382 |
if (affichageactif==1) affiche((char*)""); |
| 40 |
|
|
if (affichageactif==1) affiche((char*)"*************************"); |
| 41 |
|
|
if (affichageactif==1) affiche((char*)" Calcul par code aster"); |
| 42 |
|
|
if (affichageactif==1) affiche((char*)"*************************"); |
| 43 |
|
|
if (affichageactif==1) affiche((char*)""); |
| 44 |
|
|
if (affichageactif==1) affiche((char*)""); |
| 45 |
|
|
if (affichageactif==1) affiche((char*)""); |
| 46 |
francois |
326 |
MG_EXPORT exp; |
| 47 |
francois |
706 |
if (nomfichierparam!=NULL) exp.lire_params_aster(nomfichierparam); |
| 48 |
francois |
581 |
double version; |
| 49 |
francois |
327 |
int nb3=fem->get_nb_fem_element3(); |
| 50 |
|
|
if (nb3!=0) |
| 51 |
|
|
{ |
| 52 |
francois |
706 |
if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE)) |
| 53 |
mejrim |
597 |
{ |
| 54 |
|
|
FEM_ELEMENT3 *ele=fem->get_fem_element3(0); |
| 55 |
|
|
MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie(); |
| 56 |
|
|
char message[500]; |
| 57 |
|
|
sprintf(message," Calcul sur le volume %lu",vol->get_id()); |
| 58 |
|
|
if (affichageactif==1) affiche(message); |
| 59 |
francois |
706 |
if (code!=NULL) version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,code); |
| 60 |
|
|
else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,(char*)"00011111"); |
| 61 |
mejrim |
597 |
} |
| 62 |
|
|
if ((typeetude == MAGIC::CALCUL_ASTER::THERMIQUE)) |
| 63 |
|
|
{ |
| 64 |
|
|
FEM_ELEMENT3 *ele=fem->get_fem_element3(0); |
| 65 |
|
|
MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie(); |
| 66 |
|
|
char message[500]; |
| 67 |
|
|
sprintf(message," Calcul sur le volume %lu",vol->get_id()); |
| 68 |
|
|
if (affichageactif==1) affiche(message); |
| 69 |
|
|
version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::THERMIQUE,(char*)""); |
| 70 |
|
|
} |
| 71 |
francois |
327 |
} |
| 72 |
|
|
else |
| 73 |
|
|
{ |
| 74 |
francois |
706 |
if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE) || (typeetude == MAGIC::CALCUL_ASTER::CONTRAINTE_PLANE) || (typeetude == MAGIC::CALCUL_ASTER::DEFORMATION_PLANE) ) |
| 75 |
mejrim |
597 |
{ |
| 76 |
|
|
FEM_ELEMENT2 *ele=fem->get_fem_element2(0); |
| 77 |
|
|
MG_FACE* face=(MG_FACE*)ele->get_mg_element_maillage()->get_lien_topologie(); |
| 78 |
|
|
MG_COQUE* coq=face->get_mg_coface(0)->get_coquille()->get_mg_coque(); |
| 79 |
|
|
char message[500]; |
| 80 |
|
|
sprintf(message," Calcul sur la coque %lu",coq->get_id()); |
| 81 |
|
|
if (affichageactif==1) affiche(message); |
| 82 |
francois |
603 |
if (code!=NULL) version=exp.aster(coq,fem,nometude,typeetude,code); |
| 83 |
|
|
else version=exp.aster(coq,fem,nometude,typeetude,(char*)"11111111"); |
| 84 |
mejrim |
597 |
} |
| 85 |
francois |
327 |
} |
| 86 |
francois |
581 |
if (affichageactif==1) |
| 87 |
|
|
{ |
| 88 |
|
|
char messagever[255]; |
| 89 |
francois |
691 |
sprintf(messagever," Calcul aster \033[1;32m%.1lf\033[1;33m ",version); |
| 90 |
francois |
581 |
affiche(messagever); |
| 91 |
|
|
} |
| 92 |
francois |
326 |
char nomfichiertmp[255]; |
| 93 |
|
|
sprintf(nomfichiertmp,"%s/as_run %s.export 1>aster.log 2>&1",getenv("PATHASTER"),nometude); |
| 94 |
|
|
int codesortie=system(nomfichiertmp); |
| 95 |
|
|
if (codesortie!=0) |
| 96 |
|
|
{ |
| 97 |
|
|
sprintf(nomfichiertmp," Code de sortie aster : %d",codesortie); |
| 98 |
gervaislavoie |
382 |
if (affichageactif==1) affiche(nomfichiertmp); |
| 99 |
|
|
return codesortie; |
| 100 |
francois |
326 |
} |
| 101 |
gervaislavoie |
382 |
if (affichageactif==1) affiche((char*)" Recuperation des resultats"); |
| 102 |
francois |
326 |
MG_IMPORT imp; |
| 103 |
|
|
sprintf(nomfichiertmp,"%s.resu",nometude); |
| 104 |
francois |
581 |
std::string res=imp.aster(fem,nomfichiertmp,nometude); |
| 105 |
francois |
691 |
if (affichageactif==1) printf("%s",(char*)("\033[1;31m"+res+"\033[1;33m").c_str()); |
| 106 |
gervaislavoie |
382 |
if (avecenreg==true) |
| 107 |
|
|
{ |
| 108 |
|
|
if (affichageactif==1) affiche((char*)" Enregistrement MAGIC"); |
| 109 |
|
|
sprintf(nomfichiertmp,"%s.magic",nometude); |
| 110 |
|
|
fem->get_mg_maillage()->get_gestionnaire()->enregistrer(nomfichiertmp); |
| 111 |
|
|
if (affichageactif==1) affiche((char*)" Enregistrement GMSH"); |
| 112 |
|
|
exp.gmsh(fem,nometude); |
| 113 |
|
|
} |
| 114 |
mejrim |
597 |
} |
| 115 |
|
|
|
| 116 |
|
|
|
| 117 |
|
|
|
| 118 |
|
|
|
