aster/src/mgaster.cpp

#include "gestionversion.h"
#include "mgaster.h"
#include "fem_maillage.h"
#include "mg_gestionnaire.h"
#include "mg_export.h"
#include "mg_import.h"


MGASTER::MGASTER():affichageactif(0)
{

}

MGASTER::MGASTER(MGASTER &mdd):affichageactif(mdd.affichageactif)
{
}


MGASTER::~MGASTER()
{
}

void MGASTER::active_affichage(void (*fonc)(char*))
{
affiche=fonc;
affichageactif=1;
}

int MGASTER::calcule(FEM_MAILLAGE* fem,char *nometude,int typeetude,char *code,bool avecenreg)
{
if (affichageactif==1) affiche((char*)"");
if (affichageactif==1) affiche((char*)"*************************");
if (affichageactif==1) affiche((char*)"  Calcul par code aster");
if (affichageactif==1) affiche((char*)"*************************");
if (affichageactif==1) affiche((char*)"");
if (affichageactif==1) affiche((char*)"");
if (affichageactif==1) affiche((char*)""); 
MG_EXPORT exp;
double version;
int nb3=fem->get_nb_fem_element3();
if (nb3!=0)
  {
  if ((typeetude == MAGIC::CALCUL_ASTER::MECANIQUE))
    {
    FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
    MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
    char message[500];
    sprintf(message,"        Calcul sur le volume %lu",vol->get_id());
    if (affichageactif==1) affiche(message);
    if (code!=NULL) version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::MECANIQUE,code);
    else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::MECANIQUE,(char*)"00011111");
    }
  if ((typeetude == MAGIC::CALCUL_ASTER::THERMIQUE))
   {
      FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
      MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
      char message[500];
      sprintf(message,"        Calcul sur le volume %lu",vol->get_id());
      if (affichageactif==1) affiche(message);
      version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::THERMIQUE,(char*)"");
    }
  }
else 
  {
  if ((typeetude == MAGIC::CALCUL_ASTER::MECANIQUE) || (typeetude == MAGIC::CALCUL_ASTER::CONTRAINTE_PLANE) || (typeetude == MAGIC::CALCUL_ASTER::DEFORMATION_PLANE) )
    {
    FEM_ELEMENT2 *ele=fem->get_fem_element2(0);
    MG_FACE* face=(MG_FACE*)ele->get_mg_element_maillage()->get_lien_topologie();
    MG_COQUE* coq=face->get_mg_coface(0)->get_coquille()->get_mg_coque();
    char message[500];
    sprintf(message,"        Calcul sur la coque %lu",coq->get_id());
    if (affichageactif==1) affiche(message);
    if (code!=NULL) version=exp.aster(coq,fem,nometude,typeetude,code);
    else version=exp.aster(coq,fem,nometude,typeetude,(char*)"11111111");
    }
  }
if (affichageactif==1) 
  {
  char messagever[255];
  sprintf(messagever,"        Calcul aster \033[1;32m%.1lf\033[1;33m ",version);
  affiche(messagever);
  }
char nomfichiertmp[255];
sprintf(nomfichiertmp,"%s/as_run %s.export  1>aster.log 2>&1",getenv("PATHASTER"),nometude);
int codesortie=system(nomfichiertmp);
if (codesortie!=0)
        {
        sprintf(nomfichiertmp,"          Code de sortie aster : %d",codesortie);
        if (affichageactif==1) affiche(nomfichiertmp);
        return codesortie;
        }
if (affichageactif==1) affiche((char*)"        Recuperation des resultats");
MG_IMPORT imp;
sprintf(nomfichiertmp,"%s.resu",nometude);
std::string res=imp.aster(fem,nomfichiertmp,nometude);
if (affichageactif==1) printf("%s",(char*)("\033[1;31m"+res+"\033[1;33m").c_str());
if (avecenreg==true)
  {
  if (affichageactif==1) affiche((char*)"        Enregistrement MAGIC");
  sprintf(nomfichiertmp,"%s.magic",nometude);
  fem->get_mg_maillage()->get_gestionnaire()->enregistrer(nomfichiertmp);
  if (affichageactif==1) affiche((char*)"        Enregistrement GMSH");
  exp.gmsh(fem,nometude);
  }
}


