mtu/src/mg_pyramide.cpp

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mg_pyramide.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------


#include "gestionversion.h"
#include "mg_pyramide.h"
#include "mg_maillage.h"

MG_PYRAMIDE::MG_PYRAMIDE(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_NOEUD *mgnoeud5,class MG_TRIANGLE* mgtriangle1,MG_TRIANGLE* mgtriangle2,MG_TRIANGLE* mgtriangle3,MG_TRIANGLE* mgtriangle4,class MG_QUADRANGLE* mgquadrangle1,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4),quadrangle1(mgquadrangle1)
{
    init_pyramide();
}

MG_PYRAMIDE::MG_PYRAMIDE(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_NOEUD *mgnoeud5,class MG_TRIANGLE* mgtriangle1,MG_TRIANGLE* mgtriangle2,MG_TRIANGLE* mgtriangle3,MG_TRIANGLE* mgtriangle4,class MG_QUADRANGLE* mgquadrangle1,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4),quadrangle1(mgquadrangle1)
{
    init_pyramide();
}


MG_PYRAMIDE::MG_PYRAMIDE(MG_PYRAMIDE& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),noeud5(mdd.noeud5),triangle1(mdd.triangle1),triangle2(mdd.triangle2),triangle3(mdd.triangle3),triangle4(mdd.triangle4),quadrangle1(mdd.quadrangle1)
{
    init_pyramide();
}


MG_PYRAMIDE* MG_PYRAMIDE::dupliquer(MG_MAILLAGE *mgmai,long decalage)
{
    MG_NOEUD* nd1=mgmai->get_mg_noeudid(noeud1->get_id()+decalage);
    MG_NOEUD* nd2=mgmai->get_mg_noeudid(noeud2->get_id()+decalage);
    MG_NOEUD* nd3=mgmai->get_mg_noeudid(noeud3->get_id()+decalage);
    MG_NOEUD* nd4=mgmai->get_mg_noeudid(noeud4->get_id()+decalage);
    MG_NOEUD* nd5=mgmai->get_mg_noeudid(noeud5->get_id()+decalage);
    MG_TRIANGLE* tri1=mgmai->get_mg_triangleid(triangle1->get_id()+decalage);
    MG_TRIANGLE* tri2=mgmai->get_mg_triangleid(triangle2->get_id()+decalage);
    MG_TRIANGLE* tri3=mgmai->get_mg_triangleid(triangle3->get_id()+decalage);
    MG_TRIANGLE* tri4=mgmai->get_mg_triangleid(triangle4->get_id()+decalage);
    MG_QUADRANGLE* quad1=mgmai->get_mg_quadrangleid(quadrangle1->get_id()+decalage);
    MG_PYRAMIDE* pyr=new MG_PYRAMIDE(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,nd5,tri1,tri2,tri3,tri4,quad1,MAGIC::ORIGINE::DUPLIQUER);
    mgmai->ajouter_mg_pyramide(pyr);
    return pyr;
}

