geometrie/src/mg_tetra.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       mg_tetra.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include "mg_tetra.h"
#include "mg_maillage.h"

MG_TETRA::MG_TETRA(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_TRIANGLE* mgtriangle1,class MG_TRIANGLE* mgtriangle2,class MG_TRIANGLE* mgtriangle3,class MG_TRIANGLE* mgtriangle4,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4)
{
    init_tetra();
}

MG_TETRA::MG_TETRA(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_TRIANGLE* mgtriangle1,class MG_TRIANGLE* mgtriangle2,class MG_TRIANGLE* mgtriangle3,class MG_TRIANGLE* mgtriangle4,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4)
{
    init_tetra();
}


MG_TETRA::MG_TETRA(MG_TETRA& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),triangle1(mdd.triangle1),triangle2(mdd.triangle2),triangle3(mdd.triangle3),triangle4(mdd.triangle4)
{
    init_tetra();
}

MG_TETRA* MG_TETRA::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_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_TETRA* tet=new MG_TETRA(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,tri1,tri2,tri3,tri4,MAGIC::ORIGINE::DUPLIQUER);
    mgmai->ajouter_mg_tetra(tet);
    return tet;
}

MG_TETRA::~MG_TETRA()
{
    if (noeud1!=NULL) noeud1->get_lien_tetra()->supprimer(this);
    if (noeud2!=NULL) noeud2->get_lien_tetra()->supprimer(this);
    if (noeud3!=NULL) noeud3->get_lien_tetra()->supprimer(this);
    if (noeud4!=NULL) noeud4->get_lien_tetra()->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());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_tetra()->supprimer(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_tetra()->supprimer(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_tetra()->supprimer(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_tetra()->supprimer(this);
    if (triangle1!=NULL) {
        triangle1->decrementer();
        triangle1->get_lien_tetra()->supprimer(this);
    }
    if (triangle2!=NULL) {
        triangle2->decrementer();
        triangle2->get_lien_tetra()->supprimer(this);
    }
    if (triangle3!=NULL) {
        triangle3->decrementer();
        triangle3->get_lien_tetra()->supprimer(this);
    }
    if (triangle4!=NULL) {
        triangle4->decrementer();
        triangle4->get_lien_tetra()->supprimer(this);
    }
    if (liaison_topologique!=NULL) 
        if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
}

void MG_TETRA::init_tetra(void)
{
    if (noeud1!=NULL) noeud1->get_lien_tetra()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_tetra()->ajouter(this);
    if (noeud3!=NULL) noeud3->get_lien_tetra()->ajouter(this);
    if (noeud4!=NULL) noeud4->get_lien_tetra()->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());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_tetra()->ajouter(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_tetra()->ajouter(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_tetra()->ajouter(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_tetra()->ajouter(this);
    if (triangle1!=NULL) {
        triangle1->incrementer();
        triangle1->get_lien_tetra()->ajouter(this);
    }
    if (triangle2!=NULL) {
        triangle2->incrementer();
        triangle2->get_lien_tetra()->ajouter(this);
    }
    if (triangle3!=NULL) {
        triangle3->incrementer();
        triangle3->get_lien_tetra()->ajouter(this);
    }
    if (triangle4!=NULL) {
        triangle4->incrementer();
        triangle4->get_lien_tetra()->ajouter(this);
    }
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
}


int MG_TETRA::get_type_entite(void)
{
    return IDMG_TETRA;
}


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

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

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

MG_NOEUD* MG_TETRA::get_noeud4(void)
{
    return noeud4;
}

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

MG_TRIANGLE* MG_TETRA::get_triangle2(void)
{
    return triangle2;
}

MG_TRIANGLE* MG_TETRA::get_triangle3(void)
{
    return triangle3;
}

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

void MG_TETRA::inverse_sens(void)
{
MG_NOEUD* noeud=noeud4;
noeud4=noeud3;
noeud3=noeud;
MG_TRIANGLE* triangle=triangle1;
triangle1=triangle2;
triangle2=triangle;
}


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


BOITE_3D MG_TETRA::get_boite_3D(void)
{
    double coo1[3];
    double coo2[3];
    double coo3[3];
    double coo4[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();
    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]);
    xmax=std::max(coo1[0],coo2[0]);
    xmax=std::max(xmax,coo3[0]);
    xmax=std::max(xmax,coo4[0]);
    ymin=std::min(coo1[1],coo2[1]);
    ymin=std::min(ymin,coo3[1]);
    ymin=std::min(ymin,coo4[1]);
    ymax=std::max(coo1[1],coo2[1]);
    ymax=std::max(ymax,coo3[1]);
    ymax=std::max(ymax,coo4[1]);
    zmin=std::min(coo1[2],coo2[2]);
    zmin=std::min(zmin,coo3[2]);
    zmin=std::min(zmin,coo4[2]);
    zmax=std::max(coo1[2],coo2[2]);
    zmax=std::max(zmax,coo3[2]);
    zmax=std::max(zmax,coo4[2]);
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    return boite;
}
void MG_TETRA::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);
      }
}

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


