mstruct_analyse_fem_maillage.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_fem_maillage.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_analyse_fem_maillage.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"
#include "ot_mathematique.h"
#include <string.h>
MSTRUCT_ANALYSE_FEM_MAILLAGE::MSTRUCT_ANALYSE_FEM_MAILLAGE(void): MSTRUCT_ANALYSE()
{
 
}
 
MSTRUCT_ANALYSE_FEM_MAILLAGE::MSTRUCT_ANALYSE_FEM_MAILLAGE(MSTRUCT_VES* ves,
                                                           std::string identifiant,
                                                           long int id_fem_maillage,
                                                           double largeur_colonne_jacobien_2D_min,
                                                           double largeur_colonne_jacobien_2D_max,
                                                           double largeur_colonne_distortion_2D,
                                                           double largeur_colonne_jacobien_3D_min,
                                                           double largeur_colonne_jacobien_3D_max,
                                                           double largeur_colonne_distortion_3D,
                                                           std::string nom_groupe_forme,
                                                           BOITE_3D* boite_3d): MSTRUCT_ANALYSE(ves,identifiant,nom_groupe_forme,boite_3d)
{
  m_id_fem_maillage=id_fem_maillage;
  m_2D_jacobien_histogramme_min.fixe_largeur_colonne(largeur_colonne_jacobien_2D_min);
  m_2D_jacobien_histogramme_max.fixe_largeur_colonne(largeur_colonne_jacobien_2D_max);
  m_3D_jacobien_histogramme_min.fixe_largeur_colonne(largeur_colonne_jacobien_3D_min);
  m_3D_jacobien_histogramme_max.fixe_largeur_colonne(largeur_colonne_jacobien_3D_max);
  m_2D_distortion_histogramme.fixe_largeur_colonne(largeur_colonne_distortion_2D);
  m_3D_distortion_histogramme.fixe_largeur_colonne(largeur_colonne_distortion_3D);
  
  m_nb_element_2D=0.0;
  m_nb_element_3D=0.0;
  
  m_2D_jacobien_min_min=0.0;
  m_2D_jacobien_min_max=0.0;
  m_2D_jacobien_min_moyenne=0.0;
  m_2D_jacobien_min_ecart_type=0.0;
  
  m_2D_jacobien_max_min=0.0;
  m_2D_jacobien_max_max=0.0;
  m_2D_jacobien_max_moyenne=0.0;
  m_2D_jacobien_max_ecart_type=0.0;
  
  m_3D_jacobien_min_min=0.0;
  m_3D_jacobien_min_max=0.0;
  m_3D_jacobien_min_moyenne=0.0;
  m_3D_jacobien_min_ecart_type=0.0;
  
  m_3D_jacobien_max_min=0.0;
  m_3D_jacobien_max_max=0.0;
  m_3D_jacobien_max_moyenne=0.0;
  m_3D_jacobien_max_ecart_type=0.0;
  
  m_2D_distortion_min=0.0;
  m_2D_distortion_max=0.0;
  m_2D_distortion_moyenne=0.0;
  m_2D_distortion_ecart_type=0.0;
  m_3D_distortion_min=0.0;
  m_3D_distortion_max=0.0;
  m_3D_distortion_moyenne=0.0;
  m_3D_distortion_ecart_type=0.0;
  
  m_volume=0.0;
  m_fraction_volumique=0.0;
}
 
