mstruct_analyse_cao.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_analyse_cao.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_analyse_cao.h"
#include "mstruct_ves.h"
#include "mstruct_outils.h"
#include "mstruct_definition.h"
#include "mg_cg_forme.h"
#include "mg_cg_groupe_forme.h"
#include "mg_cg_assemblage.h"
#include "mg_cg_modele.h"
#include "magic_plot.h"
MSTRUCT_ANALYSE_CAO::MSTRUCT_ANALYSE_CAO(void): MSTRUCT_ANALYSE()
{
}
 
MSTRUCT_ANALYSE_CAO::MSTRUCT_ANALYSE_CAO(MSTRUCT_VES* ves,
                                         std::string identifiant,
                                         double largeur_colonne_volume_forme,
                                         std::string nom_groupe_forme): MSTRUCT_ANALYSE(ves,identifiant,nom_groupe_forme)
{
  m_volume_forme_histogramme.fixe_largeur_colonne(largeur_colonne_volume_forme);
}
 
MSTRUCT_ANALYSE_CAO::MSTRUCT_ANALYSE_CAO(MSTRUCT_ANALYSE_CAO& mdd,bool cpy_data): MSTRUCT_ANALYSE(mdd,cpy_data)
{
  if(cpy_data)
  {
    m_nb_forme=mdd.m_nb_forme;    
    m_nb_volume=mdd.m_nb_volume;    
    m_volume_forme_min=mdd.m_volume_forme_min;
    m_volume_forme_max=mdd.m_volume_forme_max;
    m_volume_forme_moyenne=mdd.m_volume_forme_moyenne;
    m_volume_forme_ecart_type=mdd.m_volume_forme_ecart_type;
    m_volume_forme_histogramme=OT_HISTOGRAMME(mdd.m_volume_forme_histogramme);
  }
  else
  {
    m_volume_forme_histogramme.fixe_largeur_colonne(mdd.m_volume_forme_histogramme.get_largeur_colonne());
  }
}
 
MSTRUCT_ANALYSE_CAO::~MSTRUCT_ANALYSE_CAO(void)
{
}
 
long int MSTRUCT_ANALYSE_CAO::get_nb_forme(void)
{
  return m_nb_forme;
}
 
long int MSTRUCT_ANALYSE_CAO::get_nb_volume(void)
{
  return m_nb_volume;
}
 
double MSTRUCT_ANALYSE_CAO::get_volume_forme_min(void)
{
  return m_volume_forme_min;
}
 
double MSTRUCT_ANALYSE_CAO::get_volume_forme_max(void)
{
  return m_volume_forme_max;
}
 
double MSTRUCT_ANALYSE_CAO::get_volume_forme_moyenne(void)
{
  return m_volume_forme_moyenne;
}
 
