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