Salome HOME
PR: merged from V5_1_4rc1
[modules/smesh.git] / src / StdMeshers / StdMeshers_Hexa_3D.hxx
1 //  Copyright (C) 2007-2010  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 SMESH : implementaion of SMESH idl descriptions
24 //  File   : StdMeshers_Hexa_3D.hxx
25 //           Moved here from SMESH_Hexa_3D.hxx
26 //  Author : Paul RASCLE, EDF
27 //  Module : SMESH
28 //
29 #ifndef _SMESH_HEXA_3D_HXX_
30 #define _SMESH_HEXA_3D_HXX_
31
32 #include "SMESH_StdMeshers.hxx"
33
34 #include "SMESH_3D_Algo.hxx"
35 #include "SMESH_Mesh.hxx"
36 #include "StdMeshers_Quadrangle_2D.hxx"
37 #include "Utils_SALOME_Exception.hxx"
38
39 #include "SMESH_MesherHelper.hxx"
40
41 class TopTools_IndexedMapOfShape;
42
43 typedef struct point3Dstruct
44 {
45   const SMDS_MeshNode * node;
46 } Point3DStruct;
47
48 typedef double Pt3[3];
49
50 typedef struct conv2dstruct
51 {
52   double a1; // X = a1*x + b1*y + c1 
53   double b1; // Y = a2*x + b2*y + c2
54   double c1; // a1, b1 a2, b2 in {-1,0,1}
55   double a2; // c1, c2 in {0,1}
56   double b2;
57   double c2;
58   int ia;    // I = ia*i + ib*j + ic
59   int ib;
60   int ic;
61   int ja;    // J = ja*i + jb*j + jc
62   int jb;
63   int jc;
64 } Conv2DStruct;
65
66 class STDMESHERS_EXPORT StdMeshers_Hexa_3D:
67   public SMESH_3D_Algo
68 {
69 public:
70   StdMeshers_Hexa_3D(int hypId, int studyId, SMESH_Gen* gen);
71   virtual ~StdMeshers_Hexa_3D();
72
73   virtual bool CheckHypothesis(SMESH_Mesh& aMesh,
74                                const TopoDS_Shape& aShape,
75                                SMESH_Hypothesis::Hypothesis_Status& aStatus);
76
77   virtual bool Compute(SMESH_Mesh& aMesh,  const TopoDS_Shape& aShape);
78
79   virtual bool Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHelper);
80
81   virtual bool Evaluate(SMESH_Mesh & aMesh, const TopoDS_Shape & aShape,
82                         MapShapeNbElems& aResMap);
83
84   static TopoDS_Vertex OppositeVertex(const TopoDS_Vertex& aVertex,
85                                       const TopTools_IndexedMapOfShape& aQuads0Vertices,
86                                       FaceQuadStruct* aQuads[6]);
87
88 protected:
89   TopoDS_Edge
90   EdgeNotInFace(SMESH_Mesh& aMesh,
91                 const TopoDS_Shape& aShape,
92                 const TopoDS_Face& aFace,
93                 const TopoDS_Vertex& aVertex,
94                 const TopTools_IndexedDataMapOfShapeListOfShape& MS);
95
96   int GetFaceIndex(SMESH_Mesh& aMesh,
97                    const TopoDS_Shape& aShape,
98                    const std::vector<SMESH_subMesh*>& meshFaces,
99                    const TopoDS_Vertex& V0,
100                    const TopoDS_Vertex& V1,
101                    const TopoDS_Vertex& V2,
102                    const TopoDS_Vertex& V3);
103
104   void GetConv2DCoefs(const faceQuadStruct& quad,
105                       const TopoDS_Shape& aShape,
106                       const TopoDS_Vertex& V0,
107                       const TopoDS_Vertex& V1,
108                       const TopoDS_Vertex& V2,
109                       const TopoDS_Vertex& V3,
110                       Conv2DStruct& conv);
111
112   void GetPoint(Pt3 p,
113                 int i, int j, int k,
114                 int nbx, int nby, int nbz,
115                 Point3DStruct *np,
116                 const SMESHDS_Mesh* meshDS);
117
118   bool ClearAndReturn(FaceQuadStruct* theQuads[6], const bool res);
119 };
120
121 #endif