mstruct_ves_file.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mstruct_ves_file.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_ves_file.h"
#include "mstruct_definition.h"
#include "mg_definition.h"
#include "mstruct_analyse.h"
#include "mstruct_analyse_cao.h"
#include "mstruct_analyse_propriete_massique.h"
#include "mstruct_analyse_champ.h"
#include "mstruct_analyse_orientation.h"
#include "mstruct_analyse_erosion.h"
#include "mstruct_analyse_mg_maillage.h"
#include "mstruct_analyse_fem_maillage.h"
#include "mstruct_analyse_modules_elasticite.h"
#include "mstruct_analyse_conductivite_thermique.h"
#include "mstruct_analyse_energie_hill.h"
#include "ot_boite_3d.h"
#include <string.h>
#include <map>
 
MSTRUCT_VES_FILE::MSTRUCT_VES_FILE(void)
{
  m_nom_fichier_magic=NULL;
  m_id_mgcg_modele=-1;
  m_id_mgcg_assemblage=-1;
  m_id_mg_geometrie=-1;  
  m_nom_fichier_carte=NULL;
  m_id_mg_maillage=-1;
  m_id_fem_maillage=-1;
  m_boite3D_ves = BOITE_3D(0.0,0.0,0.0,1.0,1.0,1.0);
  m_precision=1e-06;
  m_temps_geometrie=0.0;
  m_temps_materiau=0.0;
  m_temps_carte=0.0;
  m_temps_maillage=0.0;
  m_temps_fem_maillage=0.0;
  m_temps_etude=0.0;
  m_temps_calcul=0.0;
  m_temps_analyse=0.0;
}
 
MSTRUCT_VES_FILE::~MSTRUCT_VES_FILE(void)
{
  if(m_nom_fichier_magic!=NULL) delete [] m_nom_fichier_magic;
  if(m_nom_fichier_carte!=NULL) delete [] m_nom_fichier_carte;
  supprimer_tout_analyse();
}
 
char* MSTRUCT_VES_FILE::get_nom_fichier_magic(void)
{
  return m_nom_fichier_magic;
}
 
void MSTRUCT_VES_FILE::change_nom_fichier_magic(char* nom_fichier_magic)
{
  if(m_nom_fichier_magic!=NULL) delete [] m_nom_fichier_magic;
  m_nom_fichier_magic=new char[strlen(nom_fichier_magic)+1];
  strcpy(m_nom_fichier_magic,nom_fichier_magic);
}
 
long int MSTRUCT_VES_FILE::get_id_mgcg_modele(void)
{
  return m_id_mgcg_modele;
}
 
void MSTRUCT_VES_FILE::change_id_mgcg_modele(long int id)
{
   m_id_mgcg_modele=id;
}
 
long int MSTRUCT_VES_FILE::get_id_mgcg_assemblage(void)
{
  return m_id_mgcg_assemblage;
}
 
void MSTRUCT_VES_FILE::change_id_mgcg_assemblage(long int id)
{
  m_id_mgcg_assemblage=id;
}
 
long int MSTRUCT_VES_FILE::get_id_mg_geometrie(void)
{
  return m_id_mg_geometrie;
}
 
void MSTRUCT_VES_FILE::change_id_mg_geometrie(long int id)
{
  m_id_mg_geometrie=id;
}
 
char* MSTRUCT_VES_FILE::get_nom_fichier_carte(void)
{
return m_nom_fichier_carte;
}
 
void MSTRUCT_VES_FILE::change_nom_fichier_carte(char* nom_fichier_carte)
{
  if(m_nom_fichier_carte!=NULL) delete [] m_nom_fichier_carte;
  m_nom_fichier_carte=new char[strlen(nom_fichier_carte)+1];
  strcpy(m_nom_fichier_carte,nom_fichier_carte);
}
 
long int MSTRUCT_VES_FILE::get_id_mg_maillage(void)
{
  return m_id_mg_maillage;
}
 
