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