| 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 |
|
|
nvtri->change_etat_front(FRONT_PAS_A_INVERSE); |
| 31 |
|
|
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 |
|
|
nvtri->change_etat_front(FRONT_A_INVERSE); |
| 35 |
|
|
|
| 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 |
|
|
if (mtri->get_etat_front()==FRONT_PAS_A_INVERSE) |
| 57 |
|
|
ajouter_front_courant(NONFORCE,tri->get_noeud1(),tri->get_noeud2(),tri->get_noeud3(),mtri); |
| 58 |
|
|
else if (mtri->get_etat_front()==FRONT_A_INVERSE) |
| 59 |
|
|
ajouter_front_courant(NONFORCE,tri->get_noeud2(),tri->get_noeud1(),tri->get_noeud3(),mtri); |
| 60 |
|
|
} |
| 61 |
|
|
else mtri->change_etat_front(FRONT_NONACTIF); |
| 62 |
|
|
|
| 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 |
|
|
if (normtri*normext>0) tri->change_etat_front(FRONT_PAS_A_INVERSE); |
| 105 |
|
|
else tri->change_etat_front(FRONT_A_INVERSE); |
| 106 |
|
|
} |