MSTRUCT_ANALYSE_FEM_MAILLAGE::MSTRUCT_ANALYSE_FEM_MAILLAGE(MSTRUCT_ANALYSE_FEM_MAILLAGE& mdd,bool cpy_data): MSTRUCT_ANALYSE(mdd,cpy_data)
{
  m_id_fem_maillage=mdd.m_id_fem_maillage;
  if(cpy_data)
  {
    m_nb_element_2D=mdd.m_nb_element_2D;
    m_nb_element_3D=mdd.m_nb_element_3D;    
    m_2D_jacobien_min_min=mdd.m_2D_jacobien_min_min;
    m_2D_jacobien_min_max=mdd.m_2D_jacobien_min_max;
    m_2D_jacobien_min_moyenne=mdd.m_2D_jacobien_min_moyenne;
    m_2D_jacobien_min_ecart_type=mdd.m_2D_jacobien_min_ecart_type;
    m_2D_jacobien_histogramme_min=OT_HISTOGRAMME(mdd.m_2D_jacobien_histogramme_min);    
    m_2D_jacobien_max_min=mdd.m_2D_jacobien_max_min;
    m_2D_jacobien_max_max=mdd.m_2D_jacobien_max_max;
    m_2D_jacobien_max_moyenne=mdd.m_2D_jacobien_max_moyenne;
    m_2D_jacobien_max_ecart_type=mdd.m_2D_jacobien_max_ecart_type;
    m_2D_jacobien_histogramme_max=OT_HISTOGRAMME(mdd.m_2D_jacobien_histogramme_max);    
    m_3D_jacobien_min_min=mdd.m_3D_jacobien_min_min;
    m_3D_jacobien_min_max=mdd.m_3D_jacobien_min_max;
    m_3D_jacobien_min_moyenne=mdd.m_3D_jacobien_min_moyenne;
    m_3D_jacobien_min_ecart_type=mdd.m_3D_jacobien_min_ecart_type;
    m_3D_jacobien_histogramme_min=OT_HISTOGRAMME(mdd.m_3D_jacobien_histogramme_min);    
    m_3D_jacobien_max_min=mdd.m_3D_jacobien_max_min;
    m_3D_jacobien_max_max=mdd.m_3D_jacobien_max_max;
    m_3D_jacobien_max_moyenne=mdd.m_3D_jacobien_max_moyenne;
    m_3D_jacobien_max_ecart_type=mdd.m_3D_jacobien_max_ecart_type;
    m_3D_jacobien_histogramme_max=OT_HISTOGRAMME(mdd.m_3D_jacobien_histogramme_max);    
    m_2D_distortion_min=mdd.m_2D_distortion_min;
    m_2D_distortion_max=mdd.m_2D_distortion_max;
    m_2D_distortion_moyenne=mdd.m_2D_distortion_moyenne;
    m_2D_distortion_ecart_type=mdd.m_2D_distortion_ecart_type;
    m_2D_distortion_histogramme=OT_HISTOGRAMME(mdd.m_2D_distortion_histogramme);
    m_3D_distortion_min=mdd.m_3D_distortion_min;
    m_3D_distortion_max=mdd.m_3D_distortion_max;
    m_3D_distortion_moyenne=mdd.m_3D_distortion_moyenne;
    m_3D_distortion_ecart_type=mdd.m_3D_distortion_ecart_type;
    m_3D_distortion_histogramme=OT_HISTOGRAMME(mdd.m_3D_distortion_histogramme);  
    m_volume=mdd.m_volume;
    m_fraction_volumique=mdd.m_fraction_volumique;
  }
  else
  {
    m_nb_element_2D=0.0;
    m_nb_element_3D=0.0;
    
    m_2D_jacobien_min_min=0.0;
    m_2D_jacobien_min_max=0.0;
    m_2D_jacobien_min_moyenne=0.0;
    m_2D_jacobien_min_ecart_type=0.0;
    m_2D_jacobien_histogramme_min=OT_HISTOGRAMME(mdd.m_2D_jacobien_histogramme_min.get_largeur_colonne());
    
    m_2D_jacobien_max_min=0.0;
    m_2D_jacobien_max_max=0.0;
    m_2D_jacobien_max_moyenne=0.0;
    m_2D_jacobien_max_ecart_type=0.0;
    m_2D_jacobien_histogramme_max=OT_HISTOGRAMME(mdd.m_2D_jacobien_histogramme_max.get_largeur_colonne());
    
    m_3D_jacobien_min_min=0.0;
    m_3D_jacobien_min_max=0.0;
    m_3D_jacobien_min_moyenne=0.0;
    m_3D_jacobien_min_ecart_type=0.0;
    m_3D_jacobien_histogramme_min=OT_HISTOGRAMME(mdd.m_3D_jacobien_histogramme_min.get_largeur_colonne());
    
    m_3D_jacobien_max_min=0.0;
    m_3D_jacobien_max_max=0.0;
    m_3D_jacobien_max_moyenne=0.0;
    m_3D_jacobien_max_ecart_type=0.0;
    m_3D_jacobien_histogramme_max=OT_HISTOGRAMME(mdd.m_3D_jacobien_histogramme_max.get_largeur_colonne());
    
    m_2D_distortion_min=0.0;
    m_2D_distortion_max=0.0;
    m_2D_distortion_moyenne=0.0;
    m_2D_distortion_ecart_type=0.0;
    m_2D_distortion_histogramme=OT_HISTOGRAMME(mdd.m_2D_distortion_histogramme.get_largeur_colonne());
    m_3D_distortion_min=0.0;
    m_3D_distortion_max=0.0;
    m_3D_distortion_moyenne=0.0;
    m_3D_distortion_ecart_type=0.0;
    m_3D_distortion_histogramme=OT_HISTOGRAMME(mdd.m_3D_distortion_histogramme.get_largeur_colonne());
  
    m_volume=0.0;
    m_fraction_volumique=0.0;
  }
  
  
}
 
