Salome HOME
23080: [CEA 1497] Do not merge a middle node in quadratic with the extreme nodes...
[modules/smesh.git] / src / SMDS / SMDS_VolumeOfFaces.cxx
1 // Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 //  SMESH SMDS : implementaion of Salome mesh data structure
24 //  File   : SMDS_VolumeOfFaces.cxx
25 //  Author : Jean-Michel BOULCOURT
26 //  Module : SMESH
27 //
28 #ifdef _MSC_VER
29 #pragma warning(disable:4786)
30 #endif
31
32 #include "SMDS_VolumeOfFaces.hxx"
33 #include "SMDS_IteratorOfElements.hxx"
34 #include "utilities.h"
35
36 using namespace std;
37
38 //=======================================================================
39 //function : Print
40 //purpose  : 
41 //=======================================================================
42
43 void SMDS_VolumeOfFaces::Print(ostream & OS) const
44 {
45         OS << "volume <" << GetID() << "> : ";
46         int i;
47         for (i = 0; i < NbFaces()-1; ++i) OS << myFaces[i] << ",";
48         OS << myFaces[i]<< ") " << endl;
49 }
50
51
52 int SMDS_VolumeOfFaces::NbFaces() const
53 {
54         return myNbFaces;
55 }
56
57 class SMDS_VolumeOfFaces_MyIterator:public SMDS_ElemIterator
58 {
59   const SMDS_MeshFace* const *mySet;
60   int myLength;
61   int index;
62  public:
63   SMDS_VolumeOfFaces_MyIterator(const SMDS_MeshFace* const *s, int l):
64     mySet(s),myLength(l),index(0) {}
65
66   bool more()
67   {
68     return index<myLength;
69   }
70
71   const SMDS_MeshElement* next()
72   {
73     index++;
74     return mySet[index-1];
75   }
76 };
77
78 SMDS_ElemIteratorPtr SMDS_VolumeOfFaces::
79         elementsIterator(SMDSAbs_ElementType type) const
80 {
81   switch(type)
82   {
83   case SMDSAbs_Volume:
84     return SMDS_MeshElement::elementsIterator(SMDSAbs_Volume);
85   case SMDSAbs_Face:
86     return SMDS_ElemIteratorPtr(new SMDS_VolumeOfFaces_MyIterator(myFaces,myNbFaces));
87   default:
88     return SMDS_ElemIteratorPtr
89       (new SMDS_IteratorOfElements
90        (this,type,SMDS_ElemIteratorPtr
91         (new SMDS_VolumeOfFaces_MyIterator(myFaces,myNbFaces))));
92   }
93 }
94
95 SMDS_VolumeOfFaces::SMDS_VolumeOfFaces(const SMDS_MeshFace * face1,
96                                        const SMDS_MeshFace * face2,
97                                        const SMDS_MeshFace * face3,
98                                        const SMDS_MeshFace * face4)
99 {
100   //MESSAGE("****************************************************** SMDS_VolumeOfFaces");
101         myNbFaces = 4;
102         myFaces[0]=face1;
103         myFaces[1]=face2;
104         myFaces[2]=face3;
105         myFaces[3]=face4;
106         myFaces[4]=0;
107         myFaces[5]=0;
108 }
109
110 SMDS_VolumeOfFaces::SMDS_VolumeOfFaces(const SMDS_MeshFace * face1,
111                                        const SMDS_MeshFace * face2,
112                                        const SMDS_MeshFace * face3,
113                                        const SMDS_MeshFace * face4,
114                                        const SMDS_MeshFace * face5)
115 {
116   //MESSAGE("****************************************************** SMDS_VolumeOfFaces");
117         myNbFaces = 5;
118         myFaces[0]=face1;
119         myFaces[1]=face2;
120         myFaces[2]=face3;
121         myFaces[3]=face4;
122         myFaces[4]=face5;
123         myFaces[5]=0;
124 }
125
126 SMDS_VolumeOfFaces::SMDS_VolumeOfFaces(const SMDS_MeshFace * face1,
127                                        const SMDS_MeshFace * face2,
128                                        const SMDS_MeshFace * face3,
129                                        const SMDS_MeshFace * face4,
130                                        const SMDS_MeshFace * face5,
131                                        const SMDS_MeshFace * face6)
132 {
133   //MESSAGE("****************************************************** SMDS_VolumeOfFaces");
134         myNbFaces = 6;
135         myFaces[0]=face1;
136         myFaces[1]=face2;
137         myFaces[2]=face3;
138         myFaces[3]=face4;
139         myFaces[4]=face5;
140         myFaces[5]=face6;
141 }
142
143 SMDSAbs_EntityType SMDS_VolumeOfFaces::GetEntityType() const
144 {
145   SMDSAbs_EntityType aType = SMDSEntity_Tetra;
146   switch(myNbFaces)
147   {
148   case 4: aType = SMDSEntity_Tetra;   break;
149   case 5: aType = SMDSEntity_Pyramid; break;
150   case 6: aType = SMDSEntity_Penta;   break;
151   case 8:
152   default: aType = SMDSEntity_Hexa;    break;
153   }
154   return aType;
155 }
156
157 SMDSAbs_GeometryType SMDS_VolumeOfFaces::GetGeomType() const
158 {
159   SMDSAbs_GeometryType aType = SMDSGeom_NONE;
160   switch(myNbFaces)
161   {
162   case 4: aType = SMDSGeom_TETRA;   break;
163   case 5: aType = SMDSGeom_PYRAMID; break;
164   case 6: aType = SMDSGeom_PENTA;   break;
165   case 8:
166   default: aType = SMDSGeom_HEXA;   break;
167   }
168   return aType;
169 }