src/BRepBlend/BRepBlend_HCurveTool.cxx

   1 // Created on: 1995-07-17
   2 // Created by: Modelistation
   3 // Copyright (c) 1995-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17
  18 #include <Adaptor3d_HCurve.hxx>
  19 #include <BRepBlend_HCurveTool.hxx>
  20 #include <Geom_BezierCurve.hxx>
  21 #include <Geom_BSplineCurve.hxx>
  22 #include <GeomAbs_CurveType.hxx>
  23 #include <GeomAbs_Shape.hxx>
  24 #include <gp_Pnt.hxx>
  25 #include <gp_Vec.hxx>
  26 #include <Standard_DomainError.hxx>
  27 #include <Standard_NoSuchObject.hxx>
  28 #include <Standard_OutOfRange.hxx>
  29 #include <TColStd_Array1OfReal.hxx>
  30
  31 //============================================================
  32 Standard_Integer BRepBlend_HCurveTool::NbSamples (const Handle(Adaptor3d_HCurve)& C,
  33                                                        const Standard_Real U0,
  34                                                        const Standard_Real U1) {
  35   GeomAbs_CurveType typC = C->GetType();
  36   static Standard_Real nbsOther = 10.0;
  37   Standard_Real nbs = nbsOther;
  38
  39   if(typC == GeomAbs_Line)
  40     nbs = 2;
  41   else if(typC == GeomAbs_BezierCurve)
  42     nbs = 3 + C->Bezier()->NbPoles();
  43   else if(typC == GeomAbs_BSplineCurve) {
  44     Handle(Geom_BSplineCurve) BSC = C->BSpline();
  45     nbs = BSC->NbKnots();
  46     nbs*= BSC->Degree();
  47     nbs*= BSC->LastParameter()- BSC->FirstParameter();
  48     nbs/= U1-U0;
  49     if(nbs < 2.0) nbs=2;
  50   }
  51   if(nbs>50)
  52     nbs = 50;
  53   return((Standard_Integer)nbs);
  54 }