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


#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 MAGIC::TYPE_ENTITE::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:	1158
Committed:	Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size:	9370 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_tetra.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:58:53 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22	francois	283
23
24			#include "gestionversion.h"
25			#include "mg_tetra.h"
26			#include "mg_maillage.h"
27
28			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)
29			{
30			init_tetra();
31			}
32
33			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)
34			{
35			init_tetra();
36			}
37
38
39			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)
40			{
41			init_tetra();
42			}
43
44			MG_TETRA* MG_TETRA::dupliquer(MG_MAILLAGE *mgmai,long decalage)
45			{
46			MG_NOEUD* nd1=mgmai->get_mg_noeudid(noeud1->get_id()+decalage);
47			MG_NOEUD* nd2=mgmai->get_mg_noeudid(noeud2->get_id()+decalage);
48			MG_NOEUD* nd3=mgmai->get_mg_noeudid(noeud3->get_id()+decalage);
49			MG_NOEUD* nd4=mgmai->get_mg_noeudid(noeud4->get_id()+decalage);
50			MG_TRIANGLE* tri1=mgmai->get_mg_triangleid(triangle1->get_id()+decalage);
51			MG_TRIANGLE* tri2=mgmai->get_mg_triangleid(triangle2->get_id()+decalage);
52			MG_TRIANGLE* tri3=mgmai->get_mg_triangleid(triangle3->get_id()+decalage);
53			MG_TRIANGLE* tri4=mgmai->get_mg_triangleid(triangle4->get_id()+decalage);
54	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);
55	francois	283	mgmai->ajouter_mg_tetra(tet);
56			return tet;
57			}
58
59			MG_TETRA::~MG_TETRA()
60			{
61			if (noeud1!=NULL) noeud1->get_lien_tetra()->supprimer(this);
62			if (noeud2!=NULL) noeud2->get_lien_tetra()->supprimer(this);
63			if (noeud3!=NULL) noeud3->get_lien_tetra()->supprimer(this);
64			if (noeud4!=NULL) noeud4->get_lien_tetra()->supprimer(this);
65			unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
66			noeudminid=std::min(noeudminid,noeud3->get_id());
67			noeudminid=std::min(noeudminid,noeud4->get_id());
68			if (noeud1!=NULL)
69			if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_tetra()->supprimer(this);
70			if (noeud2!=NULL)
71			if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_tetra()->supprimer(this);
72			if (noeud3!=NULL)
73			if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_tetra()->supprimer(this);
74			if (noeud4!=NULL)
75			if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_tetra()->supprimer(this);
76			if (triangle1!=NULL) {
77			triangle1->decrementer();
78			triangle1->get_lien_tetra()->supprimer(this);
79			}
80			if (triangle2!=NULL) {
81			triangle2->decrementer();
82			triangle2->get_lien_tetra()->supprimer(this);
83			}
84			if (triangle3!=NULL) {
85			triangle3->decrementer();
86			triangle3->get_lien_tetra()->supprimer(this);
87			}
88			if (triangle4!=NULL) {
89			triangle4->decrementer();
90			triangle4->get_lien_tetra()->supprimer(this);
91			}
92	francois	663	if (liaison_topologique!=NULL)
93			if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
94	francois	283	}
95
96			void MG_TETRA::init_tetra(void)
97			{
98			if (noeud1!=NULL) noeud1->get_lien_tetra()->ajouter(this);
99			if (noeud2!=NULL) noeud2->get_lien_tetra()->ajouter(this);
100			if (noeud3!=NULL) noeud3->get_lien_tetra()->ajouter(this);
101			if (noeud4!=NULL) noeud4->get_lien_tetra()->ajouter(this);
102			unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
103			noeudminid=std::min(noeudminid,noeud3->get_id());
104			noeudminid=std::min(noeudminid,noeud4->get_id());
105			if (noeud1!=NULL)
106			if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_tetra()->ajouter(this);
107			if (noeud2!=NULL)
108			if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_tetra()->ajouter(this);
109			if (noeud3!=NULL)
110			if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_tetra()->ajouter(this);
111			if (noeud4!=NULL)
112			if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_tetra()->ajouter(this);
113			if (triangle1!=NULL) {
114			triangle1->incrementer();
115			triangle1->get_lien_tetra()->ajouter(this);
116			}
117			if (triangle2!=NULL) {
118			triangle2->incrementer();
119			triangle2->get_lien_tetra()->ajouter(this);
120			}
121			if (triangle3!=NULL) {
122			triangle3->incrementer();
123			triangle3->get_lien_tetra()->ajouter(this);
124			}
125			if (triangle4!=NULL) {
126			triangle4->incrementer();
127			triangle4->get_lien_tetra()->ajouter(this);
128			}
129			if (liaison_topologique==NULL) return;
130			if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
131			}
132
133
134			int MG_TETRA::get_type_entite(void)
135			{
136	francois	1150	return MAGIC::TYPE_ENTITE::IDMG_TETRA;
137	francois	283	}
138
139
140
141
142
143			MG_NOEUD* MG_TETRA::get_noeud1(void)
144			{
145			return noeud1;
146			}
147
148			MG_NOEUD* MG_TETRA::get_noeud2(void)
149			{
150			return noeud2;
151			}
152
153			MG_NOEUD* MG_TETRA::get_noeud3(void)
154			{
155			return noeud3;
156			}
157
158			MG_NOEUD* MG_TETRA::get_noeud4(void)
159			{
160			return noeud4;
161			}
162
163			MG_TRIANGLE* MG_TETRA::get_triangle1(void)
164			{
165			return triangle1;
166			}
167
168			MG_TRIANGLE* MG_TETRA::get_triangle2(void)
169			{
170			return triangle2;
171			}
172
173			MG_TRIANGLE* MG_TETRA::get_triangle3(void)
174			{
175			return triangle3;
176			}
177
178			MG_TRIANGLE* MG_TETRA::get_triangle4(void)
179			{
180			return triangle4;
181			}
182
183	francois	702	void MG_TETRA::inverse_sens(void)
184			{
185			MG_NOEUD* noeud=noeud4;
186			noeud4=noeud3;
187			noeud3=noeud;
188			MG_TRIANGLE* triangle=triangle1;
189			triangle1=triangle2;
190			triangle2=triangle;
191			}
192	francois	283
193	francois	702
194	francois	283	int MG_TETRA::get_dimension(void)
195			{
196			return 3;
197			}
198
199
200			BOITE_3D MG_TETRA::get_boite_3D(void)
201			{
202			double coo1[3];
203			double coo2[3];
204			double coo3[3];
205			double coo4[3];
206			coo1[0]=noeud1->get_x();
207			coo1[1]=noeud1->get_y();
208			coo1[2]=noeud1->get_z();
209			coo2[0]=noeud2->get_x();
210			coo2[1]=noeud2->get_y();
211			coo2[2]=noeud2->get_z();
212			coo3[0]=noeud3->get_x();
213			coo3[1]=noeud3->get_y();
214			coo3[2]=noeud3->get_z();
215			coo4[0]=noeud4->get_x();
216			coo4[1]=noeud4->get_y();
217			coo4[2]=noeud4->get_z();
218			double xmin,ymin,zmin,xmax,ymax,zmax;
219			xmin=std::min(coo1[0],coo2[0]);
220			xmin=std::min(xmin,coo3[0]);
221			xmin=std::min(xmin,coo4[0]);
222			xmax=std::max(coo1[0],coo2[0]);
223			xmax=std::max(xmax,coo3[0]);
224			xmax=std::max(xmax,coo4[0]);
225			ymin=std::min(coo1[1],coo2[1]);
226			ymin=std::min(ymin,coo3[1]);
227			ymin=std::min(ymin,coo4[1]);
228			ymax=std::max(coo1[1],coo2[1]);
229			ymax=std::max(ymax,coo3[1]);
230			ymax=std::max(ymax,coo4[1]);
231			zmin=std::min(coo1[2],coo2[2]);
232			zmin=std::min(zmin,coo3[2]);
233			zmin=std::min(zmin,coo4[2]);
234			zmax=std::max(coo1[2],coo2[2]);
235			zmax=std::max(zmax,coo3[2]);
236			zmax=std::max(zmax,coo4[2]);
237			BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
238			return boite;
239			}
240			void MG_TETRA::extrapoler_solution_noeud(void)
241			{
242	francois	377	for (int j=0;j<MAX_TYPE_SOLUTION;j++)
243			{
244			noeud1->change_solution(solution[j],j);
245			noeud2->change_solution(solution[j],j);
246			noeud3->change_solution(solution[j],j);
247			noeud4->change_solution(solution[j],j);
248			}
249	francois	283	}
250
251	francois	763	void MG_TETRA::enregistrer(std::ostream& o,double version)
252	francois	283	{
253			if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
254			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;
255			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;
256			}
257
258
259
260
261
262