ffb102ab4bcd88abad332080ff51a909d3e3ed92
[occt.git] / src / BRepToIGES / BRepToIGES_BRWire.cxx
1 // Created on: 1995-01-27
2 // Created by: Marie Jose MARTZ
3 // Copyright (c) 1995-1999 Matra Datavision
4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
5 //
6 // The content of this file is subject to the Open CASCADE Technology Public
7 // License Version 6.5 (the "License"). You may not use the content of this file
8 // except in compliance with the License. Please obtain a copy of the License
9 // at http://www.opencascade.org and read it completely before using this file.
10 //
11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
13 //
14 // The Original Code and all software distributed under the License is
15 // distributed on an "AS IS" basis, without warranty of any kind, and the
16 // Initial Developer hereby disclaims all such warranties, including without
17 // limitation, any warranties of merchantability, fitness for a particular
18 // purpose or non-infringement. Please see the License for the specific terms
19 // and conditions governing the rights and limitations under the License.
20
21 //:q4 abv 16.03.99: PRO17828 face 555: transform pcurves on SurfaceOfRevolution
22 //szv#4 S4163
23 //S4181 pdn implementing of writing IGES elementary surfaces.
24
25 #include <BRepToIGES_BRWire.ixx>
26
27 #include <BRep_Tool.hxx>
28 #include <BRepTools.hxx>
29 #include <BRepTools_WireExplorer.hxx>
30
31 #include <gp.hxx>
32 #include <gp_Ax2d.hxx>
33 #include <gp_Dir2d.hxx>
34 #include <gp_Pnt2d.hxx>
35 #include <gp_Trsf.hxx>
36 #include <gp_Trsf2d.hxx>
37
38 #include <Geom_CartesianPoint.hxx>
39 #include <Geom_ConicalSurface.hxx>
40 #include <Geom_Curve.hxx>
41 #include <Geom_CylindricalSurface.hxx>
42 #include <Geom_Plane.hxx>
43 #include <Geom_RectangularTrimmedSurface.hxx>
44 #include <Geom_SphericalSurface.hxx>
45 #include <Geom_SurfaceOfLinearExtrusion.hxx>
46 #include <Geom_SurfaceOfRevolution.hxx>
47 #include <Geom_ToroidalSurface.hxx>
48 #include <Geom_TrimmedCurve.hxx>
49
50 #include <Geom2d_Curve.hxx>
51 #include <Geom2d_Line.hxx>
52
53 #include <GeomToIGES_GeomPoint.hxx>
54 #include <GeomToIGES_GeomCurve.hxx>
55
56 #include <Geom2dToIGES_Geom2dCurve.hxx>
57
58 #include <IGESData_HArray1OfIGESEntity.hxx>
59 #include <IGESData_IGESEntity.hxx>
60 #include <IGESData_IGESModel.hxx>
61
62 #include <IGESGeom_CompositeCurve.hxx>
63 #include <IGESGeom_Point.hxx>
64
65 #include <Interface_Macros.hxx>
66
67 #include <Precision.hxx>
68
69 #include <TColStd_HSequenceOfTransient.hxx>
70
71 #include <TopoDS.hxx>
72 #include <TopoDS_Vertex.hxx>
73 #include <TopoDS_Edge.hxx>
74 #include <TopoDS_Wire.hxx>
75 #include <TopoDS_Shape.hxx>
76
77 #include <TopAbs_ShapeEnum.hxx>
78 #include <TopAbs_Orientation.hxx>
79 #include <TopExp.hxx>
80 #include <TopExp_Explorer.hxx>
81 #include <Interface_Static.hxx>
82 #include <ShapeBuild_Edge.hxx>
83 #include <Geom_OffsetSurface.hxx>
84 #include <ShapeExtend_WireData.hxx>
85 #include <ShapeFix_Wire.hxx>
86
87
88 //=============================================================================
89 // BRepToIGES_BRWire
90 //=============================================================================
91
92 BRepToIGES_BRWire::BRepToIGES_BRWire()
93 {
94 }
95
96
97 //=============================================================================
98 // BRepToIGES_BRWire
99 //=============================================================================
100
101 BRepToIGES_BRWire::BRepToIGES_BRWire
102 (const BRepToIGES_BREntity& BR)
103 : BRepToIGES_BREntity(BR)
104 {
105 }
106
107
108 //=============================================================================
109 // TransferWire
110 //=============================================================================
111
112 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferWire
113 (const TopoDS_Shape& start)
114 {
115   Handle(IGESData_IGESEntity) res;
116
117   if (start.IsNull())  return  res;
118
119   if (start.ShapeType() == TopAbs_VERTEX) {
120     TopoDS_Vertex V =  TopoDS::Vertex(start);
121     res = TransferVertex(V);
122   }  
123   else if (start.ShapeType() == TopAbs_EDGE) {
124     TopoDS_Edge E =  TopoDS::Edge(start);
125     res = TransferEdge(E, Standard_False);
126   }  
127   else if (start.ShapeType() == TopAbs_WIRE) {
128     TopoDS_Wire W =  TopoDS::Wire(start);
129     res = TransferWire(W);
130   }  
131   else {
132     // error message
133   }  
134   return res;
135 }
136
137
138 //=============================================================================
139 // TransferVertex
140 // 
141 //=============================================================================
142
143 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferVertex
144 (const TopoDS_Vertex& myvertex)
145 {
146   Handle(IGESData_IGESEntity) res;
147   if ( myvertex.IsNull()) return res;
148
149   Handle(Geom_CartesianPoint) Point;
150   Point = new Geom_CartesianPoint(BRep_Tool::Pnt(myvertex));
151   Handle(IGESData_IGESEntity) IVertex;
152   if (!Point.IsNull()) {
153     GeomToIGES_GeomPoint GP;
154     GP.SetModel(GetModel());
155     IVertex = GP.TransferPoint(Point);
156   }
157
158   if (!IVertex.IsNull()) res = IVertex;
159   return res;
160 }
161
162
163 //=============================================================================
164 // TransferVertex
165 // 
166 //=============================================================================
167
168 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferVertex
169 (const TopoDS_Vertex& myvertex,
170  const TopoDS_Edge&   myedge,
171  Standard_Real&       parameter)
172 {
173   Handle(IGESData_IGESEntity) res;
174   if ( myvertex.IsNull()) return res;
175
176   Handle(IGESData_IGESEntity) IVertex = TransferVertex(myvertex);
177
178   // returns the parameter of myvertex on myedge
179   parameter = BRep_Tool::Parameter(myvertex,myedge);
180
181   if (!IVertex.IsNull()) res = IVertex;
182   return res;
183 }
184
185
186 //=============================================================================
187 // TransferVertex
188 // 
189 //=============================================================================
190
191 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferVertex
192 (const TopoDS_Vertex& myvertex,
193  const TopoDS_Edge&   myedge,
194  const TopoDS_Face&   myface,
195  Standard_Real&       parameter)
196 {
197   Handle(IGESData_IGESEntity) res;
198   if ( myvertex.IsNull()) return res;
199
200   Handle(IGESData_IGESEntity) IVertex = TransferVertex(myvertex);
201
202   // returns the parameter of myvertex on the pcurve of myedge on myface
203   parameter = BRep_Tool::Parameter(myvertex, myedge, myface);
204
205   if (!IVertex.IsNull()) res = IVertex;
206   return res;
207 }
208
209
210 //=============================================================================
211 // TransferVertex
212 // 
213 //=============================================================================
214
215 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferVertex
216 (const TopoDS_Vertex&          myvertex,
217  const TopoDS_Edge&            myedge,
218  const Handle(Geom_Surface)&   mysurface,
219  const TopLoc_Location&        myloc,
220  Standard_Real&                parameter)
221 {
222   Handle(IGESData_IGESEntity) res;
223   if ( myvertex.IsNull()) return res;
224
225   Handle(IGESData_IGESEntity) IVertex = TransferVertex(myvertex);
226
227   // returns the parameter of myvertex on the pcurve of myedge on mysurface
228   parameter = BRep_Tool::Parameter(myvertex, myedge, mysurface, myloc);
229
230   if (!IVertex.IsNull()) res = IVertex;
231   return res;
232 }
233
234
235 //=============================================================================
236 // TransferVertex
237 // 
238 //=============================================================================
239
240 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferVertex
241 (const TopoDS_Vertex& myvertex,
242  const TopoDS_Face&   myface,
243        gp_Pnt2d&      mypoint)
244 {
245   Handle(IGESData_IGESEntity) res;
246   if ( myvertex.IsNull()) return res;
247
248   Handle(IGESData_IGESEntity) IVertex = TransferVertex(myvertex);
249
250   // returns the parameter of myvertex on myface
251   mypoint = BRep_Tool::Parameters(myvertex, myface);
252
253   if (!IVertex.IsNull()) res = IVertex;
254   return res;
255 }
256
257
258 //=============================================================================
259 // TransferEdge
260 //=============================================================================
261
262 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferEdge
263 (const TopoDS_Edge&     myedge,
264  const Standard_Boolean isBRepMode)
265 {
266   Handle(IGESData_IGESEntity) res;
267   if ( myedge.IsNull()) return res;
268
269   // returns the 3d curve of the edge and the parameter range
270   TopLoc_Location L;
271   Standard_Real First, Last, U1, U2;
272   Handle(IGESData_IGESEntity) ICurve;
273   Handle(Geom_Curve) Curve3d = BRep_Tool::Curve(myedge, L, First, Last);
274
275   //#28 rln 19.10.98 UKI60155, etc.
276   //Only Bezier will be converted into B-Spline, not Conic. This conertation
277   //will be done only inside GeomToIGES package for simplifying the code.
278
279   //#29 rln 19.10.98
280   //Unnecessary duplication of curves is removed.
281
282   if (!Curve3d.IsNull()) {
283     gp_Trsf Tr = L.Transformation();
284     if (Tr.Form() != gp_Identity)
285       Curve3d = Handle(Geom_Curve)::DownCast(Curve3d->Transformed (Tr));
286     else
287       Curve3d = Handle(Geom_Curve)::DownCast(Curve3d->Copy()); 
288
289     if (myedge.Orientation() == TopAbs_REVERSED && !isBRepMode) {
290       U1 = Curve3d->ReversedParameter(Last);
291       U2 = Curve3d->ReversedParameter(First);
292       Curve3d->Reverse();
293     }
294     else {
295       U1 = First;
296       U2 = Last;
297     }
298
299     GeomToIGES_GeomCurve GC;
300     GC.SetModel(GetModel());
301     ICurve = GC.TransferCurve(Curve3d, U1, U2);
302   }
303
304   //#31 rln 19.10.98
305   //Vertices are not translated into IGES anymore since they are not put into
306   //the model
307
308   if (!ICurve.IsNull()) res = ICurve;
309
310   SetShapeResult ( myedge, res );
311
312   return res;
313 }
314
315
316 //=============================================================================
317 // TransferEdge
318 //=============================================================================
319
320 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferEdge (const TopoDS_Edge& myedge,
321                                                               const TopoDS_Face& myface,
322                                                               const Standard_Real length,
323                                                               const Standard_Boolean isBRepMode)
324 {
325   Handle(IGESData_IGESEntity) res;
326   if ( myedge.IsNull() || GetPCurveMode() ==0 ||
327        ( ! isBRepMode && BRep_Tool::Degenerated ( myedge ) ) ) return res;
328
329   //S4181 pdn 23.04.99 adjusting length factor according to analytic mode.
330   Standard_Real myLen = length;
331   Standard_Boolean analyticMode = ( GetConvertSurfaceMode() ==0 && isBRepMode );
332   
333   // returns the 2d curve associated to myedge in the parametric space of myface
334   Standard_Real First, Last;
335
336   Handle(Geom2d_Curve) Curve2d = BRep_Tool::CurveOnSurface
337     (myedge, myface, First, Last);
338   Handle(IGESData_IGESEntity) ICurve2d;
339   //#29 rln 19.10.98
340
341   if (!Curve2d.IsNull()) {
342     // For "revolution" and "LinearExtrusion" surfaces, it is necessary
343     // to apply a translation of 2D curves to agree on the 
344     // origin (U,V) between IGES and BRep (for Cylindrical,
345     // Conical and SurfaceOfLinearExtrusion)
346     // It is necessary to invert (u,v) surfaces of revolution.
347     
348     TopLoc_Location L;
349     Handle(Geom_Surface) st = BRep_Tool::Surface(myface, L);
350     Standard_Real Ufirst, Ulast, Vfirst, Vlast;
351     BRepTools::UVBounds(myface, Ufirst, Ulast, Vfirst, Vlast);
352     Handle(Geom_Surface) Surf;
353
354     if (st->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface))) { 
355       DeclareAndCast(Geom_RectangularTrimmedSurface, rectang, st);
356       Surf = rectang->BasisSurface();
357     }
358     else 
359       Surf = st;
360
361     //:abv 19.06.02: writing (loopback) on file offseted_sphere.rle in BRep mode
362     if (st->IsKind(STANDARD_TYPE(Geom_OffsetSurface))) { 
363       DeclareAndCast(Geom_OffsetSurface, offset, Surf);
364       Surf = offset->BasisSurface();
365     }
366
367     //S4181 pdn 20.04.99 transformation of pcurves id defined by type of surface 
368     // and analytic mode.
369     // skl 18.07.2005 for OCC9490 : in following if() commented
370     // condition for SurfaceOfLinearExtrusion
371     Standard_Boolean needShift = (!analyticMode&&
372                                   ((Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface))) ||
373                                    (Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface)))/* ||
374                                    (Surf->IsKind(STANDARD_TYPE(Geom_SurfaceOfLinearExtrusion)))*/));
375     //:q4 abv 16 Mar 99: PRO17828 face 555: shift pcurves on SurfaceOfRevolution
376     if (Surf->IsKind(STANDARD_TYPE(Geom_SurfaceOfRevolution))) { 
377       Handle(Geom_SurfaceOfRevolution) rev = 
378         Handle(Geom_SurfaceOfRevolution)::DownCast ( Surf );
379       Handle(Geom_Curve) curve = rev->BasisCurve();
380       // skl 31.05.2004 for OCC6004
381       if(curve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) {
382         Handle(Geom_TrimmedCurve) tc = Handle(Geom_TrimmedCurve)::DownCast(curve);
383         curve = tc->BasisCurve();
384       }
385       if ( curve->IsKind(STANDARD_TYPE(Geom_Line)) ) needShift = Standard_True;
386     }
387     if ( needShift ) {
388       gp_Trsf2d TR;
389       TR.SetTranslation(gp_Pnt2d(0.,0.),gp_Pnt2d(0.,-Vfirst)); 
390       Curve2d = Handle(Geom2d_Curve)::DownCast(Curve2d->Transformed(TR));
391     }
392     else Curve2d = Handle(Geom2d_Curve)::DownCast(Curve2d->Copy());
393
394     if (!analyticMode&&((Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface)))  ||
395                         (Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface)))      ||
396                         (Surf->IsKind(STANDARD_TYPE(Geom_SphericalSurface))))) {
397       //#30 rln 19.10.98 transformation of pcurves for IGES Surface of Revolution
398       Curve2d->Mirror (gp_Ax2d (gp::Origin2d(), gp_Dir2d (1.,1.)));
399       Curve2d->Mirror (gp::OX2d());
400       Curve2d->Translate (gp_Vec2d (0, 2 * M_PI));
401     }
402     
403     if(Surf->IsKind(STANDARD_TYPE(Geom_SurfaceOfRevolution))||
404        (Surf->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)))){
405       Curve2d->Mirror (gp_Ax2d (gp::Origin2d(), gp_Dir2d (1.,1.)));
406       Curve2d->Mirror (gp::OX2d());
407       Curve2d->Translate (gp_Vec2d (0, 2 * M_PI));
408     }
409     
410     if (analyticMode&&(Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface)) ||
411                        Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface))))
412       myLen = M_PI/180.;
413     
414     if (analyticMode&&(Surf->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
415                        Surf->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)))) {
416       gp_Trsf2d trans;
417       trans.SetScale(gp_Pnt2d(0,0),180./M_PI);
418       Curve2d->Transform(trans);
419       First = Curve2d->TransformedParameter(First,trans);
420       Last  = Curve2d->TransformedParameter(Last, trans);
421     }
422     
423     if (analyticMode&&(Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface)))) {
424       Handle(Geom_ConicalSurface) con = Handle(Geom_ConicalSurface)::DownCast ( Surf );
425       if(con->SemiAngle() < 0) {
426         Standard_Real vApex = 2 * con->RefRadius() / Sin (con->SemiAngle());
427         Curve2d->Translate (gp_Vec2d (0, vApex));
428       }
429     }
430     //S4181 transfer functionality
431     gp_Trsf2d trans;
432     Standard_Real uFact = 1.;
433     if(isBRepMode && Surf->IsKind(STANDARD_TYPE(Geom_Plane))) {
434       trans.SetScale(gp_Pnt2d(0,0),1./GetUnit());
435     }
436     if(Surf->IsKind(STANDARD_TYPE(Geom_SurfaceOfLinearExtrusion))) {
437       Standard_Real aDiv = myLen;
438       if(aDiv < gp::Resolution())
439         aDiv = 1.;
440       //emv: changed for bug OCC22126 17.12.2010
441       trans.SetScale(gp_Pnt2d(0,0), 1. / (Vlast - Vfirst));
442       //uFact = myLen; 
443       
444       //added by skl 18.07.2005 for OCC9490
445 //      trans.SetScale(gp_Pnt2d(0,0),1./Vlast);
446
447       Standard_Real du = 1.;
448       Handle(Geom_SurfaceOfLinearExtrusion) LE =
449         Handle(Geom_SurfaceOfLinearExtrusion)::DownCast(Surf);
450       if(LE->BasisCurve()->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) {
451         Handle(Geom_TrimmedCurve) tc =
452           Handle(Geom_TrimmedCurve)::DownCast(LE->BasisCurve());
453         if(tc->BasisCurve()->IsKind(STANDARD_TYPE(Geom_Line))) {
454           Standard_Real us1,us2,vs1,vs2;
455           Surf->Bounds(us1,us2,vs1,vs2);
456           du = us2-us1;
457         }
458       }
459       //emv: changed for bug OCC22126 17.12.2010
460       uFact = (Vlast - Vfirst)/du;
461       //uFact = aDiv/du;
462     }
463     if (Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface)) ||
464         Surf->IsKind(STANDARD_TYPE(Geom_ConicalSurface)) ||
465         Surf->IsKind(STANDARD_TYPE(Geom_SurfaceOfRevolution))) { //:q4
466       uFact = 1./myLen; 
467     }
468     
469     ShapeBuild_Edge sbe;
470     Curve2d = sbe.TransformPCurve(Curve2d,trans,uFact,First,Last);
471 //      (Curve2d, Surf, First, Last, myLen, isBRepMode);
472     // if the edge is REVERSED, it is necessary to "REVERSE" the curve 2d.
473
474     // added by skl 18.07.2005 for OCC9490
475     if(Surf->IsKind(STANDARD_TYPE(Geom_SurfaceOfLinearExtrusion))) {
476       //emv: changed for bug OCC22126 17.12.2010
477       gp_Trsf2d trans1;
478       trans1.SetTranslation(gp_Vec2d(0.,-Vfirst/(Vlast-Vfirst)));
479       Curve2d = sbe.TransformPCurve(Curve2d,trans1,1.,First,Last);
480     }
481
482     if (myedge.Orientation() == TopAbs_REVERSED) {
483       Standard_Real tmpFirst = Curve2d->ReversedParameter(Last),
484                     tmpLast  = Curve2d->ReversedParameter(First);
485       Curve2d->Reverse();
486       First = tmpFirst;
487       Last  = tmpLast;
488     }
489     Geom2dToIGES_Geom2dCurve GC;
490     GC.SetModel(GetModel());
491     ICurve2d = GC.Transfer2dCurve(Curve2d, First, Last);
492   }
493
494   //#31 rln 19.10.98
495   //Vertices are not translated into IGES anymore since they are not put into
496   //the model
497
498   if (!ICurve2d.IsNull()) res = ICurve2d;
499
500   SetShapeResult ( myedge, res );
501
502   return res;
503 }
504
505
506 //=============================================================================
507 // TransferWire
508 //=============================================================================
509
510 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferWire
511 (const TopoDS_Wire& mywire)
512 {
513   Handle(IGESData_IGESEntity) res;
514   if ( mywire.IsNull()) return res;
515
516   //A composite curve is defined as an ordered list of entities
517   //consisting of a point, connect point and parametrised curve
518   //entities (excluding the CompositeCurve entity).
519
520   Handle(IGESData_IGESEntity) ent ;
521   Handle(TColStd_HSequenceOfTransient) Seq = new TColStd_HSequenceOfTransient();
522
523   TopExp_Explorer TE(mywire, TopAbs_VERTEX);
524   if ( TE.More()) {
525     BRepTools_WireExplorer WE;
526     for ( WE.Init(mywire); WE.More(); WE.Next()) {
527       TopoDS_Edge E = WE.Current();
528       if (E.IsNull()) {
529         AddWarning(mywire, "an Edge is a null entity");
530       }
531       else {
532         ent = TransferEdge(E, Standard_False);
533         if (!ent.IsNull()) Seq->Append(ent);
534       }
535     }
536   }
537   else 
538     AddWarning(mywire, " no Vertex associated to the Wire");
539   
540
541   Standard_Integer nbedges = Seq->Length();
542   Handle(IGESData_HArray1OfIGESEntity) Tab;
543   if ( nbedges == 1 ) {
544     res = ent;
545   }
546   else if ( nbedges >= 2) {
547     Tab =  new IGESData_HArray1OfIGESEntity(1,nbedges);
548     for (Standard_Integer itab = 1; itab <= nbedges; itab++) {
549       Handle(IGESData_IGESEntity) item = GetCasted(IGESData_IGESEntity, Seq->Value(itab));
550       Tab->SetValue(itab,item);
551     }
552     Handle(IGESGeom_CompositeCurve) Comp = new IGESGeom_CompositeCurve;
553     Comp->Init(Tab);
554     res = Comp;
555   }
556
557   SetShapeResult ( mywire, res );
558
559   return res;
560 }
561
562
563 //=============================================================================
564 // TransferWire
565 //=============================================================================
566
567 Handle(IGESData_IGESEntity) BRepToIGES_BRWire ::TransferWire
568 (const TopoDS_Wire &          mywire,
569  const TopoDS_Face &          myface,
570  Handle(IGESData_IGESEntity)& mycurve2d,
571  const Standard_Real          length)
572 {
573   Handle(IGESData_IGESEntity) res;
574   if ( mywire.IsNull()) return res;
575
576   Handle(IGESData_IGESEntity) ent3d ;
577   Handle(IGESData_IGESEntity) ent2d ;
578   Handle(TColStd_HSequenceOfTransient) Seq3d = new TColStd_HSequenceOfTransient();
579   Handle(TColStd_HSequenceOfTransient) Seq2d = new TColStd_HSequenceOfTransient();
580
581
582   // create a 3d CompositeCurve and a 2d CompositeCurve
583   TopExp_Explorer TE(mywire, TopAbs_VERTEX);
584   if ( TE.More()) {
585     // PTV OCC908 workaround for KAS:dev version
586     /*
587     BRepTools_WireExplorer WE;
588     for ( WE.Init(mywire,myface); WE.More(); WE.Next()) { 
589       TopoDS_Edge E = WE.Current();
590       if (E.IsNull()) {
591         AddWarning(mywire, "an Edge is a null entity");
592       }
593       else {
594         ent3d = TransferEdge(E, Standard_False);
595         if (!ent3d.IsNull()) Seq3d->Append(ent3d);
596         ent2d = TransferEdge(E, myface, length, Standard_False);
597         if (!ent2d.IsNull()) Seq2d->Append(ent2d);
598       }
599     }
600     */
601     Handle(ShapeFix_Wire) aSFW = 
602       new ShapeFix_Wire( mywire, myface, Precision::Confusion() );
603     aSFW->FixReorder();
604     Handle(ShapeExtend_WireData) aSEWD = aSFW->WireData();
605     Standard_Integer nbE = aSEWD->NbEdges();
606     for (Standard_Integer windex = 1; windex <= nbE; windex++) {
607       TopoDS_Edge E = aSEWD->Edge( windex );
608       if (E.IsNull()) {
609         AddWarning(mywire, "an Edge is a null entity");
610       }
611       else {
612         ent3d = TransferEdge(E, Standard_False);
613         if (!ent3d.IsNull()) Seq3d->Append(ent3d);
614         ent2d = TransferEdge(E, myface, length, Standard_False);
615         if (!ent2d.IsNull()) Seq2d->Append(ent2d);
616       }
617     }
618     // OCC908 end of workaround
619   }
620   else 
621     AddWarning(mywire, " no Vertex associated to the Wire");
622   
623   // Composite Curve 3D
624   Standard_Integer nb3d = Seq3d->Length();
625   Handle(IGESData_HArray1OfIGESEntity) Tab3d;
626   if ( nb3d == 1 ) {
627     res = ent3d;
628   }
629   else if (nb3d >= 2) {
630     Tab3d =  new IGESData_HArray1OfIGESEntity(1,nb3d);
631     //Standard_Integer tabval = nb3d; //szv#4:S4163:12Mar99 unused
632     for (Standard_Integer itab = 1; itab <= nb3d; itab++) {
633       Handle(IGESData_IGESEntity) item = 
634         GetCasted(IGESData_IGESEntity, Seq3d->Value(itab));
635       Tab3d->SetValue(itab,item);
636     }
637     Handle(IGESGeom_CompositeCurve) Comp = new IGESGeom_CompositeCurve;
638     Comp->Init(Tab3d);
639     res = Comp;
640   }
641
642   // Composite Curve 2D
643   Standard_Integer nb2d = Seq2d->Length();
644   Handle(IGESData_HArray1OfIGESEntity) Tab2d;
645   if ( nb2d == 1 ) {
646     mycurve2d = ent2d;
647   }
648   else if (nb2d >= 2) {
649     Tab2d =  new IGESData_HArray1OfIGESEntity(1,nb2d);
650     //Standard_Integer tabval = nb2d; //szv#4:S4163:12Mar99 unused
651     for ( Standard_Integer itab = 1; itab <= nb2d; itab++) {
652       Handle(IGESData_IGESEntity) item = 
653         GetCasted(IGESData_IGESEntity, Seq2d->Value(itab));
654       Tab2d->SetValue(itab,item);
655     }
656     Handle(IGESGeom_CompositeCurve) Comp = new IGESGeom_CompositeCurve;
657     Comp->Init(Tab2d);
658     mycurve2d = Comp;
659   }
660   
661   SetShapeResult ( mywire, res );
662
663   return res;
664 }