MG_PYRAMIDE::~MG_PYRAMIDE()
{
    if (noeud1!=NULL) noeud1->get_lien_pyramide()->supprimer(this);
    if (noeud2!=NULL) noeud2->get_lien_pyramide()->supprimer(this);
    if (noeud3!=NULL) noeud3->get_lien_pyramide()->supprimer(this);
    if (noeud4!=NULL) noeud4->get_lien_pyramide()->supprimer(this);
    if (noeud5!=NULL) noeud5->get_lien_pyramide()->supprimer(this);
    unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
    noeudminid=std::min(noeudminid,noeud3->get_id());
    noeudminid=std::min(noeudminid,noeud4->get_id());
    noeudminid=std::min(noeudminid,noeud5->get_id());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_pyramide()->supprimer(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_pyramide()->supprimer(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_pyramide()->supprimer(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_pyramide()->supprimer(this);
    if (noeud5!=NULL)
        if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_pyramide()->supprimer(this);
    if (triangle1!=NULL) {
        triangle1->decrementer();
        triangle1->get_lien_pyramide()->supprimer(this);
    }
    if (triangle2!=NULL) {
        triangle2->decrementer();
        triangle2->get_lien_pyramide()->supprimer(this);
    }
    if (triangle3!=NULL) {
        triangle3->decrementer();
        triangle3->get_lien_pyramide()->supprimer(this);
    }
    if (triangle4!=NULL) {
        triangle4->decrementer();
        triangle4->get_lien_pyramide()->supprimer(this);
    }
    if (quadrangle1!=NULL) {
        quadrangle1->decrementer();
        quadrangle1->get_lien_pyramide()->supprimer(this);
    }
    if (liaison_topologique!=NULL)
      if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
}

void MG_PYRAMIDE::init_pyramide(void)
{
    if (noeud1!=NULL) noeud1->get_lien_pyramide()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_pyramide()->ajouter(this);
    if (noeud3!=NULL) noeud3->get_lien_pyramide()->ajouter(this);
    if (noeud4!=NULL) noeud4->get_lien_pyramide()->ajouter(this);
    if (noeud5!=NULL) noeud5->get_lien_pyramide()->ajouter(this);
    unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
    noeudminid=std::min(noeudminid,noeud3->get_id());
    noeudminid=std::min(noeudminid,noeud4->get_id());
    noeudminid=std::min(noeudminid,noeud5->get_id());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_pyramide()->ajouter(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_pyramide()->ajouter(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_pyramide()->ajouter(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_pyramide()->ajouter(this);
    if (noeud5!=NULL)
        if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_pyramide()->ajouter(this);
    if (triangle1!=NULL) {
        triangle1->incrementer();
        triangle1->get_lien_pyramide()->ajouter(this);
    }
    if (triangle2!=NULL) {
        triangle2->incrementer();
        triangle2->get_lien_pyramide()->ajouter(this);
    }
    if (triangle3!=NULL) {
        triangle3->incrementer();
        triangle3->get_lien_pyramide()->ajouter(this);
    }
    if (triangle4!=NULL) {
        triangle4->incrementer();
        triangle4->get_lien_pyramide()->ajouter(this);
    }
    if (quadrangle1!=NULL) {
        quadrangle1->incrementer();
        quadrangle1->get_lien_pyramide()->ajouter(this);
    }
    
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
}


int MG_PYRAMIDE::get_type_entite(void)
{
    return MAGIC::TYPE_ENTITE::IDMG_PYRAMIDE;
}


MG_NOEUD* MG_PYRAMIDE::get_noeud1(void)
{
    return noeud1;
}

MG_NOEUD* MG_PYRAMIDE::get_noeud2(void)
{
    return noeud2;
}

MG_NOEUD* MG_PYRAMIDE::get_noeud3(void)
{
    return noeud3;
}

MG_NOEUD* MG_PYRAMIDE::get_noeud4(void)
{
    return noeud4;
}
MG_NOEUD* MG_PYRAMIDE::get_noeud5(void)
{
    return noeud5;
}


MG_QUADRANGLE* MG_PYRAMIDE::get_quadrangle1(void)
{
    return quadrangle1;
}

MG_TRIANGLE* MG_PYRAMIDE::get_triangle1(void)
{
    return triangle1;
}

MG_TRIANGLE* MG_PYRAMIDE::get_triangle2(void)
{
    return triangle2;
}
MG_TRIANGLE* MG_PYRAMIDE::get_triangle3(void)
{
    return triangle3;
}

MG_TRIANGLE* MG_PYRAMIDE::get_triangle4(void)
{
    return triangle4;
}

void MG_PYRAMIDE::inverse_sens(void)
{
MG_NOEUD* noeud=noeud4;
noeud4=noeud2;
noeud2=noeud;
MG_TRIANGLE* triangle=triangle1;
triangle1=triangle4;
triangle4=triangle;
triangle=triangle2;
triangle2=triangle3;
triangle3=triangle;
}


int MG_PYRAMIDE::get_dimension(void)
{
    return 3;
}


BOITE_3D MG_PYRAMIDE::get_boite_3D(void)
{
    double coo1[3];
    double coo2[3];
    double coo3[3];
    double coo4[3];
    double coo5[3];
    coo1[0]=noeud1->get_x();
    coo1[1]=noeud1->get_y();
    coo1[2]=noeud1->get_z();
    coo2[0]=noeud2->get_x();
    coo2[1]=noeud2->get_y();
    coo2[2]=noeud2->get_z();
    coo3[0]=noeud3->get_x();
    coo3[1]=noeud3->get_y();
    coo3[2]=noeud3->get_z();
    coo4[0]=noeud4->get_x();
    coo4[1]=noeud4->get_y();
    coo4[2]=noeud4->get_z();
    coo5[0]=noeud5->get_x();
    coo5[1]=noeud5->get_y();
    coo5[2]=noeud5->get_z();
    double xmin,ymin,zmin,xmax,ymax,zmax;
    xmin=std::min(coo1[0],coo2[0]);
    xmin=std::min(xmin,coo3[0]);
    xmin=std::min(xmin,coo4[0]);
    xmin=std::min(xmin,coo5[0]);
    xmax=std::max(coo1[0],coo2[0]);
    xmax=std::max(xmax,coo3[0]);
    xmax=std::max(xmax,coo4[0]);
    xmax=std::max(xmax,coo5[0]);
    ymin=std::min(coo1[1],coo2[1]);
    ymin=std::min(ymin,coo3[1]);
    ymin=std::min(ymin,coo4[1]);
    ymin=std::min(ymin,coo5[1]);
    ymax=std::max(coo1[1],coo2[1]);
    ymax=std::max(ymax,coo3[1]);
    ymax=std::max(ymax,coo4[1]);
    ymax=std::max(ymax,coo5[1]);
    zmin=std::min(coo1[2],coo2[2]);
    zmin=std::min(zmin,coo3[2]);
    zmin=std::min(zmin,coo4[2]);
    zmin=std::min(zmin,coo5[2]);
    zmax=std::max(coo1[2],coo2[2]);
    zmax=std::max(zmax,coo3[2]);
    zmax=std::max(zmax,coo4[2]);
    zmax=std::max(zmax,coo5[2]);
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    return boite;
}
void MG_PYRAMIDE::extrapoler_solution_noeud(void)
{
    for (int j=0;j<MAX_TYPE_SOLUTION;j++)  
      {
        noeud1->change_solution(solution[j],j);
        noeud2->change_solution(solution[j],j);
        noeud3->change_solution(solution[j],j);
        noeud4->change_solution(solution[j],j);
        noeud5->change_solution(solution[j],j);
          }
}

void MG_PYRAMIDE::enregistrer(std::ostream& o,double version)
{
    if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
    if (liaison_topologique==NULL) o << "%" << get_id() << "=PYRAMIDE(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id()  << ",$" << noeud4->get_id()  << ",$" << noeud5->get_id()  << "," << origine_creation << ");" << std::endl;
    else if (liaison_topologique->get_dimension()==3)  o << "%" << get_id() << "=PYRAMIDE($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id()  << ",$" << noeud4->get_id()  << ",$" <<noeud5->get_id() <<  "," << origine_creation << ");" << std::endl;
}


Revision:	1158
Committed:	Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size:	11256 byte(s)
Log Message:	compatibilité Ubuntu 22.04 Suppression des refeences à Windows Ajout d'une banière
#	User	Rev	Content
1	francois	1158	//####//------------------------------------------------------------
2			//####//------------------------------------------------------------
3			//####// MAGiC
4			//####// Jean Christophe Cuilliere et Vincent FRANCOIS
5			//####// Departement de Genie Mecanique - UQTR
6			//####//------------------------------------------------------------
7			//####// MAGIC est un projet de recherche de l equipe ERICCA
8			//####// du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
9			//####// http://www.uqtr.ca/ericca
10			//####// http://www.uqtr.ca/
11			//####//------------------------------------------------------------
12			//####//------------------------------------------------------------
13			//####//
14			//####// mg_pyramide.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:58:53 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22	francois	1123
23
24			#include "gestionversion.h"
25			#include "mg_pyramide.h"
26			#include "mg_maillage.h"
27
28			MG_PYRAMIDE::MG_PYRAMIDE(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD mgnoeud1,class MG_NOEUD mgnoeud2,class MG_NOEUD mgnoeud3,class MG_NOEUD mgnoeud4,class MG_NOEUD mgnoeud5,class MG_TRIANGLE mgtriangle1,MG_TRIANGLE* mgtriangle2,MG_TRIANGLE* mgtriangle3,MG_TRIANGLE* mgtriangle4,class MG_QUADRANGLE* mgquadrangle1,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4),quadrangle1(mgquadrangle1)
29			{
30			init_pyramide();
31			}
32
33			MG_PYRAMIDE::MG_PYRAMIDE(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD mgnoeud1,class MG_NOEUD mgnoeud2,class MG_NOEUD mgnoeud3,class MG_NOEUD mgnoeud4,class MG_NOEUD mgnoeud5,class MG_TRIANGLE mgtriangle1,MG_TRIANGLE* mgtriangle2,MG_TRIANGLE* mgtriangle3,MG_TRIANGLE* mgtriangle4,class MG_QUADRANGLE* mgquadrangle1,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4),quadrangle1(mgquadrangle1)
34			{
35			init_pyramide();
36			}
37
38
39			MG_PYRAMIDE::MG_PYRAMIDE(MG_PYRAMIDE& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),noeud5(mdd.noeud5),triangle1(mdd.triangle1),triangle2(mdd.triangle2),triangle3(mdd.triangle3),triangle4(mdd.triangle4),quadrangle1(mdd.quadrangle1)
40			{
41			init_pyramide();
42			}
43
44
45			MG_PYRAMIDE* MG_PYRAMIDE::dupliquer(MG_MAILLAGE *mgmai,long decalage)
46			{
47			MG_NOEUD* nd1=mgmai->get_mg_noeudid(noeud1->get_id()+decalage);
48			MG_NOEUD* nd2=mgmai->get_mg_noeudid(noeud2->get_id()+decalage);
49			MG_NOEUD* nd3=mgmai->get_mg_noeudid(noeud3->get_id()+decalage);
50			MG_NOEUD* nd4=mgmai->get_mg_noeudid(noeud4->get_id()+decalage);
51			MG_NOEUD* nd5=mgmai->get_mg_noeudid(noeud5->get_id()+decalage);
52			MG_TRIANGLE* tri1=mgmai->get_mg_triangleid(triangle1->get_id()+decalage);
53			MG_TRIANGLE* tri2=mgmai->get_mg_triangleid(triangle2->get_id()+decalage);
54			MG_TRIANGLE* tri3=mgmai->get_mg_triangleid(triangle3->get_id()+decalage);
55			MG_TRIANGLE* tri4=mgmai->get_mg_triangleid(triangle4->get_id()+decalage);
56			MG_QUADRANGLE* quad1=mgmai->get_mg_quadrangleid(quadrangle1->get_id()+decalage);
57			MG_PYRAMIDE* pyr=new MG_PYRAMIDE(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,nd5,tri1,tri2,tri3,tri4,quad1,MAGIC::ORIGINE::DUPLIQUER);
58			mgmai->ajouter_mg_pyramide(pyr);
59			return pyr;
60			}
61
62			MG_PYRAMIDE::~MG_PYRAMIDE()
63			{
64			if (noeud1!=NULL) noeud1->get_lien_pyramide()->supprimer(this);
65			if (noeud2!=NULL) noeud2->get_lien_pyramide()->supprimer(this);
66			if (noeud3!=NULL) noeud3->get_lien_pyramide()->supprimer(this);
67			if (noeud4!=NULL) noeud4->get_lien_pyramide()->supprimer(this);
68			if (noeud5!=NULL) noeud5->get_lien_pyramide()->supprimer(this);
69			unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
70			noeudminid=std::min(noeudminid,noeud3->get_id());
71			noeudminid=std::min(noeudminid,noeud4->get_id());
72			noeudminid=std::min(noeudminid,noeud5->get_id());
73			if (noeud1!=NULL)
74			if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_pyramide()->supprimer(this);
75			if (noeud2!=NULL)
76			if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_pyramide()->supprimer(this);
77			if (noeud3!=NULL)
78			if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_pyramide()->supprimer(this);
79			if (noeud4!=NULL)
80			if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_pyramide()->supprimer(this);
81			if (noeud5!=NULL)
82			if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_pyramide()->supprimer(this);
83			if (triangle1!=NULL) {
84			triangle1->decrementer();
85			triangle1->get_lien_pyramide()->supprimer(this);
86			}
87			if (triangle2!=NULL) {
88			triangle2->decrementer();
89			triangle2->get_lien_pyramide()->supprimer(this);
90			}
91			if (triangle3!=NULL) {
92			triangle3->decrementer();
93			triangle3->get_lien_pyramide()->supprimer(this);
94			}
95			if (triangle4!=NULL) {
96			triangle4->decrementer();
97			triangle4->get_lien_pyramide()->supprimer(this);
98			}
99			if (quadrangle1!=NULL) {
100			quadrangle1->decrementer();
101			quadrangle1->get_lien_pyramide()->supprimer(this);
102			}
103			if (liaison_topologique!=NULL)
104			if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
105			}
106
107			void MG_PYRAMIDE::init_pyramide(void)
108			{
109			if (noeud1!=NULL) noeud1->get_lien_pyramide()->ajouter(this);
110			if (noeud2!=NULL) noeud2->get_lien_pyramide()->ajouter(this);
111			if (noeud3!=NULL) noeud3->get_lien_pyramide()->ajouter(this);
112			if (noeud4!=NULL) noeud4->get_lien_pyramide()->ajouter(this);
113			if (noeud5!=NULL) noeud5->get_lien_pyramide()->ajouter(this);
114			unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
115			noeudminid=std::min(noeudminid,noeud3->get_id());
116			noeudminid=std::min(noeudminid,noeud4->get_id());
117			noeudminid=std::min(noeudminid,noeud5->get_id());
118			if (noeud1!=NULL)
119			if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_pyramide()->ajouter(this);
120			if (noeud2!=NULL)
121			if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_pyramide()->ajouter(this);
122			if (noeud3!=NULL)
123			if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_pyramide()->ajouter(this);
124			if (noeud4!=NULL)
125			if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_pyramide()->ajouter(this);
126			if (noeud5!=NULL)
127			if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_pyramide()->ajouter(this);
128			if (triangle1!=NULL) {
129			triangle1->incrementer();
130			triangle1->get_lien_pyramide()->ajouter(this);
131			}
132			if (triangle2!=NULL) {
133			triangle2->incrementer();
134			triangle2->get_lien_pyramide()->ajouter(this);
135			}
136			if (triangle3!=NULL) {
137			triangle3->incrementer();
138			triangle3->get_lien_pyramide()->ajouter(this);
139			}
140			if (triangle4!=NULL) {
141			triangle4->incrementer();
142			triangle4->get_lien_pyramide()->ajouter(this);
143			}
144			if (quadrangle1!=NULL) {
145			quadrangle1->incrementer();
146			quadrangle1->get_lien_pyramide()->ajouter(this);
147			}
148
149			if (liaison_topologique==NULL) return;
150			if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
151			}
152
153
154			int MG_PYRAMIDE::get_type_entite(void)
155			{
156	francois	1150	return MAGIC::TYPE_ENTITE::IDMG_PYRAMIDE;
157	francois	1123	}
158
159
160
161
162
163			MG_NOEUD* MG_PYRAMIDE::get_noeud1(void)
164			{
165			return noeud1;
166			}
167
168			MG_NOEUD* MG_PYRAMIDE::get_noeud2(void)
169			{
170			return noeud2;
171			}
172
173			MG_NOEUD* MG_PYRAMIDE::get_noeud3(void)
174			{
175			return noeud3;
176			}
177
178			MG_NOEUD* MG_PYRAMIDE::get_noeud4(void)
179			{
180			return noeud4;
181			}
182			MG_NOEUD* MG_PYRAMIDE::get_noeud5(void)
183			{
184			return noeud5;
185			}
186
187
188			MG_QUADRANGLE* MG_PYRAMIDE::get_quadrangle1(void)
189			{
190			return quadrangle1;
191			}
192
193			MG_TRIANGLE* MG_PYRAMIDE::get_triangle1(void)
194			{
195			return triangle1;
196			}
197
198			MG_TRIANGLE* MG_PYRAMIDE::get_triangle2(void)
199			{
200			return triangle2;
201			}
202			MG_TRIANGLE* MG_PYRAMIDE::get_triangle3(void)
203			{
204			return triangle3;
205			}
206
207			MG_TRIANGLE* MG_PYRAMIDE::get_triangle4(void)
208			{
209			return triangle4;
210			}
211
212			void MG_PYRAMIDE::inverse_sens(void)
213			{
214			MG_NOEUD* noeud=noeud4;
215			noeud4=noeud2;
216			noeud2=noeud;
217			MG_TRIANGLE* triangle=triangle1;
218			triangle1=triangle4;
219			triangle4=triangle;
220			triangle=triangle2;
221			triangle2=triangle3;
222			triangle3=triangle;
223			}
224
225
226			int MG_PYRAMIDE::get_dimension(void)
227			{
228			return 3;
229			}
230
231
232			BOITE_3D MG_PYRAMIDE::get_boite_3D(void)
233			{
234			double coo1[3];
235			double coo2[3];
236			double coo3[3];
237			double coo4[3];
238			double coo5[3];
239			coo1[0]=noeud1->get_x();
240			coo1[1]=noeud1->get_y();
241			coo1[2]=noeud1->get_z();
242			coo2[0]=noeud2->get_x();
243			coo2[1]=noeud2->get_y();
244			coo2[2]=noeud2->get_z();
245			coo3[0]=noeud3->get_x();
246			coo3[1]=noeud3->get_y();
247			coo3[2]=noeud3->get_z();
248			coo4[0]=noeud4->get_x();
249			coo4[1]=noeud4->get_y();
250			coo4[2]=noeud4->get_z();
251			coo5[0]=noeud5->get_x();
252			coo5[1]=noeud5->get_y();
253			coo5[2]=noeud5->get_z();
254			double xmin,ymin,zmin,xmax,ymax,zmax;
255			xmin=std::min(coo1[0],coo2[0]);
256			xmin=std::min(xmin,coo3[0]);
257			xmin=std::min(xmin,coo4[0]);
258			xmin=std::min(xmin,coo5[0]);
259			xmax=std::max(coo1[0],coo2[0]);
260			xmax=std::max(xmax,coo3[0]);
261			xmax=std::max(xmax,coo4[0]);
262			xmax=std::max(xmax,coo5[0]);
263			ymin=std::min(coo1[1],coo2[1]);
264			ymin=std::min(ymin,coo3[1]);
265			ymin=std::min(ymin,coo4[1]);
266			ymin=std::min(ymin,coo5[1]);
267			ymax=std::max(coo1[1],coo2[1]);
268			ymax=std::max(ymax,coo3[1]);
269			ymax=std::max(ymax,coo4[1]);
270			ymax=std::max(ymax,coo5[1]);
271			zmin=std::min(coo1[2],coo2[2]);
272			zmin=std::min(zmin,coo3[2]);
273			zmin=std::min(zmin,coo4[2]);
274			zmin=std::min(zmin,coo5[2]);
275			zmax=std::max(coo1[2],coo2[2]);
276			zmax=std::max(zmax,coo3[2]);
277			zmax=std::max(zmax,coo4[2]);
278			zmax=std::max(zmax,coo5[2]);
279			BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
280			return boite;
281			}
282			void MG_PYRAMIDE::extrapoler_solution_noeud(void)
283			{
284			for (int j=0;j<MAX_TYPE_SOLUTION;j++)
285			{
286			noeud1->change_solution(solution[j],j);
287			noeud2->change_solution(solution[j],j);
288			noeud3->change_solution(solution[j],j);
289			noeud4->change_solution(solution[j],j);
290			noeud5->change_solution(solution[j],j);
291			}
292			}
293
294			void MG_PYRAMIDE::enregistrer(std::ostream& o,double version)
295			{
296			if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
297			if (liaison_topologique==NULL) o << "%" << get_id() << "=PYRAMIDE(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id() << ",$" << noeud4->get_id() << ",$" << noeud5->get_id() << "," << origine_creation << ");" << std::endl;
298			else if (liaison_topologique->get_dimension()==3) o << "%" << get_id() << "=PYRAMIDE($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id() << ",$" << noeud4->get_id() << ",$" <<noeud5->get_id() << "," << origine_creation << ");" << std::endl;
299			}
300
301
302
303
304
305