2c4ea9536c2695af6bd1844995c79979f73f4e25
[occt.git] / 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 }