b311480e |
1 | // Created on: 1993-07-19 |
2 | // Created by: Remi LEQUETTE |
3 | // Copyright (c) 1993-1999 Matra Datavision |
973c2be1 |
4 | // Copyright (c) 1999-2014 OPEN CASCADE SAS |
b311480e |
5 | // |
973c2be1 |
6 | // This file is part of Open CASCADE Technology software library. |
b311480e |
7 | // |
d5f74e42 |
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 |
973c2be1 |
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. |
b311480e |
13 | // |
973c2be1 |
14 | // Alternatively, this file may be used under the terms of Open CASCADE |
15 | // commercial license or contractual agreement. |
7fd59977 |
16 | |
7fd59977 |
17 | |
42cf5bc1 |
18 | #include <Geom_BezierCurve.hxx> |
19 | #include <Geom_BSplineCurve.hxx> |
7fd59977 |
20 | #include <Geom_Circle.hxx> |
42cf5bc1 |
21 | #include <Geom_Curve.hxx> |
7fd59977 |
22 | #include <Geom_Ellipse.hxx> |
7fd59977 |
23 | #include <Geom_Hyperbola.hxx> |
42cf5bc1 |
24 | #include <Geom_Line.hxx> |
7fd59977 |
25 | #include <Geom_OffsetCurve.hxx> |
42cf5bc1 |
26 | #include <Geom_Parabola.hxx> |
27 | #include <Geom_TrimmedCurve.hxx> |
28 | #include <GeomTools.hxx> |
29 | #include <GeomTools_CurveSet.hxx> |
30 | #include <GeomTools_UndefinedTypeHandler.hxx> |
7fd59977 |
31 | #include <gp_Circ.hxx> |
32 | #include <gp_Elips.hxx> |
7fd59977 |
33 | #include <gp_Hypr.hxx> |
42cf5bc1 |
34 | #include <gp_Lin.hxx> |
35 | #include <gp_Parab.hxx> |
36 | #include <Message_ProgressIndicator.hxx> |
7fd59977 |
37 | #include <Message_ProgressSentry.hxx> |
42cf5bc1 |
38 | #include <Standard_ErrorHandler.hxx> |
39 | #include <Standard_Failure.hxx> |
40 | #include <Standard_OutOfRange.hxx> |
41 | #include <Standard_Stream.hxx> |
42 | #include <TColgp_Array1OfPnt.hxx> |
43 | #include <TColStd_Array1OfInteger.hxx> |
44 | #include <TColStd_Array1OfReal.hxx> |
7fd59977 |
45 | |
46 | #define LINE 1 |
47 | #define CIRCLE 2 |
48 | #define ELLIPSE 3 |
49 | #define PARABOLA 4 |
50 | #define HYPERBOLA 5 |
51 | #define BEZIER 6 |
52 | #define BSPLINE 7 |
53 | #define TRIMMED 8 |
54 | #define OFFSET 9 |
55 | |
56 | //======================================================================= |
57 | //function : GeomTools_CurveSet |
58 | //purpose : |
59 | //======================================================================= |
60 | |
61 | GeomTools_CurveSet::GeomTools_CurveSet() |
62 | { |
63 | } |
64 | |
65 | |
66 | //======================================================================= |
67 | //function : Clear |
68 | //purpose : |
69 | //======================================================================= |
70 | |
71 | void GeomTools_CurveSet::Clear() |
72 | { |
73 | myMap.Clear(); |
74 | } |
75 | |
76 | |
77 | //======================================================================= |
78 | //function : Add |
79 | //purpose : |
80 | //======================================================================= |
81 | |
82 | Standard_Integer GeomTools_CurveSet::Add(const Handle(Geom_Curve)& C) |
83 | { |
84 | return (C.IsNull()) ? 0 : myMap.Add(C); |
85 | } |
86 | |
87 | |
88 | //======================================================================= |
89 | //function : Curve |
90 | //purpose : |
91 | //======================================================================= |
92 | |
93 | Handle(Geom_Curve) GeomTools_CurveSet::Curve |
94 | (const Standard_Integer I)const |
95 | { |
96 | if (I <= 0 || I > myMap.Extent()) |
97 | return Handle(Geom_Curve)(); |
98 | return Handle(Geom_Curve)::DownCast(myMap(I)); |
99 | } |
100 | |
101 | |
102 | //======================================================================= |
103 | //function : Index |
104 | //purpose : |
105 | //======================================================================= |
106 | |
107 | Standard_Integer GeomTools_CurveSet::Index |
108 | (const Handle(Geom_Curve)& S)const |
109 | { |
110 | return S.IsNull() ? 0 : myMap.FindIndex(S); |
111 | } |
112 | |
113 | |
114 | //======================================================================= |
115 | //function : Print |
116 | //purpose : |
117 | //======================================================================= |
118 | |
119 | static void Print(const gp_Pnt P, |
120 | Standard_OStream& OS, |
121 | const Standard_Boolean compact) |
122 | { |
123 | OS << P.X(); |
124 | if (!compact) OS << ","; |
125 | OS << " "; |
126 | OS << P.Y(); |
127 | if (!compact) OS << ","; |
128 | OS << " "; |
129 | OS << P.Z(); |
130 | OS << " "; |
131 | } |
132 | |
133 | //======================================================================= |
134 | //function : Print |
135 | //purpose : |
136 | //======================================================================= |
137 | |
138 | static void Print(const gp_Dir D, |
139 | Standard_OStream& OS, |
140 | const Standard_Boolean compact) |
141 | { |
142 | OS << D.X(); |
143 | if (!compact) OS << ","; |
144 | OS << " "; |
145 | OS << D.Y(); |
146 | if (!compact) OS << ","; |
147 | OS << " "; |
148 | OS << D.Z(); |
149 | OS << " "; |
150 | } |
151 | |
152 | |
153 | //======================================================================= |
154 | //function : Print |
155 | //purpose : |
156 | //======================================================================= |
157 | |
158 | static void Print(const Handle(Geom_Line)& L, |
159 | Standard_OStream& OS, |
160 | const Standard_Boolean compact) |
161 | { |
162 | if (compact) |
163 | OS << LINE << " "; |
164 | else |
165 | OS << "Line"; |
166 | |
167 | gp_Lin C = L->Lin(); |
168 | if (!compact) OS << "\n Origin :"; |
169 | Print(C.Location(),OS,compact); |
170 | if (!compact) OS << "\n Axis :"; |
171 | Print(C.Direction(),OS,compact); |
172 | if (!compact) OS << "\n"; |
173 | OS << "\n"; |
174 | } |
175 | |
176 | //======================================================================= |
177 | //function : Print |
178 | //purpose : |
179 | //======================================================================= |
180 | |
181 | static void Print(const Handle(Geom_Circle)& CC, |
182 | Standard_OStream& OS, |
183 | const Standard_Boolean compact) |
184 | { |
185 | if (compact) |
186 | OS << CIRCLE << " "; |
187 | else |
188 | OS << "Circle"; |
189 | |
190 | gp_Circ C = CC->Circ(); |
191 | if (!compact) OS << "\n Center :"; |
192 | Print(C.Location(),OS,compact); |
193 | if (!compact) OS << "\n Axis :"; |
194 | Print(C.Axis().Direction(),OS,compact); |
195 | if (!compact) OS << "\n XAxis :"; |
196 | Print(C.XAxis().Direction(),OS,compact); |
197 | if (!compact) OS << "\n YAxis :"; |
198 | Print(C.YAxis().Direction(),OS,compact); |
199 | if (!compact) OS << "\n Radius :"; |
200 | OS << C.Radius(); |
201 | if (!compact) OS << "\n"; |
202 | OS << "\n"; |
203 | } |
204 | |
205 | //======================================================================= |
206 | //function : Print |
207 | //purpose : |
208 | //======================================================================= |
209 | |
210 | static void Print(const Handle(Geom_Ellipse)& E, |
211 | Standard_OStream& OS, |
212 | const Standard_Boolean compact) |
213 | { |
214 | if (compact) |
215 | OS << ELLIPSE << " "; |
216 | else |
217 | OS << "Ellipse"; |
218 | |
219 | gp_Elips C = E->Elips(); |
220 | if (!compact) OS << "\n Center :"; |
221 | Print(C.Location(),OS,compact); |
222 | if (!compact) OS << "\n Axis :"; |
223 | Print(C.Axis().Direction(),OS,compact); |
224 | if (!compact) OS << "\n XAxis :"; |
225 | Print(C.XAxis().Direction(),OS,compact); |
226 | if (!compact) OS << "\n YAxis :"; |
227 | Print(C.YAxis().Direction(),OS,compact); |
228 | if (!compact) OS << "\n Radii :"; |
229 | OS << C.MajorRadius(); |
230 | if (!compact) OS << ","; |
231 | OS << " "; |
232 | OS << C.MinorRadius(); |
233 | if (!compact) OS << "\n"; |
234 | OS << "\n"; |
235 | } |
236 | |
237 | //======================================================================= |
238 | //function : Print |
239 | //purpose : |
240 | //======================================================================= |
241 | |
242 | static void Print(const Handle(Geom_Parabola)& P, |
243 | Standard_OStream& OS, |
244 | const Standard_Boolean compact) |
245 | { |
246 | if (compact) |
247 | OS << PARABOLA << " "; |
248 | else |
249 | OS << "Parabola"; |
250 | |
251 | gp_Parab C = P->Parab(); |
252 | if (!compact) OS << "\n Center :"; |
253 | Print(C.Location(),OS,compact); |
254 | if (!compact) OS << "\n Axis :"; |
255 | Print(C.Axis().Direction(),OS,compact); |
256 | if (!compact) OS << "\n XAxis :"; |
257 | Print(C.XAxis().Direction(),OS,compact); |
258 | if (!compact) OS << "\n YAxis :"; |
259 | Print(C.YAxis().Direction(),OS,compact); |
260 | if (!compact) OS << "\n Focal :"; |
261 | OS << C.Focal(); |
262 | if (!compact) OS << "\n"; |
263 | OS << "\n"; |
264 | } |
265 | |
266 | //======================================================================= |
267 | //function : Print |
268 | //purpose : |
269 | //======================================================================= |
270 | |
271 | static void Print(const Handle(Geom_Hyperbola)& H, |
272 | Standard_OStream& OS, |
273 | const Standard_Boolean compact) |
274 | { |
275 | if (compact) |
276 | OS << HYPERBOLA << " "; |
277 | else |
278 | OS << "Hyperbola"; |
279 | |
280 | gp_Hypr C = H->Hypr(); |
281 | if (!compact) OS << "\n Center :"; |
282 | Print(C.Location(),OS,compact); |
283 | if (!compact) OS << "\n Axis :"; |
284 | Print(C.Axis().Direction(),OS,compact); |
285 | if (!compact) OS << "\n XAxis :"; |
286 | Print(C.XAxis().Direction(),OS,compact); |
287 | if (!compact) OS << "\n YAxis :"; |
288 | Print(C.YAxis().Direction(),OS,compact); |
289 | if (!compact) OS << "\n Radii :"; |
290 | OS << C.MajorRadius(); |
291 | if (!compact) OS << ","; |
292 | OS << " "; |
293 | OS << C.MinorRadius(); |
294 | if (!compact) OS << "\n"; |
295 | OS << "\n"; |
296 | } |
297 | |
298 | //======================================================================= |
299 | //function : Print |
300 | //purpose : |
301 | //======================================================================= |
302 | |
303 | static void Print(const Handle(Geom_BezierCurve)& B, |
304 | Standard_OStream& OS, |
305 | const Standard_Boolean compact) |
306 | { |
307 | if (compact) |
308 | OS << BEZIER << " "; |
309 | else |
310 | OS << "BezierCurve"; |
311 | |
312 | Standard_Boolean rational = B->IsRational(); |
313 | if (compact) |
314 | OS << (rational ? 1 : 0) << " "; |
315 | else { |
316 | if (rational) |
317 | OS << " rational"; |
318 | } |
319 | |
320 | // poles and weights |
321 | Standard_Integer i,degree = B->Degree(); |
322 | if (!compact) OS << "\n Degree :"; |
323 | OS << degree << " "; |
324 | |
325 | for (i = 1; i <= degree+1; i++) { |
04232180 |
326 | if (!compact) OS << "\n "<<std::setw(2)<<i<<" : "; |
7fd59977 |
327 | Print(B->Pole(i),OS,compact); |
328 | if (rational) |
329 | OS << " " << B->Weight(i); |
330 | if (compact) |
331 | OS << " "; |
332 | } |
333 | OS << "\n"; |
334 | if (!compact) OS << "\n"; |
335 | } |
336 | |
337 | //======================================================================= |
338 | //function : Print |
339 | //purpose : |
340 | //======================================================================= |
341 | |
342 | static void Print(const Handle(Geom_BSplineCurve)& B, |
343 | Standard_OStream& OS, |
344 | const Standard_Boolean compact) |
345 | { |
346 | if (compact) |
347 | OS << BSPLINE << " "; |
348 | else |
349 | OS << "BSplineCurve"; |
350 | |
351 | |
352 | Standard_Boolean rational = B->IsRational(); |
353 | if (compact) |
354 | OS << (rational ? 1 : 0) << " "; |
355 | else { |
356 | if (rational) |
357 | OS << " rational"; |
358 | } |
359 | |
360 | Standard_Boolean periodic = B->IsPeriodic(); |
361 | if (compact) |
362 | OS << (periodic ? 1 : 0)<< " "; |
363 | else { |
364 | if (periodic) |
365 | OS << " periodic"; |
366 | } |
367 | |
368 | // poles and weights |
369 | Standard_Integer i,degree,nbpoles,nbknots; |
370 | degree = B->Degree(); |
371 | nbpoles = B->NbPoles(); |
372 | nbknots = B->NbKnots(); |
373 | if (!compact) OS << "\n Degree "; |
374 | else OS << " "; |
375 | OS << degree; |
376 | if (!compact) OS << ","; |
377 | OS << " "; |
378 | OS << nbpoles; |
379 | if (!compact) OS << " Poles,"; |
380 | OS << " "; |
381 | OS << nbknots << " "; |
382 | if (!compact) OS << " Knots\n"; |
383 | |
384 | if (!compact) OS << "Poles :\n"; |
385 | for (i = 1; i <= nbpoles; i++) { |
04232180 |
386 | if (!compact) OS << "\n "<<std::setw(2)<<i<<" : "; |
7fd59977 |
387 | else OS << " "; |
388 | Print(B->Pole(i),OS,compact); |
389 | if (rational) |
390 | OS << " " << B->Weight(i); |
391 | } |
392 | OS << "\n"; |
393 | |
394 | if (!compact) OS << "Knots :\n"; |
395 | for (i = 1; i <= nbknots; i++) { |
04232180 |
396 | if (!compact) OS << "\n "<<std::setw(2)<<i<<" : "; |
7fd59977 |
397 | OS << " " << B->Knot(i) << " " << B->Multiplicity(i); |
398 | } |
399 | |
400 | OS << "\n"; |
401 | if (!compact) OS << "\n"; |
402 | } |
403 | |
404 | //======================================================================= |
405 | //function : Print |
406 | //purpose : |
407 | //======================================================================= |
408 | |
409 | static void Print(const Handle(Geom_TrimmedCurve)& C, |
410 | Standard_OStream& OS, |
411 | const Standard_Boolean compact) |
412 | { |
413 | if (compact) |
414 | OS << TRIMMED << " "; |
415 | else |
416 | OS << "Trimmed curve\n"; |
417 | if (!compact) OS << "Parameters : "; |
418 | OS << C->FirstParameter() << " " << C->LastParameter() << "\n"; |
419 | if (!compact) OS << "Basis curve :\n"; |
420 | GeomTools_CurveSet::PrintCurve(C->BasisCurve(),OS,compact); |
421 | } |
422 | |
423 | //======================================================================= |
424 | //function : Print |
425 | //purpose : |
426 | //======================================================================= |
427 | |
428 | static void Print(const Handle(Geom_OffsetCurve)& C, |
429 | Standard_OStream& OS, |
430 | const Standard_Boolean compact) |
431 | { |
432 | if (compact) |
433 | OS << OFFSET << " "; |
434 | else |
435 | OS << "OffsetCurve"; |
436 | if (!compact) OS << "Offset : "; |
437 | OS << C->Offset() << "\n"; |
438 | if (!compact) OS << "Direction : "; |
439 | Print(C->Direction(),OS,compact); |
440 | OS << "\n"; |
441 | if (!compact) OS << "Basis curve :\n"; |
442 | GeomTools_CurveSet::PrintCurve(C->BasisCurve(),OS,compact); |
443 | } |
444 | |
445 | //======================================================================= |
446 | //function : Print |
447 | //purpose : |
448 | //======================================================================= |
449 | |
450 | void GeomTools_CurveSet::PrintCurve(const Handle(Geom_Curve)& C, |
451 | Standard_OStream& OS, |
452 | const Standard_Boolean compact) |
453 | { |
454 | Handle(Standard_Type) TheType = C->DynamicType(); |
455 | |
456 | if ( TheType ==STANDARD_TYPE(Geom_Line)) { |
457 | Print(Handle(Geom_Line)::DownCast(C),OS,compact); |
458 | } |
459 | else if ( TheType == STANDARD_TYPE(Geom_Circle)) { |
460 | Print(Handle(Geom_Circle)::DownCast(C),OS,compact); |
461 | } |
462 | else if ( TheType == STANDARD_TYPE(Geom_Ellipse)) { |
463 | Print(Handle(Geom_Ellipse)::DownCast(C),OS,compact); |
464 | } |
465 | else if ( TheType == STANDARD_TYPE(Geom_Parabola)) { |
466 | Print(Handle(Geom_Parabola)::DownCast(C),OS,compact); |
467 | } |
468 | else if ( TheType == STANDARD_TYPE(Geom_Hyperbola)) { |
469 | Print(Handle(Geom_Hyperbola)::DownCast(C),OS,compact); |
470 | } |
471 | else if ( TheType == STANDARD_TYPE(Geom_BezierCurve)) { |
472 | Print(Handle(Geom_BezierCurve)::DownCast(C),OS,compact); |
473 | } |
474 | else if ( TheType == STANDARD_TYPE(Geom_BSplineCurve)) { |
475 | Print(Handle(Geom_BSplineCurve)::DownCast(C),OS,compact); |
476 | } |
477 | else if ( TheType == STANDARD_TYPE(Geom_TrimmedCurve)) { |
478 | Print(Handle(Geom_TrimmedCurve)::DownCast(C),OS,compact); |
479 | } |
480 | else if ( TheType == STANDARD_TYPE(Geom_OffsetCurve)) { |
481 | Print(Handle(Geom_OffsetCurve)::DownCast(C),OS,compact); |
482 | } |
483 | else { |
484 | GeomTools::GetUndefinedTypeHandler()->PrintCurve(C,OS,compact); |
485 | //if (!compact) |
486 | // OS << "****** UNKNOWN CURVE TYPE ******\n"; |
487 | //else |
04232180 |
488 | // std::cout << "****** UNKNOWN CURVE TYPE ******" << std::endl; |
7fd59977 |
489 | } |
490 | } |
491 | |
492 | //======================================================================= |
493 | //function : Dump |
494 | //purpose : |
495 | //======================================================================= |
496 | |
497 | void GeomTools_CurveSet::Dump(Standard_OStream& OS)const |
498 | { |
499 | Standard_Integer i, nbsurf = myMap.Extent(); |
500 | OS << "\n -------\n"; |
501 | OS << "Dump of "<< nbsurf << " Curves "; |
502 | OS << "\n -------\n\n"; |
503 | |
504 | for (i = 1; i <= nbsurf; i++) { |
04232180 |
505 | OS << std::setw(4) << i << " : "; |
7fd59977 |
506 | PrintCurve(Handle(Geom_Curve)::DownCast(myMap(i)),OS,Standard_False); |
507 | } |
508 | } |
509 | |
510 | |
511 | //======================================================================= |
512 | //function : Write |
513 | //purpose : |
514 | //======================================================================= |
515 | |
6d8f9f4a |
516 | void GeomTools_CurveSet::Write (Standard_OStream& OS, |
517 | const Handle(Message_ProgressIndicator)& theProgress)const |
7fd59977 |
518 | { |
60be1f9b |
519 | std::streamsize prec = OS.precision(17); |
7fd59977 |
520 | |
521 | Standard_Integer i, nbcurve = myMap.Extent(); |
522 | OS << "Curves "<< nbcurve << "\n"; |
523 | //OCC19559 |
6d8f9f4a |
524 | Message_ProgressSentry PS(theProgress, "3D Curves", 0, nbcurve, 1); |
7fd59977 |
525 | for (i = 1; i <= nbcurve && PS.More(); i++, PS.Next()) { |
7fd59977 |
526 | PrintCurve(Handle(Geom_Curve)::DownCast(myMap(i)),OS,Standard_True); |
527 | } |
7fd59977 |
528 | OS.precision(prec); |
529 | } |
530 | |
531 | |
532 | //======================================================================= |
533 | //function : ReadPnt |
534 | //purpose : |
535 | //======================================================================= |
536 | |
537 | static Standard_IStream& operator>>(Standard_IStream& IS, gp_Pnt& P) |
538 | { |
539 | Standard_Real X=0.,Y=0.,Z=0.; |
71598a83 |
540 | GeomTools::GetReal(IS, X); |
541 | GeomTools::GetReal(IS, Y); |
542 | GeomTools::GetReal(IS, Z); |
7fd59977 |
543 | P.SetCoord(X,Y,Z); |
544 | return IS; |
545 | } |
546 | |
547 | //======================================================================= |
548 | //function : ReadDir |
549 | //purpose : |
550 | //======================================================================= |
551 | |
552 | static Standard_IStream& operator>>(Standard_IStream& IS, gp_Dir& D) |
553 | { |
554 | Standard_Real X=0.,Y=0.,Z=0.; |
71598a83 |
555 | GeomTools::GetReal(IS, X); |
556 | GeomTools::GetReal(IS, Y); |
557 | GeomTools::GetReal(IS, Z); |
7fd59977 |
558 | D.SetCoord(X,Y,Z); |
559 | return IS; |
560 | } |
561 | |
562 | |
563 | //======================================================================= |
564 | //function : ReadCurve |
565 | //purpose : |
566 | //======================================================================= |
567 | |
568 | static Standard_IStream& operator>>(Standard_IStream& IS, |
569 | Handle(Geom_Line)& L) |
570 | { |
571 | gp_Pnt P(0.,0.,0.); |
572 | gp_Dir AX(1.,0.,0.); |
573 | IS >> P >> AX; |
574 | L = new Geom_Line(P,AX); |
575 | return IS; |
576 | } |
577 | |
578 | //======================================================================= |
579 | //function : ReadCurve |
580 | //purpose : |
581 | //======================================================================= |
582 | |
583 | static Standard_IStream& operator>>(Standard_IStream& IS, |
584 | Handle(Geom_Circle)& C) |
585 | { |
586 | gp_Pnt P(0.,0.,0.); |
587 | gp_Dir A(1.,0.,0.),AX(1.,0.,0.),AY(1.,0.,0.); |
588 | Standard_Real R=0.; |
71598a83 |
589 | IS >> P >> A >> AX >> AY; |
590 | GeomTools::GetReal(IS, R); |
7fd59977 |
591 | C = new Geom_Circle(gp_Ax2(P,A,AX),R); |
592 | return IS; |
593 | } |
594 | |
595 | //======================================================================= |
596 | //function : ReadCurve |
597 | //purpose : |
598 | //======================================================================= |
599 | |
600 | static Standard_IStream& operator>>(Standard_IStream& IS, |
601 | Handle(Geom_Ellipse)& E) |
602 | { |
603 | gp_Pnt P(0.,0.,0.); |
604 | gp_Dir A(1.,0.,0.),AX(1.,0.,0.),AY(1.,0.,0.); |
605 | Standard_Real R1=0.,R2=0.; |
71598a83 |
606 | IS >> P >> A >> AX >> AY; |
607 | GeomTools::GetReal(IS, R1); |
608 | GeomTools::GetReal(IS, R2); |
7fd59977 |
609 | E = new Geom_Ellipse(gp_Ax2(P,A,AX),R1,R2); |
610 | return IS; |
611 | } |
612 | |
613 | //======================================================================= |
614 | //function : ReadCurve |
615 | //purpose : |
616 | //======================================================================= |
617 | |
618 | static Standard_IStream& operator>>(Standard_IStream& IS, |
619 | Handle(Geom_Parabola)& C) |
620 | { |
621 | gp_Pnt P(0.,0.,0.); |
622 | gp_Dir A(1.,0.,0.),AX(1.,0.,0.),AY(1.,0.,0.); |
623 | Standard_Real R1=0.; |
71598a83 |
624 | IS >> P >> A >> AX >> AY; |
625 | GeomTools::GetReal(IS, R1); |
7fd59977 |
626 | C = new Geom_Parabola(gp_Ax2(P,A,AX),R1); |
627 | return IS; |
628 | } |
629 | |
630 | //======================================================================= |
631 | //function : ReadCurve |
632 | //purpose : |
633 | //======================================================================= |
634 | |
635 | static Standard_IStream& operator>>(Standard_IStream& IS, |
636 | Handle(Geom_Hyperbola)& H) |
637 | { |
638 | gp_Pnt P(0.,0.,0.); |
639 | gp_Dir A(1.,0.,0.),AX(1.,0.,0.),AY(1.,0.,0.); |
640 | Standard_Real R1=0.,R2=0.; |
71598a83 |
641 | IS >> P >> A >> AX >> AY; |
642 | GeomTools::GetReal(IS, R1); |
643 | GeomTools::GetReal(IS, R2); |
7fd59977 |
644 | H = new Geom_Hyperbola(gp_Ax2(P,A,AX),R1,R2); |
645 | return IS; |
646 | } |
647 | |
648 | //======================================================================= |
649 | //function : ReadCurve |
650 | //purpose : |
651 | //======================================================================= |
652 | |
653 | static Standard_IStream& operator>>(Standard_IStream& IS, |
654 | Handle(Geom_BezierCurve)& B) |
655 | { |
656 | Standard_Boolean rational=Standard_False; |
657 | IS >> rational; |
658 | |
659 | // poles and weights |
660 | Standard_Integer i=0,degree=0; |
661 | IS >> degree; |
662 | |
663 | TColgp_Array1OfPnt poles(1,degree+1); |
664 | TColStd_Array1OfReal weights(1,degree+1); |
665 | |
666 | for (i = 1; i <= degree+1; i++) { |
667 | IS >> poles(i); |
668 | if (rational) |
71598a83 |
669 | GeomTools::GetReal(IS, weights(i)); |
7fd59977 |
670 | } |
671 | |
672 | if (rational) |
673 | B = new Geom_BezierCurve(poles,weights); |
674 | else |
675 | B = new Geom_BezierCurve(poles); |
676 | |
677 | return IS; |
678 | } |
679 | |
680 | //======================================================================= |
681 | //function : ReadCurve |
682 | //purpose : |
683 | //======================================================================= |
684 | |
685 | static Standard_IStream& operator>>(Standard_IStream& IS, |
686 | Handle(Geom_BSplineCurve)& B) |
687 | { |
688 | |
689 | Standard_Boolean rational=Standard_False,periodic=Standard_False; |
690 | IS >> rational >> periodic; |
691 | |
692 | // poles and weights |
693 | Standard_Integer i=0,degree=0,nbpoles=0,nbknots=0; |
694 | IS >> degree >> nbpoles >> nbknots; |
695 | |
696 | TColgp_Array1OfPnt poles(1,nbpoles); |
697 | TColStd_Array1OfReal weights(1,nbpoles); |
698 | |
699 | for (i = 1; i <= nbpoles; i++) { |
700 | IS >> poles(i); |
701 | if (rational) |
71598a83 |
702 | GeomTools::GetReal(IS, weights(i)); |
7fd59977 |
703 | } |
704 | |
705 | TColStd_Array1OfReal knots(1,nbknots); |
706 | TColStd_Array1OfInteger mults(1,nbknots); |
707 | |
708 | for (i = 1; i <= nbknots; i++) { |
71598a83 |
709 | GeomTools::GetReal(IS, knots(i)); |
710 | IS >> mults(i); |
7fd59977 |
711 | } |
712 | |
713 | if (rational) |
714 | B = new Geom_BSplineCurve(poles,weights,knots,mults,degree,periodic); |
715 | else |
716 | B = new Geom_BSplineCurve(poles,knots,mults,degree,periodic); |
717 | |
718 | return IS; |
719 | } |
720 | |
721 | //======================================================================= |
722 | //function : ReadCurve |
723 | //purpose : |
724 | //======================================================================= |
725 | |
726 | static Standard_IStream& operator>>(Standard_IStream& IS, |
727 | Handle(Geom_TrimmedCurve)& C) |
728 | { |
729 | Standard_Real p1=0.,p2=0.; |
71598a83 |
730 | GeomTools::GetReal(IS, p1); |
731 | GeomTools::GetReal(IS, p2); |
aa00364d |
732 | Handle(Geom_Curve) BC = GeomTools_CurveSet::ReadCurve(IS); |
7fd59977 |
733 | C = new Geom_TrimmedCurve(BC,p1,p2); |
734 | return IS; |
735 | } |
736 | |
737 | //======================================================================= |
738 | //function : ReadCurve |
739 | //purpose : |
740 | //======================================================================= |
741 | |
742 | static Standard_IStream& operator>>(Standard_IStream& IS, |
743 | Handle(Geom_OffsetCurve)& C) |
744 | { |
745 | Standard_Real p=0.; |
71598a83 |
746 | GeomTools::GetReal(IS, p); |
7fd59977 |
747 | gp_Dir D(1.,0.,0.); |
748 | IS >> D; |
aa00364d |
749 | Handle(Geom_Curve) BC = GeomTools_CurveSet::ReadCurve(IS); |
7fd59977 |
750 | C = new Geom_OffsetCurve(BC,p,D); |
751 | return IS; |
752 | } |
753 | |
754 | //======================================================================= |
755 | //function : ReadCurve |
756 | //purpose : |
757 | //======================================================================= |
758 | |
aa00364d |
759 | Handle(Geom_Curve) GeomTools_CurveSet::ReadCurve (Standard_IStream& IS) |
7fd59977 |
760 | { |
761 | Standard_Integer ctype; |
762 | |
aa00364d |
763 | Handle(Geom_Curve) C; |
7fd59977 |
764 | try { |
765 | OCC_CATCH_SIGNALS |
766 | IS >> ctype; |
767 | switch (ctype) { |
768 | |
769 | case LINE : |
770 | { |
771 | Handle(Geom_Line) CC; |
772 | IS >> CC; |
773 | C = CC; |
774 | } |
775 | break; |
776 | |
777 | case CIRCLE : |
778 | { |
779 | Handle(Geom_Circle) CC; |
780 | IS >> CC; |
781 | C = CC; |
782 | } |
783 | break; |
784 | |
785 | case ELLIPSE : |
786 | { |
787 | Handle(Geom_Ellipse) CC; |
788 | IS >> CC; |
789 | C = CC; |
790 | } |
791 | break; |
792 | |
793 | case PARABOLA : |
794 | { |
795 | Handle(Geom_Parabola) CC; |
796 | IS >> CC; |
797 | C = CC; |
798 | } |
799 | break; |
800 | |
801 | case HYPERBOLA : |
802 | { |
803 | Handle(Geom_Hyperbola) CC; |
804 | IS >> CC; |
805 | C = CC; |
806 | } |
807 | break; |
808 | |
809 | case BEZIER : |
810 | { |
811 | Handle(Geom_BezierCurve) CC; |
812 | IS >> CC; |
813 | C = CC; |
814 | } |
815 | break; |
816 | |
817 | case BSPLINE : |
818 | { |
819 | Handle(Geom_BSplineCurve) CC; |
820 | IS >> CC; |
821 | C = CC; |
822 | } |
823 | break; |
824 | |
825 | case TRIMMED : |
826 | { |
827 | Handle(Geom_TrimmedCurve) CC; |
828 | IS >> CC; |
829 | C = CC; |
830 | } |
831 | break; |
832 | |
833 | case OFFSET : |
834 | { |
835 | Handle(Geom_OffsetCurve) CC; |
836 | IS >> CC; |
837 | C = CC; |
838 | } |
839 | break; |
840 | |
841 | default: |
842 | { |
843 | Handle(Geom_Curve) CC; |
844 | GeomTools::GetUndefinedTypeHandler()->ReadCurve(ctype,IS,CC); |
845 | C = CC; |
846 | } |
847 | } |
848 | } |
9775fa61 |
849 | catch(Standard_Failure const& anException) { |
0797d9d3 |
850 | #ifdef OCCT_DEBUG |
04232180 |
851 | std::cout <<"EXCEPTION in GeomTools_CurveSet::ReadCurve(..)!!!" << std::endl; |
852 | std::cout << anException << std::endl; |
7fd59977 |
853 | #endif |
9775fa61 |
854 | (void)anException; |
7fd59977 |
855 | } |
aa00364d |
856 | return C; |
7fd59977 |
857 | } |
858 | |
859 | //======================================================================= |
860 | //function : Read |
861 | //purpose : |
862 | //======================================================================= |
863 | |
6d8f9f4a |
864 | void GeomTools_CurveSet::Read (Standard_IStream& IS, |
865 | const Handle(Message_ProgressIndicator)& theProgress) |
7fd59977 |
866 | { |
867 | char buffer[255]; |
868 | IS >> buffer; |
869 | if (strcmp(buffer,"Curves")) { |
04232180 |
870 | std::cout << "Not a Curve table"<<std::endl; |
7fd59977 |
871 | return; |
872 | } |
873 | |
7fd59977 |
874 | Standard_Integer i, nbcurve; |
875 | IS >> nbcurve; |
876 | //OCC19559 |
6d8f9f4a |
877 | Message_ProgressSentry PS(theProgress, "3D Curves", 0, nbcurve, 1); |
7fd59977 |
878 | for (i = 1; i <= nbcurve && PS.More(); i++, PS.Next()) { |
aa00364d |
879 | Handle(Geom_Curve) C = GeomTools_CurveSet::ReadCurve (IS); |
7fd59977 |
880 | myMap.Add(C); |
881 | } |
882 | } |