]> SALOME platform Git repositories - modules/smesh.git/blob - src/StdMeshers/StdMeshers_Penta_3D.hxx
Salome HOME
speed up NotifySubMeshesHypothesisModification()
[modules/smesh.git] / src / StdMeshers / StdMeshers_Penta_3D.hxx
1 // Copyright (C) 2007-2011  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.
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 StdMeshers : implementaion of SMESH idl descriptions
24 //  File   : StdMeshers_Penta_3D.hxx
25 //  Module : SMESH
26 //
27 #ifndef StdMeshers_Penta_3D_HeaderFile
28 #define StdMeshers_Penta_3D_HeaderFile
29
30 #include "SMESH_StdMeshers.hxx"
31
32 #include <map>
33
34 ////////////////////////////////////////////////////////////////////////
35 //
36 //  class StdMeshers_SMESHBlock
37 //
38 ////////////////////////////////////////////////////////////////////////
39 #include <gp_Pnt.hxx>
40 #include <gp_XYZ.hxx>
41 #include <TopoDS_Vertex.hxx>
42 #include <TopoDS_Shell.hxx>
43 #include <TopTools_IndexedMapOfOrientedShape.hxx>
44 #include <TColStd_MapOfInteger.hxx>
45
46 #include "SMESH_Block.hxx"
47 #include "SMESH_ComputeError.hxx"
48 #include "SMESH_MesherHelper.hxx"
49 #include "SMESH_3D_Algo.hxx"
50
51 typedef std::map< double, std::vector<const SMDS_MeshNode*> > StdMeshers_IJNodeMap;
52
53 class STDMESHERS_EXPORT StdMeshers_SMESHBlock {
54  
55 public:
56   //
57   StdMeshers_SMESHBlock();
58   
59   void Load (const TopoDS_Shell& theShell);
60   
61   void Load (const TopoDS_Shell& theShell,
62              const TopoDS_Vertex& theV000,
63              const TopoDS_Vertex& theV001);
64   
65   void ComputeParameters(const gp_Pnt& thePnt, 
66                          gp_XYZ& theXYZ);
67   
68   void ComputeParameters(const gp_Pnt& thePnt,
69                          const TopoDS_Shape& theShape,
70                          gp_XYZ& theXYZ);
71   
72   void ComputeParameters(const double& theU,
73                          const TopoDS_Shape& theShape,
74                          gp_XYZ& theXYZ);
75   
76   void Point(const gp_XYZ& theParams, 
77              gp_Pnt& thePnt);
78   
79   void Point(const gp_XYZ& theParams,
80              const TopoDS_Shape& theShape, 
81              gp_Pnt& thePnt);
82   
83   int ShapeID(const TopoDS_Shape& theShape); 
84   
85   const TopoDS_Shape& Shape(const int theID);
86
87   SMESH_Block & Block() { return myTBlock; }
88
89   bool IsForwadEdge(const int theEdgeID);
90
91   int  ErrorStatus() const;
92
93   SMESH_ComputeErrorPtr GetError() const;
94
95
96 protected:
97   TopoDS_Shell                       myShell;
98   TopTools_IndexedMapOfOrientedShape myShapeIDMap;
99   SMESH_Block                        myTBlock;
100   TopoDS_Shape                       myEmptyShape;
101   std::vector<int>                   myIsEdgeForward;
102   //
103   int myErrorStatus;
104 };
105
106 ////////////////////////////////////////////////////////////////////////
107 //
108 //  class StdMeshers_TNode
109 //
110 ////////////////////////////////////////////////////////////////////////
111 #include "SMDS_MeshNode.hxx"
112
113 class STDMESHERS_EXPORT StdMeshers_TNode {
114
115 public:
116   
117   StdMeshers_TNode(){
118     myNode=NULL;
119     myXYZ.SetCoord(99., 99., 99.);
120     myShapeSupportID=-1;
121     myBaseNodeID=-1;
122   }
123   
124   void SetNode(const SMDS_MeshNode* theNode) {
125     myNode=(SMDS_MeshNode*) theNode;
126   }
127   
128   const SMDS_MeshNode* Node()const {
129     return myNode;
130   }
131
132   void SetShapeSupportID (const int theID) {
133     myShapeSupportID=theID;
134   }
135   
136   int ShapeSupportID()const {
137     return myShapeSupportID;
138   }
139   
140   void SetNormCoord (const gp_XYZ& theXYZ) {
141     myXYZ=theXYZ;
142   }
143
144   const gp_XYZ& NormCoord ()const{
145     return myXYZ;
146   }
147   
148   void SetBaseNodeID (const int theID) {
149     myBaseNodeID=theID;
150   }
151   
152   int BaseNodeID ()const{
153     return myBaseNodeID;
154   }
155
156 private:
157   SMDS_MeshNode* myNode;
158   int  myShapeSupportID;
159   gp_XYZ         myXYZ;
160   int            myBaseNodeID;
161 };
162
163 ////////////////////////////////////////////////////////////////////////
164 //
165 //  class StdMeshers_Penta_3D
166 //
167 ////////////////////////////////////////////////////////////////////////
168 #include "SMESH_Mesh.hxx"
169 #include <TopoDS_Shape.hxx>
170 //
171 class STDMESHERS_EXPORT StdMeshers_Penta_3D {
172 //
173   public: // methods
174     StdMeshers_Penta_3D();
175     
176     ~StdMeshers_Penta_3D();
177     
178     bool Compute(SMESH_Mesh& , const TopoDS_Shape& );
179     
180     int ErrorStatus() const {
181       if (myErrorStatus->IsOK())
182         return 0;
183       return myErrorStatus->myName;
184     }
185
186     SMESH_ComputeErrorPtr GetComputeError() const {
187       return myErrorStatus;
188     }
189    
190     void SetTolerance(const double theTol3D) {
191       myTol3D=theTol3D;
192     }
193     
194     double Tolerance() const {
195       return myTol3D;
196     }
197
198     bool LoadIJNodes(StdMeshers_IJNodeMap & theIJNodes,
199                      const TopoDS_Face&     theFace,
200                      const TopoDS_Edge&     theBaseEdge,
201                      SMESHDS_Mesh*          theMesh);
202     // Load nodes bound to theFace into column (vectors) and rows
203     // of theIJNodes.
204     // The value of theIJNodes map is a vector of ordered nodes so
205     // that the 0-the one lies on theBaseEdge.
206     // The key of theIJNodes map is a normalized parameter of each
207     // 0-the node on theBaseEdge.
208
209     bool Evaluate(SMESH_Mesh & aMesh, const TopoDS_Shape & aShape,
210                   MapShapeNbElems& aResMap);
211
212
213   protected: // methods
214     
215     void CheckData();
216     
217     void MakeBlock();
218
219     void MakeNodes();
220
221     double SetHorizEdgeXYZ(const gp_XYZ&                  aBNXYZ,
222                            const int                      aFaceID,
223                            std::vector<const SMDS_MeshNode*>*& aCol1,
224                            std::vector<const SMDS_MeshNode*>*& aCol2);
225
226     void ShapeSupportID(const bool theIsUpperLayer,
227                         const SMESH_Block::TShapeID theBNSSID,
228                         SMESH_Block::TShapeID& theSSID);
229
230     void FindNodeOnShape(const TopoDS_Shape& aS,
231                          const gp_XYZ& aParams,
232                          const int z,
233                          StdMeshers_TNode& aTN);
234
235     void CreateNode(const bool theIsUpperLayer,
236                     const gp_XYZ& aParams,
237                     StdMeshers_TNode& aTN);
238
239     void ClearMeshOnFxy1();
240
241     void MakeMeshOnFxy1();
242
243     void MakeConnectingMap();
244
245     int GetIndexOnLayer(const int aID);
246
247     void MakeVolumeMesh();
248   
249     void SetMesh(SMESH_Mesh& theMesh) {
250       myMesh=(void *)&theMesh;
251     }
252     
253     SMESH_Mesh* GetMesh()const {
254       return (SMESH_Mesh*)myMesh;
255     }
256     
257   protected: // fields
258     TopoDS_Shape              myShape;
259     StdMeshers_SMESHBlock     myBlock;
260     void *                    myMesh;
261     SMESH_ComputeErrorPtr     myErrorStatus;
262     //
263     std::vector <StdMeshers_TNode> myTNodes;
264     int                       myISize;
265     int                       myJSize;
266     double                    myTol3D;        // Tolerance value     
267     std::map < int, int >     myConnectingMap;
268     //
269     std::vector<StdMeshers_IJNodeMap> myWallNodesMaps; // nodes on a face
270     std::vector<gp_XYZ>               myShapeXYZ; // point on each sub-shape
271
272     bool myCreateQuadratic;
273     SMESH_MesherHelper* myTool; // tool building quadratic elements
274 };
275
276 #endif