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)
{
ini_param();  

}

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


MGASTER::~MGASTER()
{
}

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

int MGASTER::calcule(char* nomfichierparam,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;
if (nomfichierparam!=NULL) exp.lire_params_aster(nomfichierparam);
double version;
int nb3=fem->get_nb_fem_element3();
if (nb3!=0)
  {
  if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE))
    {
    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::ELASTIQUE,code);
    else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,(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*)"");
    }
  if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF))
    {
    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::ELASTIQUE_GROT_GDEF,code);
    else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF,(char*)"00011111");
    }
   if ((typeetude == MAGIC::CALCUL_ASTER::PLASTIQUE))
    {
    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::PLASTIQUE,code);
    else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::PLASTIQUE,(char*)"00011111");
    }
  }
else 
  {
  if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE) || (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];
int choix=(int)param.get_valeur("Affiche_Code_Aster");
if (choix!=0) sprintf(nomfichiertmp,"xterm -e \"%s/as_run %s.export \"",getenv("PATHASTER"),nometude);
else 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);
  }
}


void MGASTER::ini_param(void)
{
param.ajouter("Affiche_Code_Aster",0,OT_PARAMETRES::DOUBLE,"0. Affichage redirigé vers le fichier aster.log 1 affichage dans un terminal");
char fichier[500];
sprintf(fichier,"%s/.magic",getenv("HOME"));
param.lire(fichier);
param.enregistrer(fichier);
 
}