MSTRUCT_ANALYSE_FEM_MAILLAGE::~MSTRUCT_ANALYSE_FEM_MAILLAGE(void)
{
 
}
 
long int MSTRUCT_ANALYSE_FEM_MAILLAGE::get_id_fem_maillage(void)
{
  return m_id_fem_maillage;
}
 
long int MSTRUCT_ANALYSE_FEM_MAILLAGE::get_nb_element_2D(void)
{
  return m_nb_element_2D;
}
 
long int MSTRUCT_ANALYSE_FEM_MAILLAGE::get_nb_element_3D(void)
{
  return m_nb_element_3D;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_moyenne_2D(void)
{
  return m_2D_jacobien_min_moyenne;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_ecart_type_2D(void)
{
  return m_2D_jacobien_min_ecart_type;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_min_2D(void)
{
  return m_2D_jacobien_min_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_max_2D(void)
{
  return m_2D_jacobien_min_max;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distribution_jacobien_min_2D(void)
{
  return &m_2D_jacobien_histogramme_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_moyenne_2D(void)
{
  return m_2D_jacobien_max_moyenne;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_ecart_type_2D(void)
{
  return m_2D_jacobien_max_ecart_type;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_min_2D(void)
{
  return m_2D_jacobien_max_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_max_2D(void)
{
  return m_2D_jacobien_max_max;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distribution_jacobien_max_2D(void)
{
  return &m_2D_jacobien_histogramme_max;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_moyenne_3D(void)
{
  return m_3D_jacobien_min_moyenne;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_ecart_type_3D(void)
{
  return m_3D_jacobien_min_ecart_type;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_min_3D(void)
{
  return m_3D_jacobien_min_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_min_max_3D(void)
{
  return m_3D_jacobien_min_max;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distribution_jacobien_min_3D(void)
{
  return &m_3D_jacobien_histogramme_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_moyenne_3D(void)
{
  return m_3D_jacobien_max_moyenne;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_ecart_type_3D(void)
{
  return m_3D_jacobien_max_ecart_type;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_min_3D(void)
{
  return m_3D_jacobien_max_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_jacobien_max_max_3D(void)
{
  return m_3D_jacobien_max_max;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distribution_jacobien_max_3D(void)
{
  return &m_3D_jacobien_histogramme_max;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_moyenne_2D(void)
{
  return m_2D_distortion_moyenne;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_ecart_type_2D(void)
{
  return m_2D_distortion_ecart_type;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_min_2D(void)
{
  return m_2D_distortion_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_max_2D(void)
{
  return m_2D_distortion_max;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distribution_distortion_2D(void)
{
  return &m_2D_distortion_histogramme;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_moyenne_3D(void)
{
  return m_3D_distortion_moyenne;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_ecart_type_3D(void)
{
  return m_3D_distortion_ecart_type;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_min_3D(void)
{
  return m_3D_distortion_min;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distortion_max_3D(void)
{
  return m_3D_distortion_max;
}
 
OT_HISTOGRAMME* MSTRUCT_ANALYSE_FEM_MAILLAGE::get_distribution_distortion_3D(void)
{
  return &m_3D_distortion_histogramme;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_volume(void)
{
  return m_volume;
}
 
double MSTRUCT_ANALYSE_FEM_MAILLAGE::get_fraction_volumique(void)
{
  return m_fraction_volumique;
}
 
long int MSTRUCT_ANALYSE_FEM_MAILLAGE::get_type(void)
{
  return MSTRUCT::TYPE_ANALYSE::MAILLAGE_FEM;
}
 
void MSTRUCT_ANALYSE_FEM_MAILLAGE::executer(void)
{
  FEM_MAILLAGE* fem=m_ves->get_fem_maillage();
  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_fem_maillage(fem,
                                    m_boite_analyse,
                                    tpl_map_forme,
                                    m_nb_element_2D,m_nb_element_3D,
                                    m_2D_jacobien_min_moyenne, m_2D_jacobien_min_ecart_type, m_2D_jacobien_min_min, m_2D_jacobien_min_max,m_2D_jacobien_histogramme_min,
                                    m_2D_jacobien_max_moyenne, m_2D_jacobien_max_ecart_type, m_2D_jacobien_max_min, m_2D_jacobien_max_max,m_2D_jacobien_histogramme_max,
                                    m_3D_jacobien_min_moyenne, m_3D_jacobien_min_ecart_type, m_3D_jacobien_min_min, m_3D_jacobien_min_max,m_3D_jacobien_histogramme_min,
                                    m_3D_jacobien_max_moyenne, m_3D_jacobien_max_ecart_type, m_3D_jacobien_max_min, m_3D_jacobien_max_max,m_3D_jacobien_histogramme_max,
                                    m_2D_distortion_moyenne,m_2D_distortion_ecart_type,m_2D_distortion_min,m_2D_distortion_max,m_2D_distortion_histogramme,
                                    m_3D_distortion_moyenne,m_3D_distortion_ecart_type,m_3D_distortion_min,m_3D_distortion_max,m_3D_distortion_histogramme,
                                    m_volume,
                                    m_fraction_volumique);
}
 
void MSTRUCT_ANALYSE_FEM_MAILLAGE::exporter(std::ofstream& ofstrm, long i, bool avec_entete, bool avec_graph, char* prefix_graph)
{
  if(avec_entete) ofstrm << "#(1) volume(2) frac_vol(3) nb_ele_2D(4) nb_ele_3D(5) jacobien_min_2D[moy(6) ±(7) min(8) max(9)] jacobien_max_2D[moy(10) ±(11) min(12) max(13)] jacobien_min_3D[moy(14) ±(15) min(16) max(17)] jacobien_max_3D[moy(18) ±(19) min(20) max(21)] distortion_2D[moy(22) ±(23) min(24) max(25)] distortion_3D[moy(26) ±(27) min(28) max(29)]" << std::endl; 
  ofstrm << i << " "
         << m_volume << " " 
         << m_fraction_volumique << " "
         << m_nb_element_2D << " "
         << m_nb_element_3D << " "
         << m_2D_jacobien_min_moyenne << " "
         << m_2D_jacobien_min_ecart_type << " "
         << m_2D_jacobien_min_min << " "
         << m_2D_jacobien_min_max << " "
         << m_2D_jacobien_max_moyenne << " "
         << m_2D_jacobien_max_ecart_type << " "
         << m_2D_jacobien_max_min << " "
         << m_2D_jacobien_max_max << " "
         << m_3D_jacobien_min_moyenne << " "
         << m_3D_jacobien_min_ecart_type << " "
         << m_3D_jacobien_min_min << " "
         << m_3D_jacobien_min_max << " "
         << m_3D_jacobien_max_moyenne << " "
         << m_3D_jacobien_max_ecart_type << " "
         << m_3D_jacobien_max_min << " "
         << m_3D_jacobien_max_max << " "
         << m_2D_distortion_moyenne << " "
         << m_2D_distortion_ecart_type << " "
         << m_2D_distortion_min << " "
         << m_2D_distortion_max << " "
         << m_3D_distortion_moyenne << " "
         << m_3D_distortion_ecart_type << " "
         << m_3D_distortion_min << " "
         << m_3D_distortion_max << std::endl;
  if(avec_graph)
  {
    char nom_fichier[500];
    char nom_plot[500];
    MAGIC_PLOT magic_plot;
    sprintf(nom_fichier,"%s/histo_%s_jacobien_min_2D.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_javcobien_min_2d(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_javcobien_min_2d.precision(16);
    of_histo_javcobien_min_2d.setf(std::ios::showpoint);
    m_2D_jacobien_histogramme_min.exporter(of_histo_javcobien_min_2d);
    of_histo_javcobien_min_2d.close();
    sprintf(nom_plot,"%s/histo_%s_jacobien_min_2D.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="Jacobien min 2D";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_jacobien_min_2D.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);
    
    sprintf(nom_fichier,"%s/histo_%s_jacobien_max_2D.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_javcobien_max_2d(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_javcobien_max_2d.precision(16);
    of_histo_javcobien_max_2d.setf(std::ios::showpoint);
    m_2D_jacobien_histogramme_max.exporter(of_histo_javcobien_max_2d);
    of_histo_javcobien_max_2d.close();
    sprintf(nom_plot,"%s/histo_%s_jacobien_max_2D.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="Jacobien min 2D";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_jacobien_max_2D.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);    
    
    sprintf(nom_fichier,"%s/histo_%s_jacobien_min_3D.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_javcobien_min_3d(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_javcobien_min_3d.precision(16);
    of_histo_javcobien_min_3d.setf(std::ios::showpoint);
    m_3D_jacobien_histogramme_min.exporter(of_histo_javcobien_min_3d);
    of_histo_javcobien_min_3d.close();
    sprintf(nom_plot,"%s/histo_%s_jacobien_min_3D.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="Jacobien min 3D";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_jacobien_min_3D.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);
    
    sprintf(nom_fichier,"%s/histo_%s_jacobien_max_3D.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_javcobien_max_3d(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_javcobien_max_3d.precision(16);
    of_histo_javcobien_max_3d.setf(std::ios::showpoint);
    m_3D_jacobien_histogramme_max.exporter(of_histo_javcobien_max_3d);
    of_histo_javcobien_max_3d.close();
    sprintf(nom_plot,"%s/histo_%s_jacobien_max_3D.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="Jacobien max 3D";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_jacobien_max_3D.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);
    
    sprintf(nom_fichier,"%s/histo_%s_distortion_2D.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_distortion_2D(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_distortion_2D.precision(16);
    of_histo_distortion_2D.setf(std::ios::showpoint);
    m_2D_distortion_histogramme.exporter(of_histo_distortion_2D);
    of_histo_distortion_2D.close();
    sprintf(nom_plot,"%s/histo_%s_distortion_2D.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="Distortion 2D";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_distortion_2D.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);
    
    sprintf(nom_fichier,"%s/histo_%s_distortion_3D.txt",prefix_graph,m_identifiant.c_str());
    std::ofstream of_histo_distortion_3D(nom_fichier,std::ios::out|std::ios::trunc);
    of_histo_distortion_3D.precision(16);
    of_histo_distortion_3D.setf(std::ios::showpoint);
    m_3D_distortion_histogramme.exporter(of_histo_distortion_3D);
    of_histo_distortion_3D.close();
    sprintf(nom_plot,"%s/histo_%s_distortion_3D.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="Distortion 3D";
    magic_plot.histogramme_axeY_nom="";
    magic_plot.histogramme_nom_fichier_eps=std::string(m_identifiant+".eps").c_str();
    sprintf(nom_fichier,"histo_%s_distortion_3D.txt",m_identifiant.c_str());
    magic_plot.plot_script_histogramme(nom_fichier,nom_plot,true);
  }
}
 
void MSTRUCT_ANALYSE_FEM_MAILLAGE::enregistrer(std::ofstream& ofstrm)
{
  long type_analyse=get_type();
  ofstrm.write((char*)&type_analyse,sizeof(long));
  MSTRUCT_ANALYSE::enregistrer(ofstrm);
  ofstrm.write((char*)&m_id_fem_maillage,sizeof(long));
  
  ofstrm.write((char*)&m_nb_element_2D,sizeof(long));
  
  ofstrm.write((char*)&m_nb_element_3D,sizeof(long));
  
  ofstrm.write((char*)&m_2D_jacobien_min_min,sizeof(double));
  ofstrm.write((char*)&m_2D_jacobien_min_max,sizeof(double));
  ofstrm.write((char*)&m_2D_jacobien_min_moyenne,sizeof(double));
  ofstrm.write((char*)&m_2D_jacobien_min_ecart_type,sizeof(double));
  m_2D_jacobien_histogramme_min.enregistrer_bin(ofstrm);
  
  ofstrm.write((char*)&m_2D_jacobien_max_min,sizeof(double));
  ofstrm.write((char*)&m_2D_jacobien_max_max,sizeof(double));
  ofstrm.write((char*)&m_2D_jacobien_max_moyenne,sizeof(double));
  ofstrm.write((char*)&m_2D_jacobien_max_ecart_type,sizeof(double));
  m_2D_jacobien_histogramme_max.enregistrer_bin(ofstrm);
  
  ofstrm.write((char*)&m_3D_jacobien_min_min,sizeof(double));
  ofstrm.write((char*)&m_3D_jacobien_min_max,sizeof(double));
  ofstrm.write((char*)&m_3D_jacobien_min_moyenne,sizeof(double));
  ofstrm.write((char*)&m_3D_jacobien_min_ecart_type,sizeof(double));
  m_3D_jacobien_histogramme_min.enregistrer_bin(ofstrm);
  
  ofstrm.write((char*)&m_3D_jacobien_max_min,sizeof(double));
  ofstrm.write((char*)&m_3D_jacobien_max_max,sizeof(double));
  ofstrm.write((char*)&m_3D_jacobien_max_moyenne,sizeof(double));
  ofstrm.write((char*)&m_3D_jacobien_max_ecart_type,sizeof(double));
  m_3D_jacobien_histogramme_max.enregistrer_bin(ofstrm);
  
  ofstrm.write((char*)&m_2D_distortion_min,sizeof(double));
  ofstrm.write((char*)&m_2D_distortion_max,sizeof(double));
  ofstrm.write((char*)&m_2D_distortion_moyenne,sizeof(double));
  ofstrm.write((char*)&m_2D_distortion_ecart_type,sizeof(double));
  m_2D_distortion_histogramme.enregistrer_bin(ofstrm);
  
  ofstrm.write((char*)&m_3D_distortion_min,sizeof(double));
  ofstrm.write((char*)&m_3D_distortion_max,sizeof(double));
  ofstrm.write((char*)&m_3D_distortion_moyenne,sizeof(double));
  ofstrm.write((char*)&m_3D_distortion_ecart_type,sizeof(double));
  m_3D_distortion_histogramme.enregistrer_bin(ofstrm);
  
  ofstrm.write((char*)&m_volume,sizeof(double));
 
  ofstrm.write((char*)&m_fraction_volumique,sizeof(double));
}
 
void MSTRUCT_ANALYSE_FEM_MAILLAGE::ouvrir(std::ifstream& ifstrm)
{
  MSTRUCT_ANALYSE::ouvrir(ifstrm);
  ifstrm.read((char*)&m_id_fem_maillage,sizeof(long));
  
  ifstrm.read((char*)&m_nb_element_2D,sizeof(long));
  
  ifstrm.read((char*)&m_nb_element_3D,sizeof(long));
  
  ifstrm.read((char*)&m_2D_jacobien_min_min,sizeof(double));
  ifstrm.read((char*)&m_2D_jacobien_min_max,sizeof(double));
  ifstrm.read((char*)&m_2D_jacobien_min_moyenne,sizeof(double));
  ifstrm.read((char*)&m_2D_jacobien_min_ecart_type,sizeof(double));
  m_2D_jacobien_histogramme_min.ouvrir_bin(ifstrm);
  
  ifstrm.read((char*)&m_2D_jacobien_max_min,sizeof(double));
  ifstrm.read((char*)&m_2D_jacobien_max_max,sizeof(double));
  ifstrm.read((char*)&m_2D_jacobien_max_moyenne,sizeof(double));
  ifstrm.read((char*)&m_2D_jacobien_max_ecart_type,sizeof(double));
  m_2D_jacobien_histogramme_max.ouvrir_bin(ifstrm);
  
  ifstrm.read((char*)&m_3D_jacobien_min_min,sizeof(double));
  ifstrm.read((char*)&m_3D_jacobien_min_max,sizeof(double));
  ifstrm.read((char*)&m_3D_jacobien_min_moyenne,sizeof(double));
  ifstrm.read((char*)&m_3D_jacobien_min_ecart_type,sizeof(double));
  m_3D_jacobien_histogramme_min.ouvrir_bin(ifstrm);
  
  ifstrm.read((char*)&m_3D_jacobien_max_min,sizeof(double));
  ifstrm.read((char*)&m_3D_jacobien_max_max,sizeof(double));
  ifstrm.read((char*)&m_3D_jacobien_max_moyenne,sizeof(double));
  ifstrm.read((char*)&m_3D_jacobien_max_ecart_type,sizeof(double));
  m_3D_jacobien_histogramme_max.ouvrir_bin(ifstrm);
  
  ifstrm.read((char*)&m_2D_distortion_min,sizeof(double));
  ifstrm.read((char*)&m_2D_distortion_max,sizeof(double));
  ifstrm.read((char*)&m_2D_distortion_moyenne,sizeof(double));
  ifstrm.read((char*)&m_2D_distortion_ecart_type,sizeof(double));
  m_2D_distortion_histogramme.ouvrir_bin(ifstrm);
  
  ifstrm.read((char*)&m_3D_distortion_min,sizeof(double));
  ifstrm.read((char*)&m_3D_distortion_max,sizeof(double));
  ifstrm.read((char*)&m_3D_distortion_moyenne,sizeof(double));
  ifstrm.read((char*)&m_3D_distortion_ecart_type,sizeof(double));
  m_3D_distortion_histogramme.ouvrir_bin(ifstrm);
  
  ifstrm.read((char*)&m_volume,sizeof(double));
 
  ifstrm.read((char*)&m_fraction_volumique,sizeof(double));
}
 
void MSTRUCT_ANALYSE_FEM_MAILLAGE::affiche_contenu(fonction_affiche* fonc)
{
  MSTRUCT_ANALYSE::affiche_contenu(fonc);
  char ligne[5000];
  sprintf(ligne,"MSTRUCT_ANALYSE_CAO");fonc(ligne);
  sprintf(ligne,"-> ID MG_MAILLAGE : %li",m_id_fem_maillage);fonc(ligne);  
  sprintf(ligne,"-> Nb element 2D    : %li",m_nb_element_2D);fonc(ligne);
  fonc(ligne);
  sprintf(ligne,"-> Nb element 3D    : %li",m_nb_element_3D);fonc(ligne);
 
  fonc(ligne);
  
  sprintf(ligne,"-> 2D Jacobien Min moyenne :    %lf",m_2D_jacobien_min_moyenne);fonc(ligne);
  sprintf(ligne,"-> 2D Jacobien Min ecart-type : %lf",m_2D_jacobien_min_ecart_type);fonc(ligne);
  sprintf(ligne,"-> 2D Jacobien Min min :        %lf",m_2D_jacobien_min_min);fonc(ligne);
  sprintf(ligne,"-> 2D Jacobien Min max :        %lf",m_2D_jacobien_min_max);fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME 2D Jacobien Min : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_2D_jacobien_histogramme_min.get_nb_colonne(),
                                                                                                                    m_2D_jacobien_histogramme_min.get_largeur_colonne(),
                                                                                                                    m_2D_jacobien_histogramme_min.get_x_min(),
                                                                                                                    m_2D_jacobien_histogramme_min.get_x_max());
  fonc(ligne);
  
  sprintf(ligne,"-> 2D Jacobien Max moyenne :    %lf",m_2D_jacobien_max_moyenne);fonc(ligne);
  sprintf(ligne,"-> 2D Jacobien Max ecart-type : %lf",m_2D_jacobien_max_ecart_type);fonc(ligne);
  sprintf(ligne,"-> 2D Jacobien Max min :        %lf",m_2D_jacobien_max_min);fonc(ligne);
  sprintf(ligne,"-> 2D Jacobien Max max :        %lf",m_2D_jacobien_max_max);fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME 2D Jacobien Min : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_2D_jacobien_histogramme_max.get_nb_colonne(),
                                                                                                                    m_2D_jacobien_histogramme_max.get_largeur_colonne(),
                                                                                                                    m_2D_jacobien_histogramme_max.get_x_min(),
                                                                                                                    m_2D_jacobien_histogramme_max.get_x_max());
  fonc(ligne);
  
  sprintf(ligne,"-> 3D Jacobien Min moyenne :    %lf",m_3D_jacobien_min_moyenne);fonc(ligne);
  sprintf(ligne,"-> 3D Jacobien Min ecart-type : %lf",m_3D_jacobien_min_ecart_type);fonc(ligne);
  sprintf(ligne,"-> 3D Jacobien Min min :        %lf",m_3D_jacobien_min_min);fonc(ligne);
  sprintf(ligne,"-> 3D Jacobien Min max :        %lf",m_3D_jacobien_min_max);fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME 3D Jacobien Min : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_3D_jacobien_histogramme_min.get_nb_colonne(),
                                                                                                                    m_3D_jacobien_histogramme_min.get_largeur_colonne(),
                                                                                                                    m_3D_jacobien_histogramme_min.get_x_min(),
                                                                                                                    m_3D_jacobien_histogramme_min.get_x_max());
  fonc(ligne);
  
  sprintf(ligne,"-> 3D Jacobien Max moyenne :    %lf",m_3D_jacobien_max_moyenne);fonc(ligne);
  sprintf(ligne,"-> 3D Jacobien Max ecart-type : %lf",m_3D_jacobien_max_ecart_type);fonc(ligne);
  sprintf(ligne,"-> 3D Jacobien Max min :        %lf",m_3D_jacobien_max_min);fonc(ligne);
  sprintf(ligne,"-> 3D Jacobien Max max :        %lf",m_3D_jacobien_max_max);fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME 3D Jacobien Min : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_3D_jacobien_histogramme_max.get_nb_colonne(),
                                                                                                                    m_3D_jacobien_histogramme_max.get_largeur_colonne(),
                                                                                                                    m_3D_jacobien_histogramme_max.get_x_min(),
                                                                                                                    m_3D_jacobien_histogramme_max.get_x_max());
  fonc(ligne);
  
  sprintf(ligne,"-> 2D Distortion moyenne :     %lf",m_2D_distortion_moyenne);fonc(ligne);
  sprintf(ligne,"-> 2D Distortion ecart-type :  %lf",m_2D_distortion_ecart_type);fonc(ligne);
  sprintf(ligne,"-> 2D Distortion min :         %lf",m_2D_distortion_min);fonc(ligne);
  sprintf(ligne,"-> 2D Distortion max :         %lf",m_2D_distortion_max);fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME 2D Jacobien Min : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_2D_distortion_histogramme.get_nb_colonne(),
                                                                                                                    m_2D_distortion_histogramme.get_largeur_colonne(),
                                                                                                                    m_2D_distortion_histogramme.get_x_min(),
                                                                                                                    m_2D_distortion_histogramme.get_x_max());
  fonc(ligne);
  
  sprintf(ligne,"-> 3D Distortion moyenne :     %lf",m_3D_distortion_moyenne);fonc(ligne);
  sprintf(ligne,"-> 3D Distortion ecart-type :  %lf",m_3D_distortion_ecart_type);fonc(ligne);
  sprintf(ligne,"-> 3D Distortion min :         %lf",m_3D_distortion_min);fonc(ligne);
  sprintf(ligne,"-> 3D Distortion max :         %lf",m_3D_distortion_max);fonc(ligne);
  sprintf(ligne,"-> OT_HISTOGRAMME 3D Jacobien Min : nb_colonne(%li), largeur_colonne(%lf), xmin(%lf), xmax(%lf)", m_3D_distortion_histogramme.get_nb_colonne(),
                                                                                                                    m_3D_distortion_histogramme.get_largeur_colonne(),
                                                                                                                    m_3D_distortion_histogramme.get_x_min(),
                                                                                                                    m_3D_distortion_histogramme.get_x_max());
  fonc(ligne);
  sprintf(ligne,"-> Volume      : %lf",m_volume); fonc(ligne);
 
  fonc(ligne);
  sprintf(ligne,"-> Fraction volumique      : %lf",m_fraction_volumique); fonc(ligne);
 
  fonc(ligne);
}