mailleur3d_recuperation.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mailleur3d_recuperation.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "mailleur3d_recuperation.h"
#include "fct_taille.h"
#include "mg_volume.h"
#include "m3d_triangle.h"
 
 
 
 
 
MAILLEUR3D_RECUPERATION::MAILLEUR3D_RECUPERATION(MG_MAILLAGE* mgmai,MG_VOLUME* mgvol,FCT_TAILLE* fct_taille,bool save):MAILLEUR3D(mgmai,mgmai->get_mg_geometrie(),fct_taille,save,mgvol)
{
 
}
 
 
MAILLEUR3D_RECUPERATION::~MAILLEUR3D_RECUPERATION()
{
}
 
 
 
void MAILLEUR3D_RECUPERATION::cree_front(MG_VOLUME* mgvol) 
{
LISTE_MG_TRIANGLE::iterator it;
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
  {
  M3D_TRIANGLE* nvtri=(M3D_TRIANGLE*)tri;  
  nvtri->change_nouveau_numero(0);
  nvtri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_PAS_A_INVERSE);
  if (nvtri->get_lien_topologie()!=NULL)
    if (nvtri->get_lien_topologie()->get_dimension()==2)
        if (((MG_FACE*)nvtri->get_lien_topologie())->get_mg_coface(0)->get_orientation()==-1)
             nvtri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_A_INVERSE);
      
    
  }
  
  
LISTE_MG_TETRA::iterator ittet;
for (MG_TETRA* tet=mg_maillage->get_premier_tetra(ittet);tet!=NULL;tet=mg_maillage->get_suivant_tetra(ittet))
      {
      recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle1(),0);
      recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle2(),1);
      recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle3(),2);
      recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle4(),3);
      }
  
 
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
      {
        M3D_TRIANGLE* mtri=(M3D_TRIANGLE*)tri;
        //if ((mtri->get_nouveau_numero()==0) || ((mtri->get_nouveau_numero()==1) && (mtri->get_lien_topologie()->get_dimension()==3) )) 
        if ((mtri->get_nouveau_numero()==0) || ((mtri->get_nouveau_numero()==1) && (mtri->get_frontiere()!=1) )) 
          {
          if (mtri->get_etat_front()==MAGIC::MAILLEURFRONTALETAT::FRONT_PAS_A_INVERSE)
                    ajouter_front_courant(MAGIC::MAILLEURFRONTALETAT::NONFORCE,tri->get_noeud1(),tri->get_noeud2(),tri->get_noeud3(),mtri);  
          else if (mtri->get_etat_front()==MAGIC::MAILLEURFRONTALETAT::FRONT_A_INVERSE)
                    ajouter_front_courant(MAGIC::MAILLEURFRONTALETAT::NONFORCE,tri->get_noeud2(),tri->get_noeud1(),tri->get_noeud3(),mtri);  
          }
        else mtri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_NONACTIF);
          
      }
  
}
 
 
void MAILLEUR3D_RECUPERATION::recherche_front_triangle(MG_TETRA* tet,M3D_TRIANGLE* tri,int num)
{
tri->change_nouveau_numero(tri->get_nouveau_numero()+1);
MG_NOEUD* no1;
MG_NOEUD* no2;
MG_NOEUD* no3;
if (num==0)
    {
    no1=tet->get_noeud1();  
    no2=tet->get_noeud3();  
    no3=tet->get_noeud2();  
    }
if (num==1)
    {
    no1=tet->get_noeud1();  
    no2=tet->get_noeud2();  
    no3=tet->get_noeud4();  
    }
if (num==2)
    {
    no1=tet->get_noeud2();  
    no2=tet->get_noeud3();  
    no3=tet->get_noeud4();  
    }
if (num==3)
    {
    no1=tet->get_noeud1();  
    no2=tet->get_noeud4();  
    no3=tet->get_noeud3();  
    }
OT_VECTEUR_3D vec1(no1->get_coord(),no2->get_coord());
OT_VECTEUR_3D vec2(no1->get_coord(),no3->get_coord());
OT_VECTEUR_3D normext=vec1&vec2;
OT_VECTEUR_3D vec3(tri->get_noeud1()->get_coord(),tri->get_noeud2()->get_coord());
OT_VECTEUR_3D vec4(tri->get_noeud1()->get_coord(),tri->get_noeud3()->get_coord());
OT_VECTEUR_3D normtri=vec3&vec4;
if (normtri*normext>0) tri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_PAS_A_INVERSE);
else tri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_A_INVERSE);
}