]> SALOME platform Git repositories - modules/smesh.git/blob - src/Tools/MeshCut/MeshCut_Cas.cxx
Salome HOME
Compilation on Windows platform.
[modules/smesh.git] / src / Tools / MeshCut / MeshCut_Cas.cxx
1 // Copyright (C) 2006-2014  EDF R&D
2 //
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License, or (at your option) any later version.
7 //
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
11 // Lesser General Public License for more details.
12 //
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
16 //
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
18 //
19
20 #include "MeshCut_Cas.hxx"
21
22 #include "MeshCut_Globals.hxx"
23
24 using namespace MESHCUT;
25 using namespace std;
26
27 /*!
28  * Le cas 1 traduit le fait que deux des sommets du T4 initial sont dans le plan de coupe.
29  * Le point d'intersection franc trouvé est sur l'arête opposée à ces deux points du T4.
30
31  * Le T4 initial produit deux nouveaux T4.
32  */
33 void MESHCUT::cas1(int VN[6], int it4)
34 {
35   cutTetras.push_back(it4);
36   //  cout << "Cas 1 - it4=" << it4 << ", VN = " << VN[0] << " " << VN[1] << " " << VN[2] << " " << VN[3] << " " << VN[4]
37   //      << " " << VN[5] << " " << endl;
38
39   // Numéros des noeuds du TETRA4
40   int ng0 = MAILLAGE1->CNX[TETRA4][4 * it4 + 0];
41   int ng1 = MAILLAGE1->CNX[TETRA4][4 * it4 + 1];
42   int ng2 = MAILLAGE1->CNX[TETRA4][4 * it4 + 2];
43   int ng3 = MAILLAGE1->CNX[TETRA4][4 * it4 + 3];
44
45   int i1, i2;
46
47   if (VN[0] != -1)
48     {
49       // Le sommet de T4new1 servant à la détermination du groupe est le noeud 1 du T4 d'origine
50
51       // cout << "cas 1/0" << endl;
52       newCNX[TETRA4].push_back(VN[0]);
53       newCNX[TETRA4].push_back(ng1);
54       newCNX[TETRA4].push_back(ng2);
55       newCNX[TETRA4].push_back(ng3);
56       cptNouvellesMailles[TETRA4]++;
57       i1 = cptNouvellesMailles[TETRA4] - 1;
58
59       newCNX[TETRA4].push_back(ng0);
60       newCNX[TETRA4].push_back(VN[0]);
61       newCNX[TETRA4].push_back(ng2);
62       newCNX[TETRA4].push_back(ng3);
63       cptNouvellesMailles[TETRA4]++;
64       i2 = cptNouvellesMailles[TETRA4] - 1;
65
66       if (POSN[ng1 - 1] == 1)
67         {
68           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
69           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
70         }
71       else
72         {
73           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
74           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
75         }
76     }
77
78   else if (VN[1] != -1)
79     {
80       // Le sommet de T4new1 servant à la détermination du groupe est le noeud 0 du T4 d'origine
81
82       // cout << "cas 1/1" << endl;
83       newCNX[TETRA4].push_back(ng0);
84       newCNX[TETRA4].push_back(ng1);
85       newCNX[TETRA4].push_back(VN[1]);
86       newCNX[TETRA4].push_back(ng3);
87       cptNouvellesMailles[TETRA4]++;
88       i1 = cptNouvellesMailles[TETRA4] - 1;
89
90       newCNX[TETRA4].push_back(ng1);
91       newCNX[TETRA4].push_back(ng2);
92       newCNX[TETRA4].push_back(VN[1]);
93       newCNX[TETRA4].push_back(ng3);
94       cptNouvellesMailles[TETRA4]++;
95       i2 = cptNouvellesMailles[TETRA4] - 1;
96
97       if (POSN[ng0 - 1] == 1)
98         {
99           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
100           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
101         }
102       else
103         {
104           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
105           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
106         }
107     }
108
109   else if (VN[2] != -1)
110     {
111       // Le sommet de T4new1 servant à la détermination du groupe est le noeud 0 du T4 d'origine
112
113       // cout << "cas 1/2" << endl;
114       newCNX[TETRA4].push_back(ng0);
115       newCNX[TETRA4].push_back(ng1);
116       newCNX[TETRA4].push_back(ng2);
117       newCNX[TETRA4].push_back(VN[2]);
118       cptNouvellesMailles[TETRA4]++;
119       i1 = cptNouvellesMailles[TETRA4] - 1;
120
121       newCNX[TETRA4].push_back(ng1);
122       newCNX[TETRA4].push_back(ng2);
123       newCNX[TETRA4].push_back(VN[2]);
124       newCNX[TETRA4].push_back(ng3);
125       cptNouvellesMailles[TETRA4]++;
126       i2 = cptNouvellesMailles[TETRA4] - 1;
127
128       if (POSN[ng0 - 1] == 1)
129         {
130           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
131           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
132         }
133       else
134         {
135           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
136           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
137         }
138     }
139
140   else if (VN[3] != -1)
141     {
142       // Le sommet de T4new1 servant à la détermination du groupe est le noeud 1 du T4 d'origine
143
144       // cout << "cas 1/3" << endl;
145       newCNX[TETRA4].push_back(ng0);
146       newCNX[TETRA4].push_back(ng1);
147       newCNX[TETRA4].push_back(VN[3]);
148       newCNX[TETRA4].push_back(ng3);
149       cptNouvellesMailles[TETRA4]++;
150       i1 = cptNouvellesMailles[TETRA4] - 1;
151
152       newCNX[TETRA4].push_back(ng0);
153       newCNX[TETRA4].push_back(VN[3]);
154       newCNX[TETRA4].push_back(ng2);
155       newCNX[TETRA4].push_back(ng3);
156       cptNouvellesMailles[TETRA4]++;
157       i2 = cptNouvellesMailles[TETRA4] - 1;
158
159       if (POSN[ng1 - 1] == 1)
160         {
161           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
162           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
163         }
164       else
165         {
166           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
167           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
168         }
169     }
170
171   else if (VN[4] != -1)
172     {
173       // Le sommet de T4new1 servant à la détermination du groupe est le noeud 1 du T4 d'origine
174
175       // cout << "cas 1/4" << endl;
176       newCNX[TETRA4].push_back(ng0);
177       newCNX[TETRA4].push_back(ng1);
178       newCNX[TETRA4].push_back(ng2);
179       newCNX[TETRA4].push_back(VN[4]);
180       cptNouvellesMailles[TETRA4]++;
181       i1 = cptNouvellesMailles[TETRA4] - 1;
182
183       newCNX[TETRA4].push_back(ng0);
184       newCNX[TETRA4].push_back(VN[4]);
185       newCNX[TETRA4].push_back(ng2);
186       newCNX[TETRA4].push_back(ng3);
187       cptNouvellesMailles[TETRA4]++;
188       i2 = cptNouvellesMailles[TETRA4] - 1;
189
190       if (POSN[ng1 - 1] == 1)
191         {
192           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
193           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
194         }
195       else
196         {
197           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
198           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
199         }
200     }
201
202   else if (VN[5] != -1)
203     {
204       // Le sommet de T4new1 servant à la détermination du groupe est le noeud 3 du T4 d'origine
205
206       // cout << "cas 1/5" << endl;
207       newCNX[TETRA4].push_back(ng0);
208       newCNX[TETRA4].push_back(ng1);
209       newCNX[TETRA4].push_back(VN[5]);
210       newCNX[TETRA4].push_back(ng3);
211       cptNouvellesMailles[TETRA4]++;
212       i1 = cptNouvellesMailles[TETRA4] - 1;
213
214       newCNX[TETRA4].push_back(ng0);
215       newCNX[TETRA4].push_back(ng1);
216       newCNX[TETRA4].push_back(ng2);
217       newCNX[TETRA4].push_back(VN[5]);
218       cptNouvellesMailles[TETRA4]++;
219       i2 = cptNouvellesMailles[TETRA4] - 1;
220
221       if (POSN[ng3 - 1] == 1)
222         {
223           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
224           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
225         }
226       else
227         {
228           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2);
229           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1);
230         }
231     }
232
233   else
234     ERREUR(" Intersections configuration not taken into account (case cptPI=1) ");
235
236   //int nl1 = MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i1;
237   //int nl2 = MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + i2;
238   // cout << "La maille TETRA4 " << it4 << " produit les mailles TETRA4-" << nl1 << " et TETRA4-" << nl2 << endl;
239
240 }
241
242 /*!  Deux points d'intersection
243  *  Le cas 2 traduit le fait qu'un des sommets du T4 est dans le plan de coupe.
244  *  Ce sommet est celui des quatre qui n'appartient à aucune des deux arêtes sur lesquelles
245  *  un point d'intersection non -1 a été calculé.
246  *
247  *  Le T4 initial produit un nouveau T4 et un élément PYRAM5.
248  *
249  */
250 void MESHCUT::cas2(int VN[6], int it4)
251 {
252   cutTetras.push_back(it4);
253   //  cout << "Cas 2 - it4=" << it4 << ", VN = " << VN[0] << " " << VN[1] << " " << VN[2] << " " << VN[3] << " " << VN[4]
254   //      << " " << VN[5] << " " << endl;
255
256   // Numéros des noeuds du TETRA4
257   int ng0 = MAILLAGE1->CNX[TETRA4][4 * it4 + 0];
258   int ng1 = MAILLAGE1->CNX[TETRA4][4 * it4 + 1];
259   int ng2 = MAILLAGE1->CNX[TETRA4][4 * it4 + 2];
260   int ng3 = MAILLAGE1->CNX[TETRA4][4 * it4 + 3];
261
262   if (VN[0] != -1 && VN[1] != -1)
263     {
264       // Le sommet du nouveau T4 est le noeud 0 du T4 d'origine
265
266       // cout << "cas 2.01" << endl;
267       newCNX[TETRA4].push_back(VN[0]);
268       newCNX[TETRA4].push_back(ng3);
269       newCNX[TETRA4].push_back(VN[1]);
270       newCNX[TETRA4].push_back(ng0);
271       cptNouvellesMailles[TETRA4]++;
272
273       newCNX[PYRAM5].push_back(VN[0]);
274       newCNX[PYRAM5].push_back(ng1);
275       newCNX[PYRAM5].push_back(ng2);
276       newCNX[PYRAM5].push_back(VN[1]);
277       newCNX[PYRAM5].push_back(ng3);
278       cptNouvellesMailles[PYRAM5]++;
279
280       if (POSN[ng0 - 1] == 1)
281         {
282           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
283           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
284         }
285       else
286         {
287           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
288           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
289         }
290     }
291   else if (VN[0] != -1 && VN[2] != -1)
292     {
293       // Le sommet du nouveau T4 est le noeud 0 du T4 d'origine
294
295       // cout << "cas 2.02" << endl;
296       newCNX[TETRA4].push_back(VN[0]);
297       newCNX[TETRA4].push_back(VN[2]);
298       newCNX[TETRA4].push_back(ng2);
299       newCNX[TETRA4].push_back(ng0);
300       cptNouvellesMailles[TETRA4]++;
301
302       newCNX[PYRAM5].push_back(VN[0]);
303       newCNX[PYRAM5].push_back(VN[2]);
304       newCNX[PYRAM5].push_back(ng3);
305       newCNX[PYRAM5].push_back(ng1);
306       newCNX[PYRAM5].push_back(ng2);
307       cptNouvellesMailles[PYRAM5]++;
308
309       if (POSN[ng0 - 1] == 1)
310         {
311           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
312           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
313         }
314       else
315         {
316           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
317           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
318         }
319     }
320   else if (VN[0] != -1 && VN[3] != -1)
321     {
322       // Le sommet du nouveau T4 est le noeud 1 du T4 d'origine
323
324       // cout << "cas 2.03" << endl;
325       newCNX[TETRA4].push_back(VN[0]);
326       newCNX[TETRA4].push_back(VN[3]);
327       newCNX[TETRA4].push_back(ng3);
328       newCNX[TETRA4].push_back(ng1);
329       cptNouvellesMailles[TETRA4]++;
330
331       newCNX[PYRAM5].push_back(VN[0]);
332       newCNX[PYRAM5].push_back(VN[3]);
333       newCNX[PYRAM5].push_back(ng2);
334       newCNX[PYRAM5].push_back(ng0);
335       newCNX[PYRAM5].push_back(ng3);
336       cptNouvellesMailles[PYRAM5]++;
337
338       if (POSN[ng1 - 1] == 1)
339         {
340           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
341           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
342         }
343       else
344         {
345           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
346           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
347         }
348     }
349   else if (VN[0] != -1 && VN[4] != -1)
350     {
351       // Le sommet du nouveau T4 est le noeud 1 du T4 d'origine
352
353       // cout << "cas 2.04" << endl;
354       newCNX[TETRA4].push_back(VN[0]);
355       newCNX[TETRA4].push_back(ng2);
356       newCNX[TETRA4].push_back(VN[4]);
357       newCNX[TETRA4].push_back(ng1);
358       cptNouvellesMailles[TETRA4]++;
359
360       newCNX[PYRAM5].push_back(ng0);
361       newCNX[PYRAM5].push_back(ng3);
362       newCNX[PYRAM5].push_back(VN[4]);
363       newCNX[PYRAM5].push_back(VN[0]);
364       newCNX[PYRAM5].push_back(ng2);
365       cptNouvellesMailles[PYRAM5]++;
366
367       if (POSN[ng1 - 1] == 1)
368         {
369           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
370           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
371         }
372       else
373         {
374           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
375           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
376         }
377     }
378
379   else if (VN[0] != -1 && VN[5] != -1)
380     ERREUR("Case 2/05 forbidden");
381
382   else if (VN[1] != -1 && VN[2] != -1)
383     {
384       // Le sommet du nouveau T4 est le noeud 0 du T4 d'origine
385
386       // cout << "cas 2.12" << endl;
387       newCNX[TETRA4].push_back(ng1);
388       newCNX[TETRA4].push_back(VN[2]);
389       newCNX[TETRA4].push_back(VN[1]);
390       newCNX[TETRA4].push_back(ng0);
391       cptNouvellesMailles[TETRA4]++;
392
393       newCNX[PYRAM5].push_back(ng2);
394       newCNX[PYRAM5].push_back(ng3);
395       newCNX[PYRAM5].push_back(VN[2]);
396       newCNX[PYRAM5].push_back(VN[1]);
397       newCNX[PYRAM5].push_back(ng1);
398       cptNouvellesMailles[PYRAM5]++;
399
400       if (POSN[ng0 - 1] == 1)
401         {
402           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
403           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
404         }
405       else
406         {
407           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
408           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
409         }
410     }
411   else if (VN[1] != -1 && VN[3] != -1)
412     {
413       // Le sommet du nouveau T4 est le noeud 2 du T4 d'origine
414
415       // cout << "cas 2.13" << endl;
416       newCNX[TETRA4].push_back(VN[1]);
417       newCNX[TETRA4].push_back(ng3);
418       newCNX[TETRA4].push_back(VN[3]);
419       newCNX[TETRA4].push_back(ng2);
420       cptNouvellesMailles[TETRA4]++;
421
422       newCNX[PYRAM5].push_back(ng0);
423       newCNX[PYRAM5].push_back(ng1);
424       newCNX[PYRAM5].push_back(VN[3]);
425       newCNX[PYRAM5].push_back(VN[1]);
426       newCNX[PYRAM5].push_back(ng3);
427       cptNouvellesMailles[PYRAM5]++;
428
429       if (POSN[ng2 - 1] == 1)
430         {
431           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
432           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
433         }
434       else
435         {
436           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
437           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
438         }
439     }
440
441   else if (VN[1] != -1 && VN[4] != -1)
442     ERREUR("Case 2/14 excluded");
443
444   else if (VN[1] != -1 && VN[5] != -1)
445     {
446       // Le sommet du nouveau T4 est le noeud 2 du T4 d'origine
447
448       // cout << "cas 2.15" << endl;
449       newCNX[TETRA4].push_back(ng1);
450       newCNX[TETRA4].push_back(ng2);
451       newCNX[TETRA4].push_back(VN[1]);
452       newCNX[TETRA4].push_back(VN[5]);
453       cptNouvellesMailles[TETRA4]++;
454
455       newCNX[PYRAM5].push_back(VN[1]);
456       newCNX[PYRAM5].push_back(VN[5]);
457       newCNX[PYRAM5].push_back(ng3);
458       newCNX[PYRAM5].push_back(ng0);
459       newCNX[PYRAM5].push_back(ng1);
460       cptNouvellesMailles[PYRAM5]++;
461
462       if (POSN[ng2 - 1] == 1)
463         {
464           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
465           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
466         }
467       else
468         {
469           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
470           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
471         }
472     }
473
474   else if (VN[2] != -1 && VN[3] != -1)
475     ERREUR("Case 2/23 excluded");
476
477   else if (VN[2] != -1 && VN[4] != -1)
478     {
479       // Le sommet du nouveau T4 est le noeud 3 du T4 d'origine
480
481       // cout << "cas 2.24" << endl;
482       newCNX[TETRA4].push_back(VN[2]);
483       newCNX[TETRA4].push_back(VN[4]);
484       newCNX[TETRA4].push_back(ng2);
485       newCNX[TETRA4].push_back(ng3);
486       cptNouvellesMailles[TETRA4]++;
487
488       newCNX[PYRAM5].push_back(ng0);
489       newCNX[PYRAM5].push_back(VN[2]);
490       newCNX[PYRAM5].push_back(VN[4]);
491       newCNX[PYRAM5].push_back(ng1);
492       newCNX[PYRAM5].push_back(ng2);
493       cptNouvellesMailles[PYRAM5]++;
494
495       if (POSN[ng3 - 1] == 1)
496         {
497           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
498           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
499         }
500       else
501         {
502           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
503           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
504         }
505     }
506
507   else if (VN[2] != -1 && VN[5] != -1)
508     {
509       // Le sommet du nouveau T4 est le noeud 3 du T4 d'origine
510
511       // cout << "cas 2.25" << endl;
512       newCNX[TETRA4].push_back(ng1);
513       newCNX[TETRA4].push_back(VN[5]);
514       newCNX[TETRA4].push_back(VN[2]);
515       newCNX[TETRA4].push_back(ng3);
516       cptNouvellesMailles[TETRA4]++;
517
518       newCNX[PYRAM5].push_back(ng0);
519       newCNX[PYRAM5].push_back(ng2);
520       newCNX[PYRAM5].push_back(VN[5]);
521       newCNX[PYRAM5].push_back(VN[2]);
522       newCNX[PYRAM5].push_back(ng1);
523       cptNouvellesMailles[PYRAM5]++;
524
525       if (POSN[ng3 - 1] == 1)
526         {
527           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
528           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
529         }
530       else
531         {
532           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
533           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
534         }
535     }
536
537   else if (VN[3] != -1 && VN[4] != -1)
538     {
539       // Le sommet du nouveau T4 est le noeud 1 du T4 d'origine
540
541       // cout << "cas 2.34" << endl;
542       newCNX[TETRA4].push_back(ng0);
543       newCNX[TETRA4].push_back(VN[3]);
544       newCNX[TETRA4].push_back(VN[4]);
545       newCNX[TETRA4].push_back(ng1);
546       cptNouvellesMailles[TETRA4]++;
547
548       newCNX[PYRAM5].push_back(VN[3]);
549       newCNX[PYRAM5].push_back(VN[4]);
550       newCNX[PYRAM5].push_back(ng3);
551       newCNX[PYRAM5].push_back(ng2);
552       newCNX[PYRAM5].push_back(ng0);
553       cptNouvellesMailles[PYRAM5]++;
554
555       if (POSN[ng1 - 1] == 1)
556         {
557           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
558           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
559         }
560       else
561         {
562           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
563           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
564         }
565     }
566
567   else if (VN[3] != -1 && VN[5] != -1)
568     {
569       // Le sommet du nouveau T4 est le noeud 2 du T4 d'origine
570
571       // cout << "cas 2.35" << endl;
572       newCNX[TETRA4].push_back(VN[3]);
573       newCNX[TETRA4].push_back(VN[5]);
574       newCNX[TETRA4].push_back(ng2);
575       newCNX[TETRA4].push_back(ng0);
576       cptNouvellesMailles[TETRA4]++;
577
578       newCNX[PYRAM5].push_back(VN[3]);
579       newCNX[PYRAM5].push_back(ng1);
580       newCNX[PYRAM5].push_back(ng3);
581       newCNX[PYRAM5].push_back(VN[5]);
582       newCNX[PYRAM5].push_back(ng0);
583       cptNouvellesMailles[PYRAM5]++;
584
585       if (POSN[ng2 - 1] == 1)
586         {
587           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
588           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
589         }
590       else
591         {
592           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
593           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
594         }
595     }
596
597   else if (VN[4] != -1 && VN[5] != -1)
598     {
599       // Le sommet du nouveau T4 est le noeud 3 du T4 d'origine
600
601       // cout << "cas 2.35" << endl;
602       newCNX[TETRA4].push_back(ng0);
603       newCNX[TETRA4].push_back(VN[4]);
604       newCNX[TETRA4].push_back(VN[5]);
605       newCNX[TETRA4].push_back(ng3);
606       cptNouvellesMailles[TETRA4]++;
607
608       newCNX[PYRAM5].push_back(ng1);
609       newCNX[PYRAM5].push_back(VN[4]);
610       newCNX[PYRAM5].push_back(VN[5]);
611       newCNX[PYRAM5].push_back(ng2);
612       newCNX[PYRAM5].push_back(ng0);
613       cptNouvellesMailles[PYRAM5]++;
614
615       if (POSN[ng3 - 1] == 1)
616         {
617           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
618           GMmoins[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
619         }
620       else
621         {
622           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
623           GMplus[PYRAM5].push_back(MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1);
624         }
625     }
626
627   else
628     ERREUR(" Intersections configuration not taken into account (case cptPI=2) ");
629
630   // int ngT4 = MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1;
631   // int ngP5 = MAILLAGE1->EFFECTIFS_TYPES[PYRAM5] + cptNouvellesMailles[PYRAM5] - 1;
632   // cout << "La maille TETRA4 " << it4 << " produit les mailles TETRA4-" << ngT4 << " et PYRAM5-" << ngP5 << endl;
633
634 }
635
636 /*! Trois points d'intersection
637  *  ATTENTION: pour les PENTA6 on adopte la convention d'orientation SALOME :
638  *
639  *  N1 N2 N3 N4 N5 N6
640  *
641  *  où N1 N2 N3 sont les sommets du haut et N4 N5 N6 les sommets du bas
642  *  (selon l'orientation donnée par le sens des triangles)
643  */
644 void MESHCUT::cas3(int VN[6], int it4)
645 {
646   cutTetras.push_back(it4);
647   // cout << "Cas 3 - it4="<<it4<<", VN = " << VN[0] << " " << VN[1] << " " << VN[2] << " " << VN[3] << " " << VN[4] << " " << VN[5] << " " << endl;
648
649   // Numéros des noeuds du TETRA4
650   int ng0 = MAILLAGE1->CNX[TETRA4][4 * it4 + 0];
651   int ng1 = MAILLAGE1->CNX[TETRA4][4 * it4 + 1];
652   int ng2 = MAILLAGE1->CNX[TETRA4][4 * it4 + 2];
653   int ng3 = MAILLAGE1->CNX[TETRA4][4 * it4 + 3];
654
655   if (VN[0] != -1 && VN[1] != -1 && VN[2] != -1)
656     {
657
658       // Le sommet du nouveau T4 est le noeud 0 du T4 d'origine
659
660       newCNX[TETRA4].push_back(ng0);
661       newCNX[TETRA4].push_back(VN[0]);
662       newCNX[TETRA4].push_back(VN[1]);
663       newCNX[TETRA4].push_back(VN[2]);
664       cptNouvellesMailles[TETRA4]++;
665
666       newCNX[PENTA6].push_back(ng1);
667       newCNX[PENTA6].push_back(ng3);
668       newCNX[PENTA6].push_back(ng2);
669       newCNX[PENTA6].push_back(VN[0]);
670       newCNX[PENTA6].push_back(VN[2]);
671       newCNX[PENTA6].push_back(VN[1]);
672       cptNouvellesMailles[PENTA6]++;
673
674       if (POSN[ng0 - 1] == 1)
675         {
676           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
677           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
678         }
679       else
680         {
681           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
682           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
683         }
684     }
685
686   else if (VN[0] != -1 && VN[3] != -1 && VN[4] != -1)
687     {
688       // Le sommet du nouveau T4 est le noeud 1 du T4 d'origine
689
690       //cout << "cas 3.2 (noeud sommet 1)" << endl;
691       newCNX[TETRA4].push_back(VN[0]);
692       newCNX[TETRA4].push_back(ng1);
693       newCNX[TETRA4].push_back(VN[3]);
694       newCNX[TETRA4].push_back(VN[4]);
695       cptNouvellesMailles[TETRA4]++;
696
697       newCNX[PENTA6].push_back(ng0);
698       newCNX[PENTA6].push_back(ng2);
699       newCNX[PENTA6].push_back(ng3);
700       newCNX[PENTA6].push_back(VN[0]);
701       newCNX[PENTA6].push_back(VN[3]);
702       newCNX[PENTA6].push_back(VN[4]);
703       cptNouvellesMailles[PENTA6]++;
704
705       if (POSN[ng1 - 1] == 1)
706         {
707           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
708           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
709         }
710       else
711         {
712           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
713           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
714         }
715     }
716   else if (VN[1] != -1 && VN[3] != -1 && VN[5] != -1)
717     {
718       // Le sommet du nouveau T4 est le noeud 2 du T4 d'origine
719
720       //cout << "cas 3.3 (noeud sommet 2)" << endl;
721
722       newCNX[TETRA4].push_back(VN[1]);
723       newCNX[TETRA4].push_back(VN[3]);
724       newCNX[TETRA4].push_back(ng2);
725       newCNX[TETRA4].push_back(VN[5]);
726       cptNouvellesMailles[TETRA4]++;
727
728       newCNX[PENTA6].push_back(ng0);
729       newCNX[PENTA6].push_back(ng3);
730       newCNX[PENTA6].push_back(ng1);
731       newCNX[PENTA6].push_back(VN[1]);
732       newCNX[PENTA6].push_back(VN[5]);
733       newCNX[PENTA6].push_back(VN[3]);
734       cptNouvellesMailles[PENTA6]++;
735
736       if (POSN[ng2 - 1] == 1)
737         {
738           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
739           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
740         }
741       else
742         {
743           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
744           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
745         }
746     }
747   else if (VN[2] != -1 && VN[4] != -1 && VN[5] != -1)
748     {
749       // Le sommet du nouveau T4 est le noeud 3 du T4 d'origine
750
751       newCNX[TETRA4].push_back(VN[2]);
752       newCNX[TETRA4].push_back(VN[4]);
753       newCNX[TETRA4].push_back(VN[5]);
754       newCNX[TETRA4].push_back(ng3);
755       cptNouvellesMailles[TETRA4]++;
756
757       newCNX[PENTA6].push_back(ng0);
758       newCNX[PENTA6].push_back(ng1);
759       newCNX[PENTA6].push_back(ng2);
760       newCNX[PENTA6].push_back(VN[2]);
761       newCNX[PENTA6].push_back(VN[4]);
762       newCNX[PENTA6].push_back(VN[5]);
763       cptNouvellesMailles[PENTA6]++;
764
765       if (POSN[ng3 - 1] == 1)
766         {
767           GMplus[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
768           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
769         }
770       else
771         {
772           GMmoins[TETRA4].push_back(MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1);
773           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1);
774         }
775     }
776   else
777     {
778       //      cout << "Cas 3 - it4=" << it4 << ", VN = " << VN[0] << " " << VN[1] << " " << VN[2] << " " << VN[3] << " "
779       //          << VN[4] << " " << VN[5] << " " << endl;
780       //
781       //      int n0 = *(MAILLAGE1->CNX[TETRA4] + it4 * 4 + 0);
782       //      float x0 = MAILLAGE1->XX[n0 - 1];
783       //      float y0 = MAILLAGE1->YY[n0 - 1];
784       //      float z0 = MAILLAGE1->ZZ[n0 - 1];
785       //
786       //      int n1 = *(MAILLAGE1->CNX[TETRA4] + it4 * 4 + 1);
787       //      float x1 = MAILLAGE1->XX[n1 - 1];
788       //      float y1 = MAILLAGE1->YY[n1 - 1];
789       //      float z1 = MAILLAGE1->ZZ[n1 - 1];
790       //
791       //      int n2 = *(MAILLAGE1->CNX[TETRA4] + it4 * 4 + 2);
792       //      float x2 = MAILLAGE1->XX[n2 - 1];
793       //      float y2 = MAILLAGE1->YY[n2 - 1];
794       //      float z2 = MAILLAGE1->ZZ[n2 - 1];
795       //
796       //      int n3 = *(MAILLAGE1->CNX[TETRA4] + it4 * 4 + 3);
797       //      float x3 = MAILLAGE1->XX[n3 - 1];
798       //      float y3 = MAILLAGE1->YY[n3 - 1];
799       //      float z3 = MAILLAGE1->ZZ[n3 - 1];
800       //
801       //      cout << x0 << " " << y0 << " " << z0 << " " << endl;
802       //      cout << x1 << " " << y1 << " " << z1 << " " << endl;
803       //      cout << x2 << " " << y2 << " " << z2 << " " << endl;
804       //      cout << x3 << " " << y3 << " " << z3 << " " << endl;
805
806       ERREUR(" Intersections configuration not taken into account (case cptPI=3) ");
807     }
808
809   // int ngT4 = MAILLAGE1->EFFECTIFS_TYPES[TETRA4] + cptNouvellesMailles[TETRA4] - 1;
810   // int ngP6 = MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + cptNouvellesMailles[PENTA6] - 1;
811   // cout << "La maille TETRA4 " << it4 << " produit les mailles TETRA4-" << ngT4 << " et PENTA6-" << ngP6 << endl;
812 }
813
814 /*! Quatre points d'intersection
815  *
816  *  ATTENTION: pour les PENTA6 on adopte la convention d'orientation SALOME
817  *
818  *  N1 N2 N3 N4 N5 N6
819  *
820  *  où N1 N2 N3 sont les sommets du haut et N4 N5 N6 les sommets du bas
821  *  (selon l'orientation donnée par le sens des triangles)
822  */
823 void MESHCUT::cas4(int VN[6], int it4)
824 {
825   cutTetras.push_back(it4);
826   //  cout << "Cas 4 - it4=" << it4 << ", VN = " << VN[0] << " " << VN[1] << " " << VN[2] << " " << VN[3] << " " << VN[4]
827   //      << " " << VN[5] << " " << endl;
828
829   // Numéros des noeuds du TETRA4
830   int ng0 = MAILLAGE1->CNX[TETRA4][4 * it4 + 0];
831   int ng1 = MAILLAGE1->CNX[TETRA4][4 * it4 + 1];
832   int ng2 = MAILLAGE1->CNX[TETRA4][4 * it4 + 2];
833   int ng3 = MAILLAGE1->CNX[TETRA4][4 * it4 + 3];
834
835   int i1, i2; // Numéros locaux dans le type des mailles créées
836
837   if (VN[0] == -1 && VN[5] == -1)
838     {
839       // Les deux arêtes opposées sont [0,1] et [2,3]
840
841       newCNX[PENTA6].push_back(ng0);
842       newCNX[PENTA6].push_back(VN[1]);
843       newCNX[PENTA6].push_back(VN[2]);
844       newCNX[PENTA6].push_back(ng1);
845       newCNX[PENTA6].push_back(VN[3]);
846       newCNX[PENTA6].push_back(VN[4]);
847       cptNouvellesMailles[PENTA6]++;
848       i1 = cptNouvellesMailles[PENTA6] - 1;
849
850       newCNX[PENTA6].push_back(ng3);
851       newCNX[PENTA6].push_back(VN[4]);
852       newCNX[PENTA6].push_back(VN[2]);
853       newCNX[PENTA6].push_back(ng2);
854       newCNX[PENTA6].push_back(VN[3]);
855       newCNX[PENTA6].push_back(VN[1]);
856       cptNouvellesMailles[PENTA6]++;
857       i2 = cptNouvellesMailles[PENTA6] - 1;
858
859       if (POSN[ng0 - 1] == 1)
860         {
861           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1);
862           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2);
863         }
864       else
865         {
866           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2);
867           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1);
868         }
869
870     } // if ( VN[0]==-1 && VN[5]==-1 )
871
872   else if (VN[1] == -1 && VN[4] == -1)
873     {
874       // Les deux arêtes opposées sont [0,2] et [1,3]
875
876       newCNX[PENTA6].push_back(ng2);
877       newCNX[PENTA6].push_back(VN[3]);
878       newCNX[PENTA6].push_back(VN[5]);
879       newCNX[PENTA6].push_back(ng0);
880       newCNX[PENTA6].push_back(VN[0]);
881       newCNX[PENTA6].push_back(VN[2]);
882       cptNouvellesMailles[PENTA6]++;
883       i1 = cptNouvellesMailles[PENTA6] - 1;
884
885       newCNX[PENTA6].push_back(ng1);
886       newCNX[PENTA6].push_back(VN[3]);
887       newCNX[PENTA6].push_back(VN[0]);
888       newCNX[PENTA6].push_back(ng3);
889       newCNX[PENTA6].push_back(VN[5]);
890       newCNX[PENTA6].push_back(VN[2]);
891       cptNouvellesMailles[PENTA6]++;
892       i2 = cptNouvellesMailles[PENTA6] - 1;
893
894       if (POSN[ng0 - 1] == 1)
895         {
896           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1);
897           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2);
898         }
899       else
900         {
901           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2);
902           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1);
903         }
904     }
905
906   else if (VN[2] == -1 && VN[3] == -1)
907     {
908       // Les deux arêtes opposées sont [0,3] et [1,2]
909
910       newCNX[PENTA6].push_back(ng0);
911       newCNX[PENTA6].push_back(VN[0]);
912       newCNX[PENTA6].push_back(VN[1]);
913       newCNX[PENTA6].push_back(ng3);
914       newCNX[PENTA6].push_back(VN[4]);
915       newCNX[PENTA6].push_back(VN[5]);
916       cptNouvellesMailles[PENTA6]++;
917       i1 = cptNouvellesMailles[PENTA6] - 1;
918
919       newCNX[PENTA6].push_back(ng2);
920       newCNX[PENTA6].push_back(VN[5]);
921       newCNX[PENTA6].push_back(VN[1]);
922       newCNX[PENTA6].push_back(ng1);
923       newCNX[PENTA6].push_back(VN[4]);
924       newCNX[PENTA6].push_back(VN[0]);
925       cptNouvellesMailles[PENTA6]++;
926       i2 = cptNouvellesMailles[PENTA6] - 1;
927
928       if (POSN[ng0 - 1] == 1)
929         {
930           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1);
931           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2);
932         }
933       else
934         {
935           GMplus[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2);
936           GMmoins[PENTA6].push_back(MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1);
937         }
938     }
939   else
940     ERREUR(" Intersection configuration not taken into account (case cptPI=4) ");
941
942   // int nl1 = MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i1;
943   // int nl2 = MAILLAGE1->EFFECTIFS_TYPES[PENTA6] + i2;
944   // cout << "La maille TETRA4 " << it4 << " produit les mailles PENTA6-" << nl1 << " et PENTA6-" << nl2 << endl;
945 }