0031939: Coding - correction of spelling errors in comments
[occt.git] / src / BRepMesh / BRepMesh_BaseMeshAlgo.cxx
1 // Created on: 2016-04-19
2 // Copyright (c) 2016 OPEN CASCADE SAS
3 // Created by: Oleg AGASHIN
4 //
5 // This file is part of Open CASCADE Technology software library.
6 //
7 // This library is free software; you can redistribute it and/or modify it under
8 // the terms of the GNU Lesser General Public License version 2.1 as published
9 // by the Free Software Foundation, with special exception defined in the file
10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
11 // distribution for complete text of the license and disclaimer of any warranty.
12 //
13 // Alternatively, this file may be used under the terms of Open CASCADE
14 // commercial license or contractual agreement.
15
16 #include <BRepMesh_BaseMeshAlgo.hxx>
17 #include <BRepMesh_DataStructureOfDelaun.hxx>
18 #include <IMeshData_Face.hxx>
19 #include <IMeshData_Wire.hxx>
20 #include <IMeshData_Edge.hxx>
21 #include <IMeshData_PCurve.hxx>
22 #include <IMeshData_Curve.hxx>
23 #include <BRepMesh_Delaun.hxx>
24 #include <BRepMesh_ShapeTool.hxx>
25 #include <Standard_ErrorHandler.hxx>
26
27 IMPLEMENT_STANDARD_RTTIEXT(BRepMesh_BaseMeshAlgo, IMeshTools_MeshAlgo)
28
29 //=======================================================================
30 // Function: Constructor
31 // Purpose : 
32 //=======================================================================
33 BRepMesh_BaseMeshAlgo::BRepMesh_BaseMeshAlgo()
34 {
35 }
36
37 //=======================================================================
38 // Function: Destructor
39 // Purpose : 
40 //=======================================================================
41 BRepMesh_BaseMeshAlgo::~BRepMesh_BaseMeshAlgo()
42 {
43 }
44
45 //=======================================================================
46 // Function: Perform
47 // Purpose : 
48 //=======================================================================
49 void BRepMesh_BaseMeshAlgo::Perform(
50   const IMeshData::IFaceHandle& theDFace,
51   const IMeshTools_Parameters&  theParameters,
52   const Message_ProgressRange&  theRange)
53 {
54   try
55   {
56     OCC_CATCH_SIGNALS
57
58     myDFace      = theDFace;
59     myParameters = theParameters;
60     myAllocator  = new NCollection_IncAllocator(IMeshData::MEMORY_BLOCK_SIZE_HUGE);
61     myStructure  = new BRepMesh_DataStructureOfDelaun(myAllocator);
62     myNodesMap   = new VectorOfPnt(256, myAllocator);
63     myUsedNodes  = new DMapOfIntegerInteger(1, myAllocator);
64
65     if (initDataStructure())
66     {
67       if (!theRange.More())
68       {
69         return;
70       }
71       generateMesh(theRange);
72       commitSurfaceTriangulation();
73     }
74   }
75   catch (Standard_Failure const& /*theException*/)
76   {
77   }
78
79   myDFace.Nullify(); // Do not hold link to face.
80   myStructure.Nullify();
81   myNodesMap .Nullify();
82   myUsedNodes.Nullify();
83   myAllocator.Nullify();
84 }
85
86 //=======================================================================
87 //function : initDataStructure
88 //purpose  :
89 //=======================================================================
90 Standard_Boolean BRepMesh_BaseMeshAlgo::initDataStructure()
91 {
92   for (Standard_Integer aWireIt = 0; aWireIt < myDFace->WiresNb(); ++aWireIt)
93   {
94     const IMeshData::IWireHandle& aDWire = myDFace->GetWire(aWireIt);
95     if (aDWire->IsSet(IMeshData_SelfIntersectingWire))
96     {
97       // TODO: here we can add points of self-intersecting wire as fixed points
98       // in order to keep consistency of nodes with adjacent faces.
99       continue;
100     }
101
102     for (Standard_Integer aEdgeIt = 0; aEdgeIt < aDWire->EdgesNb(); ++aEdgeIt)
103     {
104       const IMeshData::IEdgeHandle    aDEdge = aDWire->GetEdge(aEdgeIt);
105       const IMeshData::ICurveHandle&  aCurve = aDEdge->GetCurve();
106       const IMeshData::IPCurveHandle& aPCurve = aDEdge->GetPCurve(
107         myDFace.get(), aDWire->GetEdgeOrientation(aEdgeIt));
108
109       const TopAbs_Orientation aOri = fixSeamEdgeOrientation(aDEdge, aPCurve);
110
111       Standard_Integer aPrevNodeIndex = -1;
112       const Standard_Integer aLastPoint = aPCurve->ParametersNb() - 1;
113       for (Standard_Integer aPointIt = 0; aPointIt <= aLastPoint; ++aPointIt)
114       {
115         const Standard_Integer aNodeIndex = registerNode(
116           aCurve ->GetPoint(aPointIt),
117           aPCurve->GetPoint(aPointIt),
118           BRepMesh_Frontier, Standard_False/*aPointIt > 0 && aPointIt < aLastPoint*/);
119
120         aPCurve->GetIndex(aPointIt) = aNodeIndex;
121         myUsedNodes->Bind(aNodeIndex, aNodeIndex);
122
123         if (aPrevNodeIndex != -1 && aPrevNodeIndex != aNodeIndex)
124         {
125           const Standard_Integer aLinksNb   = myStructure->NbLinks();
126           const Standard_Integer aLinkIndex = addLinkToMesh(aPrevNodeIndex, aNodeIndex, aOri);
127           if (aWireIt != 0 && aLinkIndex <= aLinksNb)
128           {
129             // Prevent holes around wire of zero area.
130             BRepMesh_Edge& aLink = const_cast<BRepMesh_Edge&>(myStructure->GetLink(aLinkIndex));
131             aLink.SetMovability(BRepMesh_Fixed);
132           }
133         }
134
135         aPrevNodeIndex = aNodeIndex;
136       }
137     }
138   }
139
140   return Standard_True;
141 }
142
143 //=======================================================================
144 // Function: registerNode
145 // Purpose : 
146 //=======================================================================
147 Standard_Integer BRepMesh_BaseMeshAlgo::registerNode(
148   const gp_Pnt&                  thePoint,
149   const gp_Pnt2d&                thePoint2d,
150   const BRepMesh_DegreeOfFreedom theMovability,
151   const Standard_Boolean         isForceAdd)
152 {
153   const Standard_Integer aNodeIndex = addNodeToStructure(
154     thePoint2d, myNodesMap->Size(), theMovability, isForceAdd);
155
156   if (aNodeIndex > myNodesMap->Size())
157   {
158     myNodesMap->Append(thePoint);
159   }
160
161   return aNodeIndex;
162 }
163
164 //=======================================================================
165 // Function: addNode
166 // Purpose : 
167 //=======================================================================
168 Standard_Integer BRepMesh_BaseMeshAlgo::addNodeToStructure(
169   const gp_Pnt2d&                thePoint,
170   const Standard_Integer         theLocation3d,
171   const BRepMesh_DegreeOfFreedom theMovability,
172   const Standard_Boolean         isForceAdd)
173 {
174   BRepMesh_Vertex aNode(thePoint.XY(), theLocation3d, theMovability);
175   return myStructure->AddNode(aNode, isForceAdd);
176 }
177
178 //=======================================================================
179 //function : addLinkToMesh
180 //purpose  :
181 //=======================================================================
182 Standard_Integer BRepMesh_BaseMeshAlgo::addLinkToMesh(
183   const Standard_Integer   theFirstNodeId,
184   const Standard_Integer   theLastNodeId,
185   const TopAbs_Orientation theOrientation)
186 {
187   Standard_Integer aLinkIndex;
188   if (theOrientation == TopAbs_REVERSED)
189     aLinkIndex = myStructure->AddLink(BRepMesh_Edge(theLastNodeId, theFirstNodeId, BRepMesh_Frontier));
190   else if (theOrientation == TopAbs_INTERNAL)
191     aLinkIndex = myStructure->AddLink(BRepMesh_Edge(theFirstNodeId, theLastNodeId, BRepMesh_Fixed));
192   else
193     aLinkIndex = myStructure->AddLink(BRepMesh_Edge(theFirstNodeId, theLastNodeId, BRepMesh_Frontier));
194
195   return Abs(aLinkIndex);
196 }
197
198 //=======================================================================
199 //function : fixSeamEdgeOrientation
200 //purpose  :
201 //=======================================================================
202 TopAbs_Orientation BRepMesh_BaseMeshAlgo::fixSeamEdgeOrientation(
203   const IMeshData::IEdgeHandle&   theDEdge,
204   const IMeshData::IPCurveHandle& thePCurve) const
205 {
206   for (Standard_Integer aPCurveIt = 0; aPCurveIt < theDEdge->PCurvesNb(); ++aPCurveIt)
207   {
208     const IMeshData::IPCurveHandle& aPCurve = theDEdge->GetPCurve(aPCurveIt);
209     if (aPCurve->GetFace() == myDFace && thePCurve != aPCurve)
210     {
211       // Simple check that another pcurve of seam edge does not coincide with reference one.
212       const gp_Pnt2d& aPnt1_1 = thePCurve->GetPoint(0);
213       const gp_Pnt2d& aPnt2_1 = thePCurve->GetPoint(thePCurve->ParametersNb() - 1);
214
215       const gp_Pnt2d& aPnt1_2 = aPCurve->GetPoint(0);
216       const gp_Pnt2d& aPnt2_2 = aPCurve->GetPoint(aPCurve->ParametersNb() - 1);
217
218       const Standard_Real aSqDist1 = Min(aPnt1_1.SquareDistance(aPnt1_2), aPnt1_1.SquareDistance(aPnt2_2));
219       const Standard_Real aSqDist2 = Min(aPnt2_1.SquareDistance(aPnt1_2), aPnt2_1.SquareDistance(aPnt2_2));
220       if (aSqDist1 < Precision::SquareConfusion() &&
221           aSqDist2 < Precision::SquareConfusion())
222       {
223         return TopAbs_INTERNAL;
224       }
225     }
226   }
227
228   return thePCurve->GetOrientation();
229 }
230
231 //=======================================================================
232 //function : commitSurfaceTriangulation
233 //purpose  :
234 //=======================================================================
235 void BRepMesh_BaseMeshAlgo::commitSurfaceTriangulation()
236 {
237   Handle(Poly_Triangulation) aTriangulation = collectTriangles();
238   if (aTriangulation.IsNull())
239   {
240     myDFace->SetStatus(IMeshData_Failure);
241     return;
242   }
243
244   collectNodes(aTriangulation);
245
246   aTriangulation->Deflection(myDFace->GetDeflection());
247   BRepMesh_ShapeTool::AddInFace(myDFace->GetFace(), aTriangulation);
248 }
249
250 //=======================================================================
251 //function : collectTriangles
252 //purpose  :
253 //=======================================================================
254 Handle(Poly_Triangulation) BRepMesh_BaseMeshAlgo::collectTriangles()
255 {
256   const IMeshData::MapOfInteger& aTriangles = myStructure->ElementsOfDomain();
257   if (aTriangles.IsEmpty())
258   {
259     return Handle(Poly_Triangulation)();
260   }
261
262   Poly_Array1OfTriangle aPolyTrianges(1, aTriangles.Extent());
263   IMeshData::IteratorOfMapOfInteger aTriIt(aTriangles);
264   for (Standard_Integer aTriangeId = 1; aTriIt.More(); aTriIt.Next(), ++aTriangeId)
265   {
266     const BRepMesh_Triangle& aCurElem = myStructure->GetElement(aTriIt.Key());
267
268     Standard_Integer aNode[3];
269     myStructure->ElementNodes(aCurElem, aNode);
270
271     for (Standard_Integer i = 0; i < 3; ++i)
272     {
273       if (!myUsedNodes->IsBound(aNode[i]))
274       {
275         myUsedNodes->Bind(aNode[i], myUsedNodes->Size() + 1);
276       }
277
278       aNode[i] = myUsedNodes->Find(aNode[i]);
279     }
280
281     aPolyTrianges(aTriangeId).Set(aNode[0], aNode[1], aNode[2]);
282   }
283
284   Handle(Poly_Triangulation) aTriangulation = new Poly_Triangulation(
285     myUsedNodes->Extent(), aTriangles.Extent(), Standard_True);
286
287   aTriangulation->ChangeTriangles() = aPolyTrianges;
288   return aTriangulation;
289 }
290
291 //=======================================================================
292 //function : collectNodes
293 //purpose  :
294 //=======================================================================
295 void BRepMesh_BaseMeshAlgo::collectNodes(
296   const Handle(Poly_Triangulation)& theTriangulation)
297 {
298   // Store mesh nodes
299   TColgp_Array1OfPnt&   aNodes   = theTriangulation->ChangeNodes();
300   TColgp_Array1OfPnt2d& aNodes2d = theTriangulation->ChangeUVNodes();
301
302   for (Standard_Integer i = 1; i <= myNodesMap->Size(); ++i)
303   {
304     if (myUsedNodes->IsBound(i))
305     {
306       const BRepMesh_Vertex& aVertex = myStructure->GetNode(i);
307
308       const Standard_Integer aNodeIndex = myUsedNodes->Find(i);
309       aNodes(aNodeIndex) = myNodesMap->Value(aVertex.Location3d());
310       aNodes2d(aNodeIndex) = getNodePoint2d(aVertex);
311     }
312   }
313 }
314
315 //=======================================================================
316 // Function: getNodePoint2d
317 // Purpose : 
318 //=======================================================================
319 gp_Pnt2d BRepMesh_BaseMeshAlgo::getNodePoint2d(
320   const BRepMesh_Vertex& theVertex) const
321 {
322   return theVertex.Coord();
323 }