Revision:	791
Committed:	Fri Mar 18 21:23:00 2016 UTC (9 years, 1 month ago) by francois
File size:	9291 byte(s)
Log Message:	Transferts de la definition des origines dans le namespace MAGiC
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�re et Vincent FRANCOIS
5			// D�partement de G�nie M�canique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�partement
8			// de g�nie m�canique de l'Universit� du Qu�bec �
9			// Trois Rivi�res
10			// Les librairies ne peuvent �tre utilis�es sans l'accord
11			// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// mg_tetra.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include "mg_tetra.h"
27			#include "mg_maillage.h"
28
29			MG_TETRA::MG_TETRA(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_TRIANGLE* mgtriangle1,class MG_TRIANGLE* mgtriangle2,class MG_TRIANGLE* mgtriangle3,class MG_TRIANGLE* mgtriangle4,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4)
30			{
31			init_tetra();
32			}
33
34			MG_TETRA::MG_TETRA(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD mgnoeud1,class MG_NOEUD mgnoeud2,class MG_NOEUD mgnoeud3,class MG_NOEUD mgnoeud4,class MG_TRIANGLE* mgtriangle1,class MG_TRIANGLE* mgtriangle2,class MG_TRIANGLE* mgtriangle3,class MG_TRIANGLE* mgtriangle4,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4)
35			{
36			init_tetra();
37			}
38
39
40			MG_TETRA::MG_TETRA(MG_TETRA& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),triangle1(mdd.triangle1),triangle2(mdd.triangle2),triangle3(mdd.triangle3),triangle4(mdd.triangle4)
41			{
42			init_tetra();
43			}
44
45			MG_TETRA* MG_TETRA::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_TRIANGLE* tri1=mgmai->get_mg_triangleid(triangle1->get_id()+decalage);
52			MG_TRIANGLE* tri2=mgmai->get_mg_triangleid(triangle2->get_id()+decalage);
53			MG_TRIANGLE* tri3=mgmai->get_mg_triangleid(triangle3->get_id()+decalage);
54			MG_TRIANGLE* tri4=mgmai->get_mg_triangleid(triangle4->get_id()+decalage);
55	francois	791	MG_TETRA* tet=new MG_TETRA(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,tri1,tri2,tri3,tri4,MAGIC::ORIGINE::DUPLIQUER);
56	francois	283	mgmai->ajouter_mg_tetra(tet);
57			return tet;
58			}
59
60			MG_TETRA::~MG_TETRA()
61			{
62			if (noeud1!=NULL) noeud1->get_lien_tetra()->supprimer(this);
63			if (noeud2!=NULL) noeud2->get_lien_tetra()->supprimer(this);
64			if (noeud3!=NULL) noeud3->get_lien_tetra()->supprimer(this);
65			if (noeud4!=NULL) noeud4->get_lien_tetra()->supprimer(this);
66			unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
67			noeudminid=std::min(noeudminid,noeud3->get_id());
68			noeudminid=std::min(noeudminid,noeud4->get_id());
69			if (noeud1!=NULL)
70			if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_tetra()->supprimer(this);
71			if (noeud2!=NULL)
72			if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_tetra()->supprimer(this);
73			if (noeud3!=NULL)
74			if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_tetra()->supprimer(this);
75			if (noeud4!=NULL)
76			if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_tetra()->supprimer(this);
77			if (triangle1!=NULL) {
78			triangle1->decrementer();
79			triangle1->get_lien_tetra()->supprimer(this);
80			}
81			if (triangle2!=NULL) {
82			triangle2->decrementer();
83			triangle2->get_lien_tetra()->supprimer(this);
84			}
85			if (triangle3!=NULL) {
86			triangle3->decrementer();
87			triangle3->get_lien_tetra()->supprimer(this);
88			}
89			if (triangle4!=NULL) {
90			triangle4->decrementer();
91			triangle4->get_lien_tetra()->supprimer(this);
92			}
93	francois	663	if (liaison_topologique!=NULL)
94			if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
95	francois	283	}
96
97			void MG_TETRA::init_tetra(void)
98			{
99			if (noeud1!=NULL) noeud1->get_lien_tetra()->ajouter(this);
100			if (noeud2!=NULL) noeud2->get_lien_tetra()->ajouter(this);
101			if (noeud3!=NULL) noeud3->get_lien_tetra()->ajouter(this);
102			if (noeud4!=NULL) noeud4->get_lien_tetra()->ajouter(this);
103			unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
104			noeudminid=std::min(noeudminid,noeud3->get_id());
105			noeudminid=std::min(noeudminid,noeud4->get_id());
106			if (noeud1!=NULL)
107			if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_tetra()->ajouter(this);
108			if (noeud2!=NULL)
109			if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_tetra()->ajouter(this);
110			if (noeud3!=NULL)
111			if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_tetra()->ajouter(this);
112			if (noeud4!=NULL)
113			if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_tetra()->ajouter(this);
114			if (triangle1!=NULL) {
115			triangle1->incrementer();
116			triangle1->get_lien_tetra()->ajouter(this);
117			}
118			if (triangle2!=NULL) {
119			triangle2->incrementer();
120			triangle2->get_lien_tetra()->ajouter(this);
121			}
122			if (triangle3!=NULL) {
123			triangle3->incrementer();
124			triangle3->get_lien_tetra()->ajouter(this);
125			}
126			if (triangle4!=NULL) {
127			triangle4->incrementer();
128			triangle4->get_lien_tetra()->ajouter(this);
129			}
130			if (liaison_topologique==NULL) return;
131			if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
132			}
133
134
135			int MG_TETRA::get_type_entite(void)
136			{
137			return IDMG_TETRA;
138			}
139
140
141
142
143
144			MG_NOEUD* MG_TETRA::get_noeud1(void)
145			{
146			return noeud1;
147			}
148
149			MG_NOEUD* MG_TETRA::get_noeud2(void)
150			{
151			return noeud2;
152			}
153
154			MG_NOEUD* MG_TETRA::get_noeud3(void)
155			{
156			return noeud3;
157			}
158
159			MG_NOEUD* MG_TETRA::get_noeud4(void)
160			{
161			return noeud4;
162			}
163
164			MG_TRIANGLE* MG_TETRA::get_triangle1(void)
165			{
166			return triangle1;
167			}
168
169			MG_TRIANGLE* MG_TETRA::get_triangle2(void)
170			{
171			return triangle2;
172			}
173
174			MG_TRIANGLE* MG_TETRA::get_triangle3(void)
175			{
176			return triangle3;
177			}
178
179			MG_TRIANGLE* MG_TETRA::get_triangle4(void)
180			{
181			return triangle4;
182			}
183
184	francois	702	void MG_TETRA::inverse_sens(void)
185			{
186			MG_NOEUD* noeud=noeud4;
187			noeud4=noeud3;
188			noeud3=noeud;
189			MG_TRIANGLE* triangle=triangle1;
190			triangle1=triangle2;
191			triangle2=triangle;
192			}
193	francois	283
194	francois	702
195	francois	283	int MG_TETRA::get_dimension(void)
196			{
197			return 3;
198			}
199
200
201			BOITE_3D MG_TETRA::get_boite_3D(void)
202			{
203			double coo1[3];
204			double coo2[3];
205			double coo3[3];
206			double coo4[3];
207			coo1[0]=noeud1->get_x();
208			coo1[1]=noeud1->get_y();
209			coo1[2]=noeud1->get_z();
210			coo2[0]=noeud2->get_x();
211			coo2[1]=noeud2->get_y();
212			coo2[2]=noeud2->get_z();
213			coo3[0]=noeud3->get_x();
214			coo3[1]=noeud3->get_y();
215			coo3[2]=noeud3->get_z();
216			coo4[0]=noeud4->get_x();
217			coo4[1]=noeud4->get_y();
218			coo4[2]=noeud4->get_z();
219			double xmin,ymin,zmin,xmax,ymax,zmax;
220			xmin=std::min(coo1[0],coo2[0]);
221			xmin=std::min(xmin,coo3[0]);
222			xmin=std::min(xmin,coo4[0]);
223			xmax=std::max(coo1[0],coo2[0]);
224			xmax=std::max(xmax,coo3[0]);
225			xmax=std::max(xmax,coo4[0]);
226			ymin=std::min(coo1[1],coo2[1]);
227			ymin=std::min(ymin,coo3[1]);
228			ymin=std::min(ymin,coo4[1]);
229			ymax=std::max(coo1[1],coo2[1]);
230			ymax=std::max(ymax,coo3[1]);
231			ymax=std::max(ymax,coo4[1]);
232			zmin=std::min(coo1[2],coo2[2]);
233			zmin=std::min(zmin,coo3[2]);
234			zmin=std::min(zmin,coo4[2]);
235			zmax=std::max(coo1[2],coo2[2]);
236			zmax=std::max(zmax,coo3[2]);
237			zmax=std::max(zmax,coo4[2]);
238			BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
239			return boite;
240			}
241			void MG_TETRA::extrapoler_solution_noeud(void)
242			{
243	francois	377	for (int j=0;j<MAX_TYPE_SOLUTION;j++)
244			{
245			noeud1->change_solution(solution[j],j);
246			noeud2->change_solution(solution[j],j);
247			noeud3->change_solution(solution[j],j);
248			noeud4->change_solution(solution[j],j);
249			}
250	francois	283	}
251
252	francois	763	void MG_TETRA::enregistrer(std::ostream& o,double version)
253	francois	283	{
254			if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
255			if (liaison_topologique==NULL) o << "%" << get_id() << "=TETRAEDRE(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id() << ",$" << noeud4->get_id() << "," << origine_creation << ");" << std::endl;
256			else if (liaison_topologique->get_dimension()==3) o << "%" << get_id() << "=TETRAEDRE($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id() << ",$" << noeud4->get_id() << "," << origine_creation << ");" << std::endl;
257			}
258
259
260
261
262
263