src/Geom2d/Geom2d_OffsetCurve.cxx

   1 // Created on: 1991-06-25
   2 // Created by: JCV
   3 // Copyright (c) 1991-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 //  modified by Edward AGAPOV (eap) Jan 28 2002 --- DN(), occ143(BUC60654)
  18
  19 #include <Geom2d_BezierCurve.hxx>
  20 #include <Geom2d_BSplineCurve.hxx>
  21 #include <Geom2d_Circle.hxx>
  22 #include <Geom2d_Curve.hxx>
  23 #include <Geom2d_Ellipse.hxx>
  24 #include <Geom2d_Geometry.hxx>
  25 #include <Geom2d_Hyperbola.hxx>
  26 #include <Geom2d_Line.hxx>
  27 #include <Geom2d_OffsetCurve.hxx>
  28 #include <Geom2d_Parabola.hxx>
  29 #include <Geom2d_TrimmedCurve.hxx>
  30 #include <Geom2d_UndefinedDerivative.hxx>
  31 #include <Geom2d_UndefinedValue.hxx>
  32 #include <gp.hxx>
  33 #include <gp_Pnt2d.hxx>
  34 #include <gp_Trsf2d.hxx>
  35 #include <gp_Vec2d.hxx>
  36 #include <gp_XY.hxx>
  37 #include <Precision.hxx>
  38 #include <Standard_ConstructionError.hxx>
  39 #include <Standard_NoSuchObject.hxx>
  40 #include <Standard_NotImplemented.hxx>
  41 #include <Standard_RangeError.hxx>
  42 #include <Standard_Type.hxx>
  43
  44 IMPLEMENT_STANDARD_RTTIEXT(Geom2d_OffsetCurve,Geom2d_Curve)
  45
  46 static const Standard_Real MyAngularToleranceForG1 = Precision::Angular();
  47
  48
  49 //=======================================================================
  50 //function : Copy
  51 //purpose  :
  52 //=======================================================================
  53
  54 Handle(Geom2d_Geometry) Geom2d_OffsetCurve::Copy () const
  55 {
  56   Handle(Geom2d_OffsetCurve) C;
  57   C = new Geom2d_OffsetCurve (basisCurve, offsetValue);
  58   return C;
  59 }
  60
  61
  62 //=======================================================================
  63 //function : Geom2d_OffsetCurve
  64 //purpose  : Basis curve cannot be an Offset curve or trimmed from
  65 //            offset curve.
  66 //=======================================================================
  67
  68 Geom2d_OffsetCurve::Geom2d_OffsetCurve (const Handle(Geom2d_Curve)& theCurve,
  69                                         const Standard_Real theOffset,
  70                                         const Standard_Boolean isTheNotCheckC0)
  71 : offsetValue (theOffset)
  72 {
  73   SetBasisCurve (theCurve, isTheNotCheckC0);
  74 }
  75
  76 //=======================================================================
  77 //function : Reverse
  78 //purpose  :
  79 //=======================================================================
  80
  81 void Geom2d_OffsetCurve::Reverse ()
  82 {
  83   basisCurve->Reverse();
  84   offsetValue = -offsetValue;
  85 }
  86
  87 //=======================================================================
  88 //function : ReversedParameter
  89 //purpose  :
  90 //=======================================================================
  91
  92 Standard_Real Geom2d_OffsetCurve::ReversedParameter( const Standard_Real U) const
  93 {
  94   return basisCurve->ReversedParameter( U);
  95 }
  96
  97 //=======================================================================
  98 //function : SetBasisCurve
  99 //purpose  :
 100 //=======================================================================
 101
 102 void Geom2d_OffsetCurve::SetBasisCurve (const Handle(Geom2d_Curve)& C,
 103                                         const Standard_Boolean isNotCheckC0)
 104 {
 105   const Standard_Real aUf = C->FirstParameter(),
 106                       aUl = C->LastParameter();
 107   Handle(Geom2d_Curve) aCheckingCurve = C;
 108   Standard_Boolean isTrimmed = Standard_False;
 109
 110   while(aCheckingCurve->IsKind(STANDARD_TYPE(Geom2d_TrimmedCurve)) ||
 111         aCheckingCurve->IsKind(STANDARD_TYPE(Geom2d_OffsetCurve)))
 112   {
 113     if (aCheckingCurve->IsKind(STANDARD_TYPE(Geom2d_TrimmedCurve)))
 114     {
 115       Handle(Geom2d_TrimmedCurve) aTrimC =
 116                 Handle(Geom2d_TrimmedCurve)::DownCast(aCheckingCurve);
 117       aCheckingCurve = aTrimC->BasisCurve();
 118       isTrimmed = Standard_True;
 119     }
 120
 121     if (aCheckingCurve->IsKind(STANDARD_TYPE(Geom2d_OffsetCurve)))
 122     {
 123       Handle(Geom2d_OffsetCurve) aOC =
 124                 Handle(Geom2d_OffsetCurve)::DownCast(aCheckingCurve);
 125       aCheckingCurve = aOC->BasisCurve();
 126       offsetValue += aOC->Offset();
 127     }
 128   }
 129
 130   myBasisCurveContinuity = aCheckingCurve->Continuity();
 131
 132   Standard_Boolean isC0 = !isNotCheckC0 &&
 133                           (myBasisCurveContinuity == GeomAbs_C0);
 134
 135   // Basis curve must be at least C1
 136   if (isC0 && aCheckingCurve->IsKind(STANDARD_TYPE(Geom2d_BSplineCurve)))
 137   {
 138     Handle(Geom2d_BSplineCurve) aBC = Handle(Geom2d_BSplineCurve)::DownCast(aCheckingCurve);
 139     if(aBC->IsG1(aUf, aUl, MyAngularToleranceForG1))
 140     {
 141       //Checking if basis curve has more smooth (C1, G2 and above) is not done.
 142       //It can be done in case of need.
 143       myBasisCurveContinuity = GeomAbs_G1;
 144       isC0 = Standard_False;
 145     }
 146
 147     // Raise exception if still C0
 148     if (isC0)
 149       throw Standard_ConstructionError("Offset on C0 curve");
 150   }
 151   //
 152   if(isTrimmed)
 153   {
 154     basisCurve = new Geom2d_TrimmedCurve(aCheckingCurve, aUf, aUl);
 155   }
 156   else
 157   {
 158     basisCurve = aCheckingCurve;
 159   }
 160
 161   myEvaluator = new Geom2dEvaluator_OffsetCurve(basisCurve, offsetValue);
 162 }
 163
 164 //=======================================================================
 165 //function : SetOffsetValue
 166 //purpose  :
 167 //=======================================================================
 168
 169 void Geom2d_OffsetCurve::SetOffsetValue (const Standard_Real D)
 170 {
 171   offsetValue = D;
 172   myEvaluator->SetOffsetValue(offsetValue);
 173 }
 174
 175 //=======================================================================
 176 //function : BasisCurve
 177 //purpose  :
 178 //=======================================================================
 179
 180 Handle(Geom2d_Curve) Geom2d_OffsetCurve::BasisCurve () const
 181 {
 182   return basisCurve;
 183 }
 184
 185 //=======================================================================
 186 //function : Continuity
 187 //purpose  :
 188 //=======================================================================
 189
 190 GeomAbs_Shape Geom2d_OffsetCurve::Continuity () const
 191 {
 192   GeomAbs_Shape OffsetShape=GeomAbs_C0;
 193   switch (myBasisCurveContinuity) {
 194      case GeomAbs_C0 : OffsetShape = GeomAbs_C0;   break;
 195      case GeomAbs_C1 : OffsetShape = GeomAbs_C0;   break;
 196      case GeomAbs_C2 : OffsetShape = GeomAbs_C1;   break;
 197      case GeomAbs_C3 : OffsetShape = GeomAbs_C2;   break;
 198      case GeomAbs_CN : OffsetShape = GeomAbs_CN;   break;
 199      case GeomAbs_G1 : OffsetShape = GeomAbs_G1;   break;
 200      case GeomAbs_G2 : OffsetShape = GeomAbs_G2;   break;
 201   }
 202
 203   return OffsetShape;
 204 }
 205
 206 //=======================================================================
 207 //function : D0
 208 //purpose  :
 209 //=======================================================================
 210
 211 void Geom2d_OffsetCurve::D0 (const Standard_Real theU,
 212                                    gp_Pnt2d&     theP) const
 213 {
 214   myEvaluator->D0(theU, theP);
 215 }
 216
 217 //=======================================================================
 218 //function : D1
 219 //purpose  :
 220 //=======================================================================
 221 void Geom2d_OffsetCurve::D1 (const Standard_Real theU, gp_Pnt2d& theP, gp_Vec2d& theV1) const
 222 {
 223   myEvaluator->D1(theU, theP, theV1);
 224 }
 225
 226 //=======================================================================
 227 //function : D2
 228 //purpose  :
 229 //=======================================================================
 230
 231 void Geom2d_OffsetCurve::D2 (const Standard_Real theU,
 232                                    gp_Pnt2d& theP,
 233                                    gp_Vec2d& theV1, gp_Vec2d& theV2) const
 234 {
 235   myEvaluator->D2(theU, theP, theV1, theV2);
 236 }
 237
 238
 239 //=======================================================================
 240 //function : D3
 241 //purpose  :
 242 //=======================================================================
 243
 244 void Geom2d_OffsetCurve::D3 (const Standard_Real theU,
 245                                    gp_Pnt2d& theP,
 246                                    gp_Vec2d& theV1, gp_Vec2d& theV2, gp_Vec2d& theV3) const
 247 {
 248   myEvaluator->D3(theU, theP, theV1, theV2, theV3);
 249 }
 250
 251 //=======================================================================
 252 //function : DN
 253 //purpose  :
 254 //=======================================================================
 255
 256 gp_Vec2d Geom2d_OffsetCurve::DN (const Standard_Real U,
 257                                  const Standard_Integer N) const
 258 {
 259   Standard_RangeError_Raise_if (N < 1, "Exception: Geom2d_OffsetCurve::DN(). N<1.");
 260
 261   gp_Vec2d VN, VBidon;
 262   gp_Pnt2d PBidon;
 263   switch (N) {
 264   case 1: D1( U, PBidon, VN); break;
 265   case 2: D2( U, PBidon, VBidon, VN); break;
 266   case 3: D3( U, PBidon, VBidon, VBidon, VN); break;
 267   default:
 268     throw Standard_NotImplemented("Exception: Derivative order is greater than 3. "
 269       "Cannot compute of derivative.");
 270   }
 271
 272   return VN;
 273 }
 274
 275 //=======================================================================
 276 //function : FirstParameter
 277 //purpose  :
 278 //=======================================================================
 279
 280 Standard_Real Geom2d_OffsetCurve::FirstParameter () const
 281 {
 282   return basisCurve->FirstParameter();
 283 }
 284
 285 //=======================================================================
 286 //function : LastParameter
 287 //purpose  :
 288 //=======================================================================
 289
 290 Standard_Real Geom2d_OffsetCurve::LastParameter () const
 291 {
 292   return basisCurve->LastParameter();
 293 }
 294
 295
 296 //=======================================================================
 297 //function : Offset
 298 //purpose  :
 299 //=======================================================================
 300
 301 Standard_Real Geom2d_OffsetCurve::Offset () const
 302 { return offsetValue; }
 303
 304 //=======================================================================
 305 //function : IsClosed
 306 //purpose  :
 307 //=======================================================================
 308
 309 Standard_Boolean Geom2d_OffsetCurve::IsClosed () const
 310 {
 311   gp_Pnt2d PF, PL;
 312   D0(FirstParameter(),PF);
 313   D0(LastParameter(),PL);
 314   return ( PF.Distance(PL) <= gp::Resolution());
 315 }
 316
 317 //=======================================================================
 318 //function : IsCN
 319 //purpose  :
 320 //=======================================================================
 321
 322 Standard_Boolean Geom2d_OffsetCurve::IsCN (const Standard_Integer N) const
 323 {
 324   Standard_RangeError_Raise_if (N < 0, " " );
 325   return basisCurve->IsCN (N + 1);
 326 }
 327
 328 //=======================================================================
 329 //function : IsPeriodic
 330 //purpose  :
 331 //=======================================================================
 332
 333 Standard_Boolean Geom2d_OffsetCurve::IsPeriodic () const
 334 {
 335   return basisCurve->IsPeriodic();
 336 }
 337
 338 //=======================================================================
 339 //function : Period
 340 //purpose  :
 341 //=======================================================================
 342
 343 Standard_Real Geom2d_OffsetCurve::Period() const
 344 {
 345   return basisCurve->Period();
 346 }
 347
 348 //=======================================================================
 349 //function : Transform
 350 //purpose  :
 351 //=======================================================================
 352
 353 void Geom2d_OffsetCurve::Transform (const gp_Trsf2d& T)
 354 {
 355   basisCurve->Transform (T);
 356   offsetValue *= Abs(T.ScaleFactor());
 357
 358   myEvaluator->SetOffsetValue(offsetValue);
 359 }
 360
 361
 362 //=======================================================================
 363 //function : TransformedParameter
 364 //purpose  :
 365 //=======================================================================
 366
 367 Standard_Real Geom2d_OffsetCurve::TransformedParameter(const Standard_Real U,
 368                                                         const gp_Trsf2d& T) const
 369 {
 370   return basisCurve->TransformedParameter(U,T);
 371 }
 372
 373 //=======================================================================
 374 //function : ParametricTransformation
 375 //purpose  :
 376 //=======================================================================
 377
 378 Standard_Real Geom2d_OffsetCurve::ParametricTransformation(const gp_Trsf2d& T) const
 379 {
 380   return basisCurve->ParametricTransformation(T);
 381 }
 382
 383 //=======================================================================
 384 //function : GetBasisCurveContinuity
 385 //purpose  :
 386 //=======================================================================
 387 GeomAbs_Shape Geom2d_OffsetCurve::GetBasisCurveContinuity() const
 388 {
 389   return myBasisCurveContinuity;
 390 }