fct_taille_fem_solution_generateur_microstructure.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fct_taille_fem_solution_generateur_microstructure.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "fct_taille_fem_solution_generateur_microstructure.h"
#include "ot_geometrie.h"
#include "tpl_grille.h"
#include <stack>
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE(MG_GESTIONNAIRE* gesttmp,
                                                                                                     MG_GEOMETRIE* geotmp,
                                                                                                     double dg,
                                                                                                     double dg_eps_particule,
                                                                                                     int fechantillon,
                                                                                                     int nb_cellule_direction,
                                                                                                     char* nomsolution,
                                                                                                     TPL_MAP_ENTITE<MG_VOLUME*>& tpl_map_volume_matrice,
                                                                                                     long nb_couche_min,
                                                                                                     long nb_pas,
                                                                                                     double facteur_augmentation): FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(gesttmp,
                                                                                                                                                                            geotmp,
                                                                                                                                                                            dg,
                                                                                                                                                                            fechantillon,
                                                                                                                                                                            nb_cellule_direction,
                                                                                                                                                                            nb_cellule_direction,
                                                                                                                                                                            nb_cellule_direction,
                                                                                                                                                                            nomsolution)
{
  m_nb_couche_min=nb_couche_min;
  m_nb_pas=nb_pas;
  m_nb_cellule_direction=nb_cellule_direction;
  m_facteur_augmentation=facteur_augmentation;
  m_dg_eps_particule=dg_eps_particule;
  TPL_MAP_ENTITE<MG_VOLUME*>::ITERATEUR it_volume;
  for(MG_VOLUME* volume=tpl_map_volume_matrice.get_premier(it_volume);volume!=NULL;volume=tpl_map_volume_matrice.get_suivant(it_volume))
  {
    m_tpl_map_volume_matrice.ajouter(volume);
  }
  
}
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE(FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE& mdd): FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(mdd)
{
  m_nb_couche_min=mdd.m_nb_couche_min;
  m_nb_pas=mdd.m_nb_pas;
  m_nb_cellule_direction=mdd.m_nb_cellule_direction;
  m_facteur_augmentation=mdd.m_facteur_augmentation;
  m_dg_eps_particule=mdd.m_dg_eps_particule;
  TPL_MAP_ENTITE<MG_VOLUME*>::ITERATEUR it_volume;
  for(MG_VOLUME* volume=mdd.m_tpl_map_volume_matrice.get_premier(it_volume);volume!=NULL;volume=mdd.m_tpl_map_volume_matrice.get_suivant(it_volume))
  {
    m_tpl_map_volume_matrice.ajouter(volume);
  }
}
 
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::~FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE()
{
 
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::construit_particulier(void)
{
  echantillonnage();
}
 
bool FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::lien_topologique_entre(MG_ELEMENT_TOPOLOGIQUE* mg_ele_topo1,MG_ELEMENT_TOPOLOGIQUE* mg_ele_topo2)
{
  TPL_MAP_ENTITE<MG_SOMMET*> map_sommet_volume1;
  OT_GEOMETRIE::get_map_mg_sommet_sous_jacent(mg_ele_topo1,map_sommet_volume1);
  TPL_MAP_ENTITE<MG_SOMMET*> map_sommet_volume2;
  OT_GEOMETRIE::get_map_mg_sommet_sous_jacent(mg_ele_topo2,map_sommet_volume2);
  TPL_MAP_ENTITE<MG_SOMMET*>::ITERATEUR it_sommet;
  for(MG_SOMMET* sommet=map_sommet_volume1.get_premier(it_sommet);sommet!=NULL;sommet=map_sommet_volume1.get_suivant(it_sommet))
  {
    if(map_sommet_volume2.existe(sommet)) return true;
  }
  return false;
}
 
 
int FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage(void)
{
  char message[500];
  sprintf(message,"Creation des TPL_GRILLE");affiche(message);
  cree_grille();
  LISTE_FEM_ELEMENT3::iterator it_ele3;
  for(FEM_ELEMENT3* ele3=fem->get_premier_element3(it_ele3);ele3!=NULL;ele3=fem->get_suivant_element3(it_ele3))
  {
    ele3->change_solution(ecart_nodal);
  }
  BOITE_3D boite3d_microstructure;
  TPL_MAP_ENTITE<MG_FACE*> tpl_map_face_matrice;
  TPL_MAP_ENTITE<MG_VOLUME*>::ITERATEUR it_tpl_map_volume;
  MG_VOLUME* volume_matrice=m_tpl_map_volume_matrice.get_premier(it_tpl_map_volume);
  boite3d_microstructure=volume_matrice->get_boite_3D();
  for(;volume_matrice!=NULL;volume_matrice=m_tpl_map_volume_matrice.get_suivant(it_tpl_map_volume))
  {
    OT_GEOMETRIE::get_map_mg_face_sous_jacent(volume_matrice,tpl_map_face_matrice);
    boite3d_microstructure=boite3d_microstructure+volume_matrice->get_boite_3D();
  }
  TPL_GRILLE<ELE_TOPO_GRILLE*> grille_ele_topo;
  grille_ele_topo.initialiser(boite3d_microstructure.get_xmin(),boite3d_microstructure.get_ymin(),boite3d_microstructure.get_zmin(),
                              boite3d_microstructure.get_xmax(),boite3d_microstructure.get_ymax(),boite3d_microstructure.get_zmax(),
                              nx,ny,nz);
  long id_ele_topo=0;
  std::stack<ELE_TOPO_GRILLE*> stack_ele_topo_grille;
  LISTE_MG_FACE::iterator it_face;
  for(MG_FACE* face=geo->get_premier_face(it_face);face!=NULL;face=geo->get_suivant_face(it_face))
  {
    ELE_TOPO_GRILLE* ele_topo_grille=new ELE_TOPO_GRILLE(face,id_ele_topo); id_ele_topo++;
    grille_ele_topo.inserer(ele_topo_grille);
    stack_ele_topo_grille.push(ele_topo_grille);
  }
  LISTE_MG_ARETE::iterator it_arete;
  for(MG_ARETE* arete=geo->get_premier_arete(it_arete);arete!=NULL;arete=geo->get_suivant_arete(it_arete))
  {
    ELE_TOPO_GRILLE* ele_topo_grille=new ELE_TOPO_GRILLE(arete,id_ele_topo); id_ele_topo++;
    grille_ele_topo.inserer(ele_topo_grille);
    stack_ele_topo_grille.push(ele_topo_grille);
  }
  LISTE_MG_SOMMET::iterator it_sommet;
  for(MG_SOMMET* sommet=geo->get_premier_sommet(it_sommet);sommet!=NULL;sommet=geo->get_suivant_sommet(it_sommet))
  {
    ELE_TOPO_GRILLE* ele_topo_grille=new ELE_TOPO_GRILLE(sommet,id_ele_topo); id_ele_topo++;
    grille_ele_topo.inserer(ele_topo_grille);
    stack_ele_topo_grille.push(ele_topo_grille);
  }
  sprintf(message,"Echantionnage : Topo - Topo");affiche(message);
  double nb_ele_topo_grille=stack_ele_topo_grille.size();
  while(!stack_ele_topo_grille.empty())
  {
    sprintf(message,"\033[1A\033[KEchantionnage : Topo - Topo [%.2lf%%]",((nb_ele_topo_grille-stack_ele_topo_grille.size())/nb_ele_topo_grille)*100.0);affiche(message);
    ELE_TOPO_GRILLE* ele_topo_grille=stack_ele_topo_grille.top(); stack_ele_topo_grille.pop();
    TPL_MAP_ENTITE<ELE_TOPO_GRILLE*> tpl_map_topo_grille_trouve;
    BOITE_3D boite_recherche(ele_topo_grille->get_boite_3D().get_xmin()-m_nb_couche_min*ecart_nodal,
                             ele_topo_grille->get_boite_3D().get_ymin()-m_nb_couche_min*ecart_nodal,
                             ele_topo_grille->get_boite_3D().get_zmin()-m_nb_couche_min*ecart_nodal,
                             ele_topo_grille->get_boite_3D().get_xmax()+m_nb_couche_min*ecart_nodal,
                             ele_topo_grille->get_boite_3D().get_ymax()+m_nb_couche_min*ecart_nodal,
                             ele_topo_grille->get_boite_3D().get_zmax()+m_nb_couche_min*ecart_nodal);
    grille_ele_topo.rechercher(boite_recherche,tpl_map_topo_grille_trouve);
    TPL_MAP_ENTITE<ELE_TOPO_GRILLE*>::ITERATEUR it_trouve;
    for(ELE_TOPO_GRILLE* ele_topo_grille_trouve=tpl_map_topo_grille_trouve.get_premier(it_trouve);ele_topo_grille_trouve!=NULL;ele_topo_grille_trouve=tpl_map_topo_grille_trouve.get_suivant(it_trouve))
    {
      if(ele_topo_grille==ele_topo_grille_trouve) continue;
      MG_ELEMENT_TOPOLOGIQUE* ele_topo=ele_topo_grille->get_mgeletopo();
      MG_ELEMENT_TOPOLOGIQUE* ele_topo_trouve=ele_topo_grille_trouve->get_mgeletopo();
      if(OT_GEOMETRIE::get_lien_topologique(ele_topo,ele_topo_trouve)==OT_GEOMETRIE::TYPE_LIEN_TOPOLOGIQUE::AUCUN)
      {
        echantillonnage_ele_topo_ele_topo(ele_topo,ele_topo_trouve);
      }
    }
    grille_ele_topo.supprimer(ele_topo_grille);
    delete ele_topo_grille;    
  }
  sprintf(message,"\033[1A\033[KEchantionnage : Topo - Topo [%.2lf%%]",((nb_ele_topo_grille-stack_ele_topo_grille.size())/nb_ele_topo_grille)*100.0);affiche(message);
  //if(fabs(ecart_nodal-m_dg_face_particule)>0.00001)
  {
    sprintf(message,"Echantionnage : Courbure");affiche(message);
    long nb_volume=geo->get_nb_mg_volume()-m_tpl_map_volume_matrice.get_nb();
    long i_volume=0;
    LISTE_MG_VOLUME::iterator it_volume;
    for(MG_VOLUME* volume=geo->get_premier_volume(it_volume);volume!=NULL;volume=geo->get_suivant_volume(it_volume))
    {
      if(!m_tpl_map_volume_matrice.existe(volume))
      {
        sprintf(message,"\033[1A\033[KEchantionnage : Courbure [%.2lf%%]",100.0*i_volume/nb_volume);affiche(message);  
        TPL_MAP_ENTITE<MG_ARETE*> tpl_map_arete;
        TPL_MAP_ENTITE<MG_ARETE*>::ITERATEUR it_a;
        OT_GEOMETRIE::get_map_mg_arete_sous_jacent(volume,tpl_map_arete);
        for(MG_ARETE* arete=tpl_map_arete.get_premier(it_a);arete!=NULL;arete=tpl_map_arete.get_suivant(it_a))
        {
          echantillonnage_arete(arete);
        }
        TPL_MAP_ENTITE<MG_FACE*> tpl_map_face;
        TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_f;
        OT_GEOMETRIE::get_map_mg_face_sous_jacent(volume,tpl_map_face);
        for(MG_FACE* face=tpl_map_face.get_premier(it_f);face!=NULL;face=tpl_map_face.get_suivant(it_f))
        {
          echantillonnage_face(face);
        }
        i_volume++;
      }
    }
    sprintf(message,"\033[1A\033[KEchantionnage : Courbure [%.2lf%%]",100.0*i_volume/nb_volume);affiche(message);  
  }
  sprintf(message,"Lissage de la solution");affiche(message);
  lisser_solution();
  sprintf(message,"Extrapolation aux noeuds");affiche(message);
  extrapoler_solution_noeud();
  int i=0; 
  LISTE_FEM_NOEUD::iterator it;
  for (FEM_NOEUD* nod=fem->get_premier_noeud(it);nod!=NULL;nod=fem->get_suivant_noeud(it))
  {   
    double ecartt=nod->get_solution(0); 
    sol->ecrire(1.0/ecartt/ecartt,i,0,0);
    sol->ecrire(0.0,i,0,1);
    sol->ecrire(0.0,i,0,2);
    sol->ecrire(0.0,i,0,3);
    sol->ecrire(1./ecartt/ecartt,i,0,4);
    sol->ecrire(0.0,i,0,5);
    sol->ecrire(0.0,i,0,6);
    sol->ecrire(0.0,i,0,7);
    sol->ecrire(1./ecartt/ecartt,i,0,8);
    i++;
   }
   return OK;
}
 
 
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_ele_topo_ele_topo(MG_ELEMENT_TOPOLOGIQUE* eletopo1, MG_ELEMENT_TOPOLOGIQUE* eletopo2)
{
  if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
  {
    echantillonnage_face_face((MG_FACE*)eletopo1,(MG_FACE*)eletopo2);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
  {
    echantillonnage_arete_face((MG_ARETE*)eletopo1,(MG_FACE*)eletopo2);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
  {
    echantillonnage_arete_face((MG_ARETE*)eletopo2,(MG_FACE*)eletopo1);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
  {
    echantillonnage_sommet_face((MG_SOMMET*)eletopo1,(MG_FACE*)eletopo2);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
  {
    echantillonnage_sommet_face((MG_SOMMET*)eletopo2,(MG_FACE*)eletopo1);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
  {
    echantillonnage_arete_arete((MG_ARETE*)eletopo1,(MG_ARETE*)eletopo2);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
  {
    echantillonnage_sommet_arete((MG_SOMMET*)eletopo1,(MG_ARETE*)eletopo2);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
  {
    echantillonnage_sommet_arete((MG_SOMMET*)eletopo2,(MG_ARETE*)eletopo1);
  }
  else if(eletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET && eletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
  {
    echantillonnage_sommet_sommet((MG_SOMMET*)eletopo1,(MG_SOMMET*)eletopo2);
  }
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_face_face(MG_FACE* mgface1, MG_FACE* mgface2)
{
  double distance_min;
  if(OT_GEOMETRIE::get_distance_min_mg_face_mg_face(mgface1,mgface2,distance_min)==OK)
  {
    if((distance_min/ecart_nodal)>m_nb_couche_min) return;
  }  
  double umin,umax,vmin,vmax;
  OT_GEOMETRIE::get_param_face(mgface1,umin,umax,vmin,vmax);
  if(mgface1->get_surface()->est_periodique_u())
  {
    umin=0;
    umax=mgface1->get_surface()->get_periode_u();
  }
  if(mgface1->get_surface()->est_periodique_v())
  {
    vmin=0;
    vmax=mgface1->get_surface()->get_periode_v();
  }
  double pas_u = (umax-umin)/m_nb_pas;
  double pas_v = (vmax-vmin)/m_nb_pas;
  double uv1[2];
  /*for(long i=0;i<m_nb_pas;i++)
  {
    uv1[0]=umin+i*pas_u;
    for(long j=0;j<m_nb_pas;j++)
    {
      uv1[1] = vmin+j*pas_v;
      double xyz1[3];      
      mgface1->evaluer(uv1,xyz1);
      if(OT_GEOMETRIE::est_dans_mg_face(xyz1,mgface1)==FAIL) continue;*/
  std::vector<double> lstechan; 
  mgface1->get_echantillonnage(m_nb_pas*m_nb_pas,lstechan);   
  for(long i=0;i<lstechan.size()/5;i++)
  {
      double xyz1[3]={lstechan[5*i],lstechan[5*i+1],lstechan[5*i+2]};
      double xyz2[3];
      if(OT_GEOMETRIE::projection_au_plus_pres_sur_mg_face(xyz1,mgface1,mgface2,xyz2)==OK)
      {
        double dist = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
        if((dist/ecart_nodal)<m_nb_couche_min)
        {
          double En_local = dist/m_nb_couche_min;
          TPL_MAP_ENTITE<FEM_ELEMENT3*> tpl_map_ele3;
          get_fem_ele3_droite(xyz1,xyz2,tpl_map_ele3,8);
          TPL_MAP_ENTITE<FEM_ELEMENT3*>::ITERATEUR it_ele3;
          for(FEM_ELEMENT3* ele3=tpl_map_ele3.get_premier(it_ele3);ele3!=NULL;ele3=tpl_map_ele3.get_suivant(it_ele3))
          {
            if(ele3->get_solution(0)>En_local)
            {
              ele3->change_solution(En_local);
              m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
            }
          }
        }
      }
    }
  //}
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_arete_face(MG_ARETE* mgarete, MG_FACE* mgface)
{
  double distance_min;
  if(OT_GEOMETRIE::get_distance_min_mg_eletopo_mg_eletopo(mgarete,mgface,distance_min)==OK)
  {
    if((distance_min/ecart_nodal)>m_nb_couche_min) return;
  }
  double umin,umax,vmin,vmax;
  OT_GEOMETRIE::get_param_face(mgface,umin,umax,vmin,vmax);
  if(mgface->get_surface()->est_periodique_u())
  {
    umin=0;
    umax=mgface->get_surface()->get_periode_u();
  }
  if(mgface->get_surface()->est_periodique_v())
  {
    vmin=0;
    vmax=mgface->get_surface()->get_periode_v();
  }
  double pas_u = (umax-umin)/m_nb_pas;
  double pas_v = (vmax-vmin)/m_nb_pas;
  double uv1[2];
  /*for(long i=0;i<m_nb_pas;i++)
  {
    uv1[0]=umin+i*pas_u;
    for(long j=0;j<m_nb_pas;j++)
    {
      uv1[1] = vmin+j*pas_v;
      double xyz1[3];      
      mgface->evaluer(uv1,xyz1);
      if(OT_GEOMETRIE::est_dans_mg_face(xyz1,mgface)==FAIL) continue;*/
    std::vector<double> lstechan; 
    mgface->get_echantillonnage(m_nb_pas*m_nb_pas,lstechan);   
      for(long i=0;i<lstechan.size()/5;i++)
      {
      double xyz1[3]={lstechan[5*i],lstechan[5*i+1],lstechan[5*i+2]};
      double xyz2[3];
      if(OT_GEOMETRIE::projection_au_plus_pres_sur_mg_arete(xyz1,mgarete,xyz2)==OK)
      {
        double dist = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
        if((dist/ecart_nodal)<m_nb_couche_min)
        {
          double En_local = dist/m_nb_couche_min;
          TPL_MAP_ENTITE<FEM_ELEMENT3*> tpl_map_ele3;
          get_fem_ele3_droite(xyz1,xyz2,tpl_map_ele3,8);
          TPL_MAP_ENTITE<FEM_ELEMENT3*>::ITERATEUR it_ele3;
          for(FEM_ELEMENT3* ele3=tpl_map_ele3.get_premier(it_ele3);ele3!=NULL;ele3=tpl_map_ele3.get_suivant(it_ele3))
          {
            if(ele3->get_solution(0)>En_local)
            {
              ele3->change_solution(En_local);
              m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
            }
          }
        }
      }
    }
  //}
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_sommet_face(MG_SOMMET* mgsommet, MG_FACE* mgface)
{
  double distance_min;
  if(OT_GEOMETRIE::get_distance_min_mg_eletopo_mg_eletopo(mgsommet,mgface,distance_min)==OK)
  {
    if((distance_min/ecart_nodal)>m_nb_couche_min) return;
  }
  double umin,umax,vmin,vmax;
  OT_GEOMETRIE::get_param_face(mgface,umin,umax,vmin,vmax);
  if(mgface->get_surface()->est_periodique_u())
  {
    umin=0;
    umax=mgface->get_surface()->get_periode_u();
  }
  if(mgface->get_surface()->est_periodique_v())
  {
    vmin=0;
    vmax=mgface->get_surface()->get_periode_v();
  }
  double pas_u = (umax-umin)/m_nb_pas;
  double pas_v = (vmax-vmin)/m_nb_pas;
  double uv1[2];
  double xyz2[3];
  mgsommet->get_point()->evaluer(xyz2);
 /* for(long i=0;i<m_nb_pas;i++)
  {
    uv1[0]=umin+i*pas_u;
    for(long j=0;j<m_nb_pas;j++)
    {
      uv1[1] = vmin+j*pas_v;
      double xyz1[3];      
      mgface->evaluer(uv1,xyz1);
      if(OT_GEOMETRIE::est_dans_mg_face(xyz1,mgface)==FAIL) continue;*/
   std::vector<double> lstechan; 
   mgface->get_echantillonnage(m_nb_pas*m_nb_pas,lstechan);   
   for(long i=0;i<lstechan.size()/5;i++)
     {
      double xyz1[3]={lstechan[5*i],lstechan[5*i+1],lstechan[5*i+2]}; 
      double dist = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
      if((dist/ecart_nodal)<m_nb_couche_min || (mgsommet->est_pole() && (dist/ecart_nodal)<2.0))
      {
        double En_local;
        if(mgsommet->est_pole())
        {
          En_local=dist/2.5;
          if((dist/m_nb_couche_min)<En_local) En_local=dist/m_nb_couche_min;
        }
        else En_local=dist/m_nb_couche_min;
        TPL_MAP_ENTITE<FEM_ELEMENT3*> tpl_map_ele3;
        get_fem_ele3_droite(xyz1,xyz2,tpl_map_ele3,8);
        TPL_MAP_ENTITE<FEM_ELEMENT3*>::ITERATEUR it_ele3;
        for(FEM_ELEMENT3* ele3=tpl_map_ele3.get_premier(it_ele3);ele3!=NULL;ele3=tpl_map_ele3.get_suivant(it_ele3))
        {
          if(ele3->get_solution(0)>En_local)
          {
            ele3->change_solution(En_local);
            m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
          }
        }
      }      
    }
  //}
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_arete_arete(MG_ARETE* mgarete1, MG_ARETE* mgarete2)
{
  double distance_min;
  if(OT_GEOMETRIE::get_distance_min_mg_eletopo_mg_eletopo(mgarete1,mgarete2,distance_min)==OK)
  {
    if((distance_min/ecart_nodal)>m_nb_couche_min) return;
  }
  double tmin=mgarete1->get_tmin();
  double tmax=mgarete1->get_tmax();
  double pas_t = (tmax-tmin)/m_nb_pas;
  for(long i=0;i<m_nb_pas;i++)
  {
    double t = tmin+i*pas_t;
    double xyz1[3];
    mgarete1->evaluer(t,xyz1);
    double xyz2[3];
    if(OT_GEOMETRIE::projection_au_plus_pres_sur_mg_arete(xyz1,mgarete2,xyz2)==OK)
    {
      double dist = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
      if((dist/ecart_nodal)<m_nb_couche_min)
      {
        double En_local = dist/m_nb_couche_min;
        TPL_MAP_ENTITE<FEM_ELEMENT3*> tpl_map_ele3;
        get_fem_ele3_droite(xyz1,xyz2,tpl_map_ele3,8);
        TPL_MAP_ENTITE<FEM_ELEMENT3*>::ITERATEUR it_ele3;
        for(FEM_ELEMENT3* ele3=tpl_map_ele3.get_premier(it_ele3);ele3!=NULL;ele3=tpl_map_ele3.get_suivant(it_ele3))
        {
          if(ele3->get_solution(0)>En_local)
          {
            ele3->change_solution(En_local);
            m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
          }
        }
      }
    }
  }
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_sommet_arete(MG_SOMMET* mgsommet, MG_ARETE* mgarete)
{
  double distance_min;
  if(OT_GEOMETRIE::get_distance_min_mg_eletopo_mg_eletopo(mgsommet,mgarete,distance_min)==OK)
  {
    if((distance_min/ecart_nodal)>m_nb_couche_min) return;
  }
  double tmin=mgarete->get_tmin();
  double tmax=mgarete->get_tmax();
  double pas_t = (tmax-tmin)/m_nb_pas;
  double xyz2[3];
  mgsommet->get_point()->evaluer(xyz2);
  for(long i=0;i<m_nb_pas;i++)
  {
    double t = tmin+i*pas_t;
    double xyz1[3];
    mgarete->evaluer(t,xyz1);
    double dist = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
    if((dist/ecart_nodal)<m_nb_couche_min || (mgsommet->est_pole() && (dist/ecart_nodal)<2.0))
    {
      double En_local;
      if(mgsommet->est_pole())
      {
        En_local=dist/2.5;
        if((dist/m_nb_couche_min)<En_local) En_local=dist/m_nb_couche_min;
      }
      else En_local=dist/m_nb_couche_min;
      TPL_MAP_ENTITE<FEM_ELEMENT3*> tpl_map_ele3;
      get_fem_ele3_droite(xyz1,xyz2,tpl_map_ele3,8);
      TPL_MAP_ENTITE<FEM_ELEMENT3*>::ITERATEUR it_ele3;
      for(FEM_ELEMENT3* ele3=tpl_map_ele3.get_premier(it_ele3);ele3!=NULL;ele3=tpl_map_ele3.get_suivant(it_ele3))
      {
        if(ele3->get_solution(0)>En_local)
        {
          ele3->change_solution(En_local);
          m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
        }
      }
    }    
  }
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_sommet_sommet(MG_SOMMET* mgsommet1, MG_SOMMET* mgsommet2)
{
  double xyz1[3];
  double xyz2[3];
  mgsommet1->get_point()->evaluer(xyz1);
  mgsommet2->get_point()->evaluer(xyz2);
  double dist = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
  if((dist/ecart_nodal)<m_nb_couche_min || ((mgsommet1->est_pole() || mgsommet2->est_pole())&& (dist/ecart_nodal)<2.0))
    {
    double En_local;
    if(mgsommet1->est_pole() || mgsommet2->est_pole())
    {
      En_local=dist/3.0;
      if((dist/m_nb_couche_min)<En_local) En_local=dist/m_nb_couche_min;
    }
    else En_local=dist/m_nb_couche_min;
    TPL_MAP_ENTITE<FEM_ELEMENT3*> tpl_map_ele3;
    get_fem_ele3_droite(xyz1,xyz2,tpl_map_ele3,8);
    TPL_MAP_ENTITE<FEM_ELEMENT3*>::ITERATEUR it_ele3;
    for(FEM_ELEMENT3* ele3=tpl_map_ele3.get_premier(it_ele3);ele3!=NULL;ele3=tpl_map_ele3.get_suivant(it_ele3))
    {
      if(ele3->get_solution(0)>En_local)
      {
        ele3->change_solution(En_local);
        m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
      }
    }
  }
}
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_arete(MG_ARETE* mgarete)
{
  double tmin=mgarete->get_tmin();
  double tmax=mgarete->get_tmax();
    
    
  for(long i=0;i<m_nb_pas+1;i++)
  {
      double t=tmin+i*1.0/m_nb_pas*(tmax-tmin);
      double xyz[3];
      mgarete->evaluer(t,xyz);
      TPL_CELLULE_GRILLE<FEM_ELEMENT3*> *cellule=m_grille_ele3.get_cellule(xyz[0],xyz[1],xyz[2]);
      if(cellule!=NULL)
      {
        FEM_ELEMENT3* ele3 = cellule->get_entite(0);
        double cmax=mgarete->get_M(t);
        double m_dg_face_particule=sqrt(8.*m_dg_eps_particule/cmax/cmax);
        if(ele3->get_solution(0)>m_dg_face_particule)
            {
            ele3->change_solution(m_dg_face_particule);
            m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
            }
        }
  }
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::echantillonnage_face(MG_FACE* mgface)
{
  std::vector<double> lstechan; 
  mgface->get_echantillonnage(m_nb_pas*m_nb_pas,lstechan);   
  for(long i=0;i<lstechan.size()/5;i++)
  {
      double xyz[3]={lstechan[5*i],lstechan[5*i+1],lstechan[5*i+2]};
      TPL_CELLULE_GRILLE<FEM_ELEMENT3*> *cellule=m_grille_ele3.get_cellule(xyz[0],xyz[1],xyz[2]);
      if(cellule!=NULL)
      {
        FEM_ELEMENT3* ele3 = cellule->get_entite(0);
    double uv[2]={lstechan[5*i+3],lstechan[5*i+4]};
    double cmax,cmin;
    mgface->get_courbure(uv,cmax,cmin);
        double m_dg_face_particule=sqrt(8.*m_dg_eps_particule/cmax/cmax);
    if(ele3->get_solution(0)>m_dg_face_particule)
        {
          ele3->change_solution(m_dg_face_particule);
          m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3,ele3));
        }
      }
    }
  //}  
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::lisser_solution(void)
{
  LISTE_FEM_NOEUD::iterator itn;
  double xmin=1e308,ymin=1e308,zmin=1e308,xmax=-1e308,ymax=-1e308,zmax=-1e308;
  for (FEM_NOEUD* no=fem->get_premier_noeud(itn);no!=NULL;no=fem->get_suivant_noeud(itn))
  {
    if (no->get_x()>xmax) xmax=no->get_x();
    if (no->get_y()>ymax) ymax=no->get_y();
    if (no->get_z()>zmax) zmax=no->get_z();
    if (no->get_x()<xmin) xmin=no->get_x();
    if (no->get_y()<ymin) ymin=no->get_y();
    if (no->get_z()<zmin) zmin=no->get_z();
  }
  double pas = (xmax-xmin)/nx;
  double facteur_augmentation_distance = m_facteur_augmentation/nx;
 
  while(!m_map_fem_ele3.empty())
  {
    std::map<FEM_ELEMENT3*,FEM_ELEMENT3*>::iterator it_map_ele3;
    it_map_ele3 = m_map_fem_ele3.begin();
    FEM_ELEMENT3* ele3=it_map_ele3->second;
    double ecart_recherche=ecart_nodal+1000;    
    for(;it_map_ele3!=m_map_fem_ele3.end();it_map_ele3++)
    {
      if(it_map_ele3->second->get_solution(0)<ecart_recherche)
      {
        ele3=it_map_ele3->second;
        ecart_recherche=ele3->get_solution(0);
      }
    }
    double centre[3];
    ele3->get_boite_3D().get_centre(centre);
    double en_centre=ele3->get_solution(0);
    TPL_CELLULE_GRILLE<FEM_ELEMENT3*> *tab_cell[26];
    for(int i=0;i<26;i++)tab_cell[i]=NULL;
    
    tab_cell[0] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1]-pas,centre[2]-pas);
    tab_cell[1] = m_grille_ele3.get_cellule(centre[0],centre[1]-pas,centre[2]-pas);
    tab_cell[2] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1]-pas,centre[2]-pas);
    tab_cell[3] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1],centre[2]-pas);
    tab_cell[4] = m_grille_ele3.get_cellule(centre[0],centre[1],centre[2]-pas);
    tab_cell[5] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1],centre[2]-pas);
    tab_cell[6] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1]+pas,centre[2]-pas);
    tab_cell[7] = m_grille_ele3.get_cellule(centre[0],centre[1]+pas,centre[2]-pas);
    tab_cell[8] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1]+pas,centre[2]-pas);
    
    tab_cell[9] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1]-pas,centre[2]);
    tab_cell[10] = m_grille_ele3.get_cellule(centre[0],centre[1]-pas,centre[2]);
    tab_cell[11] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1]-pas,centre[2]);
    tab_cell[12] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1],centre[2]);
    tab_cell[13] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1],centre[2]);
    tab_cell[14] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1]+pas,centre[2]);
    tab_cell[15] = m_grille_ele3.get_cellule(centre[0],centre[1]+pas,centre[2]);
    tab_cell[16] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1]+pas,centre[2]);
    
    tab_cell[17] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1]-pas,centre[2]+pas);
    tab_cell[18] = m_grille_ele3.get_cellule(centre[0],centre[1]-pas,centre[2]+pas);
    tab_cell[19] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1]-pas,centre[2]+pas);
    tab_cell[20] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1],centre[2]+pas);
    tab_cell[21] = m_grille_ele3.get_cellule(centre[0],centre[1],centre[2]+pas);
    tab_cell[22] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1],centre[2]+pas);
    tab_cell[23] = m_grille_ele3.get_cellule(centre[0]-pas,centre[1]+pas,centre[2]+pas);
    tab_cell[24] = m_grille_ele3.get_cellule(centre[0],centre[1]+pas,centre[2]+pas);
    tab_cell[25] = m_grille_ele3.get_cellule(centre[0]+pas,centre[1]+pas,centre[2]+pas);
    for(int i=0;i<26;i++)
    {
      if(tab_cell[i]!=NULL)
      if(tab_cell[i]->get_nb_entite()==1)
      {
        FEM_ELEMENT3* ele3_voisin = tab_cell[i]->get_entite(0);
        if(ele3_voisin->get_solution(0)>en_centre*(1.+facteur_augmentation_distance))
        {
          ele3_voisin->change_solution(en_centre*(1.+facteur_augmentation_distance));
          m_map_fem_ele3.insert(std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*>(ele3_voisin,ele3_voisin));
        }
      }
    }
   m_map_fem_ele3.erase(ele3);
  }  
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::extrapoler_solution_noeud(void)
{
  LISTE_FEM_NOEUD::iterator it_fem_nd;
  for(FEM_NOEUD* femnd=fem->get_premier_noeud(it_fem_nd);femnd!=NULL;femnd=fem->get_suivant_noeud(it_fem_nd))
  {
    double ecart_nodal_noeud=0.0;
    TPL_LISTE_ENTITE<FEM_ELEMENT3*> *liste_fem_element3;
    liste_fem_element3 = femnd->get_lien_element3();    
    for(int i=0;i<liste_fem_element3->get_nb();i++)
    {
      ecart_nodal_noeud+=liste_fem_element3->get(i)->get_solution(0);
    }
    ecart_nodal_noeud=ecart_nodal_noeud/liste_fem_element3->get_nb();
    femnd->change_solution(ecart_nodal_noeud);
  }
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::cree_grille(void)
{
  LISTE_FEM_NOEUD::iterator itn;
  double xmin=1e308,ymin=1e308,zmin=1e308,xmax=-1e308,ymax=-1e308,zmax=-1e308;
  for (FEM_NOEUD* no=fem->get_premier_noeud(itn);no!=NULL;no=fem->get_suivant_noeud(itn))
  {
    if (no->get_x()>xmax) xmax=no->get_x();
    if (no->get_y()>ymax) ymax=no->get_y();
    if (no->get_z()>zmax) zmax=no->get_z();
    if (no->get_x()<xmin) xmin=no->get_x();
    if (no->get_y()<ymin) ymin=no->get_y();
    if (no->get_z()<zmin) zmin=no->get_z();
  }
  m_grille_ele3.initialiser(xmin,
                          ymin,
                          zmin,
                          xmax,
                          ymax,
                          zmax,
                          nx,ny,nz);
  double pas = (xmax-xmin)/nx;
  LISTE_FEM_ELEMENT3::iterator it_ele3;
  for(FEM_ELEMENT3* ele3=fem->get_premier_element3(it_ele3);ele3!=NULL;ele3=fem->get_suivant_element3(it_ele3))
  {
    double centre[3];
    ele3->get_boite_3D().get_centre(centre);
    TPL_CELLULE_GRILLE<FEM_ELEMENT3*>* cellule = m_grille_ele3.get_cellule(centre[0],centre[1],centre[2]);
    cellule->lst_entite_A.ajouter(ele3);
  }  
}
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_MICROSTRUCTURE::get_fem_ele3_droite(double* xyz1, double* xyz2, TPL_MAP_ENTITE< FEM_ELEMENT3* >& tpl_map_ele3, int nb_pas)
{
  double t=0;
  double dt=1.0/nb_pas;
  double xyz[3];
  double dxyz[3];
  dxyz[0]=xyz2[0]-xyz1[0];
  dxyz[1]=xyz2[1]-xyz1[1];
  dxyz[2]=xyz2[2]-xyz1[2];
  for(int i=0;i<nb_pas+1;i++)
  {
    xyz[0]=xyz1[0]+t*dxyz[0];
    xyz[1]=xyz1[1]+t*dxyz[1];
    xyz[2]=xyz1[2]+t*dxyz[2];
    TPL_CELLULE_GRILLE<FEM_ELEMENT3*>* cellule = m_grille_ele3.get_cellule(xyz[0],xyz[1],xyz[2]);
    if(cellule!=NULL) tpl_map_ele3.ajouter(cellule->get_entite(0));
    t+=dt;
  }
}