OCCT

// 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>

//=================================================================================================

Handle(IMeshData::ListOfPnt2d) BRepMesh_SphereRangeSplitter::GenerateSurfaceNodes(
  const IMeshTools_Parameters& theParameters) const
{
  // Calculate parameters for iteration in V direction
  double aStep =
    0.7
    * GCPnts_TangentialDeflection::ArcAngularStep(GetDFace()->GetSurface()->Sphere().Radius(),
                                                  GetDFace()->GetDeflection(),
                                                  theParameters.Angle,
                                                  theParameters.MinSize);

  const std::pair<double, double>* aRange[2] = {&GetRangeV(), &GetRangeU()};

  std::pair<double, double> aStepAndOffset[2];
  computeStep(*aRange[0], aStep, aStepAndOffset[0]);
  computeStep(*aRange[1], aStep, aStepAndOffset[1]);

  const occ::handle<NCollection_IncAllocator> aTmpAlloc =
    new NCollection_IncAllocator(IMeshData::MEMORY_BLOCK_SIZE_HUGE);
  Handle(IMeshData::ListOfPnt2d) aNodes = new IMeshData::ListOfPnt2d(aTmpAlloc);

  const double aHalfDu = aStepAndOffset[1].first * 0.5;
  bool         Shift   = false;
  double       aPasV   = aRange[0]->first + aStepAndOffset[0].first;
  for (; aPasV < aStepAndOffset[0].second; aPasV += aStepAndOffset[0].first)
  {
    Shift              = !Shift;
    const double d     = (Shift) ? aHalfDu : 0.;
    double       aPasU = aRange[1]->first + d;
    for (; aPasU < aStepAndOffset[1].second; aPasU += aStepAndOffset[1].first)
    {
      aNodes->Append(gp_Pnt2d(aPasU, aPasV));
    }
  }

  return aNodes;
}