mailleur/src/mailleur3d_recuperation.cpp

#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);
}
Revision:	1150
Committed:	Tue May 28 12:57:28 2024 UTC (11 months, 2 weeks ago) by francois
Original Path:	*magic/lib/mailleur_auto/src/mailleur3d_recuperation.cpp*
File size:	3585 byte(s)
Log Message:	deplacement des contantes dans le namespace MAGIC
#	User	Rev	Content
1	francois	558	#include "gestionversion.h"
2			#include "mailleur3d_recuperation.h"
3			#include "fct_taille.h"
4			#include "mg_volume.h"
5			#include "m3d_triangle.h"
6
7
8
9
10
11			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)
12			{
13
14			}
15
16
17			MAILLEUR3D_RECUPERATION::~MAILLEUR3D_RECUPERATION()
18			{
19			}
20
21
22
23			void MAILLEUR3D_RECUPERATION::cree_front(MG_VOLUME* mgvol)
24			{
25			LISTE_MG_TRIANGLE::iterator it;
26			for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
27			{
28			M3D_TRIANGLE* nvtri=(M3D_TRIANGLE*)tri;
29			nvtri->change_nouveau_numero(0);
30	francois	1150	nvtri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_PAS_A_INVERSE);
31	francois	558	if (nvtri->get_lien_topologie()!=NULL)
32			if (nvtri->get_lien_topologie()->get_dimension()==2)
33			if (((MG_FACE*)nvtri->get_lien_topologie())->get_mg_coface(0)->get_orientation()==-1)
34	francois	1150	nvtri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_A_INVERSE);
35	francois	558
36
37			}
38
39
40			LISTE_MG_TETRA::iterator ittet;
41			for (MG_TETRA* tet=mg_maillage->get_premier_tetra(ittet);tet!=NULL;tet=mg_maillage->get_suivant_tetra(ittet))
42			{
43			recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle1(),0);
44			recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle2(),1);
45			recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle3(),2);
46			recherche_front_triangle(tet,(M3D_TRIANGLE*)tet->get_triangle4(),3);
47			}
48
49
50			for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
51			{
52			M3D_TRIANGLE* mtri=(M3D_TRIANGLE*)tri;
53			//if ((mtri->get_nouveau_numero()==0) \|\| ((mtri->get_nouveau_numero()==1) && (mtri->get_lien_topologie()->get_dimension()==3) ))
54			if ((mtri->get_nouveau_numero()==0) \|\| ((mtri->get_nouveau_numero()==1) && (mtri->get_frontiere()!=1) ))
55			{
56	francois	1150	if (mtri->get_etat_front()==MAGIC::MAILLEURFRONTALETAT::FRONT_PAS_A_INVERSE)
57			ajouter_front_courant(MAGIC::MAILLEURFRONTALETAT::NONFORCE,tri->get_noeud1(),tri->get_noeud2(),tri->get_noeud3(),mtri);
58			else if (mtri->get_etat_front()==MAGIC::MAILLEURFRONTALETAT::FRONT_A_INVERSE)
59			ajouter_front_courant(MAGIC::MAILLEURFRONTALETAT::NONFORCE,tri->get_noeud2(),tri->get_noeud1(),tri->get_noeud3(),mtri);
60	francois	558	}
61	francois	1150	else mtri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_NONACTIF);
62	francois	558
63			}
64
65			}
66
67
68			void MAILLEUR3D_RECUPERATION::recherche_front_triangle(MG_TETRA* tet,M3D_TRIANGLE* tri,int num)
69			{
70			tri->change_nouveau_numero(tri->get_nouveau_numero()+1);
71			MG_NOEUD* no1;
72			MG_NOEUD* no2;
73			MG_NOEUD* no3;
74			if (num==0)
75			{
76			no1=tet->get_noeud1();
77			no2=tet->get_noeud3();
78			no3=tet->get_noeud2();
79			}
80			if (num==1)
81			{
82			no1=tet->get_noeud1();
83			no2=tet->get_noeud2();
84			no3=tet->get_noeud4();
85			}
86			if (num==2)
87			{
88			no1=tet->get_noeud2();
89			no2=tet->get_noeud3();
90			no3=tet->get_noeud4();
91			}
92			if (num==3)
93			{
94			no1=tet->get_noeud1();
95			no2=tet->get_noeud4();
96			no3=tet->get_noeud3();
97			}
98			OT_VECTEUR_3D vec1(no1->get_coord(),no2->get_coord());
99			OT_VECTEUR_3D vec2(no1->get_coord(),no3->get_coord());
100			OT_VECTEUR_3D normext=vec1&vec2;
101			OT_VECTEUR_3D vec3(tri->get_noeud1()->get_coord(),tri->get_noeud2()->get_coord());
102			OT_VECTEUR_3D vec4(tri->get_noeud1()->get_coord(),tri->get_noeud3()->get_coord());
103			OT_VECTEUR_3D normtri=vec3&vec4;
104	francois	1150	if (normtri*normext>0) tri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_PAS_A_INVERSE);
105			else tri->change_etat_front(MAGIC::MAILLEURFRONTALETAT::FRONT_A_INVERSE);
106	francois	558	}