ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/mtu/src/mg_penta.cpp
Revision: 1158
Committed: Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size: 12266 byte(s)
Log Message:
compatibilité Ubuntu 22.04
Suppression des refeences à Windows
Ajout d'une banière

File Contents

# 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_penta.cpp
15     //####//
16     //####//------------------------------------------------------------
17     //####//------------------------------------------------------------
18     //####// COPYRIGHT 2000-2024
19     //####// jeu 13 jun 2024 11:58:55 EDT
20     //####//------------------------------------------------------------
21     //####//------------------------------------------------------------
22 francois 876
23    
24     #include "gestionversion.h"
25     #include "mg_penta.h"
26     #include "mg_maillage.h"
27    
28     MG_PENTA::MG_PENTA(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_NOEUD *mgnoeud6,class MG_TRIANGLE* mgtriangle1,MG_TRIANGLE* mgtriangle2,class MG_QUADRANGLE* mgquadrangle1,class MG_QUADRANGLE* mgquadrangle2,class MG_QUADRANGLE* mgquadrangle3,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),noeud6(mgnoeud6),triangle1(mgtriangle1),triangle2(mgtriangle2),quadrangle1(mgquadrangle1),quadrangle2(mgquadrangle2),quadrangle3(mgquadrangle3)
29     {
30     init_penta();
31     }
32    
33     MG_PENTA::MG_PENTA(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_NOEUD *mgnoeud6,class MG_TRIANGLE* mgtriangle1,MG_TRIANGLE* mgtriangle2,class MG_QUADRANGLE* mgquadrangle1,class MG_QUADRANGLE* mgquadrangle2,class MG_QUADRANGLE* mgquadrangle3,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),noeud6(mgnoeud6),triangle1(mgtriangle1),triangle2(mgtriangle2),quadrangle1(mgquadrangle1),quadrangle2(mgquadrangle2),quadrangle3(mgquadrangle3)
34     {
35     init_penta();
36     }
37    
38    
39     MG_PENTA::MG_PENTA(MG_PENTA& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),noeud5(mdd.noeud5),noeud6(mdd.noeud6),triangle1(mdd.triangle1),triangle2(mdd.triangle2),quadrangle1(mdd.quadrangle1),quadrangle2(mdd.quadrangle2),quadrangle3(mdd.quadrangle3)
40     {
41     init_penta();
42     }
43    
44    
45     MG_PENTA* MG_PENTA::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_NOEUD* nd6=mgmai->get_mg_noeudid(noeud6->get_id()+decalage);
53     MG_TRIANGLE* tri1=mgmai->get_mg_triangleid(triangle1->get_id()+decalage);
54     MG_TRIANGLE* tri2=mgmai->get_mg_triangleid(triangle2->get_id()+decalage);
55     MG_QUADRANGLE* quad1=mgmai->get_mg_quadrangleid(quadrangle1->get_id()+decalage);
56     MG_QUADRANGLE* quad2=mgmai->get_mg_quadrangleid(quadrangle2->get_id()+decalage);
57     MG_QUADRANGLE* quad3=mgmai->get_mg_quadrangleid(quadrangle3->get_id()+decalage);
58     MG_PENTA* pen=new MG_PENTA(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,nd5,nd6,tri1,tri2,quad1,quad2,quad3,MAGIC::ORIGINE::DUPLIQUER);
59     mgmai->ajouter_mg_penta(pen);
60     return pen;
61     }
62    
63     MG_PENTA::~MG_PENTA()
64     {
65     if (noeud1!=NULL) noeud1->get_lien_penta()->supprimer(this);
66     if (noeud2!=NULL) noeud2->get_lien_penta()->supprimer(this);
67     if (noeud3!=NULL) noeud3->get_lien_penta()->supprimer(this);
68     if (noeud4!=NULL) noeud4->get_lien_penta()->supprimer(this);
69     if (noeud5!=NULL) noeud5->get_lien_penta()->supprimer(this);
70     if (noeud6!=NULL) noeud6->get_lien_penta()->supprimer(this);
71     unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
72     noeudminid=std::min(noeudminid,noeud3->get_id());
73     noeudminid=std::min(noeudminid,noeud4->get_id());
74     noeudminid=std::min(noeudminid,noeud5->get_id());
75     noeudminid=std::min(noeudminid,noeud6->get_id());
76     if (noeud1!=NULL)
77     if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_penta()->supprimer(this);
78     if (noeud2!=NULL)
79     if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_penta()->supprimer(this);
80     if (noeud3!=NULL)
81     if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_penta()->supprimer(this);
82     if (noeud4!=NULL)
83     if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_penta()->supprimer(this);
84     if (noeud5!=NULL)
85     if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_penta()->supprimer(this);
86     if (noeud6!=NULL)
87     if (noeud6->get_id()==noeudminid) noeud6->get_lien_petit_penta()->supprimer(this);
88     if (triangle1!=NULL) {
89     triangle1->decrementer();
90     triangle1->get_lien_penta()->supprimer(this);
91     }
92     if (triangle2!=NULL) {
93     triangle2->decrementer();
94     triangle2->get_lien_penta()->supprimer(this);
95     }
96     if (quadrangle1!=NULL) {
97     quadrangle1->decrementer();
98     quadrangle1->get_lien_penta()->supprimer(this);
99     }
100     if (quadrangle2!=NULL) {
101     quadrangle2->decrementer();
102     quadrangle2->get_lien_penta()->supprimer(this);
103     }
104     if (quadrangle3!=NULL) {
105     quadrangle3->decrementer();
106     quadrangle3->get_lien_penta()->supprimer(this);
107     }
108     if (liaison_topologique!=NULL)
109     if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
110     }
111    
112     void MG_PENTA::init_penta(void)
113     {
114     if (noeud1!=NULL) noeud1->get_lien_penta()->ajouter(this);
115     if (noeud2!=NULL) noeud2->get_lien_penta()->ajouter(this);
116     if (noeud3!=NULL) noeud3->get_lien_penta()->ajouter(this);
117     if (noeud4!=NULL) noeud4->get_lien_penta()->ajouter(this);
118     if (noeud5!=NULL) noeud5->get_lien_penta()->ajouter(this);
119     if (noeud6!=NULL) noeud6->get_lien_penta()->ajouter(this);
120     unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
121     noeudminid=std::min(noeudminid,noeud3->get_id());
122     noeudminid=std::min(noeudminid,noeud4->get_id());
123     noeudminid=std::min(noeudminid,noeud5->get_id());
124     noeudminid=std::min(noeudminid,noeud6->get_id());
125     if (noeud1!=NULL)
126     if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_penta()->ajouter(this);
127     if (noeud2!=NULL)
128     if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_penta()->ajouter(this);
129     if (noeud3!=NULL)
130     if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_penta()->ajouter(this);
131     if (noeud4!=NULL)
132     if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_penta()->ajouter(this);
133     if (noeud5!=NULL)
134     if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_penta()->ajouter(this);
135     if (noeud6!=NULL)
136     if (noeud6->get_id()==noeudminid) noeud6->get_lien_petit_penta()->ajouter(this);
137     if (triangle1!=NULL) {
138     triangle1->incrementer();
139     triangle1->get_lien_penta()->ajouter(this);
140     }
141     if (triangle2!=NULL) {
142     triangle2->incrementer();
143     triangle2->get_lien_penta()->ajouter(this);
144     }
145     if (quadrangle1!=NULL) {
146     quadrangle1->incrementer();
147     quadrangle1->get_lien_penta()->ajouter(this);
148     }
149     if (quadrangle2!=NULL) {
150     quadrangle2->incrementer();
151     quadrangle2->get_lien_penta()->ajouter(this);
152     }
153     if (quadrangle3!=NULL) {
154     quadrangle3->incrementer();
155     quadrangle3->get_lien_penta()->ajouter(this);
156     }
157    
158     if (liaison_topologique==NULL) return;
159     if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
160     }
161    
162    
163     int MG_PENTA::get_type_entite(void)
164     {
165 francois 1150 return MAGIC::TYPE_ENTITE::IDMG_PENTA;
166 francois 876 }
167    
168    
169    
170    
171    
172     MG_NOEUD* MG_PENTA::get_noeud1(void)
173     {
174     return noeud1;
175     }
176    
177     MG_NOEUD* MG_PENTA::get_noeud2(void)
178     {
179     return noeud2;
180     }
181    
182     MG_NOEUD* MG_PENTA::get_noeud3(void)
183     {
184     return noeud3;
185     }
186    
187     MG_NOEUD* MG_PENTA::get_noeud4(void)
188     {
189     return noeud4;
190     }
191     MG_NOEUD* MG_PENTA::get_noeud5(void)
192     {
193     return noeud5;
194     }
195     MG_NOEUD* MG_PENTA::get_noeud6(void)
196     {
197     return noeud6;
198     }
199    
200     MG_QUADRANGLE* MG_PENTA::get_quadrangle1(void)
201     {
202     return quadrangle1;
203     }
204    
205     MG_QUADRANGLE* MG_PENTA::get_quadrangle2(void)
206     {
207     return quadrangle2;
208     }
209    
210     MG_QUADRANGLE* MG_PENTA::get_quadrangle3(void)
211     {
212     return quadrangle3;
213     }
214    
215     MG_TRIANGLE* MG_PENTA::get_triangle1(void)
216     {
217     return triangle1;
218     }
219    
220     MG_TRIANGLE* MG_PENTA::get_triangle2(void)
221     {
222     return triangle2;
223     }
224    
225     void MG_PENTA::inverse_sens(void)
226     {
227     MG_NOEUD* noeud=noeud3;
228     noeud3=noeud1;
229     noeud1=noeud;
230     noeud=noeud6;
231     noeud6=noeud4;
232     noeud4=noeud;
233     MG_QUADRANGLE* quadrangle=quadrangle1;
234     quadrangle1=quadrangle2;
235     quadrangle2=quadrangle;
236     }
237    
238    
239     int MG_PENTA::get_dimension(void)
240     {
241     return 3;
242     }
243    
244    
245     BOITE_3D MG_PENTA::get_boite_3D(void)
246     {
247     double coo1[3];
248     double coo2[3];
249     double coo3[3];
250     double coo4[3];
251     double coo5[3];
252     double coo6[3];
253     coo1[0]=noeud1->get_x();
254     coo1[1]=noeud1->get_y();
255     coo1[2]=noeud1->get_z();
256     coo2[0]=noeud2->get_x();
257     coo2[1]=noeud2->get_y();
258     coo2[2]=noeud2->get_z();
259     coo3[0]=noeud3->get_x();
260     coo3[1]=noeud3->get_y();
261     coo3[2]=noeud3->get_z();
262     coo4[0]=noeud4->get_x();
263     coo4[1]=noeud4->get_y();
264     coo4[2]=noeud4->get_z();
265     coo5[0]=noeud5->get_x();
266     coo5[1]=noeud5->get_y();
267     coo5[2]=noeud5->get_z();
268     coo6[0]=noeud6->get_x();
269     coo6[1]=noeud6->get_y();
270     coo6[2]=noeud6->get_z();
271     double xmin,ymin,zmin,xmax,ymax,zmax;
272     xmin=std::min(coo1[0],coo2[0]);
273     xmin=std::min(xmin,coo3[0]);
274     xmin=std::min(xmin,coo4[0]);
275     xmin=std::min(xmin,coo5[0]);
276     xmin=std::min(xmin,coo6[0]);
277     xmax=std::max(coo1[0],coo2[0]);
278     xmax=std::max(xmax,coo3[0]);
279     xmax=std::max(xmax,coo4[0]);
280     xmax=std::max(xmax,coo5[0]);
281     xmax=std::max(xmax,coo6[0]);
282     ymin=std::min(coo1[1],coo2[1]);
283     ymin=std::min(ymin,coo3[1]);
284     ymin=std::min(ymin,coo4[1]);
285     ymin=std::min(ymin,coo5[1]);
286     ymin=std::min(ymin,coo6[1]);
287     ymax=std::max(coo1[1],coo2[1]);
288     ymax=std::max(ymax,coo3[1]);
289     ymax=std::max(ymax,coo4[1]);
290     ymax=std::max(ymax,coo5[1]);
291     ymax=std::max(ymax,coo6[1]);
292     zmin=std::min(coo1[2],coo2[2]);
293     zmin=std::min(zmin,coo3[2]);
294     zmin=std::min(zmin,coo4[2]);
295     zmin=std::min(zmin,coo5[2]);
296     zmin=std::min(zmin,coo6[2]);
297     zmax=std::max(coo1[2],coo2[2]);
298     zmax=std::max(zmax,coo3[2]);
299     zmax=std::max(zmax,coo4[2]);
300     zmax=std::max(zmax,coo5[2]);
301     zmax=std::max(zmax,coo6[2]);
302     BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
303     return boite;
304     }
305     void MG_PENTA::extrapoler_solution_noeud(void)
306     {
307     for (int j=0;j<MAX_TYPE_SOLUTION;j++)
308     {
309     noeud1->change_solution(solution[j],j);
310     noeud2->change_solution(solution[j],j);
311     noeud3->change_solution(solution[j],j);
312     noeud4->change_solution(solution[j],j);
313     noeud5->change_solution(solution[j],j);
314     noeud6->change_solution(solution[j],j);
315     }
316     }
317    
318     void MG_PENTA::enregistrer(std::ostream& o,double version)
319     {
320     if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
321 francois 1123 if (liaison_topologique==NULL) o << "%" << get_id() << "=PENTAEDRE(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id() << ",$" << noeud4->get_id() << ",$" << noeud5->get_id() << ",$" << noeud6->get_id() << "," << origine_creation << ");" << std::endl;
322 francois 876 else if (liaison_topologique->get_dimension()==3) o << "%" << get_id() << "=PENTAEDRE($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id() << ",$" << noeud4->get_id() << ",$" <<noeud5->get_id() << ",$" << noeud6->get_id() << "," << origine_creation << ");" << std::endl;
323     }
324    
325    
326    
327    
328    
329