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 // 24-Aug-95 : xab removed C1 and C2 test : appeller D1 et D2
18 // avec discernement !
19 // 19-09-97 : JPI correction derivee seconde
20
21 #include <Geom_BezierCurve.hxx>
22 #include <Geom_BSplineCurve.hxx>
23 #include <Geom_Circle.hxx>
24 #include <Geom_Curve.hxx>
25 #include <Geom_Ellipse.hxx>
26 #include <Geom_Geometry.hxx>
27 #include <Geom_Hyperbola.hxx>
28 #include <Geom_Line.hxx>
29 #include <Geom_OffsetCurve.hxx>
30 #include <Geom_Parabola.hxx>
31 #include <Geom_TrimmedCurve.hxx>
32 #include <Geom_UndefinedDerivative.hxx>
33 #include <Geom_UndefinedValue.hxx>
34 #include <gp.hxx>
35 #include <gp_Dir.hxx>
36 #include <gp_Pnt.hxx>
37 #include <gp_Trsf.hxx>
38 #include <gp_Vec.hxx>
39 #include <gp_XYZ.hxx>
40 #include <Standard_ConstructionError.hxx>
41 #include <Standard_NoSuchObject.hxx>
42 #include <Standard_NotImplemented.hxx>
43 #include <Standard_RangeError.hxx>
44 #include <Standard_Type.hxx>
45
46 IMPLEMENT_STANDARD_RTTIEXT(Geom_OffsetCurve,Geom_Curve)
47
48 static const Standard_Real MyAngularToleranceForG1 = Precision::Angular();
49
50
51 //=======================================================================
52 //function : Copy
53 //purpose :
54 //=======================================================================
55
Handle(Geom_Geometry)56 Handle(Geom_Geometry) Geom_OffsetCurve::Copy () const {
57
58 Handle(Geom_OffsetCurve) C;
59 C = new Geom_OffsetCurve (basisCurve, offsetValue, direction);
60 return C;
61 }
62
63
64
65 //=======================================================================
66 //function : Geom_OffsetCurve
67 //purpose : Basis curve cannot be an Offset curve or trimmed from
68 // offset curve.
69 //=======================================================================
70
Geom_OffsetCurve(const Handle (Geom_Curve)& theCurve,const Standard_Real theOffset,const gp_Dir & theDir,const Standard_Boolean isTheNotCheckC0)71 Geom_OffsetCurve::Geom_OffsetCurve (const Handle(Geom_Curve)& theCurve,
72 const Standard_Real theOffset,
73 const gp_Dir& theDir,
74 const Standard_Boolean isTheNotCheckC0)
75 : direction(theDir), offsetValue(theOffset)
76 {
77 SetBasisCurve (theCurve, isTheNotCheckC0);
78 }
79
80
81 //=======================================================================
82 //function : Reverse
83 //purpose :
84 //=======================================================================
85
Reverse()86 void Geom_OffsetCurve::Reverse ()
87 {
88 basisCurve->Reverse();
89 offsetValue = -offsetValue;
90 myEvaluator->SetOffsetValue(offsetValue);
91 }
92
93
94 //=======================================================================
95 //function : ReversedParameter
96 //purpose :
97 //=======================================================================
98
ReversedParameter(const Standard_Real U) const99 Standard_Real Geom_OffsetCurve::ReversedParameter( const Standard_Real U) const
100 {
101 return basisCurve->ReversedParameter( U);
102 }
103
104 //=======================================================================
105 //function : Direction
106 //purpose :
107 //=======================================================================
108
Direction() const109 const gp_Dir& Geom_OffsetCurve::Direction () const
110 { return direction; }
111
112 //=======================================================================
113 //function : SetDirection
114 //purpose :
115 //=======================================================================
116
SetDirection(const gp_Dir & V)117 void Geom_OffsetCurve::SetDirection (const gp_Dir& V)
118 {
119 direction = V;
120 myEvaluator->SetOffsetDirection(direction);
121 }
122
123 //=======================================================================
124 //function : SetOffsetValue
125 //purpose :
126 //=======================================================================
127
SetOffsetValue(const Standard_Real D)128 void Geom_OffsetCurve::SetOffsetValue (const Standard_Real D)
129 {
130 offsetValue = D;
131 myEvaluator->SetOffsetValue(offsetValue);
132 }
133
134
135 //=======================================================================
136 //function : IsPeriodic
137 //purpose :
138 //=======================================================================
139
IsPeriodic() const140 Standard_Boolean Geom_OffsetCurve::IsPeriodic () const
141 {
142 return basisCurve->IsPeriodic();
143 }
144
145 //=======================================================================
146 //function : Period
147 //purpose :
148 //=======================================================================
149
Period() const150 Standard_Real Geom_OffsetCurve::Period () const
151 {
152 return basisCurve->Period();
153 }
154
155 //=======================================================================
156 //function : SetBasisCurve
157 //purpose :
158 //=======================================================================
159
SetBasisCurve(const Handle (Geom_Curve)& C,const Standard_Boolean isNotCheckC0)160 void Geom_OffsetCurve::SetBasisCurve (const Handle(Geom_Curve)& C,
161 const Standard_Boolean isNotCheckC0)
162 {
163 const Standard_Real aUf = C->FirstParameter(),
164 aUl = C->LastParameter();
165 Handle(Geom_Curve) aCheckingCurve = Handle(Geom_Curve)::DownCast(C->Copy());
166 Standard_Boolean isTrimmed = Standard_False;
167
168 while(aCheckingCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve)) ||
169 aCheckingCurve->IsKind(STANDARD_TYPE(Geom_OffsetCurve)))
170 {
171 if (aCheckingCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve)))
172 {
173 Handle(Geom_TrimmedCurve) aTrimC =
174 Handle(Geom_TrimmedCurve)::DownCast(aCheckingCurve);
175 aCheckingCurve = aTrimC->BasisCurve();
176 isTrimmed = Standard_True;
177 }
178
179 if (aCheckingCurve->IsKind(STANDARD_TYPE(Geom_OffsetCurve)))
180 {
181 Handle(Geom_OffsetCurve) aOC =
182 Handle(Geom_OffsetCurve)::DownCast(aCheckingCurve);
183 aCheckingCurve = aOC->BasisCurve();
184 Standard_Real PrevOff = aOC->Offset();
185 gp_Vec V1(aOC->Direction());
186 gp_Vec V2(direction);
187 gp_Vec Vdir(PrevOff*V1 + offsetValue*V2);
188
189 if (offsetValue >= 0.)
190 {
191 offsetValue = Vdir.Magnitude();
192 direction.SetXYZ(Vdir.XYZ());
193 }
194 else
195 {
196 offsetValue = -Vdir.Magnitude();
197 direction.SetXYZ((-Vdir).XYZ());
198 }
199 }
200 }
201
202 myBasisCurveContinuity = aCheckingCurve->Continuity();
203
204 Standard_Boolean isC0 = !isNotCheckC0 &&
205 (myBasisCurveContinuity == GeomAbs_C0);
206
207 // Basis curve must be at least C1
208 if (isC0 && aCheckingCurve->IsKind(STANDARD_TYPE(Geom_BSplineCurve)))
209 {
210 Handle(Geom_BSplineCurve) aBC = Handle(Geom_BSplineCurve)::DownCast(aCheckingCurve);
211 if(aBC->IsG1(aUf, aUl, MyAngularToleranceForG1))
212 {
213 //Checking if basis curve has more smooth (C1, G2 and above) is not done.
214 //It can be done in case of need.
215 myBasisCurveContinuity = GeomAbs_G1;
216 isC0 = Standard_False;
217 }
218
219 // Raise exception if still C0
220 if (isC0)
221 throw Standard_ConstructionError("Offset on C0 curve");
222 }
223 //
224 if(isTrimmed)
225 {
226 basisCurve = new Geom_TrimmedCurve(aCheckingCurve, aUf, aUl);
227 }
228 else
229 {
230 basisCurve = aCheckingCurve;
231 }
232
233 myEvaluator = new GeomEvaluator_OffsetCurve(basisCurve, offsetValue, direction);
234 }
235
236
237
238 //=======================================================================
239 //function : BasisCurve
240 //purpose :
241 //=======================================================================
242
Handle(Geom_Curve)243 Handle(Geom_Curve) Geom_OffsetCurve::BasisCurve () const
244 {
245 return basisCurve;
246 }
247
248
249 //=======================================================================
250 //function : Continuity
251 //purpose :
252 //=======================================================================
253
Continuity() const254 GeomAbs_Shape Geom_OffsetCurve::Continuity () const {
255
256 GeomAbs_Shape OffsetShape=GeomAbs_C0;
257 switch (myBasisCurveContinuity) {
258 case GeomAbs_C0 : OffsetShape = GeomAbs_C0; break;
259 case GeomAbs_C1 : OffsetShape = GeomAbs_C0; break;
260 case GeomAbs_C2 : OffsetShape = GeomAbs_C1; break;
261 case GeomAbs_C3 : OffsetShape = GeomAbs_C2; break;
262 case GeomAbs_CN : OffsetShape = GeomAbs_CN; break;
263 case GeomAbs_G1 : OffsetShape = GeomAbs_G1; break;
264 case GeomAbs_G2 : OffsetShape = GeomAbs_G2; break;
265 }
266 return OffsetShape;
267 }
268
269
270 //=======================================================================
271 //function : D0
272 //purpose :
273 //=======================================================================
274
D0(const Standard_Real U,gp_Pnt & P) const275 void Geom_OffsetCurve::D0 (const Standard_Real U, gp_Pnt& P) const
276 {
277 myEvaluator->D0(U, P);
278 }
279
280 //=======================================================================
281 //function : D1
282 //purpose :
283 //=======================================================================
284
D1(const Standard_Real U,gp_Pnt & P,gp_Vec & V1) const285 void Geom_OffsetCurve::D1 (const Standard_Real U, gp_Pnt& P, gp_Vec& V1) const
286 {
287 myEvaluator->D1(U, P, V1);
288 }
289
290 //=======================================================================
291 //function : D2
292 //purpose :
293 //=======================================================================
294
D2(const Standard_Real U,gp_Pnt & P,gp_Vec & V1,gp_Vec & V2) const295 void Geom_OffsetCurve::D2 (const Standard_Real U, gp_Pnt& P, gp_Vec& V1, gp_Vec& V2) const
296 {
297 myEvaluator->D2(U, P, V1, V2);
298 }
299
300 //=======================================================================
301 //function : D3
302 //purpose :
303 //=======================================================================
304
D3(const Standard_Real theU,gp_Pnt & theP,gp_Vec & theV1,gp_Vec & theV2,gp_Vec & theV3) const305 void Geom_OffsetCurve::D3 (const Standard_Real theU, gp_Pnt& theP,
306 gp_Vec& theV1, gp_Vec& theV2, gp_Vec& theV3) const
307 {
308 myEvaluator->D3(theU, theP, theV1, theV2, theV3);
309 }
310
311
312 //=======================================================================
313 //function : DN
314 //purpose :
315 //=======================================================================
316
DN(const Standard_Real U,const Standard_Integer N) const317 gp_Vec Geom_OffsetCurve::DN (const Standard_Real U, const Standard_Integer N) const
318 {
319 Standard_RangeError_Raise_if (N < 1, "Exception: "
320 "Geom_OffsetCurve::DN(...). N<1.");
321
322 gp_Vec VN, Vtemp;
323 gp_Pnt Ptemp;
324 switch (N)
325 {
326 case 1:
327 D1( U, Ptemp, VN);
328 break;
329 case 2:
330 D2( U, Ptemp, Vtemp, VN);
331 break;
332 case 3:
333 D3( U, Ptemp, Vtemp, Vtemp, VN);
334 break;
335 default:
336 throw Standard_NotImplemented("Exception: "
337 "Derivative order is greater than 3. Cannot compute of derivative.");
338 }
339
340 return VN;
341 }
342
343
344 //=======================================================================
345 //function : FirstParameter
346 //purpose :
347 //=======================================================================
348
FirstParameter() const349 Standard_Real Geom_OffsetCurve::FirstParameter () const
350 {
351 return basisCurve->FirstParameter();
352 }
353
354
355 //=======================================================================
356 //function : LastParameter
357 //purpose :
358 //=======================================================================
359
LastParameter() const360 Standard_Real Geom_OffsetCurve::LastParameter () const
361 {
362 return basisCurve->LastParameter();
363 }
364
365
366 //=======================================================================
367 //function : Offset
368 //purpose :
369 //=======================================================================
370
Offset() const371 Standard_Real Geom_OffsetCurve::Offset () const
372 { return offsetValue; }
373
374
375 //=======================================================================
376 //function : IsClosed
377 //purpose :
378 //=======================================================================
379
IsClosed() const380 Standard_Boolean Geom_OffsetCurve::IsClosed () const
381 {
382 gp_Pnt PF,PL;
383 D0(FirstParameter(),PF);
384 D0(LastParameter(),PL);
385 return ( PF.Distance(PL) <= gp::Resolution());
386 }
387
388
389
390 //=======================================================================
391 //function : IsCN
392 //purpose :
393 //=======================================================================
394
IsCN(const Standard_Integer N) const395 Standard_Boolean Geom_OffsetCurve::IsCN (const Standard_Integer N) const {
396
397 Standard_RangeError_Raise_if (N < 0, " ");
398 return basisCurve->IsCN (N + 1);
399 }
400
401
402 //=======================================================================
403 //function : Transform
404 //purpose :
405 //=======================================================================
406
Transform(const gp_Trsf & T)407 void Geom_OffsetCurve::Transform (const gp_Trsf& T)
408 {
409 basisCurve->Transform (T);
410 direction.Transform(T);
411 offsetValue *= T.ScaleFactor();
412
413 myEvaluator->SetOffsetValue(offsetValue);
414 myEvaluator->SetOffsetDirection(direction);
415 }
416
417 //=======================================================================
418 //function : TransformedParameter
419 //purpose :
420 //=======================================================================
421
TransformedParameter(const Standard_Real U,const gp_Trsf & T) const422 Standard_Real Geom_OffsetCurve::TransformedParameter(const Standard_Real U,
423 const gp_Trsf& T) const
424 {
425 return basisCurve->TransformedParameter(U,T);
426 }
427
428 //=======================================================================
429 //function : ParametricTransformation
430 //purpose :
431 //=======================================================================
432
ParametricTransformation(const gp_Trsf & T) const433 Standard_Real Geom_OffsetCurve::ParametricTransformation(const gp_Trsf& T)
434 const
435 {
436 return basisCurve->ParametricTransformation(T);
437 }
438
439 //=======================================================================
440 //function : GetBasisCurveContinuity
441 //purpose :
442 //=======================================================================
GetBasisCurveContinuity() const443 GeomAbs_Shape Geom_OffsetCurve::GetBasisCurveContinuity() const
444 {
445 return myBasisCurveContinuity;
446 }
447
448 //=======================================================================
449 //function : DumpJson
450 //purpose :
451 //=======================================================================
DumpJson(Standard_OStream & theOStream,Standard_Integer theDepth) const452 void Geom_OffsetCurve::DumpJson (Standard_OStream& theOStream, Standard_Integer theDepth) const
453 {
454 OCCT_DUMP_TRANSIENT_CLASS_BEGIN (theOStream)
455
456 OCCT_DUMP_BASE_CLASS (theOStream, theDepth, Geom_Curve)
457
458 OCCT_DUMP_FIELD_VALUES_DUMPED (theOStream, theDepth, basisCurve.get())
459 OCCT_DUMP_FIELD_VALUES_DUMPED (theOStream, theDepth, &direction)
460
461 OCCT_DUMP_FIELD_VALUE_NUMERICAL (theOStream, offsetValue)
462 OCCT_DUMP_FIELD_VALUE_NUMERICAL (theOStream, myBasisCurveContinuity)
463 }
464