0026937: Eliminate NO_CXX_EXCEPTION macro support
[occt.git] / src / BRepFill / BRepFill_Draft.cxx
1 // Created on: 1998-06-08
2 // Created by: Stephanie HUMEAU
3 // Copyright (c) 1998-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17
18 #include <Adaptor3d_HCurve.hxx>
19 #include <Adaptor3d_Surface.hxx>
20 #include <Bnd_Box.hxx>
21 #include <BndLib_Add3dCurve.hxx>
22 #include <BndLib_AddSurface.hxx>
23 #include <BRep_Builder.hxx>
24 #include <BRep_Tool.hxx>
25 #include <BRepAdaptor_Curve.hxx>
26 #include <BRepAdaptor_Surface.hxx>
27 #include <BRepAlgo_DSAccess.hxx>
28 #include <BRepBuilderAPI_Sewing.hxx>
29 #include <BRepClass3d_SolidClassifier.hxx>
30 #include <BRepExtrema_DistShapeShape.hxx>
31 #include <BRepFill_DataMapOfShapeHArray2OfShape.hxx>
32 #include <BRepFill_Draft.hxx>
33 #include <BRepFill_DraftLaw.hxx>
34 #include <BRepFill_SectionLaw.hxx>
35 #include <BRepFill_ShapeLaw.hxx>
36 #include <BRepFill_Sweep.hxx>
37 #include <BRepLib_FindSurface.hxx>
38 #include <BRepLib_MakeEdge.hxx>
39 #include <BRepLib_MakeFace.hxx>
40 #include <BRepLib_MakeWire.hxx>
41 #include <BRepTools.hxx>
42 #include <Geom_Geometry.hxx>
43 #include <Geom_Line.hxx>
44 #include <Geom_Plane.hxx>
45 #include <Geom_RectangularTrimmedSurface.hxx>
46 #include <Geom_Surface.hxx>
47 #include <Geom_TrimmedCurve.hxx>
48 #include <GeomAdaptor_HSurface.hxx>
49 #include <GeomAdaptor_Surface.hxx>
50 #include <GeomFill_LocationDraft.hxx>
51 #include <GeomLProp_SLProps.hxx>
52 #include <gp_Ax3.hxx>
53 #include <gp_Dir.hxx>
54 #include <gp_Lin.hxx>
55 #include <gp_Mat.hxx>
56 #include <gp_Pnt.hxx>
57 #include <gp_Trsf.hxx>
58 #include <Precision.hxx>
59 #include <Standard_NoSuchObject.hxx>
60 #include <StdFail_NotDone.hxx>
61 #include <TColgp_Array1OfPnt.hxx>
62 #include <TColStd_Array1OfReal.hxx>
63 #include <TopAbs.hxx>
64 #include <TopExp.hxx>
65 #include <TopExp_Explorer.hxx>
66 #include <TopoDS.hxx>
67 #include <TopoDS_Edge.hxx>
68 #include <TopoDS_Iterator.hxx>
69 #include <TopoDS_Shape.hxx>
70 #include <TopoDS_Shell.hxx>
71 #include <TopoDS_Solid.hxx>
72 #include <TopoDS_Wire.hxx>
73 #include <TopTools_ListIteratorOfListOfShape.hxx>
74 #include <TopTools_ListOfShape.hxx>
75
76 #ifdef DRAW
77 #include <Geom_Circle.hxx>
78 #include <gp.hxx>
79 #include <DBRep.hxx>
80 #include <DrawTrSurf.hxx>
81 static Standard_Boolean Affich = 0;
82 #endif
83
84 //=======================================================================
85 //function : Trsf
86 //purpose  : 
87 //======================================================================
88 static void ComputeTrsf(const TopoDS_Wire& W,
89                         const gp_Dir& D,
90                         Bnd_Box& Box,
91                         gp_Trsf& Tf)
92 {
93   // Calculate approximate barycenter
94   BRepTools_WireExplorer Exp(W);
95 // Class BRep_Tool without fields and without Constructor :
96 //  BRep_Tool BT;
97   gp_XYZ Bary(0.,0.,0.);
98   Standard_Integer nb;
99
100   for (nb=0; Exp.More(); Exp.Next()) {
101 //    Bary += BT.Pnt(Exp.CurrentVertex()).XYZ();
102     Bary += BRep_Tool::Pnt(Exp.CurrentVertex()).XYZ();
103     nb++;
104   }
105   Bary /= nb;
106
107   // Calculate the Transformation  
108   gp_Ax3 N(Bary, D);
109   Tf.SetTransformation(N);
110   BRepAdaptor_Curve AC;
111 //  BndLib_Add3dCurve BC;  
112
113   // transformation to the wire
114   TopoDS_Wire TheW = W;
115   TopLoc_Location Loc(Tf);
116   TheW.Location(Loc);
117
118
119   // Calculate the box
120   Box.SetVoid();
121   for (Exp.Init(TheW); Exp.More(); Exp.Next()) {
122     AC.Initialize(Exp.Current());
123 //    BC.Add(AC, 0.1, Box);
124     BndLib_Add3dCurve::Add(AC, 0.1, Box);
125   }
126 }
127
128 //=======================================================================
129 //function : Length
130 //purpose  : 
131 //======================================================================
132 static Standard_Real Longueur(const Bnd_Box& WBox,
133                               const Bnd_Box& SBox,
134                               gp_Dir& D,
135                               gp_Pnt& P)
136 {
137   // face of the box most remoted from the face input in 
138   // the direction of skin
139   Standard_Real Xmin,Ymin,Zmin,Xmax,Ymax,Zmax,WZmin,WZmax,L;
140
141   //"coord" of the box
142   WBox.Get(Xmin,Ymin,Zmin,Xmax,Ymax,Zmax);
143   WZmin = Zmin;
144   WZmax = Zmax;   
145
146   SBox.Get(Xmin,Ymin,Zmin,Xmax,Ymax,Zmax);
147   P.SetCoord( (Xmin+Xmax)/2, (Ymin+Ymax)/2, Zmax);
148
149   if (Zmax < WZmin) {
150     // Skin in the wrong direction. Invert...
151     D.Reverse();
152     L = WZmax - Zmin;
153     P.SetZ(Zmin);
154   }
155   else {
156     L = Zmax - WZmin;
157   }
158   return L;
159 }
160
161 //=======================================================================
162 //function : GoodOrientation
163 //purpose  : Check if the law is oriented to have an exterior skin
164 //======================================================================
165 static Standard_Boolean GoodOrientation(const Bnd_Box& B,
166                                         const Handle(BRepFill_LocationLaw)& Law,
167                                         const gp_Dir& D)
168 {
169   Standard_Real f, l, r, t;
170   Standard_Real aXmin, aYmin, aZmin, aXmax, aYmax, aZmax;
171
172   B.Get(aXmin, aYmin, aZmin, aXmax, aYmax, aZmax);
173   gp_Pnt P1(aXmin, aYmin, aZmin), P2(aXmax, aYmax, aZmax);
174   gp_Vec V(P1, P2);
175
176   Law->CurvilinearBounds(Law->NbLaw(), f, l);
177   r = V.Magnitude()/l;
178  
179   Standard_Integer ii, Ind;
180 //#ifndef OCCT_DEBUG
181   Standard_Integer Nb = (Standard_Integer) (4+(10*r));
182 //#else
183 //  Standard_Integer Nb = 4+(10*r);
184 //#endif
185   r = l/Nb;
186
187   Nb++; // Number of points
188
189   TColgp_Array1OfPnt Pnts(1, Nb);
190   Handle(Adaptor3d_HCurve) AC;
191   gp_XYZ Bary(0.,0.,0.);
192   
193   for (ii=1; ii<=Nb; ii++) {
194     Law->Parameter((ii-1)*r, Ind, t);
195     AC = Law->Law(Ind)->GetCurve();
196     AC->D0(t,  Pnts(ii));
197     Bary+= Pnts(ii).XYZ();
198   }
199
200   Bary /= Nb;
201   gp_Pnt Centre(Bary);
202   gp_Vec Normal(D.XYZ());
203   Standard_Real Angle = 0;
204   gp_Vec Ref(Centre,  Pnts(1));
205
206   for (ii=2; ii<=Nb; ii++) {
207     gp_Vec R(Centre, Pnts(ii));
208     Angle += Ref.AngleWithRef(R, Normal);
209     Ref = R;
210   }
211
212   return (Angle >= 0);
213 }
214
215 //=======================================================================
216 //function : Constructeur
217 //purpose  : 
218 //======================================================================
219  BRepFill_Draft::BRepFill_Draft(const TopoDS_Shape& S,
220                                 const gp_Dir& Dir,
221                                 const Standard_Real Angle)
222 {
223   myLoc.Nullify();
224   mySec.Nullify();
225   myFaces.Nullify();
226   mySections.Nullify(); 
227    
228   switch (S.ShapeType()) {
229   case TopAbs_WIRE :
230     {
231       myWire = TopoDS::Wire(S);
232       break;
233     }
234   case TopAbs_FACE :
235     {  
236       TopoDS_Iterator Exp (S);
237       myWire = TopoDS::Wire(Exp.Value());
238       break;
239     }
240   case TopAbs_SHELL :
241     {  
242       TopTools_ListOfShape List;
243       TopTools_IndexedDataMapOfShapeListOfShape edgemap;
244       TopExp::MapShapesAndAncestors(S,TopAbs_EDGE,TopAbs_FACE,edgemap);
245       Standard_Integer iedge, nbf;
246       for (iedge = 1; iedge <= edgemap.Extent(); iedge++) {   
247         const TopoDS_Edge& theEdge = TopoDS::Edge(edgemap.FindKey(iedge));    
248         // skip degenerated edges
249         if (!BRep_Tool::Degenerated(theEdge)) {
250           nbf = edgemap(iedge).Extent();
251           if (nbf==1) List.Append(theEdge);
252         }
253       }
254
255       if( List.Extent()>0) {
256         BRepLib_MakeWire MW;
257         MW.Add(List);
258         BRepLib_WireError Err = MW.Error();
259         if (Err == BRepLib_WireDone) {
260           myWire = MW.Wire();
261         }
262         else {
263 #ifdef OCCT_DEBUG
264           cout << "Error in MakeWire" << endl;
265 #endif 
266           throw Standard_ConstructionError("BRepFill_Draft");
267         }
268       }
269       else {
270 #ifdef OCCT_DEBUG
271           cout << "No Free Borders !" << endl;
272 #endif 
273           throw Standard_ConstructionError("BRepFill_Draft");
274       }
275       break;
276     }
277     default :
278       throw Standard_ConstructionError("BRepFill_Draft");
279   }
280
281   // Attention to closed non declared wires !
282   if (!myWire.Closed()) {
283     TopoDS_Vertex Vf, Vl;
284     TopExp::Vertices(myWire, Vf, Vl);
285     if (Vf.IsSame(Vl)) myWire.Closed(Standard_True);
286   }
287 #if DRAW
288   if (Affich) {
289     DBRep::Set("TheWire", myWire);
290   }
291 #endif
292
293   myAngle = Abs(Angle);
294   myDir = Dir;
295   myTop = S;
296   myDone = Standard_False;
297   myTol = 1.e-4;
298   myCont = GeomAbs_C1;
299   SetOptions();
300   SetDraft();
301 }
302
303 //=======================================================================
304 //function :SetOptions
305 //purpose  : Defines the style
306 //======================================================================
307  void BRepFill_Draft::SetOptions(const BRepFill_TransitionStyle Style,
308                                  const Standard_Real Min,
309                                  const Standard_Real Max)
310 {
311   myStyle =  Style;
312   angmin = Min; 
313   angmax = Max;
314 }
315
316 //=======================================================================
317 //function :SetDraft
318 //purpose  :
319 //======================================================================
320  void BRepFill_Draft::SetDraft(const Standard_Boolean Internal)
321 {
322   IsInternal = Internal;
323 }
324
325
326 //=======================================================================
327 //function :Perform
328 //purpose  : calculate a surface of skinning
329 //======================================================================
330  void BRepFill_Draft::Perform(const Standard_Real LengthMax)
331 {
332   Handle(Geom_Surface) S;
333   S.Nullify();
334   Bnd_Box WBox;//, SBox;
335   gp_Trsf Trsf;
336
337
338   ComputeTrsf(myWire, myDir, WBox, Trsf); 
339   Init(S, LengthMax, WBox);
340   BuildShell(S); 
341   Sewing(); 
342 }
343
344 //=======================================================================
345 //function :Perform
346 //purpose  : calculate a surface of skinning
347 //======================================================================
348  void BRepFill_Draft::Perform(const Handle(Geom_Surface)& Surface,
349                               const Standard_Boolean KeepInsideSurface)
350 {
351   Bnd_Box WBox, SBox;
352   gp_Trsf Trsf;
353   gp_Pnt Pt;
354   Standard_Real L;
355
356   ComputeTrsf(myWire, myDir, WBox, Trsf);
357
358   // box with bounds of the stop surface  
359   Handle(Geom_Surface) Surf;
360   Surf =   Handle(Geom_Surface)::DownCast(Surface->Transformed(Trsf));
361   GeomAdaptor_Surface S1 (Surf);   
362 //  BndLib_AddSurface AS; 
363 //  AS.Add(S1, 0.1, SBox);
364   BndLib_AddSurface::Add(S1, 0.1, SBox);
365
366   // calculate the maximum length of the rule.
367   L = Longueur(WBox, SBox, myDir, Pt);
368   L /= Abs(Cos(myAngle));
369
370   // Construction
371   Init(Surface, L, WBox);
372   BuildShell(Surface, !KeepInsideSurface);
373   Sewing(); 
374 }
375
376 //================================================================
377 //function :Perform
378 //purpose  : calculate the surface of skinning, stopped by a shape
379 //================================================================
380  void BRepFill_Draft::Perform(const TopoDS_Shape& StopShape,
381                               const Standard_Boolean KeepOutSide)
382 {
383   Bnd_Box WBox, SBox;
384   gp_Trsf Trsf;
385   gp_Pnt Pt;
386   Standard_Real L;
387
388   ComputeTrsf(myWire, myDir, WBox, Trsf); 
389
390 // bounding box of the stop shape
391   Bnd_Box BSurf;//, TheBox;
392   Standard_Real Umin, Umax, Vmin, Vmax;
393 //  BRepTools B;
394 //  BRep_Tool BT;
395   Handle(Geom_Surface) Surf;
396  
397 //  BndLib_AddSurface AS;
398   
399   TopExp_Explorer Ex (StopShape, TopAbs_FACE);
400
401   SBox.SetVoid();
402   while (Ex.More()) { // parse faces of the stop shape
403 //    B.UVBounds(TopoDS::Face(Ex.Current()), Umin,Umax,Vmin,Vmax); 
404     BRepTools::UVBounds(TopoDS::Face(Ex.Current()), Umin,Umax,Vmin,Vmax); 
405     Surf = Handle(Geom_Surface)::DownCast(
406 //     BT.Surface(TopoDS::Face(Ex.Current()))->Transformed(Trsf) ); 
407      BRep_Tool::Surface(TopoDS::Face(Ex.Current()))->Transformed(Trsf) ); 
408     GeomAdaptor_Surface S1 (Surf);
409 // bounding box of the current face
410 //    AS.Add(S1, Umin, Umax, Vmin, Vmax, 0.1, BSurf); 
411     BndLib_AddSurface::Add(S1, Umin, Umax, Vmin, Vmax, 0.1, BSurf);
412     SBox.Add(BSurf);    // group boxes
413     Ex.Next();
414   }// while_Ex
415
416   // calculate the maximum length of the rule.
417   L = Longueur(WBox, SBox, myDir, Pt);
418   L /= Abs(Cos(myAngle));
419
420 // surface of stop
421   gp_Trsf Inv;
422   Inv = Trsf.Inverted(); // inverted transformation
423   Pt.Transform(Inv); // coordinate in the absolute reference
424   Handle(Geom_Plane) Plan = new (Geom_Plane)(Pt, myDir);
425   Surf = new (Geom_RectangularTrimmedSurface) (Plan,-L, L, -L, L); 
426
427 #if DRAW
428   if (Affich) {
429     char* Temp = "ThePlan" ;
430     DrawTrSurf::Set(Temp, Surf);
431 //    DrawTrSurf::Set("ThePlan", Surf);
432   }
433 #endif
434
435 // Sweeping and restriction
436   Init(Plan,  L*1.01, WBox);
437   BuildShell(Surf, Standard_False);
438   Fuse(StopShape,  KeepOutSide);
439   Sewing();
440 }
441
442 //=======================================================================
443 //function : Init
444 //purpose  : Construction of laws.
445 //======================================================================
446  void BRepFill_Draft::Init(const Handle(Geom_Surface)& ,
447                            const Standard_Real Length,
448                            const Bnd_Box&  Box)
449 {
450   Standard_Boolean B;
451
452 // law of positioning   
453   Handle(GeomFill_LocationDraft) Loc
454     = new (GeomFill_LocationDraft) (myDir, myAngle); 
455   myLoc = new (BRepFill_DraftLaw) (myWire, Loc);
456
457   B = GoodOrientation(Box, myLoc, myDir);
458
459   if (IsInternal ^ (!B) )  {
460     myAngle = - myAngle;
461     Loc->SetAngle(myAngle);
462     myLoc = new (BRepFill_DraftLaw) (myWire, Loc);
463   }
464
465   myLoc->CleanLaw(angmin); // Clean small discontinuities.
466
467 // law of section
468 // generating line is straight and parallel to binormal.
469   gp_Pnt P(0, 0, 0);
470   gp_Vec D (0., 1., 0.);
471
472 // Control of the orientation
473   Standard_Real f,l;
474   myLoc->Law(1)->GetDomain(f,l);
475   gp_Mat M;
476
477   gp_Vec Bid; 
478   myLoc->Law(1)->D0( (f+l)/2, M, Bid);
479   gp_Dir BN(M.Column(2));
480   
481   Standard_Real ang = myDir.Angle(BN);
482   if (ang > M_PI/2) D.Reverse();    
483   Handle(Geom_Line) L = new (Geom_Line) (P, D);
484
485   Handle(Geom_Curve) TC = new (Geom_TrimmedCurve) (L, 0, Length);
486
487
488 #if DRAW
489   if (Affich > 2) {
490      TC = new (Geom_Circle) (gp::XOY(), Length);
491   }
492 #endif 
493
494   BRepLib_MakeEdge ME(TC);
495   TopoDS_Edge EG = ME.Edge();
496   
497   BRepLib_MakeWire MW(EG);
498   TopoDS_Wire G = MW.Wire();
499
500   mySec = new (BRepFill_ShapeLaw) (G, Standard_True);
501 }
502
503
504 //=======================================================================
505 //function : BuildShell
506 //purpose  : Construction of the skinning surface 
507 //======================================================================
508  void BRepFill_Draft::BuildShell(const Handle(Geom_Surface)& Surf,
509                                  const Standard_Boolean  KeepOutSide) 
510 {
511 // construction of the surface
512   BRepFill_Sweep Sweep(mySec, myLoc, Standard_True);
513   Sweep.SetTolerance(myTol);
514   Sweep.SetAngularControl(angmin, angmax);
515   TopTools_MapOfShape Dummy;
516   BRepFill_DataMapOfShapeHArray2OfShape Dummy2;
517   BRepFill_DataMapOfShapeHArray2OfShape Dummy3;
518   Sweep.Build(Dummy, Dummy2, Dummy3, myStyle, myCont);
519   if (Sweep.IsDone()) {
520     myShape = Sweep.Shape();
521     myShell = TopoDS::Shell(myShape);
522     myFaces    = Sweep.SubShape();
523     mySections = Sweep.Sections();    
524     myDone = Standard_True;
525     // Control of the orientation
526     Standard_Boolean out=Standard_True;
527     TopExp_Explorer ex(myShell,TopAbs_FACE);
528     TopoDS_Face F;
529     F = TopoDS::Face(ex.Current());
530     BRepAdaptor_Surface SF(F);
531     Standard_Real u, v;
532     gp_Pnt P;
533     gp_Vec V1, V2, V;
534     u = SF.FirstUParameter();
535     v = SF.FirstVParameter();
536     SF.D1(u,v,P,V1,V2);
537     V = V1.Crossed(V2);
538     if (F.Orientation() == TopAbs_REVERSED) V.Reverse();
539     if (V.Magnitude() > 1.e-10) {
540       out = myDir.Angle(V) > M_PI/2;
541     }
542     if (out == IsInternal) {
543       myShell.Reverse();
544       myShape.Reverse();
545     }    
546   }
547   else {
548     myDone = Standard_False;
549     return;
550   }
551
552   if (!Surf.IsNull()) { // Add the face at end
553
554   // Waiting the use of traces & retriction in BRepFill_Sweep
555   // Make Fuse.
556     BRepLib_MakeFace MkF;
557     MkF.Init(Surf, Standard_True, Precision::Confusion());
558     Fuse(MkF.Face(), KeepOutSide);
559   }
560 }
561
562
563 //=======================================================================
564 //function : Fuse
565 //purpose  : Boolean operation between the skin and the  
566 //           stop shape 
567 //======================================================================
568  Standard_Boolean BRepFill_Draft::Fuse(const TopoDS_Shape& StopShape,
569                            const Standard_Boolean KeepOutSide)
570 {
571   BRep_Builder B;
572   Standard_Boolean issolid = Standard_False;
573   TopoDS_Solid Sol1, Sol2;
574   TopAbs_State State1 = TopAbs_OUT,  State2 = TopAbs_OUT;
575
576
577   if (myShape.ShapeType()==TopAbs_SOLID) {
578     Sol1 = TopoDS::Solid(myShape);
579     issolid = Standard_True;
580   }
581   else {
582     B.MakeSolid(Sol1);
583     B.Add(Sol1, myShape); // shell => solid (for fusion)    
584   }
585
586
587   switch (StopShape.ShapeType()) {
588   case TopAbs_COMPOUND :
589     {
590       TopoDS_Iterator It(StopShape);
591       return Fuse(It.Value(), KeepOutSide);
592     }
593   case TopAbs_SOLID : 
594     {
595       Sol2 = TopoDS::Solid(StopShape);
596       break;
597     }
598  case TopAbs_SHELL :
599     {
600       B.MakeSolid(Sol2);
601       B.Add(Sol2, StopShape); // shell => solid (for fusion) 
602       break;
603     }
604
605  case TopAbs_FACE :
606     {
607       TopoDS_Shell S;
608       B.MakeShell(S);
609       B.Add(S, StopShape);
610       S.Closed (BRep_Tool::IsClosed (S));
611       B.MakeSolid(Sol2);
612       B.Add(Sol2, S); // shell => solid (for fusion)
613       break;
614     }
615
616   default :
617     {
618       return Standard_False; // Impossible to do
619     }
620   }
621       
622   BRepAlgo_DSAccess DSA;
623   DSA.Load(Sol1, Sol2);
624   DSA.Intersect(Sol1, Sol2); // intersection of 2 solids
625  
626 // removal of edges corresponding to "unused" intersections
627   Standard_Integer NbPaquet;
628 // gp_Pnt P1,P2;
629   TopoDS_Vertex V,V1;
630   TopTools_ListOfShape List;
631   List  = DSA.GetSectionEdgeSet();// list of edges
632   
633   NbPaquet = List.Extent();
634
635   if (NbPaquet == 0) {
636 #if DRAW
637     cout << "No fusion" << endl;
638     DBRep::Set("DepPart", Sol1);
639     DBRep::Set("StopPart", Sol2);
640 #endif
641     return Standard_False;
642   }
643
644   if (NbPaquet > 1) {
645     // It is required to select packs.
646     TColStd_Array1OfReal Dist(1, NbPaquet);
647     TopTools_ListIteratorOfListOfShape it(List);
648     Standard_Real D, Dmin = 1.e10;
649     Standard_Integer ii;
650  
651     //Classify the packs by distance.
652     BRepExtrema_DistShapeShape Dist2;
653     Dist2.LoadS1( myWire );
654     for (ii=1; it.More();it.Next(),ii++){
655       Dist2.LoadS2( it.Value() );
656       Dist2.Perform();
657       if (Dist2.IsDone()) {
658         D = Dist2.Value();
659         Dist(ii) = D;
660         if (D < Dmin) Dmin = D;
661       }
662       else
663         Dist(ii) = 1.e10;
664     }
665
666     // remove edges "farther" than Dmin
667     for (ii=1, it.Initialize(List); it.More();it.Next(), ii++){
668       if (Dist(ii) > Dmin) {
669         DSA.SuppressEdgeSet(it.Value());
670       }
671 #if DRAW
672      else if (Affich) {
673         DBRep::Set("KeepEdges", it.Value());
674       }
675 #endif
676     }
677   }
678
679   if (StopShape.ShapeType() != TopAbs_SOLID) {
680     // It is required to choose the state by the geometry
681
682     //(1) Return an edge of section
683     List  = DSA.GetSectionEdgeSet();// list of edges
684     TopTools_ListIteratorOfListOfShape it(List);
685     TopoDS_Iterator iter(it.Value());
686     TopoDS_Edge E = TopoDS::Edge(iter.Value());
687
688     //(2) Return geometry on StopShape
689 // Class BRep_Tool without fields and without Constructor :
690 //    BRep_Tool BT;
691     Handle(Geom_Surface) S;
692     Handle(Geom2d_Curve) C2d;
693     gp_Pnt2d P2d;
694     Standard_Real f,l;  
695     TopLoc_Location L;
696 //    BT.CurveOnSurface(E, C2d, S, L, f, l, 2);
697     BRep_Tool::CurveOnSurface(E, C2d, S, L, f, l, 2);
698
699     // Find a normal.
700     C2d->D0((f+l)/2,P2d); 
701     GeomLProp_SLProps SP(S, P2d.X(), P2d.Y(), 1, 1.e-12);
702     if (! SP.IsNormalDefined()) {
703       C2d->D0((3*f+l)/4,P2d);
704       SP.SetParameters(P2d.X(), P2d.Y());
705       if ( !SP.IsNormalDefined()) {
706         C2d->D0((f+3*l)/4,P2d);
707         SP.SetParameters(P2d.X(), P2d.Y());
708       }
709     }
710
711     // Subtract State1
712     if (myDir.Angle(SP.Normal()) < M_PI/2)  State1 = TopAbs_IN;
713     else  State1 = TopAbs_OUT;
714   }
715
716   if (! KeepOutSide) { // Invert State2;
717     if (State2 == TopAbs_IN) State2 = TopAbs_OUT;
718     else State2 = TopAbs_IN;
719   }
720  
721 //recalculate the final shape
722   TopoDS_Shape result = DSA.Merge(State1, State2); 
723
724   if (issolid) myShape =  result;
725   else {
726     TopExp_Explorer Exp;
727     Exp.Init(result, TopAbs_SHELL);
728     if (Exp.More()) myShape = Exp.Current();
729   }
730
731 // Update the History 
732   Standard_Integer ii;
733   for (ii=1; ii<=myLoc->NbLaw(); ii++) {
734     const TopTools_ListOfShape& L = DSA.Modified(myFaces->Value(1,ii));
735     if (L.Extent()>0) 
736       myFaces->SetValue(1, ii, L.First());
737   }
738   for (ii=1; ii<=myLoc->NbLaw()+1; ii++) {
739     const TopTools_ListOfShape& L = DSA.Modified(mySections->Value(1,ii));
740     if (L.Extent()>0) 
741       mySections->SetValue(1, ii, L.First());
742   } 
743  
744   return Standard_True;
745 }
746
747 //=======================================================================
748 //function : Sewing 
749 //purpose  : Assemble the skin with the above face
750 //======================================================================
751  Standard_Boolean BRepFill_Draft::Sewing()
752 {
753   Standard_Boolean ToAss;
754   Standard_Boolean Ok = Standard_False;
755   ToAss = (myTop.ShapeType() != TopAbs_WIRE);
756
757   if ((!ToAss) || (!myDone)) return Standard_False;
758
759    // Assembly make a shell from the faces of the shape + the input shape
760   Handle(BRepBuilderAPI_Sewing) Ass =  new BRepBuilderAPI_Sewing(5*myTol, Standard_True, 
761                                                      Standard_True, Standard_False);
762   Ass->Add(myShape);
763   Ass->Add(myTop);
764   ToAss = Standard_True;
765  
766
767   Standard_Integer NbCE;
768   
769   Ass->Perform();
770   // Check if the assembly is real.
771   NbCE = Ass->NbContigousEdges();
772   
773   if (NbCE > 0) {
774     TopoDS_Shape res;
775     res = Ass->SewedShape();   
776     if ((res.ShapeType() == TopAbs_SHELL)||
777         (res.ShapeType() == TopAbs_SOLID))  {
778       myShape = res;
779       Ok = Standard_True;
780     }
781     else if (res.ShapeType() == TopAbs_COMPOUND) {
782       TopoDS_Iterator It(res);
783       res = It.Value();
784       It.Next();
785       if (!It.More()) {//Only one part => this is correct
786         myShape = res;
787         Ok = Standard_True;
788       }
789     }
790   }
791
792   if (Ok) {
793     // Update the History
794     Standard_Integer ii;
795     for (ii=1; ii<=myLoc->NbLaw(); ii++) {
796       if (Ass->IsModified(myFaces->Value(1,ii)))     
797         myFaces->SetValue(1, ii, 
798                           Ass->Modified(myFaces->Value(1,ii)));
799     }
800     for (ii=1; ii<=myLoc->NbLaw()+1; ii++) {
801       if (Ass->IsModified(mySections->Value(1,ii))) 
802         mySections->SetValue(1, ii, 
803                              Ass->Modified(mySections->Value(1,ii)));
804     }   
805
806     if (myShape.Closed()) { // Make a Solid
807       TopoDS_Solid solid;
808       BRep_Builder BS;
809       BS.MakeSolid(solid); 
810       BS.Add(solid,TopoDS::Shell(myShape));
811       
812       BRepClass3d_SolidClassifier SC(solid);
813       SC.PerformInfinitePoint(Precision::Confusion());
814       if ( SC.State() == TopAbs_IN) {
815         BS.MakeSolid(solid);
816         myShape.Reverse();
817         BS.Add(solid,TopoDS::Shell(myShape));
818       }
819       myShape = solid;
820     }    
821   }
822 #ifdef OCCT_DEBUG
823   else cout << "Draft : No assembly !" << endl;
824 #endif
825   return Ok;
826 }
827
828 //=======================================================================
829 //function : Generated
830 //purpose  : return a sub-part generated by sweeping
831 //======================================================================
832  const TopTools_ListOfShape& 
833  BRepFill_Draft::Generated(const TopoDS_Shape& S) 
834 {
835   myGenerated.Clear();
836   TopoDS_Edge E;
837   Standard_Integer ii;
838   E = TopoDS::Edge(S);
839   if (E.IsNull()) {
840    for (ii=0; ii<=myLoc->NbLaw(); ii++)
841       if (E.IsSame(myLoc->Vertex(ii))) {
842         myGenerated.Append(mySections->Value(1, ii+1));
843         break;
844       } 
845   }
846   else {
847     for (ii=1; ii<=myLoc->NbLaw(); ii++)
848       if (E.IsSame(myLoc->Edge(ii))) {
849         myGenerated.Append(myFaces->Value(1, ii));
850         break;
851       }
852   }
853   
854   return myGenerated;
855 }
856
857 //=======================================================================
858 //function : Shape
859 //purpose  : return the complete shape
860 //======================================================================
861  TopoDS_Shape BRepFill_Draft::Shape()const
862 {
863   return myShape;
864 }
865
866 //=====================================================================
867 //function : Shell
868 //purpose  : surface of skinning with the input face (=>shell)
869 //=====================================================================
870  TopoDS_Shell BRepFill_Draft::Shell()const
871 {
872   return myShell;
873 }
874
875 //=======================================================================
876 //function : IsDone
877 //purpose  : 
878 //======================================================================
879  Standard_Boolean BRepFill_Draft::IsDone()const
880 {
881   return myDone;
882 }