projects / occt.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
0031939: Coding - correction of spelling errors in comments
[occt.git] / src / IGESToBRep / IGESToBRep_TopoSurface.cxx
CommitLineData
973c2be1 1// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e 2//
973c2be1 3// This file is part of Open CASCADE Technology software library.
b311480e 4//
d5f74e42 5// This library is free software; you can redistribute it and/or modify it under
6// the terms of the GNU Lesser General Public License version 2.1 as published
973c2be1 7// by the Free Software Foundation, with special exception defined in the file
8// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
9// distribution for complete text of the license and disclaimer of any warranty.
b311480e 10//
973c2be1 11// Alternatively, this file may be used under the terms of Open CASCADE
12// commercial license or contractual agreement.
b311480e 13
7fd59977 14//=======================================================================
7fd59977 15//modified:
16// Integration to ensure SCCS base integrity
17// 21.12.98 rln, gka S4054
18// 28.12.98 dce S3767 New messaging system
19//#55,#56 rln 25.12.98 UKI60878
20//:l1 abv 10.01.99: USA60022 7289: fix missing seam
21//#63 rln 19.01.99 UKI60878 no offset if C0 surface is converted into the grid of C1 surfaces
22//%13 pdn 15.02.99 USA60293 entities 792, 8604 .. handling of C0 ruled surfaces, tabulated cylindres,
23// and surfaces of revolution.
24//:p4 abv, pdn 23.02.99: PRO9234 #15720: call BRepTools::Update() for faces
25//%14 pdn 24.02.99 implementing of ShapeFix_Face on IGES
26// pdn 17.04.99 S4181: Implementing of reading IGES elementary surfaces.
27// pdn 10.05.99 S4137: Using modified ShapeDivide tools
28//#11 smh 22.12.99 BUC60625 Transform axis.
29//#12 smh 12.12.99 FRA62468 - Using conversion to B-Spline for Offset surface
30//=======================================================================
31
42cf5bc1 32#include <BRep_Builder.hxx>
33#include <BRep_Tool.hxx>
7fd59977 34#include <BRepAdaptor_Curve.hxx>
153fee01 35#include <BRepBuilderAPI_GTransform.hxx>
42cf5bc1 36#include <BRepBuilderAPI_MakeFace.hxx>
7fd59977 37#include <BRepFill.hxx>
42cf5bc1 38#include <BRepGProp.hxx>
7fd59977 39#include <BRepLib_MakeFace.hxx>
42cf5bc1 40#include <BRepPrimAPI_MakePrism.hxx>
41#include <BRepPrimAPI_MakeRevol.hxx>
42#include <BRepTools.hxx>
43#include <BSplCLib.hxx>
44#include <ElCLib.hxx>
7fd59977 45#include <ElSLib.hxx>
7fd59977 46#include <Geom2d_Curve.hxx>
42cf5bc1 47#include <Geom2d_Line.hxx>
7fd59977 48#include <Geom_BezierCurve.hxx>
42cf5bc1 49#include <Geom_BSplineCurve.hxx>
50#include <Geom_BSplineSurface.hxx>
7fd59977 51#include <Geom_ConicalSurface.hxx>
42cf5bc1 52#include <Geom_Curve.hxx>
7fd59977 53#include <Geom_CylindricalSurface.hxx>
42cf5bc1 54#include <Geom_Line.hxx>
7fd59977 55#include <Geom_OffsetSurface.hxx>
56#include <Geom_Plane.hxx>
7fd59977 57#include <Geom_RectangularTrimmedSurface.hxx>
58#include <Geom_SphericalSurface.hxx>
42cf5bc1 59#include <Geom_Surface.hxx>
60#include <Geom_SurfaceOfLinearExtrusion.hxx>
61#include <Geom_SurfaceOfRevolution.hxx>
7fd59977 62#include <Geom_ToroidalSurface.hxx>
42cf5bc1 63#include <Geom_TrimmedCurve.hxx>
64#include <GeomAbs_Shape.hxx>
65#include <GeomConvert.hxx>
66#include <GeomLib.hxx>
67#include <gp.hxx>
68#include <gp_Ax1.hxx>
69#include <gp_Cone.hxx>
70#include <gp_Cylinder.hxx>
71#include <gp_Dir.hxx>
72#include <gp_Dir2d.hxx>
73#include <gp_GTrsf.hxx>
74#include <gp_Pln.hxx>
75#include <gp_Pnt.hxx>
76#include <gp_Sphere.hxx>
77#include <gp_Torus.hxx>
78#include <gp_Trsf.hxx>
79#include <gp_Trsf2d.hxx>
80#include <gp_Vec.hxx>
81#include <GProp.hxx>
82#include <GProp_GProps.hxx>
83#include <IGESBasic_SingleParent.hxx>
7fd59977 84#include <IGESData_IGESEntity.hxx>
42cf5bc1 85#include <IGESData_IGESModel.hxx>
7fd59977 86#include <IGESData_ToolLocation.hxx>
153fee01 87#include <IGESData_TransfEntity.hxx>
7fd59977 88#include <IGESGeom_BoundedSurface.hxx>
42cf5bc1 89#include <IGESGeom_BSplineSurface.hxx>
90#include <IGESGeom_CircularArc.hxx>
7fd59977 91#include <IGESGeom_CurveOnSurface.hxx>
42cf5bc1 92#include <IGESGeom_Direction.hxx>
7fd59977 93#include <IGESGeom_Line.hxx>
42cf5bc1 94#include <IGESGeom_OffsetSurface.hxx>
95#include <IGESGeom_Plane.hxx>
7fd59977 96#include <IGESGeom_Point.hxx>
42cf5bc1 97#include <IGESGeom_RuledSurface.hxx>
98#include <IGESGeom_SurfaceOfRevolution.hxx>
99#include <IGESGeom_TabulatedCylinder.hxx>
100#include <IGESGeom_TrimmedSurface.hxx>
7fd59977 101#include <IGESSolid_ConicalSurface.hxx>
42cf5bc1 102#include <IGESSolid_CylindricalSurface.hxx>
7fd59977 103#include <IGESSolid_PlaneSurface.hxx>
42cf5bc1 104#include <IGESSolid_SphericalSurface.hxx>
105#include <IGESSolid_ToroidalSurface.hxx>
106#include <IGESToBRep.hxx>
107#include <IGESToBRep_BasicCurve.hxx>
108#include <IGESToBRep_BasicSurface.hxx>
109#include <IGESToBRep_CurveAndSurface.hxx>
110#include <IGESToBRep_TopoCurve.hxx>
111#include <IGESToBRep_TopoSurface.hxx>
7fd59977 112#include <Interface_Macros.hxx>
42cf5bc1 113#include <Message_Msg.hxx>
7fd59977 114#include <Precision.hxx>
42cf5bc1 115#include <ShapeAlgo.hxx>
116#include <ShapeAlgo_AlgoContainer.hxx>
117#include <ShapeAnalysis.hxx>
118#include <ShapeExtend_WireData.hxx>
119#include <ShapeFix_Wire.hxx>
120#include <Standard_ErrorHandler.hxx>
121#include <TColgp_Array1OfPnt.hxx>
122#include <TColStd_Array1OfInteger.hxx>
123#include <TColStd_Array1OfReal.hxx>
7fd59977 124#include <TopAbs.hxx>
42cf5bc1 125#include <TopExp.hxx>
126#include <TopExp_Explorer.hxx>
127#include <TopLoc_Location.hxx>
7fd59977 128#include <TopoDS.hxx>
129#include <TopoDS_Edge.hxx>
130#include <TopoDS_Face.hxx>
131#include <TopoDS_Iterator.hxx>
132#include <TopoDS_Shape.hxx>
133#include <TopoDS_Shell.hxx>
134#include <TopoDS_Vertex.hxx>
135#include <TopoDS_Wire.hxx>
136
7fd59977 137#include <stdio.h>
42cf5bc1 138//S4054: ShapeTool_MakeWire -> ShapeExtend_WireData //:g8: BRepLib_MakeWire -> ShapeTool_MakeWire
7fd59977 139//:e3
42cf5bc1 140//#16
7fd59977 141//S4054
7fd59977 142//S3767
7fd59977 143//=======================================================================
144//function : IGESToBRep_TopoSurface
145//purpose :
146//=======================================================================
7fd59977 147IGESToBRep_TopoSurface::IGESToBRep_TopoSurface()
148 :IGESToBRep_CurveAndSurface()
149{
150}
151
152
153//=======================================================================
154//function : IGESToBRep_TopoSurface
155//purpose :
156//=======================================================================
157
158IGESToBRep_TopoSurface::IGESToBRep_TopoSurface
159 (const IGESToBRep_CurveAndSurface& CS)
160 :IGESToBRep_CurveAndSurface(CS)
161{
162}
163
164
165//=======================================================================
166//function : IGESToBRep_TopoSurface
167//purpose :
168//=======================================================================
169
170IGESToBRep_TopoSurface::IGESToBRep_TopoSurface
171 (const Standard_Real eps,
172 const Standard_Real epsCoeff,
173 const Standard_Real epsGeom,
174 const Standard_Boolean mode,
175 const Standard_Boolean modeapprox,
176 const Standard_Boolean optimized)
177 :IGESToBRep_CurveAndSurface(eps, epsCoeff, epsGeom, mode,
178 modeapprox, optimized)
179{
180}
181
182static Standard_Boolean extractCurve3d (const TopoDS_Shape& theEdges,
183 Handle(Geom_Curve)& theCurve)
184{
185 TopExp_Explorer anExp(theEdges, TopAbs_EDGE);
186 Standard_Integer howMuch = 0;
1d47d8d0 187 Standard_Real f = 0., l = 0.;
7fd59977 188 for (; anExp.More(); anExp.Next()) {
189 TopoDS_Edge anEdge = TopoDS::Edge(anExp.Current());
190 if (anEdge.IsNull())
191 continue;
192 howMuch++;
193 theCurve = BRep_Tool::Curve(anEdge, f, l);
194 }
195 if ( howMuch != 1 || theCurve.IsNull() )
196 return Standard_False;
197
198 if ( f != theCurve->FirstParameter() || l != theCurve->LastParameter() )
199 theCurve = new Geom_TrimmedCurve ( theCurve, f, l );
200 return Standard_True;
201}
202
203
204//=======================================================================
205//function : TransferTopoSurface
206//purpose :
207//=======================================================================
208
209TopoDS_Shape IGESToBRep_TopoSurface::TransferTopoSurface
210 (const Handle(IGESData_IGESEntity)& st)
211{ // Declaration of messages//
212 // DCE 22/12/98
213 //Message_Msg msg1005("IGES_1005");
214 ////////////////////////////
215
216 TopoDS_Shape res;
217 TheULength = 1.;
218 //S4054
219 if (st.IsNull()) {
220 Message_Msg msg1005("IGES_1005");
221 SendFail(st, msg1005);
222 }
223 ////modified by jgv, 20.11.2009 for OCC21487///
224 else if (HasShapeResult(st))
225 {
226 res = GetShapeResult(st);
227 return res;
228 }
229 ///////////////////////////////////////////////
230 else if (IGESToBRep::IsBasicSurface(st)) {
231 res = TransferTopoBasicSurface(st);
232 }
233 else if (st->IsKind(STANDARD_TYPE(IGESGeom_TrimmedSurface))) {
234 DeclareAndCast(IGESGeom_TrimmedSurface, st144, st);
235 res = TransferTrimmedSurface(st144);
236 }
237 else if (st->IsKind(STANDARD_TYPE(IGESGeom_SurfaceOfRevolution))) {
238 DeclareAndCast(IGESGeom_SurfaceOfRevolution, st120, st);
239 res = TransferSurfaceOfRevolution(st120);
240 }
241 else if (st->IsKind(STANDARD_TYPE(IGESGeom_TabulatedCylinder))) {
242 DeclareAndCast(IGESGeom_TabulatedCylinder, st122, st);
243 res = TransferTabulatedCylinder(st122);
244 }
245 else if (st->IsKind(STANDARD_TYPE(IGESGeom_RuledSurface))) {
246 DeclareAndCast(IGESGeom_RuledSurface, st118, st);
247 res = TransferRuledSurface(st118);
248 }
249 else if (st->IsKind(STANDARD_TYPE(IGESGeom_Plane))) {
250 DeclareAndCast(IGESGeom_Plane, st108, st);
251 res = TransferPlane(st108);
252 }
253 else if (st->IsKind(STANDARD_TYPE(IGESGeom_BoundedSurface))) {
254 DeclareAndCast(IGESGeom_BoundedSurface, st143, st);
255 res = TransferBoundedSurface(st143);
256 }
257 else if (st->IsKind(STANDARD_TYPE(IGESGeom_OffsetSurface))) {
258 DeclareAndCast(IGESGeom_OffsetSurface, st140, st);
259 res = TransferOffsetSurface(st140);
260 }
261 //S4181 pdn IGESSolid_PlaneSurface recognized as basic surface
262 else if (st->IsKind(STANDARD_TYPE(IGESBasic_SingleParent))) {
263 DeclareAndCast(IGESBasic_SingleParent,st402_9,st);
264 res = TransferPerforate(st402_9); // limite : Planes seulement
265 }
266 else {
267 // AddFail(st, "The IGESEntity is not a Topologic Surface.");
268 }
269 SetShapeResult (st, res);
270 return res;
271}
272
273//=======================================================================
274//function : TransferTopoBasicSurface
275//purpose :
276//=======================================================================
277
278TopoDS_Shape IGESToBRep_TopoSurface::TransferTopoBasicSurface
279 (const Handle(IGESData_IGESEntity)& st)
280{ // Declaration of messages//
281 // DCE 22/12/98
282 //Message_Msg msg1005("IGES_1005");
283 ////////////////////////////
284
285 TopoDS_Shape res;
286
287 if (st.IsNull()) {
288 Message_Msg msg1005("IGES_1005");
289 SendFail(st, msg1005);
290 return res;
291 }
292 if (!IGESToBRep::IsBasicSurface(st)) {
293 // AddFail(st, "BasicSurface Transfer Error : Not Allowed IGESEntity"); This message can not occur.
294 return res;
295 }
296
297 IGESToBRep_BasicSurface BS(*this);
298
299 Handle(Geom_Surface) surf = BS.TransferBasicSurface(st);
300 if (surf.IsNull()) {
301 // AddFail(st, "Surface Conversion Error"); Messages have ever been Added in the called function.
302 return res;
303 }
304
305 //#9 rln 26/02/98 UKI60106
306 if (surf->Continuity() < GeomAbs_C1) {
307 Message_Msg msg1250("IGES_1250");
308 SendWarning(st, msg1250);
309 }
310 if(surf->IsKind(STANDARD_TYPE(Geom_Plane))){
311 BRep_Builder B;
312 TopoDS_Face plane;
313 B.MakeFace(plane);
314 B.UpdateFace(plane, surf, TopLoc_Location(), Precision::Confusion());
315 res = plane;
316 }
317 else {
1c72dff6 318 BRepLib_MakeFace makeFace(surf, Precision::Confusion());
7fd59977 319 res = makeFace.Face();
320 }
321
322 if (st->HasTransf()) {
323 gp_Trsf trsf;
324 SetEpsilon(1.E-04);
325 if (IGESData_ToolLocation::ConvertLocation
326 (GetEpsilon(),st->CompoundLocation(),trsf,GetUnitFactor())) {
327 TopLoc_Location locFace(trsf);
328 res.Move(locFace);
329 }
330 else {
331 Message_Msg msg1035("IGES_1035");
332 SendWarning(st, msg1035);
333 }
334 }
335 return res;
336}
337
338
339//=======================================================================
340//function : TransferRuledSurface
341//purpose :
342//=======================================================================
343static void reparamBSpline(Handle(Geom_Curve)& curve,
344 const Standard_Real First,
345 const Standard_Real Last)
346{
347 Handle (Geom_BSplineCurve) bscurve;
348 if (!curve->IsKind (STANDARD_TYPE (Geom_BSplineCurve))) {
349 if (curve->FirstParameter() < First || curve->LastParameter() > Last)
350 curve = new Geom_TrimmedCurve (curve, First, Last);
351 bscurve = GeomConvert::CurveToBSplineCurve (curve, Convert_RationalC1);
352 }
353 else {
354 bscurve = Handle (Geom_BSplineCurve)::DownCast (curve);
355 bscurve->Segment (First, Last);
356 }
357
358 if (bscurve.IsNull())
359 return;
360
361 TColStd_Array1OfReal Knots(1, bscurve->NbKnots());
362 bscurve->Knots(Knots);
363 BSplCLib::Reparametrize (0., 1., Knots);
364 bscurve->SetKnots(Knots);
365 curve = bscurve;
366}
367
368static void ReparamCurve(TopoDS_Edge& edge)
369{
370 TopLoc_Location L;
371 Standard_Real First, Last;
372
373 Handle (Geom_Curve) curve = Handle (Geom_Curve)::DownCast (BRep_Tool::Curve ( edge, L, First, Last )->Copy());
374 //if ( Abs (First) <= Precision::PConfusion() && Abs (Last - 1.) <= Precision::PConfusion() ) return;
375 if(!curve->IsKind(STANDARD_TYPE(Geom_Line))) return;
376
377 reparamBSpline( curve, First, Last );
378
379 BRep_Builder B;
380 B.UpdateEdge ( edge, curve, L, Precision::Confusion() );
381 B.Range ( edge, 0., 1 );
382}
383
384
385//=======================================================================
386//function : TransferRuledSurface
387//purpose :
388//=======================================================================
389
390TopoDS_Shape IGESToBRep_TopoSurface::TransferRuledSurface
391 (const Handle(IGESGeom_RuledSurface)& st)
392{ // Declaration of messages//
393 // DCE 22/12/98
394 //Message_Msg msg1005("IGES_1005");
395 ////////////////////////////////
396 TopoDS_Shape res;
397
398 if (st.IsNull()) {
399 Message_Msg msg1005("IGES_1005");
400 SendFail(st, msg1005);
401 return res;
402 }
403
404 IGESToBRep_TopoCurve TC(*this);
405 //%13 pdn 12.02.99
406 TC.SetContinuity (0);
407 Handle(IGESData_IGESEntity) igesCurve1 = st->FirstCurve();
408 Handle(IGESData_IGESEntity) igesCurve2 = st->SecondCurve();
409
410 if (igesCurve1.IsNull()) {
411 Message_Msg msg148("XSTEP_148");
412 SendFail(st, msg148); // Curve Reading Error : Null IGESEntity
413 return res;
414 }
415 if (igesCurve2.IsNull()) {
416 Message_Msg msg149("XSTEP_149");
417 SendFail(st, msg149); // Curve Reading Error : Null IGESEntity
418 return res;
419 }
420
421 Standard_Integer nbEdges1, nbEdges2;
422 TopoDS_Shape shape1, shape2;
423 TopoDS_Wire wire1, wire2;
424 TopoDS_Wire newWire1, newWire2;
425 //TopoDS_Edge edge1, edge2; // skl
426
427 if (IGESToBRep::IsTopoCurve(igesCurve1)) {
428 shape1 = TC.TransferTopoCurve(igesCurve1);
429 if (shape1.IsNull()) {
430 Message_Msg msg1156("IGES_1156");
431 const Standard_CString typeName(igesCurve1->DynamicType()->Name());
432 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesCurve1);
433 msg1156.Arg(typeName);
434 msg1156.Arg(label);
435 SendFail(st, msg1156);
436 return res;
437 }
438
439 //%13 pdn 15.02.99
440 //added by rln on 03/12/97
441 //if shape1 is a wire it means that the curve1 in file was of continuity C0
442 //in order to get a face instead of shell when to BRepFill shape1
21c7c457 443 //should be retransferred with contionuity C0 (to get an edge). Once shape1
7fd59977 444 //has been built with C0, it is useless to require C1 from shape2 because
445 //anyway resulting surface was of continuity C0. Thus shape2 is built with C0
446// if (shape1.ShapeType() != TopAbs_EDGE) {
447// TC.SetContinuity (0);
448// shape1 = TC.TransferTopoCurve(igesCurve1);
449// if (shape1.IsNull()) {
450// Message_Msg msg1156("IGES_1156");
451// const Standard_CString typeName(igesCurve1->DynamicType()->Name());
452// Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesCurve1);
453// msg1156.Arg(typeName);
454// msg1156.Arg(label);
455// SendFail(st, msg1156);
456// return res;
457// }
458// else {
459// Message_Msg msg1250("IGES_1250");
460// SendWarning (st, msg1250); // RuledSurface was built with continuity C0
461// }
462// }
463
464 TopAbs_ShapeEnum shapeEnum1 = shape1.ShapeType();
465 switch (shapeEnum1) {
466 case TopAbs_EDGE :
467 {
468 TopoDS_Edge edge1 = TopoDS::Edge(shape1);
469 ReparamCurve(edge1);
470 nbEdges1 = 1;
471 }
472 break;
473 case TopAbs_WIRE :
474 {
475 wire1 = TopoDS::Wire(shape1);
476 nbEdges1 = 0;
477 for (TopoDS_Iterator hulot(wire1); hulot.More(); hulot.Next()) {
478 TopoDS_Edge edge1 = TopoDS::Edge(hulot.Value());
479 ReparamCurve(edge1);
480 nbEdges1++;
481 }
482 }
483 break;
484 default:
485 {
486 // AddFail(st, "Curve Conversion Error."); This message can not occur.
487 return res;
488 }
489 //break; //szv#4:S4163:12Mar99 unreachable
490 }
491 }
492 else {
493 Message_Msg msg148("XSTEP_148");
494 SendFail(st, msg148);
495 // Curve Type not Allowed.
496 return res;
497 }
498
499 if (IGESToBRep::IsTopoCurve(igesCurve2)) {
500 shape2 = TC.TransferTopoCurve(igesCurve2);
501 // dirflg = 0 join first to first, last to last
502 // dirflg = 1 join first to last, last to first
503
504 if (shape2.IsNull()) {
505 Message_Msg msg1156("IGES_1156");
506 const Standard_CString typeName(igesCurve2->DynamicType()->Name());
507 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesCurve2);
508 msg1156.Arg(typeName);
509 msg1156.Arg(label);
510 SendFail(st, msg1156);
511 // Curve Conversion Error.
512 return res;
513 }
514 Standard_Integer dirflag = st->DirectionFlag ();
515 // if (dirflag == 1){ // gka BUC60685
516
517 // shape2.Reverse();
518 //}
519
520 TopAbs_ShapeEnum shapeEnum2 = shape2.ShapeType();
521 switch (shapeEnum2) {
522 case TopAbs_EDGE :
523 {
524 TopoDS_Edge edge2 = TopoDS::Edge(shape2);
525 ReparamCurve(edge2);
526 if (dirflag == 1)
527 shape2.Reverse();
528 nbEdges2 = 1;
529 }
530 break;
531 case TopAbs_WIRE :
532 {
533 wire2 = TopoDS::Wire(shape2);
534 nbEdges2 = 0;
535 for (TopoDS_Iterator cousto(wire2); cousto.More(); cousto.Next()) {
536 TopoDS_Edge edge2 = TopoDS::Edge(cousto.Value());
537 ReparamCurve(edge2);
538 nbEdges2++;
539 }
540 if (dirflag == 1) { //gka BUC60685
541 Handle(ShapeExtend_WireData) sewd2 = new ShapeExtend_WireData;
542 sewd2->Add(shape2);
543 sewd2->Reverse();
544 wire2 = sewd2->Wire();
545 }
546 }
547 break;
548 default:
549 {
550 // AddFail(st, "Curve Conversion Error.");
551 return res;
552 }
553 //break; //szv#4:S4163:12Mar99 unreachable
554 }
555 }
556 else {
557 Message_Msg msg149("XSTEP_149");
558 SendFail(st, msg149);
559 // Curve Type not Allowed
560 return res;
561 }
562
563
564 if (nbEdges1 != nbEdges2) {
565 if (nbEdges1 == 1) {
566 Handle(ShapeExtend_WireData) sewd1 = new ShapeExtend_WireData;
567 sewd1->Add(shape1);
568 wire1 = sewd1->Wire();
569 }
570 else if (nbEdges2 == 1) {
571 Handle(ShapeExtend_WireData) sewd2 = new ShapeExtend_WireData;
572 sewd2->Add(shape2);
573 wire2 = sewd2->Wire();
574 }
575
576 if (!st->IsRuledByParameter()) {
577 // AddWarning (st,"Compute by parametric constant ratio");
578 }
579 if (!ShapeAlgo::AlgoContainer()->HomoWires
580 (wire1, wire2, newWire1, newWire2, st->IsRuledByParameter())) {
581 Message_Msg msg1255("IGES_1255");// "Ruled Surface Construction Error");
582 SendFail(st, msg1255);
583 return res;
584 }
585 nbEdges1 = 2; // a number > 1
586 }
587
588 else if (nbEdges1 != 1) {
589 newWire1 = wire1;
590 newWire2 = wire2;
591 }
592
593
594 if (nbEdges1 == 1) {
595
596 //:e3 abv 31 Mar 98: UK4.igs 3170: ruled surface with directixes - line
597 // In IGES, line is parametrised [0;1] - this should have been kept !
598 // Let us detect the case and remake curve as bspline [0;1]
599 for ( Standard_Integer i=1; i <=2; i++ ) {
600 //#43 rln 20.11.98 S4054 BUC50047 entity D463 (circles as generatrices [0, 2*PI])
601 //reparameterisation should be for all curves not with range [0, 1] (see IGES)
602 TopoDS_Edge edge = TopoDS::Edge ( i==1 ? shape1 : shape2 );
603 //ReparamCurve(edge);
604 TopLoc_Location L;
605 Standard_Real First, Last;
606 Handle (Geom_Curve) curve = Handle (Geom_Curve)::DownCast (BRep_Tool::Curve ( edge, L, First, Last )->Copy());
607 if ( Abs (First) <= Precision::PConfusion() && Abs (Last - 1.) <= Precision::PConfusion() ) continue;
608
609 Handle (Geom_BSplineCurve) bscurve;
610 if (!curve->IsKind (STANDARD_TYPE (Geom_BSplineCurve))) {
611 if (curve->FirstParameter() < First || curve->LastParameter() > Last)
612 curve = new Geom_TrimmedCurve (curve, First, Last);
613 bscurve = GeomConvert::CurveToBSplineCurve (curve, Convert_RationalC1);
614 }
615 else {
616 bscurve = Handle (Geom_BSplineCurve)::DownCast (curve);
617 bscurve->Segment (First, Last);
618 }
619 TColStd_Array1OfReal Knots(1, bscurve->NbKnots());
620 bscurve->Knots(Knots);
621 BSplCLib::Reparametrize (0., 1., Knots);
622 bscurve->SetKnots(Knots);
623
624 BRep_Builder B;
625 B.UpdateEdge ( edge, bscurve, L, Precision::Confusion() );
626 B.Range ( edge, 0., 1 );
627 if ( i ==1 ) shape1 = edge;
628 else shape2 = edge;
629 }
630
631 res = BRepFill::Face(TopoDS::Edge(shape1), TopoDS::Edge(shape2));
632 Handle(Geom_Surface) surf = BRep_Tool::Surface(TopoDS::Face(res));
633 if(surf->Continuity()==GeomAbs_C0) {
634 Message_Msg msg1250("IGES_1250");
635 SendWarning (st, msg1250);
636 }
637 }
638 else {
639 res = BRepFill::Shell(newWire1, newWire2);
640 }
641 if (res.IsNull()) {
642 Message_Msg msg1255("IGES_1255");// "Ruled Surface Construction Error");
643 SendFail(st, msg1255);
644 return res;
645 }
646
647
648 if (st->HasTransf()) {
649 gp_Trsf trsf;
650 SetEpsilon(1.E-04);
651 if (IGESData_ToolLocation::ConvertLocation
652 (GetEpsilon(),st->CompoundLocation(), trsf,GetUnitFactor())) {
653 TopLoc_Location shapeLoc(trsf);
654 res.Move(shapeLoc);
655 }
656 else {
657 Message_Msg msg1035("IGES_1035");
658 SendWarning(st,msg1035); // Transformation : not a similarity
659 }
660 }
661 return res;
662}
663
664
665//=======================================================================
666//function : TransferSurfaceOfRevolution
667//purpose :
668//=======================================================================
669
670TopoDS_Shape IGESToBRep_TopoSurface::TransferSurfaceOfRevolution
671 (const Handle(IGESGeom_SurfaceOfRevolution)& st)
672{ // Declaration of messages//
673 // DCE 22/12/98
674 //Message_Msg msg1005("IGES_1005");
675 ////////////////////////////////
676 TopoDS_Shape res;
677 if (st.IsNull()) {
678 Message_Msg msg1005("IGES_1005");
679 SendFail(st, msg1005);
680 return res;
681 }
682
683 IGESToBRep_TopoCurve TC(*this);
684 IGESToBRep_BasicCurve BC(*this);
685 Handle(IGESData_IGESEntity) igesGeneratrix = st->Generatrix();
686 Handle(IGESGeom_Line) igesAxis = st->AxisOfRevolution();
687
688 if (igesGeneratrix.IsNull() || !IGESToBRep::IsTopoCurve(igesGeneratrix) ) {
689 Message_Msg msg153("XSTEP_153");
690 SendFail(st, msg153);
691 // Generatrix Reading Error : Null IGESEntity
692 // Generatrix : Not Allowed IGESEntity.
693 return res;
694 }
695
696 DeclareAndCast(IGESGeom_Line,srgen,st->Generatrix());
697 if (!srgen.IsNull()) {
698 gp_Pnt gen1 = srgen->StartPoint();
699 gp_Pnt gen2 = srgen->EndPoint();
700 TheULength = gen1.Distance(gen2)*GetUnitFactor();
701 }
702
703 if (igesAxis.IsNull()) {
704 Message_Msg msg152("XSTEP_152");
705 SendFail(st, msg152);
706 return res;
707 }
708
709 //%13 pdn 15.02.99
710 TC.SetContinuity(0);
711 TopoDS_Shape generatrix = TC.TransferTopoCurve(igesGeneratrix);
712 if (generatrix.IsNull()) {
713 Message_Msg msg1156("IGES_1156");
714 const Standard_CString typeName("generatrix");
715 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesGeneratrix);
716 msg1156.Arg(typeName);
717 msg1156.Arg(label);
718 SendFail(st, msg1156);
719 // Generatrix Conversion Error.
720 return res;
721 }
722
723 gp_Trsf startLoc;
724 gp_Pnt pt1 = igesAxis->TransformedStartPoint(); //smh#11
725 gp_Pnt pt2 = igesAxis->TransformedEndPoint(); //smh#11
726 pt1.Scale(gp_Pnt(0,0,0),GetUnitFactor());
727 pt2.Scale(gp_Pnt(0,0,0),GetUnitFactor());
728 //#30 rln 19.10.98 To keep IGES surface normal CAS.CADE axis = reversed IGES axis
729 //CAS.CADE SA = 2*PI - IGES TA
730 //CAS.CADE TA = 2*PI - IGES SA
731 //gp_Ax1 revolAxis(pt1, gp_Dir(gp_Vec(pt1, pt2)));
732 //Standard_Real startAngle = st->StartAngle();
733 //Standard_Real endAngle = st->EndAngle();
734 gp_Ax1 revolAxis(pt1, gp_Dir( gp_Vec (pt2, pt1)));
c6541a0c
D		 735 Standard_Real startAngle = 2 * M_PI - st->EndAngle();
736 Standard_Real endAngle = 2 * M_PI - st->StartAngle();
7fd59977 737 Standard_Real deltaAngle = endAngle - startAngle;
c6541a0c
D		 738 Standard_Boolean IsFullAngle = ( deltaAngle > 2.*M_PI-Precision::PConfusion() );
739 if (IsFullAngle) deltaAngle = 2.*M_PI; // ** CKY 18-SEP-1996
7fd59977 740 // il faudra translater les courbes 2d de startAngle pour
741 // etre en phase IGES et BRep
742 startLoc.SetRotation(revolAxis, startAngle);
743 generatrix.Move(startLoc);
744
745 // PTV file D44-11325-6.igs OCC660 depends on OCC450
746 // PTV 29.05.2002 OCC450 create Surface of Revolution by native API
747 // file NIC_file5.igs
748 // (BRepPrimAPI_MakeRevol replace surface of revolution by plane then 3D and 2D curves are inconsistent;
749 // and 3D is ignored. As result shape is rectangle instead circle shape.
750 Handle(Geom_Curve) aBasisCurve;
9d109e39 751
7fd59977 752 {
9d109e39 753 try
754 {
755 OCC_CATCH_SIGNALS
756 if (extractCurve3d(generatrix, aBasisCurve))
757 {
758 BRepBuilderAPI_MakeFace aMakeF;
759 Handle(Geom_Surface) aResultSurf =
760 new Geom_SurfaceOfRevolution(aBasisCurve, revolAxis);
761
762 if ( !aResultSurf.IsNull())
763 {
764 if (!IsFullAngle)
765 {
766 const Standard_Real VF = aBasisCurve->FirstParameter();
767 const Standard_Real VL = aBasisCurve->LastParameter();
768
769 // PTV 29.08.2002 begin of OCC663 Trim surface by correct parameters
770 const Standard_Real UF = 0;
771 const Standard_Real UL = endAngle - startAngle;
772 // PTV 29.08.2002 end of OCC663
773
440906ef 774 aMakeF.Init(aResultSurf, UF, UL, VF, VL, Precision::Confusion());
9d109e39 775 }//if (!IsFullAngle)
776 else
777 {
440906ef 778 aMakeF.Init(aResultSurf, Standard_True, Precision::Confusion());
9d109e39 779 }
780
781 if (aMakeF.IsDone())
782 res = aMakeF.Face();
783 }//if ( !aResultSurf.IsNull())
784 }//if (extractCurve3d(generatrix, aBasisCurve))
7fd59977 785 }
9775fa61 786 catch (Standard_Failure const& anException) {
0797d9d3 787#ifdef OCCT_DEBUG
04232180 788 std::cout << "Warning: IgesToBRep_TopoSurface::"
9d109e39 789 "TransferSurfaceOfRevolution(): exception by Geom: ";
04232180 790 anException.Print ( std::cout ); std::cout << std::endl;
7fd59977 791#endif
9775fa61 792 (void)anException;
9d109e39 793 }//catch (Standard_Failure)
7fd59977 794 }
795
796 if ( res.IsNull() ) {
797 // do as usual.
798
799 BRepPrimAPI_MakeRevol revol(generatrix, revolAxis, deltaAngle);
800//mjm: si debug IsDone()est fait :
801// if (!revol.IsDone()) {
802// AddFail(st, "Revol Construction Error.");
803// return res;
804// }
805 res = revol.Shape();
806 }
807 //%13 pdn 15.02.99
808 if (res.ShapeType() == TopAbs_FACE) {
809 Handle(Geom_Surface) surf = BRep_Tool::Surface(TopoDS::Face(res));
810 if(surf->Continuity()==GeomAbs_C0) {
811 Message_Msg msg1250("IGES_1250");
812 SendWarning (st, msg1250);
813 }
814 }
815
816 if (st->HasTransf()) {
817 gp_Trsf trsf;
818 SetEpsilon(1.E-04);
819 if (IGESData_ToolLocation::ConvertLocation
820 (GetEpsilon(), st->CompoundLocation(), trsf, GetUnitFactor())) {
821 TopLoc_Location shapeLoc(trsf);
822 res.Move(shapeLoc);
823 }
824 else {
825 Message_Msg msg1035("IGES_1035");
826 SendWarning(st,msg1035); // Transformation : not a similarity
827 }
828 }
829
830 return res;
831}
832
833
834//=======================================================================
835//function : TransferTabulatedCylinder
836//purpose :
837//=======================================================================
838
839TopoDS_Shape IGESToBRep_TopoSurface::TransferTabulatedCylinder
840 (const Handle(IGESGeom_TabulatedCylinder)& st)
841{ // Declaration of messages//
842 // DCE 22/12/98
843 //Message_Msg msg1005("IGES_1005");
844 ////////////////////////////////
845 TopoDS_Shape res;
846 if (st.IsNull()) {
847 Message_Msg msg1005("IGES_1005");
848 SendFail(st, msg1005);
849 return res;
850 }
851
852 IGESToBRep_TopoCurve TC(*this);
853// TopoDS_Edge firstEdge;//commented by rln on 02/12/97
854
855 Handle(IGESData_IGESEntity) igesDirectrix = st->Directrix();
856 if (igesDirectrix.IsNull() || !IGESToBRep::IsTopoCurve(igesDirectrix) ) {
857 Message_Msg msg153("XSTEP_153");
858 SendFail(st, msg153);
859 // Directrix Reading Error : Null IGESEntity
860 //Directrix, not allowed IGESEntity
861 return res;
862 }
863
864 //%13 pdn 15.02.99
865 TC.SetContinuity(0);
866 TopoDS_Shape directrix = TC.TransferTopoCurve(igesDirectrix);
867 if (directrix.IsNull()) {
868 Message_Msg msg1156("IGES_1156");
869 const Standard_CString typeName("directrix");
870 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesDirectrix);
871 msg1156.Arg(typeName);
872 msg1156.Arg(label);
873 SendFail(st, msg1156);
874 // Directrix Conversion Error.
875 return res;
876 }
877
878 //modified by rln on 03/12/97
879 //TopoDS_Vertex firstVertex = TopExp::FirstVertex(firstEdge);
880 TopoDS_Vertex firstVertex, lastVertex;
881 ShapeAnalysis::FindBounds (directrix, firstVertex, lastVertex);
882 gp_Pnt pt1 = BRep_Tool::Pnt(firstVertex);
883 gp_Pnt pt2 = st->EndPoint();
884 pt2.Scale(gp_Pnt(0,0,0),GetUnitFactor());
885
886 TheULength = pt1.Distance(pt2);
887 if(TheULength < Precision::Confusion()) {
888 Message_Msg msg("Tabulated cylinder with zero length");
889 SendFail (st, msg); // TabulatedCylinder was built with continuity C0
890 return res;
891 }
892
893 // PTV file D44-11325-6.igs OCC660 depends on OCC450
894 // PTV 29.05.2002 OCC450 create Surface of LinearExtrusion by native API
895 // see description about problem in Surface of Revolution
896 Handle(Geom_Curve) aBasisCurve;
897 {
898 try {
899 OCC_CATCH_SIGNALS
900 if (extractCurve3d(directrix, aBasisCurve)) {
7fd59977 901 gp_Vec dir (pt1, pt2);
902 Handle(Geom_Surface) aResultSurf =
903 new Geom_SurfaceOfLinearExtrusion(aBasisCurve, dir);
904 if (!aResultSurf.IsNull()) {
9d109e39 905 //aResultSurf =
906 // new Geom_RectangularTrimmedSurface(aResultSurf,
907 // aBasisCurve->FirstParameter(),
908 // aBasisCurve->LastParameter(),
909 // 0., dir.Magnitude() );
910 BRepBuilderAPI_MakeFace aMakeF(aResultSurf, aBasisCurve->FirstParameter(),
7fd59977 911 aBasisCurve->LastParameter(),
9d109e39 912 0., dir.Magnitude(),
913 Precision::Confusion());
914 if (aMakeF.IsDone())
915 res = aMakeF.Face();
7fd59977 916 }
917 }
918 }
9775fa61 919 catch (Standard_Failure const& anException) {
0797d9d3 920#ifdef OCCT_DEBUG
04232180 921 std::cout << "Warning: IgesToBRep_TopoSurface::TransferTabulatedCylinder(): exception by Geom: ";
922 anException.Print ( std::cout ); std::cout << std::endl;
7fd59977 923#endif
9775fa61 924 (void)anException;
7fd59977 925 }
926 }
927
928 if ( res.IsNull() ) {
929 // do as usual.
930 BRepPrimAPI_MakePrism prism(directrix, gp_Vec(pt1, pt2));
931//mjm: si debug IsDone() est fait
932// if (!prism.IsDone()) {
933// AddFail(st, "Prism Construction Error.");
934// return res;
935// }
936 res = prism.Shape();
937 }
938 //#16 rln 08/04/98 coq-inf-support.igs entity 2105
939 //CAS.CADE can parametrize SurfaceOfLinearExtrusion with generatrix opposite to Vec(pt1, pt2)
940 //and with parametrization V > 0, while in IGES TabulatedCylinder is parametrized with positive V
941 //direction exactly in the direction Vec(pt1, pt2)
942 if (res.ShapeType() == TopAbs_FACE) {
943 Standard_Real UMin, UMax, VMin, VMax;
944 BRepTools::UVBounds (TopoDS::Face (res), UMin, UMax, VMin, VMax);
945 if (VMax <= Precision::PConfusion() && VMin < -Precision::PConfusion()) {
946 TheULength *= -1;
947 res.Reverse();
948 }
949 Handle(Geom_Surface) surf = BRep_Tool::Surface(TopoDS::Face(res));
950 if(surf->Continuity()==GeomAbs_C0) {
951 Message_Msg msg1250("IGES_1250");
952 SendWarning (st, msg1250);
953 }
954 }
955
956 if (st->HasTransf()) {
957 gp_Trsf trsf;
958 SetEpsilon(1.E-04);
959 if (IGESData_ToolLocation::ConvertLocation
960 (GetEpsilon(),st->CompoundLocation(), trsf, GetUnitFactor())) {
961 TopLoc_Location shapeLoc(trsf);
962 res.Move(shapeLoc);
963 }
964 else {
965 Message_Msg msg1035("IGES_1035");
966 SendWarning(st,msg1035); // Transformation : not a similarity
967 }
968 }
969 return res;
970}
971
972
973//=======================================================================
974//function : TransferOffsetSurface
975//purpose :
976//=======================================================================
977
978TopoDS_Shape IGESToBRep_TopoSurface::TransferOffsetSurface
979 (const Handle(IGESGeom_OffsetSurface)& st)
980{ // Declaration of messages//
981 // DCE 22/12/98
982 //Message_Msg msg1005("IGES_1005");
983 ////////////////////////////////
984 TopoDS_Shape res;
985 if (st.IsNull()) {
986 Message_Msg msg1005("IGES_1005");
987 SendFail(st, msg1005);
988 return res;
989 }
990
991 TopoDS_Shape igesShape;
992 TopoDS_Face face;
993 TopLoc_Location basisLoc;
994
995 Handle (IGESData_IGESEntity) igesSrf = st->Surface();
996 if (igesSrf.IsNull() || !IGESToBRep::IsTopoSurface(igesSrf) ) {
997 Message_Msg msg164("XSTEP_164");
998 SendFail(st, msg164);
999 // Basis Surface Reading Error : Null IGESEntity
1000 // Basis Surface Transfer Error : Not Allowed IGESEntity
1001 return res;
1002 }
1003
1004 igesShape = TransferTopoSurface(igesSrf);
1005 if (igesShape.IsNull()) {
1006 Message_Msg msg1156("IGES_1156");
1007 const Standard_CString typeName("basis surface");
1008 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesSrf);
1009 msg1156.Arg(typeName);
1010 msg1156.Arg(label);
1011 SendFail(st, msg1156); // Basis Surface Conversion Error.
1012 return res;
1013 }
1014
1015 TopAbs_ShapeEnum shapeEnum = igesShape.ShapeType();
1016 switch (shapeEnum) {
1017 case TopAbs_FACE :
1018 {
1019 face = TopoDS::Face(igesShape);
1020 break;
1021 }
1022 case TopAbs_SHELL :
1023 {
7fd59977 1024 TopoDS_Iterator dabovil(igesShape);
1025 if (dabovil.More()) {
21c7c457 1026 SendWarning(st, "The First Surface only will be transferred.");
b1811c1d 1027 face = TopoDS::Face(dabovil.Value());
1028 break;
7fd59977 1029 }
1030 /* else AddF("... */
1031 }
b1811c1d 1032 Standard_FALLTHROUGH
7fd59977 1033 default:
1034 {
1035 Message_Msg msg1156("IGES_1156");
1036 const Standard_CString typeName("basis surface");
1037 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesSrf);
1038 msg1156.Arg(typeName);
1039 msg1156.Arg(label);
1040 SendFail(st, msg1156); // Basis Surface Conversion Error.
1041 return res;
1042 }
1043 }
1044
1045
1046 //Handle(Geom_Surface) geomSupport = BRep_Tool::Surface(face, basisLoc);
1047 // attention on ne peut construire une Geom_OffsetSurface que
1048 // si la surface de base est au moins C1, sinon on plante !
1049 //#56 rln 25.12.98 UKI60878 entity D593 (Offset surface on C0 B-Spline)
1050 //Trying to eliminate previous limitation on processing only C1 surfaces
1051 Handle(Geom_Surface) geomSupport = BRep_Tool::Surface(face);
1052 Handle(Geom_OffsetSurface) basisSrf;
1053
1054 if (geomSupport->IsKind(STANDARD_TYPE(Geom_OffsetSurface))) {
1055 DeclareAndCast(Geom_OffsetSurface, geom140, geomSupport);
1056 geom140->SetOffsetValue(basisSrf->Offset() +
1057 st->Distance()*GetUnitFactor());
1058 basisSrf = geom140;
1059 }
1060 else {
1061 if (geomSupport->Continuity() == GeomAbs_C0) {
1062 res = ShapeAlgo::AlgoContainer()->C0ShapeToC1Shape (face, Abs (st->Distance()) * GetUnitFactor());
1063 if(res.ShapeType()!=TopAbs_FACE) {
1064 Message_Msg msg1266("IGES_1266");
1065 SendFail(st, msg1266);//Basis surface is C0-continuous and cannot be corrected to C1-continuous.
1066 return res;
1067 }
1068 else {
1069 geomSupport = BRep_Tool::Surface (TopoDS::Face(res));
1070 if (geomSupport->Continuity() == GeomAbs_C0) {
1071 Message_Msg msg1266("IGES_1266");
1072 SendFail(st, msg1266);//Basis surface is C0-continuous and cannot be corrected to C1-continuous.
1073 res.Nullify();
1074 return res;
1075 }
1076 }
1077 Message_Msg msg1267("IGES_1267");
1078 SendWarning(st, msg1267);//Basis surface is C0-continuous but was corrected to C1-continuous
1079 }
1080 //smh#12
1081 if (res.IsNull()) res = face;
1082 geomSupport = BRep_Tool::Surface (TopoDS::Face(res));
1083 Standard_Real umin, umax, vmin, vmax;
1084 geomSupport->Bounds (umin, umax, vmin, vmax);
1085 if (Precision::IsInfinite (umin) || Precision::IsInfinite (umax) ||
1086 Precision::IsInfinite (vmin) || Precision::IsInfinite (vmax)) {
1087 // convert to C1 B-Spline
1088 BRepTools::UVBounds (face, umin, umax, vmin, vmax);
1089 Handle(Geom_RectangularTrimmedSurface) TS = new Geom_RectangularTrimmedSurface (geomSupport, umin, umax, vmin, vmax);
1090 Handle (Geom_BSplineSurface) BS = ShapeAlgo::AlgoContainer()->ConvertSurfaceToBSpline(TS, umin, umax, vmin, vmax);
1091 if (BS.IsNull() || BS->Continuity() == GeomAbs_C0) {
1092 Message_Msg msg1265("IGES_1265");
1093 SendFail(st, msg1265); // OffsetSurface Construction Error.
1094 return res;
1095 }
1096 else {
1097 geomSupport = BS;
1098 }
1099 }
1100 basisSrf = new Geom_OffsetSurface(geomSupport, st->Distance()*GetUnitFactor());
1101 }
1102
1c72dff6 1103 BRepLib_MakeFace MF(basisSrf, Precision::Confusion());
7fd59977 1104 if(!MF.IsDone()) {
1105 Message_Msg msg1265("IGES_1265");
1106 SendFail(st, msg1265); // OffsetSurface Construction Error.
1107 return res;
1108 }
1109
1110 res = MF.Face();
1111
1112 if (st->HasTransf()) {
1113 gp_Trsf trsf;
1114 SetEpsilon(1.E-04);
1115 if (IGESData_ToolLocation::ConvertLocation
1116 (GetEpsilon(),st->CompoundLocation(),trsf, GetUnitFactor())) {
1117 TopLoc_Location loc2(trsf);
1118 res.Move(loc2);
1119 }
1120 else {
1121 Message_Msg msg1035("IGES_1035");
1122 SendWarning(st,msg1035); // Transformation : not a similarity
1123 }
1124 }
1125 return res;
1126}
1127
1128
1129//=======================================================================
1130//function : TransferTrimmedSurface
1131//purpose :
1132//=======================================================================
1133
1134TopoDS_Shape IGESToBRep_TopoSurface::TransferTrimmedSurface
1135 (const Handle(IGESGeom_TrimmedSurface)& st)
1136{ // Declaration of messages//
1137 // DCE 22/12/98
1138 //Message_Msg msg1005("IGES_1005");
1139 ////////////////////////////////
1140 TopoDS_Shape res;
1141 if (st.IsNull()) {
1142 Message_Msg msg1005("IGES_1005");
1143 SendFail(st, msg1005);
1144 return res;
1145 }
1146
1147 TopAbs_ShapeEnum shapeEnum;
1148 IGESToBRep_TopoCurve TC(*this);
1149
1150 Handle (IGESData_IGESEntity) igesSurface = st->Surface();
1151 if (igesSurface.IsNull() || !IGESToBRep::IsTopoSurface(igesSurface) ) {
1152 Message_Msg msg169("XSTEP_169");
1153 SendFail(st, msg169);
1154 // BasicSurface Transfer Error : Null IGESEntity
1155 // Basis Surface, not Allowed IGESEntity.
1156 return res;
1157 }
1158 gp_Trsf2d trans;
1159 Standard_Real uFact;
1160 TopoDS_Face face, faceres;
1161
1162 TopoDS_Shape myshape = ParamSurface(igesSurface, trans, uFact);
1163
1164 if (!myshape.IsNull()) {
1165 shapeEnum = myshape.ShapeType();
1166 switch (shapeEnum) {
1167 case TopAbs_FACE :
1168 {
1169 face = TopoDS::Face(myshape);
1170 faceres = face;
1171 break;
1172 }
1173 case TopAbs_SHELL :
1174 {
1175 TopoDS_Iterator IT(myshape);
1176 Standard_Integer nbfaces = 0;
1177 for (; IT.More(); IT.Next()) {
1178 nbfaces++;
1179 face = TopoDS::Face(IT.Value());
1180 faceres = face;
1181 }
1182 //szv#4:S4163:12Mar99 optimized
1183 if (nbfaces != 1) {
1184 Message_Msg msg1156("IGES_1156");
1185 const Standard_CString typeName("basis surface");
1186 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesSurface);
1187 msg1156.Arg(typeName);
1188 msg1156.Arg(label);
1189 SendFail(st, msg1156); // Not Implemented Trimmed Composite Surface.
1190 return myshape;
1191 }
1192 }
1193 break;
1194 default:
1195 {
1196 Message_Msg msg1156("IGES_1156");
1197 const Standard_CString typeName("basis surface");
1198 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesSurface);
1199 msg1156.Arg(typeName);
1200 msg1156.Arg(label);
1201 SendFail(st, msg1156); // Basis Surface Conversion Error.
1202 return res;
1203 }
1204 }
1205 }
1206 else {
1207 return res;
1208 }
1209
1210 //obtaining a surface
1211 TopLoc_Location L;
1212 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(face, L);
1213 TC.SetSurface(aSurf);
1214
1215 if (st->HasOuterContour()) {
1216 face.EmptyCopy();
1217 TopoDS_Shape myshape1 = TC.TransferCurveOnFace (face, st->OuterContour(), trans, uFact, Standard_False);
1218 // si ca se passe mal , on recupere au moins la face avec NaturalRestriction
fdabc211 1219 if (myshape1 .IsNull()) {
1220 face = faceres;
1221 BRep_Builder B;
1222 B.NaturalRestriction(face,Standard_False);
1223 }
7fd59977 1224 }
1225 for (Standard_Integer i = 1; i <= st->NbInnerContours(); i++) {
1226 TopoDS_Shape myshape2 = TC.TransferCurveOnFace (face, st->InnerContour(i), trans, uFact, Standard_False);
1227 }
153fee01 1228
1229 Handle(IGESData_TransfEntity) aTransf = st->Transf();
1230 if (!aTransf.IsNull()) {
1231 // make transformation
1232 gp_GTrsf aGT = aTransf->Value();
1233 gp_XYZ aTrans = aGT.TranslationPart();
1234 gp_Mat aMat = aGT.VectorialPart();
1235 Standard_Real s1 = aMat.Value(1, 1)*aMat.Value(1, 1) + aMat.Value(2, 1)*aMat.Value(2, 1) + aMat.Value(3, 1)*aMat.Value(3, 1);
1236 Standard_Real s2 = aMat.Value(1, 2)*aMat.Value(1, 2) + aMat.Value(2, 2)*aMat.Value(2, 2) + aMat.Value(3, 2)*aMat.Value(3, 2);
1237 Standard_Real s3 = aMat.Value(1, 3)*aMat.Value(1, 3) + aMat.Value(2, 3)*aMat.Value(2, 3) + aMat.Value(3, 3)*aMat.Value(3, 3);
1238 if (fabs(s1 - s2) > Precision::Confusion() || fabs(s1 - s3) > Precision::Confusion()) {
1239 BRepBuilderAPI_GTransform aTransform(aGT);
1240 aTransform.Perform(face, Standard_True);
1241 if (aTransform.IsDone()) {
1242 if (aTransform.Shape().ShapeType() == TopAbs_FACE) {
1243 face = TopoDS::Face(aTransform.Shape());
1244 }
1245 }
1246 }
1247 else {
1248 Standard_Real tmpVal = fabs(aMat.Value(1, 1) - 1.) + fabs(aMat.Value(1, 2)) + fabs(aMat.Value(1, 3)) +
1249 fabs(aMat.Value(2, 1)) + fabs(aMat.Value(2, 2) - 1.) + fabs(aMat.Value(2, 3)) +
1250 fabs(aMat.Value(3, 1)) + fabs(aMat.Value(3, 2)) + fabs(aMat.Value(3, 3) - 1.);
1251 if ((tmpVal + aTrans.Modulus()) > Precision::Confusion()) {
1252 // not Identity
1253 gp_Trsf aT;
1254 aT.SetValues(
1255 aMat.Value(1, 1), aMat.Value(1, 2), aMat.Value(1, 3), aTrans.X(),
1256 aMat.Value(2, 1), aMat.Value(2, 2), aMat.Value(2, 3), aTrans.Y(),
1257 aMat.Value(3, 1), aMat.Value(3, 2), aMat.Value(3, 3), aTrans.Z());
1258 TopLoc_Location aLoc(aT);
1259 face.Move(aLoc);
1260 }
1261 }
1262 }
1263
7fd59977 1264 BRepTools::Update ( face ); //:p4
1265 //%16 pdn 08.04.99
1266 return face;
1267}
1268
1269
1270//=======================================================================
1271//function : TransferBoundedSurface
1272//purpose :
1273//=======================================================================
1274
1275TopoDS_Shape IGESToBRep_TopoSurface::TransferBoundedSurface
1276 (const Handle(IGESGeom_BoundedSurface)& st)
1277{ // Declaration of messages//
1278 // DCE 22/12/98
1279 //Message_Msg msg1005("IGES_1005");
1280 ////////////////////////////////
1281 TopoDS_Shape res;
1282 if (st.IsNull()) {
1283 Message_Msg msg1005("IGES_1005");
1284 SendFail(st, msg1005);
1285 return res;
1286 }
1287
1288 if (st->RepresentationType()==0) {
1289 Message_Msg msg1275("IGES_1275");
1290 SendWarning(st, msg1275);
1291 // Model Space Representation Not Implemented : the result will be the basis surface
1292 }
1293
1294 TopAbs_ShapeEnum shapeEnum;
1295 IGESToBRep_TopoCurve TC(*this);
1296 Handle (IGESData_IGESEntity) igesSrf = st->Surface();
1297 if (igesSrf.IsNull() || !IGESToBRep::IsTopoSurface(igesSrf) ) {
1298 Message_Msg msg166("XSTEP_166");
1299 SendFail( st, msg166);
1300 // Basis Surface Transfer Error : Null IGESEntity.
1301 // Basis Surface Transfer Error : Not Allowed IGESEntity.
1302 return res;
1303 }
1304 gp_Trsf2d trans;
1305 Standard_Real uFact;
1306 TopoDS_Face face;
1307
1308 TopoDS_Shape myshape = ParamSurface(igesSrf, trans, uFact);
1309
1310 if (myshape.IsNull()) {
1311 //#55 rln 24.12.98 UKI60878 entity D593
0797d9d3 1312#ifdef OCCT_DEBUG
04232180 1313 std::cout << "Fail: IGESToBRep_TopoSurface::TransferBoundedSurface UntrimmedSurface is translated into Null" << std::endl;
7fd59977 1314#endif
1315 return res;
1316 }
1317 else {
1318 shapeEnum = myshape.ShapeType();
1319 switch (shapeEnum) {
1320 case TopAbs_FACE :
1321 {
1322 face = TopoDS::Face(myshape);
1323 }
1324 break;
1325 case TopAbs_SHELL :
1326 {
1327 TopoDS_Iterator IT(myshape);
1328 Standard_Integer nbfaces = 0;
1329 for (; IT.More(); IT.Next()) {
1330 nbfaces++;
1331 face = TopoDS::Face(IT.Value());
1332 }
1333 //szv#4:S4163:12Mar99 optimized
1334 if (nbfaces != 1) {
1335 Message_Msg msg1156("IGES_1156");
1336 const Standard_CString typeName("basis surface");
1337 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesSrf);
1338 msg1156.Arg(typeName);
1339 msg1156.Arg(label);
1340 SendFail(st, msg1156);
1341 // Not Implemented Trimmed Composite Surface.
1342 return myshape;
1343 }
1344 }
1345 break;
1346 default:
1347 {
1348 Message_Msg msg1156("IGES_1156");
1349 const Standard_CString typeName("basis surface");
1350 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(igesSrf);
1351 msg1156.Arg(typeName);
1352 msg1156.Arg(label);
1353 SendFail(st, msg1156);
1354 return res;
1355 }
1356 }
1357 }
1358
1359 face.EmptyCopy();
1360 for (Standard_Integer i = 1; i <= st->NbBoundaries(); i++)
1361 TC.TransferBoundaryOnFace(face, st->Boundary(i), trans, uFact);
1362
1363 BRepTools::Update ( face ); //:p4
1364 //#22 rln 01.06.98 UK3.igs entity 1279
1365// ShapeFix_Face sff ( face );
1366// sff.FixMissingSeam(); //:l1 abv 10 Jan 99: USA60022 7289: fix missing seam
1367// if(sff.FixSmallAreaWire()) { //%14 pdn 24.02,99: USA60293: fix small area wires.
1368// AddFail(st, "Small area wire detected, dropped");
1369// }
1370// sff.FixOrientation();
1371// face = sff.Face();
1372 //%16 pdn 08.04.99
1373 return face;
1374}
1375
1376
1377//=======================================================================
1378//function : TransferPlane
1379//purpose :
1380//=======================================================================
1381
1382TopoDS_Shape IGESToBRep_TopoSurface::TransferPlane
1383 (const Handle(IGESGeom_Plane)& st)
1384{ // Declaration of messages//
1385 // DCE 22/12/98
1386 //Message_Msg msg1005("IGES_1005");
1387 ////////////////////////////////
1388 TopoDS_Shape res;
1389 if (st.IsNull()) {
1390 Message_Msg msg1005("IGES_1005");
1391 SendFail(st, msg1005);
1392 return res;
1393 }
1394
1395 gp_Pln pln;
1396 gp_Trsf trsf;
1397 res = TransferPlaneParts (st, pln,trsf,Standard_True);
1398// res contient (en principe ...) une Face avec eventuellement un Wire
1399// il reste a la mettre en position
1400 if (trsf.Form() != gp_Identity) {
1401 TopLoc_Location loc(trsf);
1402 res.Location(loc);
1403 }
1404 return res;
1405}
1406
1407
1408//=======================================================================
1409//function : TransferPlaneSurface
1410//purpose : this function transferred into IGESToBRep_BasicSurface
1411//=======================================================================
1412
1413
1414//=======================================================================
1415//function : TransferPerforate
1416//purpose :
1417//=======================================================================
1418
1419TopoDS_Shape IGESToBRep_TopoSurface::TransferPerforate
1420 (const Handle(IGESBasic_SingleParent)& st)
1421{ // Declaration of messages//
1422 // DCE 22/12/98
1423 //Message_Msg msg1005("IGES_1005");
1424 ////////////////////////////////
1425 TopoDS_Shape res;
1426 if (st.IsNull()) {
1427 Message_Msg msg1005("IGES_1005");
1428 SendFail(st, msg1005);
1429 return res;
1430 }
1431
1432 //char mess[100];
1433 gp_Pln pln;
1434 gp_Trsf trsf;
1435 DeclareAndCast(IGESGeom_Plane,p0,st->SingleParent());
1436 BRep_Builder B;
1437 if (p0.IsNull()) {
1438 Message_Msg msg206("XSTEP_206");
1439 SendFail(st, msg206);
1440 // SingleParent does not describe a holed face
1441 return res;
1442 }
1443 res = TransferPlaneParts (p0,pln,trsf,Standard_True);
1444//res demarre avec la face et son contour externe
1445 Standard_Integer nb = st->NbChildren();
1446 for (Standard_Integer i = 1; i <= nb; i ++) {
1447 DeclareAndCast(IGESGeom_Plane,pi,st->Child(i));
1448 if (pi.IsNull()) {
1449 Message_Msg msg1285("IGES_1285");
1450 msg1285.Arg(i);
1451 // A child is not a Plane, skipped, n0 %d
1452 SendWarning(st,msg1285);
1453 continue;
1454 }
1455 gp_Pln pli;
1456 gp_Trsf trsi;
1457 TopoDS_Shape wire = TransferPlaneParts (pi,pli,trsi,Standard_False);
1458// si ce n est pas un Wire, sauter
1459 if (wire.ShapeType() != TopAbs_WIRE) {
1460 Message_Msg msg1156("IGES_1156");
1461 const Standard_CString typeName("hole");
1462 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(pi);
1463 msg1156.Arg(typeName);
1464 msg1156.Arg(label);
1465 SendWarning(st, msg1156);
1466 // A hole could not be transferred, skipped, n0 %d
1467 continue;
1468 }
1469// coplanaires ? verifier
1470 if (!pln.Position().IsCoplanar(pli.Position(),GetEpsGeom(),GetEpsilon())) {
1471 Message_Msg msg1295("IGES_1295");
1472 msg1295.Arg(i);
1473 SendWarning(st,msg1295);
1474 // "A hole is not well coplanar to the face, n0 %d",i);
1475 }
1476// Ne pas oublier de composer la transformation locale a ce Wire
1477 if (trsi.Form() != gp_Identity) {
1478 TopLoc_Location locw(trsi);
1479 wire.Location(locw);
1480 }
1481 B.Add (res,wire);
1482 }
1483// Enfin, appliquer la trsf globale
1484 if (trsf.Form() != gp_Identity) {
1485 TopLoc_Location loc(trsf);
1486 res.Location(loc);
1487 }
1488 return res;
1489}
1490
1491
1492//=======================================================================
1493//function : TransferPlaneParts
1494//purpose :
1495//=======================================================================
1496TopoDS_Shape IGESToBRep_TopoSurface::TransferPlaneParts(const Handle(IGESGeom_Plane)& st,
1497 gp_Pln& pln,
1498 gp_Trsf& trsf,
1499 const Standard_Boolean first)
1500{ // Declaration of messages//
1501 // DCE 22/12/98
1502 //Message_Msg msg1005("IGES_1005");
1503 ////////////////////////////////
1504 TopoDS_Shape res;
1505 if (st.IsNull()) {
1506 Message_Msg msg1005("IGES_1005");
1507 SendFail(st, msg1005);
1508 return res;
1509 }
1510
1511 Standard_Real a, b, c, d;
1512 // equation de Geom : ax + by + cz + d = 0.0;
1513 // equation de IGES : ax + by + cz = d;
1514 st->Equation(a, b, c, d);
1515 pln = gp_Pln(a, b, c, -d);
1516
1517 TopoDS_Face plane;
1518 TopoDS_Wire wire;
1519 BRep_Builder B;
1520 if (first) {
1521 B.MakeFace(plane); // Just to create a empty Plane with a Tshape.
1522 Handle (Geom_Plane) geomPln = new Geom_Plane(pln);
1523 geomPln->Scale(gp_Pnt(0,0,0),GetUnitFactor());
1524// ATTENTION, ici on CALCULE la trsf, on ne l`applique pas ...
1525 //S4054: B.UpdateFace (plane, geomPln, TopLoc_Location(),
1526 //GetEpsGeom()*GetUnitFactor());
1527 B.UpdateFace (plane, geomPln, TopLoc_Location(), Precision::Confusion());
1528 //:3 by ABV 5 Nov 97: set Infinite() flag (see below for unsetting)
1529 plane.Infinite ( Standard_True ); //:3
1530 }
1531
1532// ATTENTION, ici on CALCULE la trsf, on ne l'appliquera qu'a la fin !
1533 if (st->HasTransf()) {
1534 SetEpsilon(1.E-04);
1535 if (!IGESData_ToolLocation::ConvertLocation
1536 (GetEpsilon(), st->CompoundLocation(),trsf,GetUnitFactor())) {
1537 Message_Msg msg1035("IGES_1035");
1538 SendWarning(st,msg1035); // Transformation : not a similarity
1539 }
1540 }
1541
1542 if (st->HasBoundingCurve()) {
1543 IGESToBRep_TopoCurve TC(*this);
1544 Handle(IGESData_IGESEntity) crv = st->BoundingCurve();
1545
1546 if (crv.IsNull()) {
1547 Message_Msg msg1300("IGES_1300");
1548 SendWarning(st,msg1300);
1549 //:4 by ABV 5 Nov 97: plane cannot be trimmed - let it be infinite
1550 //:4 NOTE: NB "else"
1551 //:4 return res;
1552 }
1553 else //:4
1554
1555 if (IGESToBRep::IsTopoCurve(crv)) {
1556 gp_Trsf trans;
1557 if (crv->IsKind(STANDARD_TYPE(IGESGeom_CurveOnSurface))) {
1558 DeclareAndCast(IGESGeom_CurveOnSurface, crv142, crv);
45c0ba36 1559 TopoDS_Shape myshape = TC.TransferCurveOnFace (plane, crv142, trans, TheULength, Standard_False);
1560
7fd59977 1561 //:3 by ABV 5 Nov 97: set plane to be finite
1562 if ( first ) {
1563 TopExp_Explorer ws ( plane, TopAbs_WIRE );
1564 if ( ws.More() ) plane.Infinite ( Standard_False );
1565 }
1566 }
1567 else {
1568 TopoDS_Shape shape = TC.TransferTopoCurve(crv);
1569 TopAbs_ShapeEnum shapeEnum = shape.ShapeType();
1570 switch (shapeEnum) {
1571 case TopAbs_EDGE :
1572 {
1573 TopoDS_Edge edge = TopoDS::Edge(shape);
1574 Handle(ShapeExtend_WireData) sewd = new ShapeExtend_WireData;
1575 sewd->Add(edge);
1576 wire = sewd->Wire();
1577 }
1578 break;
1579 case TopAbs_WIRE :
1580 {
1581 wire = TopoDS::Wire(shape);
1582 }
1583 break;
1584 default:
1585 {
1586 Message_Msg msg1156("IGES_1156");
1587 const Standard_CString typeName("Bounding curve");
1588 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(crv);
1589 msg1156.Arg(typeName);
1590 msg1156.Arg(label);
1591 SendWarning(st, msg1156);
1592 if (first) res = plane;
1593 else res = wire;
1594 return res;
1595 }
1596 }
1597 //S4054 CTS18953 entity 14
1598 Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire;
1599 sfw->Load (wire);
1600 sfw->FixConnected (GetMaxTol());
1601 wire = sfw->Wire();
1602 BRepLib_MakeFace MF(pln, wire, Standard_False);
1603 if (!MF.IsDone()) {
1604 // AddFail(st, "Plane Construction Error.");
1605 return res;
1606 }
1607
1608 TopoDS_Face F = MF.Face();
1609 GProp_GProps G;
1610 BRepGProp::SurfaceProperties(F,G);
1611 if (G.Mass() < 0) {
1612 if(!st->HasBoundingCurveHole())
1613 wire.Reverse();
1614 }
1615 else
1616 if( st->HasBoundingCurveHole())
1617 wire.Reverse();
1618 //:3 by ABV 5 Nov 97: set plane to be finite
1619 //:3 if (first) B.Add (plane,wire);
1620 if ( first ) {
1621 B.Add ( plane, wire );
1622 plane.Infinite ( Standard_False );
1623 }
1624 //BRepLib_MakeFace MP(pln, wire);
1625 //plane = MP.Face();
1626 }
1627 }
1628 else {
1629 Message_Msg msg1156("IGES_1156");
1630 const Standard_CString typeName("Bounding curve");
1631 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(crv);
1632 msg1156.Arg(typeName);
1633 msg1156.Arg(label);
1634 SendWarning(st, msg1156);
1635 // Plane Cannot Be Trimmed.
1636 }
1637 }
1638
1639 if (first) res = plane;
1640 else res = wire;
1641 return res;
1642}
1643
1644
1645//=======================================================================
1646//function : ParamSurface
1647//purpose :
1648//=======================================================================
1649TopoDS_Shape IGESToBRep_TopoSurface::ParamSurface(const Handle(IGESData_IGESEntity)& st,
1650 gp_Trsf2d& trans,
1651 Standard_Real& uFact)
1652{ // Declaration of messages//
1653 // DCE 22/12/98
1654 //Message_Msg msg1005("IGES_1005");
1655 ////////////////////////////////
1656
1657 TopoDS_Shape res;
1658
1659 TopoDS_Shape basisSurface = TransferTopoSurface(st);
1660 Standard_Real uscale = 1.;
1661 Standard_Real cscale = TheULength;
1662 if (basisSurface.IsNull()) {
1663 Message_Msg msg1005("IGES_1005");
1664 SendFail(st, msg1005);
1665 return res;
1666 }
1667
1668 TopAbs_ShapeEnum shapeEnum;
1669 shapeEnum = basisSurface.ShapeType();
1670 TopoDS_Face face;
1671 switch (shapeEnum) {
1672 case TopAbs_FACE :
1673 {
1674 face = TopoDS::Face(basisSurface);
1675 break;
1676 }
1677 case TopAbs_SHELL :
1678 {
1679 TopoDS_Iterator IT(basisSurface);
1680 Standard_Integer nbfaces = 0;
1681 for (; IT.More(); IT.Next()) {
1682 nbfaces++;
1683 face = TopoDS::Face(IT.Value());
1684 }
1685 //szv#4:S4163:12Mar99 optimized
1686 if (nbfaces != 1) {
1687 Message_Msg msg1156("IGES_1156");
1688 const Standard_CString typeName("basis surface");
1689 Handle(TCollection_HAsciiString) label = GetModel()->StringLabel(st);
1690 msg1156.Arg(typeName);
1691 msg1156.Arg(label);
1692 SendWarning(st, msg1156);
1693 return basisSurface;
1694 }
1695 }
1696 break;
1697 default:
1698 {
1699 //AddFail(st, "Basis Surface Transfer Error.");
1700 return res;
1701 }
1702 }
1703
1704 //S4181 pdn 19.04.99 defining shift of parametric space on base
1705 // of CAS.CADE type of surface
1706 Standard_Real paramu = 0., paramv = 0.;
1707 TopLoc_Location L;
1708 TopoDS_Edge theedge;
1709 Handle(Geom_Surface) Surf = BRep_Tool::Surface(face);
1710
1711 if (Surf->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface))) {
1712 DeclareAndCast(Geom_RectangularTrimmedSurface, rectang, Surf);
1713 Surf = rectang->BasisSurface();
1714 }
1715
1716 if ((Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface))) ||
1717 (Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface))) ||
1718 (Surf->IsKind(STANDARD_TYPE(Geom_ToroidalSurface))) ||
1719 (Surf->IsKind(STANDARD_TYPE(Geom_SphericalSurface)))) {
1720 TopExp_Explorer TE;
1721 for ( TE.Init(face,TopAbs_EDGE); TE.More(); TE.Next()){
1722 TopoDS_Edge myedge = TopoDS::Edge(TE.Current());
1723 Standard_Real First, Last;
1724 Handle(Geom2d_Curve) Curve2d = BRep_Tool::CurveOnSurface
1725 (myedge, face, First, Last);
1726 if ( Curve2d->IsKind(STANDARD_TYPE(Geom2d_Line))) {
1727 DeclareAndCast(Geom2d_Line, Line2d, Curve2d);
1728 if (Line2d->Direction().IsParallel(gp::DY2d(),Precision::Angular())){
1729 theedge = myedge;
1730 break;
1731 }
1732 }
1733 }
1734
1735 Standard_Real First, Last;
1736 Handle(Geom_Curve) Curve3d = BRep_Tool::Curve(theedge, First, Last);
1737 if (Precision::IsNegativeInfinite(First)) First = 0.;
1738
1739 if (Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface))){
1740 DeclareAndCast(Geom_CylindricalSurface, Cyl, Surf);
1741 gp_Cylinder TheCyl = Cyl->Cylinder();
1742 ElSLib::CylinderParameters(TheCyl.Position(),
1743 TheCyl.Radius(),
1744 Curve3d->Value(First), paramu , paramv );
1745 }
1746
1747 else if (Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface))){
1748 DeclareAndCast(Geom_ConicalSurface, Cone, Surf);
1749 gp_Cone TheCone = Cone->Cone();
1750 ElSLib::ConeParameters(TheCone.Position(),
1751 TheCone.RefRadius(),
1752 TheCone.SemiAngle(),
1753 Curve3d->Value(First), paramu , paramv );
1754 }
1755
1756 else if (Surf->IsKind(STANDARD_TYPE(Geom_ToroidalSurface))){
1757 DeclareAndCast(Geom_ToroidalSurface, Tore, Surf);
1758 gp_Torus TheTore = Tore->Torus();
1759 ElSLib::TorusParameters(TheTore.Position(),
1760 TheTore.MajorRadius(),
1761 TheTore.MinorRadius(),
1762 Curve3d->Value(First), paramu , paramv );
1763 }
1764 else if (Surf->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) &&
1765 //: abv 18.06.02: loopback on s1.brep BRep mode, face 32 - the spherical surface (192)
1766 //: in IGES is from [-90,90] by V, i.e. similar to OCC, thus only scaling is enough
1767 ! st->IsKind (STANDARD_TYPE(IGESSolid_SphericalSurface)) ) {
1768 DeclareAndCast(Geom_SphericalSurface, Sphere, Surf);
1769 gp_Sphere TheSphere = Sphere->Sphere();
1770 ElSLib::SphereParameters(TheSphere.Position(),
1771 TheSphere.Radius(),
1772 Curve3d->Value(First), paramu , paramv );
1773 }
1774
1775 //#88 rln 06.04.99 CTS60168, BLEND.IGS entity 68
1776 //Generatrix is Transformed Circular Arc. When creating CAS.CADE surface parameterization
1777 //has changed.
1778 if (st->IsKind (STANDARD_TYPE (IGESGeom_SurfaceOfRevolution))) {
1779 DeclareAndCast (IGESGeom_SurfaceOfRevolution, revol, st);
1780 Handle(IGESData_IGESEntity) generatrix = revol->Generatrix();
1781 if (generatrix->IsKind (STANDARD_TYPE (IGESGeom_CircularArc))) {
1782 DeclareAndCast (IGESGeom_CircularArc, circ, generatrix);
1783 gp_Pnt2d startpoint = circ->StartPoint();
1784 paramv -= ElCLib::Parameter (gp_Circ2d (gp_Ax2d (circ->Center(), gp_Dir2d(1,0)), circ->Radius()), startpoint);
1785 if (Surf->IsKind (STANDARD_TYPE(Geom_SphericalSurface)))
c6541a0c 1786 paramv += ShapeAnalysis::AdjustToPeriod(paramv, - M_PI, M_PI);
7fd59977 1787 else if (Surf->IsKind (STANDARD_TYPE(Geom_ToroidalSurface)))
c6541a0c 1788 paramv += ShapeAnalysis::AdjustToPeriod(paramv, 0, M_PI * 2);
7fd59977 1789 }
1790 }
1791 else if (st->IsKind (STANDARD_TYPE (IGESGeom_TabulatedCylinder))) {
1792 DeclareAndCast (IGESGeom_TabulatedCylinder, cylinder, st);
1793 Handle(IGESData_IGESEntity) directrix = cylinder->Directrix();
1794 if (directrix->IsKind (STANDARD_TYPE (IGESGeom_CircularArc))) {
1795 DeclareAndCast (IGESGeom_CircularArc, circ, directrix);
1796 gp_Pnt2d startpoint = circ->StartPoint();
1797 paramu -= ElCLib::Parameter (gp_Circ2d (gp_Ax2d (circ->Center(), gp_Dir2d(1,0)), circ->Radius()), startpoint);
c6541a0c 1798 paramu += ShapeAnalysis::AdjustToPeriod(paramu, 0, M_PI * 2);
7fd59977 1799 }
1800 }
1801
1802 }
1803
1804 if ( Abs(paramu) <= Precision::Confusion()) paramu = 0.;
1805 if ( Abs(paramv) <= Precision::Confusion()) paramv = 0.;
1806
1807 //S4181 pdn 16.04.99 computation of transformation depending on
1808 //IGES Type of surface
1809 Handle(IGESData_IGESEntity) isrf = st;
1810 if (isrf->IsKind(STANDARD_TYPE(IGESGeom_OffsetSurface))){
1811 DeclareAndCast(IGESGeom_OffsetSurface, offsurf, isrf);
1812 isrf = offsurf->Surface();
1813 }
1814 if (isrf->IsKind(STANDARD_TYPE(IGESGeom_SurfaceOfRevolution))) {
1815 DeclareAndCast(IGESGeom_SurfaceOfRevolution, st120, isrf);
1816 //S4181 pdn 19.04.99 defining transformation matrix
1817 gp_Trsf2d tmp;
c6541a0c 1818 tmp.SetTranslation(gp_Vec2d (0, -2 * M_PI));
7fd59977 1819 trans.PreMultiply(tmp);
1820 tmp.SetMirror(gp::OX2d());
1821 trans.PreMultiply(tmp);
1822 tmp.SetMirror(gp_Ax2d (gp::Origin2d(), gp_Dir2d (1.,1.)));
1823 trans.PreMultiply(tmp);
1824 uscale = 1./cscale;
1825 //#30 rln 19.10.98
1826 //CAS.CADE SA = 2*PI - IGES TA
1827 //paramu = st120->StartAngle();
c6541a0c 1828 paramu = -(2 * M_PI - st120->EndAngle());
7fd59977 1829 }
1830 else
1831 paramu = 0.;
1832
1833 if (isrf->IsKind(STANDARD_TYPE(IGESGeom_RuledSurface))) {
1834 uscale = 1./cscale;
1835 }
1836
1837 // corrected skl 13.11.2001 for BUC61054
1838 if (isrf->IsKind(STANDARD_TYPE(IGESGeom_TabulatedCylinder))) {
1839 Handle(IGESGeom_TabulatedCylinder) igtc = Handle(IGESGeom_TabulatedCylinder)::DownCast(isrf);
1840 Handle(IGESData_IGESEntity) idie = igtc->Directrix();
1841 Standard_Real uln=1;
a9da59ff 1842 Standard_Real Umin,Umax,Vmin,Vmax;
1843 //scaling parameterization from [0,1]
1844 Surf->Bounds(Umin,Umax,Vmin,Vmax);
1845 uln = Abs(Umax-Umin);
1846 //computing shift of pcurves
7fd59977 1847 uscale = uln/cscale;
a9da59ff 1848 paramu = Umin/uln;
7fd59977 1849 }
1850
1851 if (isrf->IsKind(STANDARD_TYPE(IGESSolid_CylindricalSurface))||
1852 isrf->IsKind(STANDARD_TYPE(IGESSolid_ConicalSurface))) {
c6541a0c 1853 uscale = M_PI/180.;
7fd59977 1854 }
1855
1856 if (isrf->IsKind(STANDARD_TYPE(IGESSolid_SphericalSurface))) {
c6541a0c 1857 cscale = M_PI/180.;
7fd59977 1858 uscale = 1.;
1859 }
1860
1861 if (isrf->IsKind(STANDARD_TYPE(IGESSolid_ToroidalSurface))) {
1862 gp_Trsf2d tmp;
1863 tmp.SetTranslation(gp_Vec2d (0, -360.)); // in IGES terms
1864 trans.PreMultiply(tmp);
1865 tmp.SetMirror(gp::OX2d());
1866 trans.PreMultiply(tmp);
1867 tmp.SetMirror(gp_Ax2d (gp::Origin2d(), gp_Dir2d (1.,1.)));
1868 trans.PreMultiply(tmp);
1869 if(paramv > 0)
c6541a0c
D		 1870 paramv = paramv*180./M_PI;
1871 cscale = M_PI/180.;
7fd59977 1872 uscale = 1.;
1873 }
1874
1875 gp_Trsf2d tmp;
1876 tmp.SetTranslation(gp_Pnt2d(0.,0.), gp_Pnt2d(paramu,paramv));
1877 trans.PreMultiply(tmp);
1878
1879 tmp.SetScale(gp_Pnt2d(0,0),cscale);
1880 trans.PreMultiply(tmp);
1881 uFact = uscale;
1882 return face;
1883}
1884
Open CASCADE Technology repository