0026106: BRepMesh - revision of data model
[occt.git] / src / BRepMesh / BRepMesh_SphereRangeSplitter.cxx
diff --git a/src/BRepMesh/BRepMesh_SphereRangeSplitter.cxx b/src/BRepMesh/BRepMesh_SphereRangeSplitter.cxx
new file mode 100644 (file)
index 0000000..4a3112e
--- /dev/null
@@ -0,0 +1,59 @@
+// Created on: 2016-07-07
+// Copyright (c) 2016 OPEN CASCADE SAS
+// Created by: Oleg AGASHIN
+//
+// This file is part of Open CASCADE Technology software library.
+//
+// This library is free software; you can redistribute it and/or modify it under
+// the terms of the GNU Lesser General Public License version 2.1 as published
+// by the Free Software Foundation, with special exception defined in the file
+// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
+// distribution for complete text of the license and disclaimer of any warranty.
+//
+// Alternatively, this file may be used under the terms of Open CASCADE
+// commercial license or contractual agreement.
+
+#include <BRepMesh_SphereRangeSplitter.hxx>
+#include <GCPnts_TangentialDeflection.hxx>
+
+//=======================================================================
+// Function: GenerateSurfaceNodes
+// Purpose : 
+//=======================================================================
+Handle(IMeshData::ListOfPnt2d) BRepMesh_SphereRangeSplitter::GenerateSurfaceNodes(
+  const IMeshTools_Parameters& theParameters) const
+{
+  // Calculate parameters for iteration in V direction
+  Standard_Real aStep = 0.7 * GCPnts_TangentialDeflection::ArcAngularStep(
+    GetDFace()->GetSurface()->Sphere().Radius(), GetDFace()->GetDeflection(),
+    theParameters.Angle, theParameters.MinSize);
+
+  const std::pair<Standard_Real, Standard_Real>* aRange[2] = {
+    &GetRangeV(),
+    &GetRangeU()
+  };
+
+  std::pair<Standard_Real, Standard_Real> aStepAndOffset[2];
+  computeStep(*aRange[0], aStep, aStepAndOffset[0]);
+  computeStep(*aRange[1], aStep, aStepAndOffset[1]);
+
+  const Handle(NCollection_IncAllocator) aTmpAlloc =
+    new NCollection_IncAllocator(IMeshData::MEMORY_BLOCK_SIZE_HUGE);
+  Handle(IMeshData::ListOfPnt2d) aNodes = new IMeshData::ListOfPnt2d(aTmpAlloc);
+
+  const Standard_Real aHalfDu = aStepAndOffset[1].first * 0.5;
+  Standard_Boolean Shift = Standard_False;
+  Standard_Real aPasV = aRange[0]->first + aStepAndOffset[0].first;
+  for (; aPasV < aStepAndOffset[0].second; aPasV += aStepAndOffset[0].first)
+  {
+    Shift = !Shift;
+    const Standard_Real d = (Shift) ? aHalfDu : 0.;
+    Standard_Real aPasU = aRange[1]->first + d;
+    for (; aPasU < aStepAndOffset[1].second; aPasU += aStepAndOffset[1].first)
+    {
+      aNodes->Append(gp_Pnt2d(aPasU, aPasV));
+    }
+  }
+
+  return aNodes;
+}