void MSTRUCT_VES_FILE::change_id_mg_maillage(long int id)
{
  m_id_mg_maillage=id;
}
 
long int MSTRUCT_VES_FILE::get_id_fem_maillage(void)
{
  return m_id_fem_maillage;
}
 
void MSTRUCT_VES_FILE::change_id_fem_maillage(long int id)
{
  m_id_fem_maillage=id;
}
 
BOITE_3D MSTRUCT_VES_FILE::get_boite_3D_ves(void)
{
  return m_boite3D_ves;
}
 
void MSTRUCT_VES_FILE::change_boite_3D_ves(BOITE_3D boite3D)
{
  m_boite3D_ves=boite3D;
}
 
double MSTRUCT_VES_FILE::get_precision(void)
{
  return m_precision;
}
 
void MSTRUCT_VES_FILE::change_precision(double precision)
{
  m_precision=precision;
}
 
double MSTRUCT_VES_FILE::get_temps_geometrie(void)
{
  return m_temps_geometrie;
}
 
void MSTRUCT_VES_FILE::change_temps_geometrie(double temps)
{
  m_temps_geometrie=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_materiau(void)
{
  return m_temps_materiau;
}
 
void MSTRUCT_VES_FILE::change_temps_materiau(double temps)
{
  m_temps_materiau=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_carte(void)
{
  return m_temps_carte;
}
 
void MSTRUCT_VES_FILE::change_temps_carte(double temps)
{
  m_temps_carte=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_maillage(void)
{
  return m_temps_maillage;
}
 
void MSTRUCT_VES_FILE::change_temps_maillage(double temps)
{
  m_temps_maillage=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_fem_maillage(void)
{
  return m_temps_fem_maillage;
}
 
void MSTRUCT_VES_FILE::change_temps_fem_maillage(double temps)
{
  m_temps_fem_maillage=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_etude(void)
{
  return m_temps_etude;
}
 
void MSTRUCT_VES_FILE::change_temps_etude(double temps)
{
  m_temps_etude=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_calcul(void)
{
  return m_temps_calcul;
}
 
void MSTRUCT_VES_FILE::change_temps_calcul(double temps)
{
  m_temps_calcul=temps;
}
 
double MSTRUCT_VES_FILE::get_temps_analyse(void)
{
  return m_temps_analyse;
}
 
void MSTRUCT_VES_FILE::change_temps_analyse(double temps)
{
  m_temps_analyse=temps;
}
 
long int MSTRUCT_VES_FILE::get_nb_analyse(void)
{
  return m_map_analyse_ves.size();
}
 
int MSTRUCT_VES_FILE::ajouter_analyse(MSTRUCT_ANALYSE* analyse_ves)
{
  std::pair<std::map<std::string,MSTRUCT_ANALYSE*>::iterator,bool> p=m_map_analyse_ves.insert(std::pair<std::string,MSTRUCT_ANALYSE*>(analyse_ves->get_identifiant(),analyse_ves));
  if (!p.second)  
    return FAIL;
  return OK;
}
 
MSTRUCT_ANALYSE* MSTRUCT_VES_FILE::get_premiere_analyse(std::map< std::string, MSTRUCT_ANALYSE* >::iterator& it)
{
  it=m_map_analyse_ves.begin();
  if(it==m_map_analyse_ves.end())  
    return NULL;
  return it->second; 
}
 
MSTRUCT_ANALYSE* MSTRUCT_VES_FILE::get_suivante_analyse(std::map< std::string, MSTRUCT_ANALYSE* >::iterator& it)
{
  it++;
  if(it==m_map_analyse_ves.end())  
    return NULL;
  return it->second;
}
 
MSTRUCT_ANALYSE* MSTRUCT_VES_FILE::get_analyse(std::string identifiant)
{
  std::map<std::string,MSTRUCT_ANALYSE*>::iterator it=m_map_analyse_ves.find(identifiant);
  if (it==m_map_analyse_ves.end())
  {
    return NULL;
  }
  return ((*it).second);
}
 
int MSTRUCT_VES_FILE::supprimer_analyse(std::string identifiant)
{
  std::map<std::string,MSTRUCT_ANALYSE*>::iterator it=m_map_analyse_ves.find(identifiant);
  if (it==m_map_analyse_ves.end())
  {
    m_map_analyse_ves.erase(it);
    delete it->second;
    return OK;
  }
  return FAIL;
}
 
int MSTRUCT_VES_FILE::supprimer_tout_analyse(void)
{
  std::map<std::string,MSTRUCT_ANALYSE*>::iterator it;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it);analyse!=NULL;analyse=get_suivante_analyse(it))
  {
    delete analyse;
  }
  m_map_analyse_ves.clear();
  return OK;
}
 
int MSTRUCT_VES_FILE::enregistrer(char* nom_fichier_ves)
{
  std::ofstream ofstrm(nom_fichier_ves, std::ios::out | std::ios::binary);
  size_t len_fichier_magic;
  if(m_nom_fichier_magic!=NULL)
  {
    len_fichier_magic = strlen(m_nom_fichier_magic);
    ofstrm.write((char*)&len_fichier_magic,sizeof(size_t));
    ofstrm.write((char*)m_nom_fichier_magic,len_fichier_magic);
  }
  else
  {
    len_fichier_magic=0;
    ofstrm.write((char*)&len_fichier_magic,sizeof(size_t));
  }
  ofstrm.write((char*)&m_id_mgcg_modele,sizeof(long));
  ofstrm.write((char*)&m_id_mgcg_assemblage,sizeof(long));
  ofstrm.write((char*)&m_id_mg_geometrie,sizeof(long));
  size_t len_fichier_carte;
  if(m_nom_fichier_carte!=NULL)
  {
    len_fichier_carte=strlen(m_nom_fichier_carte);
    ofstrm.write((char*)&len_fichier_carte,sizeof(size_t));
    ofstrm.write((char*)m_nom_fichier_carte,len_fichier_carte);
  }  
  else 
  {
    len_fichier_carte=0;
    ofstrm.write((char*)&len_fichier_carte,sizeof(size_t)); 
  }   
  ofstrm.write((char*)&m_id_mg_maillage,sizeof(long));
  ofstrm.write((char*)&m_id_fem_maillage,sizeof(long));
  double xmin=m_boite3D_ves.get_xmin();
  ofstrm.write((char*)&xmin,sizeof(double));
  double ymin=m_boite3D_ves.get_ymin();
  ofstrm.write((char*)&ymin,sizeof(double));
  double zmin=m_boite3D_ves.get_zmin();
  ofstrm.write((char*)&zmin,sizeof(double));
  double xmax=m_boite3D_ves.get_xmax();
  ofstrm.write((char*)&xmax,sizeof(double));
  double ymax=m_boite3D_ves.get_ymax();
  ofstrm.write((char*)&ymax,sizeof(double));
  double zmax=m_boite3D_ves.get_zmax();
  ofstrm.write((char*)&zmax,sizeof(double));
  ofstrm.write((char*)&m_precision,sizeof(double));
  ofstrm.write((char*)&m_temps_geometrie,sizeof(double));
  ofstrm.write((char*)&m_temps_materiau,sizeof(double));
  ofstrm.write((char*)&m_temps_carte,sizeof(double));
  ofstrm.write((char*)&m_temps_maillage,sizeof(double));
  ofstrm.write((char*)&m_temps_fem_maillage,sizeof(double));
  ofstrm.write((char*)&m_temps_etude,sizeof(double));
  ofstrm.write((char*)&m_temps_calcul,sizeof(double));
  ofstrm.write((char*)&m_temps_analyse,sizeof(double));
  long nb_analyse = get_nb_analyse();
  ofstrm.write((char*)&nb_analyse,sizeof(long));
  std::map<std::string,MSTRUCT_ANALYSE*>::iterator it;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it);analyse!=NULL;analyse=get_suivante_analyse(it))
  {
    analyse->enregistrer(ofstrm);
  }
  ofstrm.close();
return 0;
}
 
int MSTRUCT_VES_FILE::ouvrir(char* nom_fichier_ves)
{
  std::ifstream ifstrm(nom_fichier_ves, std::ios::in | std::ios::binary);
  if(!ifstrm.is_open()) 
  {
    std::cerr << "*** MSTRUCT_VES_FILE::ouvrir : ERREUR ***" << std::endl;
    return FAIL;
  }
  size_t len_fichier_magic;
  ifstrm.read((char*)&len_fichier_magic,sizeof(size_t));
  if(len_fichier_magic>0)
  {
    m_nom_fichier_magic = new char[len_fichier_magic+1];
    ifstrm.read((char*)m_nom_fichier_magic,len_fichier_magic);
    m_nom_fichier_magic[len_fichier_magic]='\0';
  }
  
  ifstrm.read((char*)&m_id_mgcg_modele,sizeof(long));
  ifstrm.read((char*)&m_id_mgcg_assemblage,sizeof(long));
  ifstrm.read((char*)&m_id_mg_geometrie,sizeof(long));
  size_t len_fichier_carte;
  ifstrm.read((char*)&len_fichier_carte,sizeof(size_t));
  if(len_fichier_carte>0)
  {
    m_nom_fichier_carte = new char[len_fichier_carte+1];  
    ifstrm.read((char*)m_nom_fichier_carte,len_fichier_carte);
    m_nom_fichier_carte[len_fichier_carte]='\0';
  }
  else m_nom_fichier_carte=NULL;
  ifstrm.read((char*)&m_id_mg_maillage,sizeof(long));
  ifstrm.read((char*)&m_id_fem_maillage,sizeof(long));
  double xmin;
  ifstrm.read((char*)&xmin,sizeof(double));
  double ymin;
  ifstrm.read((char*)&ymin,sizeof(double));
  double zmin;
  ifstrm.read((char*)&zmin,sizeof(double));
  double xmax;
  ifstrm.read((char*)&xmax,sizeof(double));
  double ymax;
  ifstrm.read((char*)&ymax,sizeof(double));
  double zmax;
  ifstrm.read((char*)&zmax,sizeof(double));
  m_boite3D_ves = BOITE_3D(xmin,ymin,zmin,xmax,ymax,zmax);
  ifstrm.read((char*)&m_precision,sizeof(double));
  ifstrm.read((char*)&m_temps_geometrie,sizeof(double));
  ifstrm.read((char*)&m_temps_materiau,sizeof(double));
  ifstrm.read((char*)&m_temps_carte,sizeof(double));
  ifstrm.read((char*)&m_temps_maillage,sizeof(double));
  ifstrm.read((char*)&m_temps_fem_maillage,sizeof(double));
  ifstrm.read((char*)&m_temps_etude,sizeof(double));
  ifstrm.read((char*)&m_temps_calcul,sizeof(double));
  ifstrm.read((char*)&m_temps_analyse,sizeof(double));
  long nb_analyse;
  ifstrm.read((char*)&nb_analyse,sizeof(long));
  for(long i=0;i<nb_analyse;i++)
  {
    long type_analyse;
    ifstrm.read((char*)&type_analyse,sizeof(long));
    MSTRUCT_ANALYSE *analyse;
    if(type_analyse==MSTRUCT::TYPE_ANALYSE::CHAMP)
    {
      analyse = new MSTRUCT_ANALYSE_CHAMP;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::ORIENTATION)
    {
      analyse = new MSTRUCT_ANALYSE_ORIENTATION;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::CAO)
    {
      analyse = new MSTRUCT_ANALYSE_CAO;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::PROPRIETE_MASSIQUE)
    {
      analyse = new MSTRUCT_ANALYSE_PROPRIETE_MASSIQUE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MAILLAGE_MG)
    {
      analyse = new MSTRUCT_ANALYSE_MG_MAILLAGE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MAILLAGE_FEM)
    {
      analyse = new MSTRUCT_ANALYSE_FEM_MAILLAGE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::EROSION)
    {
      analyse = new MSTRUCT_ANALYSE_EROSION;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MODULES_ELASTICITE)
    {
      analyse = new MSTRUCT_ANALYSE_MODULES_ELASTICITE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::CONDUCTIVITE_THERMIQUE)
    {
      analyse = new MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::ENERGIE_HILL)
    {
      analyse = new MSTRUCT_ANALYSE_ENERGIE_HILL;
      analyse->ouvrir(ifstrm);      
    }
    ajouter_analyse(analyse);
  }
  ifstrm.close();
return OK;
    
}
 
void MSTRUCT_VES_FILE::affiche_contenu(fonction_affiche* fonc)
{
  char ligne[5000];
  sprintf(ligne,"MSTRUCT_VES_FILE"); fonc(ligne);
  sprintf(ligne,"-> Nom fichier MAGiC : %s",m_nom_fichier_magic); fonc(ligne);
  sprintf(ligne,"-> ID MG_CG_MODELE : %li",m_id_mgcg_modele); fonc(ligne);
  sprintf(ligne,"-> ID MG_CG_ASSEMBLAGE : %li",m_id_mgcg_assemblage); fonc(ligne);
  sprintf(ligne,"-> ID MG_GEOMETRIE : %li",m_id_mg_geometrie); fonc(ligne);
  sprintf(ligne,"-> Nom fichier carte : %s",m_nom_fichier_carte); fonc(ligne);
  sprintf(ligne,"-> ID MG_MAILLAGE : %li",m_id_mg_maillage); fonc(ligne);
  sprintf(ligne,"-> ID FEM_MAILLAGE : %li",m_id_fem_maillage); fonc(ligne);
  sprintf(ligne,"-> Boite 3D VES : [%lf,%lf,%lf,%lf,%lf,%lf]", m_boite3D_ves.get_xmin(),
                                                              m_boite3D_ves.get_ymin(),
                                                              m_boite3D_ves.get_zmin(),
                                                              m_boite3D_ves.get_xmax(),
                                                              m_boite3D_ves.get_ymax(),
                                                              m_boite3D_ves.get_zmax()); fonc(ligne);
  sprintf(ligne,"-> Precision : %lf",m_precision); fonc(ligne);
  sprintf(ligne,"-> Temps geometrie : %lf",m_temps_geometrie); fonc(ligne);
  sprintf(ligne,"-> Temps materiau : %lf",m_temps_materiau); fonc(ligne);
  sprintf(ligne,"-> Temps carte : %lf",m_temps_carte); fonc(ligne);
  sprintf(ligne,"-> Temps maillage : %lf",m_temps_maillage); fonc(ligne);
  sprintf(ligne,"-> Temps fem_maillage : %lf",m_temps_fem_maillage); fonc(ligne);
  sprintf(ligne,"-> Temps etude : %lf",m_temps_etude); fonc(ligne);
  sprintf(ligne,"-> Temps calcul : %lf",m_temps_calcul); fonc(ligne);
  sprintf(ligne,"-> Temps analyse : %lf",m_temps_analyse); fonc(ligne);
  sprintf(ligne," "); fonc(ligne);
  std::map<std::string,MSTRUCT_ANALYSE*>::iterator it;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it);analyse!=NULL;analyse=get_suivante_analyse(it))
  {
    analyse->affiche_contenu(fonc);
    sprintf(ligne," "); fonc(ligne);
  }
}