OCCT

// Created on: 1992-10-22
// Created by: Laurent BUCHARD
// Copyright (c) 1992-1999 Matra Datavision
// Copyright (c) 1999-2014 OPEN CASCADE SAS
//
// 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 <Geom2d_BezierCurve.hxx>
#include <Geom2dInt_Geom2dCurveTool.hxx>
#include <GeomAbs_CurveType.hxx>
#include <Precision.hxx>

//============================================================
int Geom2dInt_Geom2dCurveTool::NbSamples(const Adaptor2d_Curve2d& C,
                                         const double             U0,
                                         const double             U1)
{
  GeomAbs_CurveType typC = C.GetType();
  int               nbs  = C.NbSamples();

  if (typC == GeomAbs_BSplineCurve)
  {
    double t = C.LastParameter() - C.FirstParameter();
    if (t > Precision::PConfusion())
    {
      double t1 = U1 - U0;
      if (t1 < 0.0)
      {
        t1 = -t1;
      }
      nbs = C.NbKnots();
      nbs *= C.Degree();
      double anb = t1 / t * nbs;
      nbs        = (int)anb;

      int aMinPntNb = std::max(C.Degree() + 1, 4);
      if (nbs < aMinPntNb)
      {
        nbs = aMinPntNb;
      }
    }
  }
  else if (typC == GeomAbs_Circle)
  {
    // Try to reach deflection = eps*R, eps = 0.01
    const double minR = 1.; // eps = 0.01
    double       R    = C.Circle().Radius();
    if (R > minR)
    {
      double angl = 0.283079; // 2.*std::acos(1. - eps);
      int    n    = RealToInt(std::abs(U1 - U0) / angl);
      nbs         = std::max(n, nbs);
    }
  }

  if (nbs > 300)
  {
    nbs = 300;
  }
  return nbs;
}

//============================================================
int Geom2dInt_Geom2dCurveTool::NbSamples(const Adaptor2d_Curve2d& C)
{
  int               nbs  = C.NbSamples();
  GeomAbs_CurveType typC = C.GetType();
  if (typC == GeomAbs_Circle)
  {
    // Try to reach deflection = eps*R, eps = 0.01
    const double minR = 1.; // eps = 0.01
    double       R    = C.Circle().Radius();
    if (R > minR)
    {
      double angl = 0.283079; // 2.*std::acos(1. - eps);
      int    n    = RealToInt((C.LastParameter() - C.FirstParameter()) / angl);
      nbs         = std::max(n, nbs);
    }
  }

  return nbs;
}