double MSTRUCT_ANALYSE_CAO::get_volume_forme_ecart_type(void)
{
  return m_volume_forme_ecart_type;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_CAO::get_distribution_volume_forme(void)
{
  return &m_volume_forme_histogramme;
}
 
double MSTRUCT_ANALYSE_CAO::get_volume(void)
{
  return m_volume;
}
 
double MSTRUCT_ANALYSE_CAO::get_fraction_volumique(void)
{
  return m_fraction_volumique;
}
 
long int MSTRUCT_ANALYSE_CAO::get_type(void)
{
  return MSTRUCT::TYPE_ANALYSE::CAO;
}
 
void MSTRUCT_ANALYSE_CAO::executer(void)
{
  MG_CG_GROUPE_FORME* groupe_forme=NULL;
  if(m_nom_groupe_forme!="ALL") groupe_forme=m_ves->get_mgcg_modele()->get_mgcg_groupe_forme(m_nom_groupe_forme);
  TPL_MAP_ENTITE<MG_CG_FORME*> tpl_map_forme;
  if(groupe_forme!=NULL)
  {
    std::map<long,MG_CG_FORME*>::iterator it_forme;
    for(MG_CG_FORME* forme=groupe_forme->get_premiere_mgcg_forme(it_forme);forme!=NULL;forme=groupe_forme->get_suivante_mgcg_forme(it_forme))
    {
      tpl_map_forme.ajouter(forme);
    }
  }
  else
  {
    std::map<long,MG_CG_FORME*>::iterator it_forme;
    MG_CG_ASSEMBLAGE* mgcg_ass=m_ves->get_mgcg_assemblage();
    for(MG_CG_FORME*forme=mgcg_ass->get_premiere_mgcg_forme(it_forme);forme!=NULL;forme=mgcg_ass->get_suivante_mgcg_forme(it_forme))
    {
      tpl_map_forme.ajouter(forme);
    }
  }
  MSTRUCT_OUTILS::statistiques_CAO(m_ves->get_boite3d_ves(),
                                    tpl_map_forme,
                                    m_nb_forme,
                                    m_nb_volume,
                                    m_volume_forme_moyenne,
                                    m_volume_forme_ecart_type,
                                    m_volume_forme_min,
                                    m_volume_forme_max,
                                    m_volume_forme_histogramme,
                                    m_volume,
                                    m_fraction_volumique);
}
 
void MSTRUCT_ANALYSE_CAO::exporter(std::ofstream& ofstrm, long i, bool avec_entete, bool avec_graph, char* prefix_graph)
{
  if(avec_entete)
    ofstrm << "#(1) volume(1) frac_vol(2) nb_forme(3) nb_volume(4) volume_forme[moy(5) ±(6) min(7) max(8)]" << std::endl;
  ofstrm << i << " " 
         << m_volume << " "
         << m_fraction_volumique << " "
         << m_nb_forme << " "
         << m_nb_volume << " "
         << m_volume_forme_moyenne << " "
         << m_volume_forme_ecart_type << " "
         << m_volume_forme_min << " "
         << m_volume_forme_max << std::endl;
  if(avec_graph)
  {
    char nom_fichier[500];
    char nom_plot[500];
    MAGIC_PLOT magic_plot;    
    sprintf(nom_fichier,"%s/histo_%s_volume_forme.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_volume_forme(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_volume_forme.precision(16);
    of_histo_volume_forme.setf(std::ios::showpoint);
    m_volume_forme_histogramme.exporter(of_histo_volume_forme);
    of_histo_volume_forme.close();    
    sprintf(nom_plot,"%s/histo_%s_volume_forme.plt",prefix_graph,m_identifiant.c_str());
    magic_plot.histogramme_couleur=MAGIC_PLOT::COULEUR::BLEU;
    magic_plot.histogramme_titre=m_identifiant.c_str();
    magic_plot.histogramme_axeX_nom="Volume forme";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_volume_forme.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);
  }
  
}
 
void MSTRUCT_ANALYSE_CAO::enregistrer(std::ofstream& ofstrm)
{
  long type_analyse=get_type();
  ofstrm.write((char*)&type_analyse,sizeof(long));
  MSTRUCT_ANALYSE::enregistrer(ofstrm);
  ofstrm.write((char*)&m_nb_forme,sizeof(long));
  ofstrm.write((char*)&m_nb_volume,sizeof(long));
  ofstrm.write((char*)&m_volume_forme_moyenne,sizeof(double));
  ofstrm.write((char*)&m_volume_forme_ecart_type,sizeof(double));
  ofstrm.write((char*)&m_volume_forme_min,sizeof(double));
  ofstrm.write((char*)&m_volume_forme_max,sizeof(double));
  m_volume_forme_histogramme.enregistrer_bin(ofstrm);
  ofstrm.write((char*)&m_volume,sizeof(double));
  ofstrm.write((char*)&m_fraction_volumique,sizeof(double));
}
 
void MSTRUCT_ANALYSE_CAO::ouvrir(std::ifstream& ifstrm)
{
  MSTRUCT_ANALYSE::ouvrir(ifstrm);
  ifstrm.read((char*)&m_nb_forme,sizeof(long));
  ifstrm.read((char*)&m_nb_volume,sizeof(long));
  ifstrm.read((char*)&m_volume_forme_moyenne,sizeof(double));
  ifstrm.read((char*)&m_volume_forme_ecart_type,sizeof(double));
  ifstrm.read((char*)&m_volume_forme_min,sizeof(double));
  ifstrm.read((char*)&m_volume_forme_max,sizeof(double));
  m_volume_forme_histogramme.ouvrir_bin(ifstrm);
  ifstrm.read((char*)&m_volume,sizeof(double));
  ifstrm.read((char*)&m_fraction_volumique,sizeof(double));
}
 
void MSTRUCT_ANALYSE_CAO::affiche_contenu(fonction_affiche* fonc)
{
  MSTRUCT_ANALYSE::affiche_contenu(fonc);
  char ligne[5000];
  sprintf(ligne,"MSTRUCT_ANALYSE_CAO");fonc(ligne);
  sprintf(ligne,"-> Nb forme    : %li",m_nb_forme); fonc(ligne);
  sprintf(ligne,"-> Nb volume    : %li",m_nb_volume); fonc(ligne);
  sprintf(ligne,"-> Volume      : %lf",m_volume); fonc(ligne);
  sprintf(ligne,"-> Fraction volumique      : %lf",m_fraction_volumique); fonc(ligne);
  sprintf(ligne,"-> Moyenne volume forme    : %lf",m_volume_forme_moyenne); fonc(ligne);
  sprintf(ligne,"-> Ecart-type volume forme : %lf",m_volume_forme_ecart_type); fonc(ligne);
  sprintf(ligne,"-> Min volume forme        : %lf",m_volume_forme_min); fonc(ligne);
  sprintf(ligne,"-> Max volume forme        : %lf",m_volume_forme_max); fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME volume forme : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_volume_forme_histogramme.get_nb_colonne(),
                                                                                                                m_volume_forme_histogramme.get_largeur_colonne(),
                                                                                                                m_volume_forme_histogramme.get_x_min(),
                                                                                                                m_volume_forme_histogramme.get_x_max());
  fonc(ligne);
}