Revision:	691
Committed:	Mon Jul 13 22:17:05 2015 UTC (9 years, 10 months ago) by francois
File size:	3667 byte(s)
Log Message:	Sortir des resultats directement en time step sous gmsh. Pour sela faire des fem_solutions sous MAGIC avec un nom de solution pour cahque pas de calcul selon le format suivant : nomsolution%numpas. Les pas doivent commencer 0 et etre incrementé de 1. Il n epeut pas y avoir de trou dans les pas. Faire renumérotation si nécessaire
#	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	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	mejrim	597	int MGASTER::calcule(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	581	double version;
48	francois	327	int nb3=fem->get_nb_fem_element3();
49			if (nb3!=0)
50			{
51	mejrim	597	if ((typeetude == MAGIC::CALCUL_ASTER::MECANIQUE))
52			{
53			FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
54			MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
55			char message[500];
56			sprintf(message," Calcul sur le volume %lu",vol->get_id());
57			if (affichageactif==1) affiche(message);
58			if (code!=NULL) version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::MECANIQUE,code);
59			else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::MECANIQUE,(char*)"00011111");
60			}
61			if ((typeetude == MAGIC::CALCUL_ASTER::THERMIQUE))
62			{
63			FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
64			MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
65			char message[500];
66			sprintf(message," Calcul sur le volume %lu",vol->get_id());
67			if (affichageactif==1) affiche(message);
68			version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::THERMIQUE,(char*)"");
69			}
70	francois	327	}
71			else
72			{
73	francois	603	if ((typeetude == MAGIC::CALCUL_ASTER::MECANIQUE) \|\| (typeetude == MAGIC::CALCUL_ASTER::CONTRAINTE_PLANE) \|\| (typeetude == MAGIC::CALCUL_ASTER::DEFORMATION_PLANE) )
74	mejrim	597	{
75			FEM_ELEMENT2 *ele=fem->get_fem_element2(0);
76			MG_FACE* face=(MG_FACE*)ele->get_mg_element_maillage()->get_lien_topologie();
77			MG_COQUE* coq=face->get_mg_coface(0)->get_coquille()->get_mg_coque();
78			char message[500];
79			sprintf(message," Calcul sur la coque %lu",coq->get_id());
80			if (affichageactif==1) affiche(message);
81	francois	603	if (code!=NULL) version=exp.aster(coq,fem,nometude,typeetude,code);
82			else version=exp.aster(coq,fem,nometude,typeetude,(char*)"11111111");
83	mejrim	597	}
84	francois	327	}
85	francois	581	if (affichageactif==1)
86			{
87			char messagever[255];
88	francois	691	sprintf(messagever," Calcul aster \033[1;32m%.1lf\033[1;33m ",version);
89	francois	581	affiche(messagever);
90			}
91	francois	326	char nomfichiertmp[255];
92			sprintf(nomfichiertmp,"%s/as_run %s.export 1>aster.log 2>&1",getenv("PATHASTER"),nometude);
93			int codesortie=system(nomfichiertmp);
94			if (codesortie!=0)
95			{
96			sprintf(nomfichiertmp," Code de sortie aster : %d",codesortie);
97	gervaislavoie	382	if (affichageactif==1) affiche(nomfichiertmp);
98			return codesortie;
99	francois	326	}
100	gervaislavoie	382	if (affichageactif==1) affiche((char*)" Recuperation des resultats");
101	francois	326	MG_IMPORT imp;
102			sprintf(nomfichiertmp,"%s.resu",nometude);
103	francois	581	std::string res=imp.aster(fem,nomfichiertmp,nometude);
104	francois	691	if (affichageactif==1) printf("%s",(char*)("\033[1;31m"+res+"\033[1;33m").c_str());
105	gervaislavoie	382	if (avecenreg==true)
106			{
107			if (affichageactif==1) affiche((char*)" Enregistrement MAGIC");
108			sprintf(nomfichiertmp,"%s.magic",nometude);
109			fem->get_mg_maillage()->get_gestionnaire()->enregistrer(nomfichiertmp);
110			if (affichageactif==1) affiche((char*)" Enregistrement GMSH");
111			exp.gmsh(fem,nometude);
112			}
113	mejrim	597	}
114
115
116
117