0026106: BRepMesh - revision of data model
[occt.git] / src / BRepMesh / BRepMesh_SphereRangeSplitter.cxx
1 // Created on: 2016-07-07
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_SphereRangeSplitter.hxx>
17 #include <GCPnts_TangentialDeflection.hxx>
18
19 //=======================================================================
20 // Function: GenerateSurfaceNodes
21 // Purpose : 
22 //=======================================================================
23 Handle(IMeshData::ListOfPnt2d) BRepMesh_SphereRangeSplitter::GenerateSurfaceNodes(
24   const IMeshTools_Parameters& theParameters) const
25 {
26   // Calculate parameters for iteration in V direction
27   Standard_Real aStep = 0.7 * GCPnts_TangentialDeflection::ArcAngularStep(
28     GetDFace()->GetSurface()->Sphere().Radius(), GetDFace()->GetDeflection(),
29     theParameters.Angle, theParameters.MinSize);
30
31   const std::pair<Standard_Real, Standard_Real>* aRange[2] = {
32     &GetRangeV(),
33     &GetRangeU()
34   };
35
36   std::pair<Standard_Real, Standard_Real> aStepAndOffset[2];
37   computeStep(*aRange[0], aStep, aStepAndOffset[0]);
38   computeStep(*aRange[1], aStep, aStepAndOffset[1]);
39
40   const Handle(NCollection_IncAllocator) aTmpAlloc =
41     new NCollection_IncAllocator(IMeshData::MEMORY_BLOCK_SIZE_HUGE);
42   Handle(IMeshData::ListOfPnt2d) aNodes = new IMeshData::ListOfPnt2d(aTmpAlloc);
43
44   const Standard_Real aHalfDu = aStepAndOffset[1].first * 0.5;
45   Standard_Boolean Shift = Standard_False;
46   Standard_Real aPasV = aRange[0]->first + aStepAndOffset[0].first;
47   for (; aPasV < aStepAndOffset[0].second; aPasV += aStepAndOffset[0].first)
48   {
49     Shift = !Shift;
50     const Standard_Real d = (Shift) ? aHalfDu : 0.;
51     Standard_Real aPasU = aRange[1]->first + d;
52     for (; aPasU < aStepAndOffset[1].second; aPasU += aStepAndOffset[1].first)
53     {
54       aNodes->Append(gp_Pnt2d(aPasU, aPasV));
55     }
56   }
57
58   return aNodes;
59 }