Revision:	725
Committed:	Thu Sep 10 22:13:51 2015 UTC (9 years, 8 months ago) by francois
File size:	5319 byte(s)
Log Message:	parametrisation de la sortie de code aster via le fichier de parametre ~/.magic
#	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	francois	725	ini_param();
20	francois	326
21			}
22
23	francois	470	MGASTER::MGASTER(MGASTER &mdd):affichageactif(mdd.affichageactif)
24	francois	326	{
25	francois	725	ini_param();
26	francois	326	}
27
28
29			MGASTER::~MGASTER()
30			{
31			}
32
33			void MGASTER::active_affichage(void (fonc)(char))
34			{
35			affiche=fonc;
36			affichageactif=1;
37			}
38
39	francois	706	int MGASTER::calcule(char* nomfichierparam,FEM_MAILLAGE* fem,char nometude,int typeetude,char code,bool avecenreg)
40	francois	326	{
41	gervaislavoie	382	if (affichageactif==1) affiche((char*)"");
42			if (affichageactif==1) affiche((char)"************************");
43			if (affichageactif==1) affiche((char*)" Calcul par code aster");
44			if (affichageactif==1) affiche((char)"************************");
45			if (affichageactif==1) affiche((char*)"");
46			if (affichageactif==1) affiche((char*)"");
47			if (affichageactif==1) affiche((char*)"");
48	francois	326	MG_EXPORT exp;
49	francois	706	if (nomfichierparam!=NULL) exp.lire_params_aster(nomfichierparam);
50	francois	581	double version;
51	francois	327	int nb3=fem->get_nb_fem_element3();
52			if (nb3!=0)
53			{
54	francois	706	if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE))
55	mejrim	597	{
56			FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
57			MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
58			char message[500];
59			sprintf(message," Calcul sur le volume %lu",vol->get_id());
60			if (affichageactif==1) affiche(message);
61	francois	706	if (code!=NULL) version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,code);
62			else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,(char*)"00011111");
63	mejrim	597	}
64			if ((typeetude == MAGIC::CALCUL_ASTER::THERMIQUE))
65			{
66			FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
67			MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
68			char message[500];
69			sprintf(message," Calcul sur le volume %lu",vol->get_id());
70			if (affichageactif==1) affiche(message);
71			version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::THERMIQUE,(char*)"");
72			}
73	couturad	719	if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF))
74			{
75			FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
76			MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
77			char message[500];
78			sprintf(message," Calcul sur le volume %lu",vol->get_id());
79			if (affichageactif==1) affiche(message);
80			if (code!=NULL) version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF,code);
81			else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF,(char*)"00011111");
82			}
83			if ((typeetude == MAGIC::CALCUL_ASTER::PLASTIQUE))
84			{
85			FEM_ELEMENT3 *ele=fem->get_fem_element3(0);
86			MG_VOLUME* vol=(MG_VOLUME*)ele->get_mg_element_maillage()->get_lien_topologie();
87			char message[500];
88			sprintf(message," Calcul sur le volume %lu",vol->get_id());
89			if (affichageactif==1) affiche(message);
90			if (code!=NULL) version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::PLASTIQUE,code);
91			else version=exp.aster(vol,fem,nometude,MAGIC::CALCUL_ASTER::PLASTIQUE,(char*)"00011111");
92			}
93	francois	327	}
94			else
95			{
96	francois	706	if ((typeetude == MAGIC::CALCUL_ASTER::ELASTIQUE) \|\| (typeetude == MAGIC::CALCUL_ASTER::CONTRAINTE_PLANE) \|\| (typeetude == MAGIC::CALCUL_ASTER::DEFORMATION_PLANE) )
97	mejrim	597	{
98			FEM_ELEMENT2 *ele=fem->get_fem_element2(0);
99			MG_FACE* face=(MG_FACE*)ele->get_mg_element_maillage()->get_lien_topologie();
100			MG_COQUE* coq=face->get_mg_coface(0)->get_coquille()->get_mg_coque();
101			char message[500];
102			sprintf(message," Calcul sur la coque %lu",coq->get_id());
103			if (affichageactif==1) affiche(message);
104	francois	603	if (code!=NULL) version=exp.aster(coq,fem,nometude,typeetude,code);
105			else version=exp.aster(coq,fem,nometude,typeetude,(char*)"11111111");
106	mejrim	597	}
107	francois	327	}
108	francois	581	if (affichageactif==1)
109			{
110			char messagever[255];
111	francois	691	sprintf(messagever," Calcul aster \033[1;32m%.1lf\033[1;33m ",version);
112	francois	581	affiche(messagever);
113			}
114	francois	326	char nomfichiertmp[255];
115	francois	725	int choix=(int)param.get_valeur("Affiche_Code_Aster");
116			if (choix!=0) sprintf(nomfichiertmp,"xterm -e \"%s/as_run %s.export \"",getenv("PATHASTER"),nometude);
117			else sprintf(nomfichiertmp,"%s/as_run %s.export 1>aster.log 2>&1",getenv("PATHASTER"),nometude);
118	francois	326	int codesortie=system(nomfichiertmp);
119			if (codesortie!=0)
120			{
121			sprintf(nomfichiertmp," Code de sortie aster : %d",codesortie);
122	gervaislavoie	382	if (affichageactif==1) affiche(nomfichiertmp);
123			return codesortie;
124	francois	326	}
125	gervaislavoie	382	if (affichageactif==1) affiche((char*)" Recuperation des resultats");
126	francois	326	MG_IMPORT imp;
127			sprintf(nomfichiertmp,"%s.resu",nometude);
128	francois	581	std::string res=imp.aster(fem,nomfichiertmp,nometude);
129	francois	691	if (affichageactif==1) printf("%s",(char*)("\033[1;31m"+res+"\033[1;33m").c_str());
130	gervaislavoie	382	if (avecenreg==true)
131			{
132			if (affichageactif==1) affiche((char*)" Enregistrement MAGIC");
133			sprintf(nomfichiertmp,"%s.magic",nometude);
134			fem->get_mg_maillage()->get_gestionnaire()->enregistrer(nomfichiertmp);
135			if (affichageactif==1) affiche((char*)" Enregistrement GMSH");
136			exp.gmsh(fem,nometude);
137			}
138	mejrim	597	}
139
140
141	francois	725	void MGASTER::ini_param(void)
142			{
143			param.ajouter("Affiche_Code_Aster",0,OT_PARAMETRES::DOUBLE,"0. Affichage redirigé vers le fichier aster.log 1 affichage dans un terminal");
144			char fichier[500];
145			sprintf(fichier,"%s/.magic",getenv("HOME"));
146			param.lire(fichier);
147			param.enregistrer(fichier);
148
149			}
150	mejrim	597