0024023: Revamp the OCCT Handle -- downcast (automatic)
[occt.git] / src / BRepFeat / BRepFeat_MakeRevol.cxx
1 // Created on: 1996-02-13
2 // Created by: Jacques GOUSSARD
3 // Copyright (c) 1996-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 #include <BRepFeat_MakeRevol.ixx>
18
19 #include <BRepFeat.hxx>
20 #include <LocOpe.hxx>
21 #include <LocOpe_Revol.hxx>
22 #include <LocOpe_Gluer.hxx>
23 #include <LocOpe_FindEdges.hxx>
24 #include <LocOpe_SequenceOfCirc.hxx>
25 #include <LocOpe_BuildShape.hxx>
26 #include <LocOpe_CSIntersector.hxx>
27 #include <LocOpe_PntFace.hxx>
28
29 #include <gp_Vec.hxx>
30 #include <gp_Pnt.hxx>
31 #include <gp_Ax1.hxx>
32 #include <gp_Pnt2d.hxx>
33 #include <TColgp_SequenceOfPnt.hxx>
34 #include <Geom_Surface.hxx>
35 #include <Geom_Curve.hxx>
36 #include <Geom_Circle.hxx>
37 #include <Geom2d_Curve.hxx>
38 #include <Geom_RectangularTrimmedSurface.hxx>
39 #include <Bnd_Box.hxx>
40
41 #include <BRepSweep_Revol.hxx>
42 #include <BRep_Builder.hxx>
43 #include <BRep_Tool.hxx>
44 #include <TopExp_Explorer.hxx>
45 #include <TopTools_MapOfShape.hxx>
46 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
47 #include <TopTools_ListIteratorOfListOfShape.hxx>
48 #include <TopTools_MapIteratorOfMapOfShape.hxx>
49 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
50
51 #include <BRepAlgoAPI_Cut.hxx>
52 #include <BRepAlgoAPI_Fuse.hxx>
53
54 #include <BRepLib_MakeFace.hxx>
55 #include <BRepTools_Modifier.hxx>
56 #include <BRepTools_TrsfModification.hxx>
57
58 #include <Standard_ConstructionError.hxx>
59
60 #include <TopoDS_Solid.hxx>
61 #include <TopoDS_Compound.hxx>
62 #include <TopoDS_Shell.hxx>
63
64 #include <gp_Pln.hxx>
65 #include <Geom_Plane.hxx>
66
67
68 #include <TopExp.hxx>
69 #include <Precision.hxx>
70 #include <BRepBndLib.hxx>
71 #include <BRepBuilderAPI.hxx>
72 #include <BRepFeat.hxx>
73
74 #include <ElCLib.hxx>
75
76 #ifdef OCCT_DEBUG
77 extern Standard_Boolean BRepFeat_GettraceFEAT();
78 #endif
79
80 static void MajMap(const TopoDS_Shape&, // base
81   const LocOpe_Revol&,
82   TopTools_DataMapOfShapeListOfShape&, // myMap
83   TopoDS_Shape&,  // myFShape
84   TopoDS_Shape&); // myLShape
85
86
87 static void VerifGluedFaces(const TopoDS_Face& theSkface,
88   const TopoDS_Shape& thePbase,
89   Handle(Geom_Curve)& theBCurve,
90   TColGeom_SequenceOfCurve& theCurves,
91   LocOpe_Revol& theRevol,
92   TopTools_DataMapOfShapeShape& theMap);
93
94
95 static Standard_Boolean ToFuse(const TopoDS_Face& ,
96   const TopoDS_Face&);
97
98
99
100
101 //=======================================================================
102 //function : Init
103 //purpose  : 
104 //=======================================================================
105
106 void BRepFeat_MakeRevol::Init(const TopoDS_Shape& Sbase,
107   const TopoDS_Shape& Pbase,
108   const TopoDS_Face& Skface,
109   const gp_Ax1& Axis,
110   const Standard_Integer Mode,
111   const Standard_Boolean Modify)
112 {
113 #ifdef OCCT_DEBUG
114   Standard_Boolean trc = BRepFeat_GettraceFEAT();
115   if (trc) cout << "BRepFeat_MakeRevol::Init" << endl;
116 #endif
117   myAxis   = Axis;
118   myPbase  = Pbase;
119   mySbase  = Sbase;
120   BasisShapeValid();
121   mySkface = Skface;
122   SketchFaceValid();
123   myPbase  = Pbase;
124   mySlface.Clear();
125   if(Mode == 0) {
126     myFuse   = Standard_False;
127     myJustFeat = Standard_False;
128   }
129   else if(Mode == 1) {
130     myFuse   = Standard_True;
131     myJustFeat = Standard_False;
132   }
133   else if(Mode == 2) {
134     myFuse   = Standard_True;
135     myJustFeat = Standard_True;
136   }
137   else {    
138   }
139   myModify = Modify;
140   myJustGluer = Standard_False;
141
142   //-------------- ifv
143   //  mySkface.Nullify();
144   //-------------- ifv
145
146
147   myShape.Nullify();
148   myMap.Clear();
149   myFShape.Nullify();
150   myLShape.Nullify();
151   TopExp_Explorer exp;
152   for (exp.Init(mySbase,TopAbs_FACE);exp.More();exp.Next()) {
153     TopTools_ListOfShape thelist;
154     myMap.Bind(exp.Current(), thelist);
155     myMap(exp.Current()).Append(exp.Current());
156   }
157 #ifdef OCCT_DEBUG
158   if (trc) {
159     if (myJustFeat)  cout << " Just Feature" << endl;
160     if (myFuse)  cout << " Fuse" << endl;
161     if (!myFuse)  cout << " Cut" << endl;
162     if (!myModify) cout << " Modify = 0" << endl;
163   }
164 #endif 
165 }
166
167
168 //=======================================================================
169 //function : Add
170 //purpose  : add faces add edges of sliding
171 //=======================================================================
172
173 void BRepFeat_MakeRevol::Add(const TopoDS_Edge& E,
174   const TopoDS_Face& F)
175 {
176 #ifdef OCCT_DEBUG
177   Standard_Boolean trc = BRepFeat_GettraceFEAT();
178   if (trc) cout << "BRepFeat_MakeRevol::Add(Edge,face)" << endl;
179 #endif
180   TopExp_Explorer exp;
181   for (exp.Init(mySbase,TopAbs_FACE);exp.More();exp.Next()) {
182     if (exp.Current().IsSame(F)) {
183       break;
184     }
185   }
186   if (!exp.More()) {
187     Standard_ConstructionError::Raise();
188   }
189
190   for (exp.Init(myPbase,TopAbs_EDGE);exp.More();exp.Next()) {
191     if (exp.Current().IsSame(E)) {
192       break;
193     }
194   }
195   if (!exp.More()) {
196     Standard_ConstructionError::Raise();
197   }
198
199   if (!mySlface.IsBound(F)) {
200     TopTools_ListOfShape thelist;
201     mySlface.Bind(F, thelist);
202   }
203   TopTools_ListIteratorOfListOfShape itl(mySlface(F));
204   for (; itl.More();itl.Next()) {
205     if (itl.Value().IsSame(E)) {
206       break;
207     }
208   }
209   if (!itl.More()) {
210     mySlface(F).Append(E);
211   }
212 }
213
214
215 //=======================================================================
216 //function : Perform
217 //purpose  : 
218 //=======================================================================
219
220 void BRepFeat_MakeRevol::Perform(const Standard_Real Angle)
221 {
222 #ifdef OCCT_DEBUG
223   Standard_Boolean trc = BRepFeat_GettraceFEAT();
224   if (trc) cout << "BRepFeat_MakeRevol::Perform(Angle)" << endl;
225 #endif
226   mySFrom.Nullify();
227   ShapeFromValid();
228   mySUntil.Nullify();
229   ShapeUntilValid();
230   myGluedF.Clear();
231   myPerfSelection = BRepFeat_NoSelection;
232   PerfSelectionValid();
233   Standard_Boolean RevolComp = (2*M_PI-Abs(Angle) <= Precision::Angular());
234   LocOpe_Revol theRevol;
235   Standard_Real angledec = 0.;
236   TopExp_Explorer exp;
237   if(RevolComp) {
238     /*
239     if (!mySkface.IsNull() || !mySlface.IsEmpty()) {
240     for (exp.Init(mySbase,TopAbs_FACE); exp.More(); exp.Next()) {
241     if (exp.Current().IsSame(mySkface)) {
242     angledec = M_PI/5; // pourquoi pas
243     if (myFuse) angledec = -angledec;
244     break;
245     }
246     }
247     }
248     */
249     mySkface.Nullify();
250   }
251   if(angledec == 0.) theRevol.Perform(myPbase, myAxis, Angle);
252   else theRevol.Perform(myPbase, myAxis, Angle, angledec);
253
254   TopoDS_Shape VraiRevol = theRevol.Shape();
255
256   MajMap(myPbase,theRevol,myMap,myFShape,myLShape);
257
258   myGShape = VraiRevol;
259   GeneratedShapeValid();
260   TopoDS_Shape Base = theRevol.FirstShape();
261   exp.Init(Base, TopAbs_FACE);
262   TopoDS_Face theBase = TopoDS::Face(exp.Current());
263   exp.Next();
264   if(exp.More()) {
265     NotDone();
266     myStatusError = BRepFeat_InvFirstShape;
267     return;
268   }
269
270   TopoDS_Face FFace;
271
272   Standard_Boolean found = Standard_False;
273
274   if(!mySkface.IsNull() || !mySlface.IsEmpty()) {
275     if(myLShape.ShapeType() == TopAbs_WIRE) {
276       TopExp_Explorer ex1(VraiRevol, TopAbs_FACE);
277       for(; ex1.More(); ex1.Next()) {
278         TopExp_Explorer ex2(ex1.Current(), TopAbs_WIRE);
279         for(; ex2.More(); ex2.Next()) {
280           if(ex2.Current().IsSame(myLShape)) {
281             FFace = TopoDS::Face(ex1.Current());
282             found = Standard_True;
283             break;
284           }
285         }
286         if(found) break;
287       }
288     }
289
290     TopExp_Explorer exp(mySbase, TopAbs_FACE);
291     for(; exp.More(); exp.Next()) {
292       const TopoDS_Face& ff = TopoDS::Face(exp.Current());
293       if(ToFuse(ff, FFace)) {
294         TopTools_DataMapOfShapeListOfShape sl;
295         if(!FFace.IsSame(myPbase) && BRepFeat::IsInside(ff, FFace)) 
296           break;
297       }
298     }
299   }
300   GluedFacesValid();
301   if (!mySkface.IsNull()) {
302     VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF);
303   }
304
305   if(myGluedF.IsEmpty()) {
306     if(myFuse == 1) {
307       //modified by NIZNHY-PKV Thu Mar 21 18:15:06 2002 f
308       //BRepAlgo_Fuse f(mySbase, myGShape);
309       //myShape = f.Shape();
310       //UpdateDescendants(f.Builder(), myShape, Standard_False);
311       BRepAlgoAPI_Fuse f(mySbase, myGShape);
312       myShape = f.Shape();
313       UpdateDescendants(f, myShape, Standard_False);
314       //modified by NIZNHY-PKV Thu Mar 21 18:15:11 2002 t
315       Done();
316     }
317     else if(myFuse == 0) {
318       //modified by NIZNHY-PKV Thu Mar 21 18:15:37 2002 f
319       //BRepAlgo_Cut c(mySbase, myGShape);
320       //myShape = c.Shape();
321       //UpdateDescendants(c.Builder(), myShape, Standard_False);
322       BRepAlgoAPI_Cut c(mySbase, myGShape);
323       myShape = c.Shape();
324       UpdateDescendants(c, myShape, Standard_False);
325       //modified by NIZNHY-PKV Thu Mar 21 18:15:47 2002 t
326       Done();
327     }
328     else {
329       myShape = myGShape;
330       Done();
331     }
332   }
333   else {
334     theRevol.Curves(myCurves);
335     myBCurve = theRevol.BarycCurve();
336     GlobalPerform();
337   }
338 }
339
340
341 //=======================================================================
342 //function : Perform
343 //purpose  : feature till shape Until
344 //=======================================================================
345
346 void BRepFeat_MakeRevol::Perform(const TopoDS_Shape& Until)
347 {
348 #ifdef OCCT_DEBUG
349   Standard_Boolean trc = BRepFeat_GettraceFEAT();
350   if (trc) cout << "BRepFeat_MakeRevol::Perform(Until)" << endl;
351 #endif
352   Standard_Real Angle = 0.;
353   Standard_Boolean TourComplet = Standard_False;
354
355   if (Until.IsNull()) {
356     Standard_ConstructionError::Raise();
357   }
358   TopExp_Explorer exp(Until, TopAbs_FACE);
359   if (!exp.More()) {
360     Standard_ConstructionError::Raise();
361   }
362   if (!mySkface.IsNull() && Until.IsSame(mySkface)) {
363     Angle = 2*M_PI;
364     TourComplet = Standard_True;
365   }
366   myGluedF.Clear();
367   myPerfSelection = BRepFeat_SelectionU;
368   PerfSelectionValid();
369   mySFrom.Nullify();
370   ShapeFromValid();
371   mySUntil = Until;
372   Standard_Boolean Trf = TransformShapeFU(1);
373   ShapeUntilValid();
374
375   // Do systematically almost complete revolution
376   // BRepSweep_Revol theRevol(myPbase,myAxis,2.*M_PI-10.*Precision::Angular());
377   LocOpe_Revol theRevol;
378   if(!TourComplet) {
379     Angle = 2.*M_PI- 3*M_PI/180.;
380 #ifdef OCCT_DEBUG
381     if (trc) cout << " No complete Revolution" << endl;
382 #endif
383   }
384   theRevol.Perform(myPbase, myAxis, Angle);
385   TopoDS_Shape VraiRevol = theRevol.Shape();
386   MajMap(myPbase,theRevol,myMap,myFShape,myLShape);
387
388
389   if(!Trf) {
390
391     myGShape = VraiRevol;
392     GeneratedShapeValid();
393
394     TopoDS_Shape Base = theRevol.FirstShape();
395     exp.Init(Base, TopAbs_FACE);
396     TopoDS_Face theBase = TopoDS::Face(exp.Current());
397     exp.Next();
398     if(exp.More()) {
399       NotDone();
400       myStatusError = BRepFeat_InvFirstShape;
401       return;
402     }
403     GluedFacesValid();
404     //VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF);
405
406     theRevol.Curves(myCurves);
407     myBCurve = theRevol.BarycCurve();
408     GlobalPerform();
409   }
410   else {
411     TColGeom_SequenceOfCurve scur;
412     theRevol.Curves(myCurves);
413     myBCurve = theRevol.BarycCurve();    
414     scur.Clear();    
415     scur.Append(myBCurve);
416     LocOpe_CSIntersector ASI(mySUntil);
417     ASI.Perform(scur);
418     if (ASI.IsDone() && ASI.NbPoints(1) >=1) {
419       TopAbs_Orientation Or = ASI.Point(1,1).Orientation();
420       TopoDS_Face FUntil = ASI.Point(1,1).Face();
421       TopoDS_Shape Comp;
422       BRep_Builder B;
423       B.MakeCompound(TopoDS::Compound(Comp));
424       TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, Or);
425       if (!S.IsNull()) B.Add(Comp,S);
426       //modified by NIZNHY-PKV Thu Mar 21 18:17:31 2002 f
427       //BRepAlgo_Cut trP(VraiRevol,Comp);
428       BRepAlgoAPI_Cut trP(VraiRevol,Comp);
429       //modified by NIZNHY-PKV Thu Mar 21 18:17:37 2002 t
430       TopoDS_Shape Cutsh = trP.Shape();
431       TopExp_Explorer ex(Cutsh, TopAbs_SOLID);
432       for(; ex.More(); ex.Next()) {
433         TopExp_Explorer ex1(ex.Current(), TopAbs_FACE);
434         for(; ex1.More(); ex1.Next()) {
435           const TopoDS_Face& fac = TopoDS::Face(ex1.Current());
436           if(fac.IsSame(myPbase)) {
437             VraiRevol = ex.Current();
438             break;
439           }
440         }
441       }
442       if(myFuse == 1) {
443         //modified by NIZNHY-PKV Thu Mar 21 18:17:53 2002 f
444         //BRepAlgo_Fuse f(mySbase, VraiRevol);
445         //myShape = f.Shape();
446         //UpdateDescendants(f.Builder(), myShape, Standard_False);
447         BRepAlgoAPI_Fuse f(mySbase, VraiRevol);
448         myShape = f.Shape();
449         UpdateDescendants(f, myShape, Standard_False);
450         //modified by NIZNHY-PKV Thu Mar 21 18:17:57 2002 t
451         Done();
452       }
453       else if(myFuse == 0) {
454         //modified by NIZNHY-PKV Thu Mar 21 18:18:23 2002 f
455         //BRepAlgo_Cut c(mySbase, VraiRevol);
456         //myShape = c.Shape();
457         //UpdateDescendants(c.Builder(), myShape, Standard_False);
458         BRepAlgoAPI_Cut c(mySbase, VraiRevol);
459         myShape = c.Shape();
460         UpdateDescendants(c, myShape, Standard_False);
461         //modified by NIZNHY-PKV Thu Mar 21 18:18:28 2002 t
462         Done();
463       }
464       else {
465         myShape = VraiRevol;
466         Done(); 
467       }
468     }         
469   }
470   // Loop of control of descendance
471   /*
472   TopExp_Explorer expr(mySbase, TopAbs_FACE);
473   char nom1[20], nom2[20];
474   Standard_Integer ii = 0;
475   for(; expr.More(); expr.Next()) {
476   ii++;
477   sprintf(nom1, "faceinitial_%d", ii);
478   DBRep::Set(nom1, expr.Current());
479   Standard_Integer jj = 0;
480   const TopTools_ListOfShape& list = Modified(expr.Current());
481   TopTools_ListIteratorOfListOfShape ite(list);
482   for(; ite.More(); ite.Next()) {
483   jj++;
484   sprintf(nom2, "facemodifie_%d_%d", ii, jj);
485   DBRep::Set(nom2, ite.Value());
486   }
487   }
488
489   expr.Init(myPbase, TopAbs_EDGE); 
490   ii=0; 
491   for(; expr.More(); expr.Next()) {
492   ii++;
493   sprintf(nom1, "edgeinitial_%d", ii); 
494   DBRep::Set(nom1, expr.Current()); 
495   Standard_Integer jj = 0;
496   const TopTools_ListOfShape& list = Generated(expr.Current());
497   TopTools_ListIteratorOfListOfShape ite(list);
498   for(; ite.More(); ite.Next()) {
499   jj++;
500   sprintf(nom2, "facegeneree_%d_%d", ii, jj);
501   DBRep::Set(nom2, ite.Value());
502   }
503   }
504   */
505 }
506
507
508 //=======================================================================
509 //function : Perform
510 //purpose  : feature limited by two shapes
511 //=======================================================================
512
513 void BRepFeat_MakeRevol::Perform(const TopoDS_Shape& From,
514   const TopoDS_Shape& Until)
515 {
516 #ifdef OCCT_DEBUG
517   Standard_Boolean trc = BRepFeat_GettraceFEAT();
518   if (trc) cout << "BRepFeat_MakeRevol::Perform(From,Until)" << endl;
519 #endif
520   if (From.IsNull() || Until.IsNull()) {
521     Standard_ConstructionError::Raise();
522   }
523   if (!mySkface.IsNull()) {
524     if (From.IsSame(mySkface)) {
525       myJustGluer = Standard_True;
526       Perform(Until);
527       if (myJustGluer) return;
528     }
529     else if (Until.IsSame(mySkface)) {
530       myJustGluer = Standard_True;
531       myAxis.Reverse();
532       Perform(From);
533       if (myJustGluer) return;
534     }
535   }
536
537   myGluedF.Clear();
538   myPerfSelection = BRepFeat_SelectionFU;
539   PerfSelectionValid();
540
541   TopExp_Explorer exp(From, TopAbs_FACE);
542   if (!exp.More()) {
543     Standard_ConstructionError::Raise();
544   }
545   exp.Init(Until, TopAbs_FACE);
546   if (!exp.More()) {
547     Standard_ConstructionError::Raise();
548   }
549
550   mySFrom = From;
551   Standard_Boolean Trff = TransformShapeFU(0);
552   ShapeFromValid();
553   mySUntil = Until;
554   Standard_Boolean Trfu = TransformShapeFU(1);
555   ShapeUntilValid();  
556
557   if(Trfu != Trff) {
558     NotDone();
559     myStatusError = BRepFeat_IncTypes;
560     return;
561   }
562
563   LocOpe_Revol theRevol;
564   theRevol.Perform(myPbase, myAxis, 2*M_PI);
565   TopoDS_Shape VraiRevol = theRevol.Shape();
566
567   MajMap(myPbase,theRevol,myMap,myFShape,myLShape);
568
569   if(!Trff) {    
570     myGShape = VraiRevol;
571     GeneratedShapeValid();
572     GluedFacesValid();
573     //    VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF);
574
575     theRevol.Curves(myCurves);
576     myBCurve = theRevol.BarycCurve();
577     GlobalPerform();
578   }
579   else {
580     theRevol.Curves(myCurves);
581     myBCurve = theRevol.BarycCurve();    
582     TColGeom_SequenceOfCurve scur;
583     scur.Clear();    
584     scur.Append(myBCurve);
585     LocOpe_CSIntersector ASI1(mySUntil);
586     LocOpe_CSIntersector ASI2(mySFrom);
587     ASI1.Perform(scur);
588     ASI2.Perform(scur);
589     TopAbs_Orientation OrU, OrF;
590     TopoDS_Face FFrom, FUntil;
591     Standard_Real PrF, PrU;
592     if (ASI1.IsDone() && ASI1.NbPoints(1) >=1) {
593       OrU = ASI1.Point(1,1).Orientation();
594       FUntil = ASI1.Point(1,1).Face();
595       PrU = ASI1.Point(1,1).Parameter();
596     }
597     else {
598       NotDone();
599       myStatusError = BRepFeat_NoIntersectU;
600       return;
601     }
602     if (ASI2.IsDone() && ASI2.NbPoints(1) >=1) {
603       Standard_Real pr1 = ASI2.Point(1,1).Parameter();
604       pr1 = ElCLib::InPeriod(pr1,PrU-2*M_PI,PrU);
605       Standard_Real pr2 = ASI2.Point(1,ASI2.NbPoints(1)).Parameter();
606       pr2 = ElCLib::InPeriod(pr2,PrU-2*M_PI,PrU);
607       //OrF = OrU;
608       OrF = TopAbs::Reverse(OrU);
609       FFrom = ASI2.Point(1,1).Face();
610       PrF = Max(pr1, pr2);
611     }
612     else {
613       NotDone();
614       myStatusError = BRepFeat_NoIntersectF;
615       return;
616     }
617     if(!(PrU > PrF)) {
618       NotDone();
619       myStatusError = BRepFeat_IncParameter;
620       return;
621     }
622     TopoDS_Shape Comp;
623     BRep_Builder B;
624     B.MakeCompound(TopoDS::Compound(Comp));
625     TopoDS_Solid SF = BRepFeat::Tool(mySFrom, FFrom, OrF);
626     if (!SF.IsNull()) B.Add(Comp,SF);
627     TopoDS_Solid SU = BRepFeat::Tool(mySUntil, FUntil, OrU);
628     if (!SU.IsNull()) B.Add(Comp,SU);
629     //modified by NIZNHY-PKV Thu Mar 21 18:18:54 2002 f
630     //BRepAlgo_Cut trP(VraiRevol,Comp);
631     BRepAlgoAPI_Cut trP(VraiRevol,Comp);
632     //modified by NIZNHY-PKV Thu Mar 21 18:18:57 2002 t
633     TopoDS_Shape Cutsh = trP.Shape();
634     TopExp_Explorer ex(Cutsh, TopAbs_SOLID);
635     //    Standard_Real PrF = BRepFeat::ParametricBarycenter(mySFrom, myBCurve);
636     //    Standard_Real PrU = BRepFeat::ParametricBarycenter(mySUntil, myBCurve);
637     VraiRevol = ex.Current();
638     for(; ex.More(); ex.Next()) {
639       Standard_Real PrCur = BRepFeat::
640         ParametricBarycenter(ex.Current(), myBCurve);
641       if(PrF <= PrCur && PrU >= PrCur) {
642         VraiRevol = ex.Current();
643         break;
644       }
645     }
646     if(myFuse == 1 && !myJustFeat) {
647       //modified by NIZNHY-PKV Thu Mar 21 18:19:14 2002 f
648       //BRepAlgo_Fuse f(mySbase, VraiRevol);
649       //myShape = f.Shape();
650       //UpdateDescendants(f.Builder(), myShape, Standard_False);
651       BRepAlgoAPI_Fuse f(mySbase, VraiRevol);
652       myShape = f.Shape();
653       UpdateDescendants(f, myShape, Standard_False);
654       //modified by NIZNHY-PKV Thu Mar 21 18:19:18 2002 t
655       Done();
656     }
657     else if(myFuse == 0 && !myJustFeat) {
658       //modified by NIZNHY-PKV Thu Mar 21 18:19:46 2002 f
659       //BRepAlgo_Cut c(mySbase, VraiRevol);
660       //myShape = c.Shape();
661       //UpdateDescendants(c.Builder(), myShape, Standard_False);
662       BRepAlgoAPI_Cut c(mySbase, VraiRevol);
663       myShape = c.Shape();
664       UpdateDescendants(c, myShape, Standard_False);
665       //modified by NIZNHY-PKV Thu Mar 21 18:19:50 2002 t
666       Done();
667     }
668     else {
669       myShape = VraiRevol;
670       Done();   
671     }
672   }
673 }
674
675
676 //=======================================================================
677 //function : PerformThruAll
678 //purpose  : feature throughout the initial shape
679 //=======================================================================
680
681 void BRepFeat_MakeRevol::PerformThruAll()
682 {
683 #ifdef OCCT_DEBUG
684   Standard_Boolean trc = BRepFeat_GettraceFEAT();
685   if (trc) cout << "BRepFeat_MakeRevol::PerformThruAll()" << endl;
686 #endif
687   Perform(2.*M_PI);
688 }
689
690 //=======================================================================
691 //function : PerformUntilAngle
692 //purpose  : feature till shape Until defined with the angle
693 //=======================================================================
694
695 void BRepFeat_MakeRevol::PerformUntilAngle(const TopoDS_Shape& Until,
696   const Standard_Real Angle)
697 {
698 #ifdef OCCT_DEBUG
699   Standard_Boolean trc = BRepFeat_GettraceFEAT();
700   if (trc) cout << "BRepFeat_MakeRevol::PerformUntilAngle(Until,Angle)" << endl;
701 #endif
702   if (Until.IsNull()) {
703     Perform(Angle);
704   }
705   if(Angle == 0) {
706     Perform(Until);
707   }
708   TopExp_Explorer exp(Until, TopAbs_FACE);
709   if (!exp.More()) {
710     Standard_ConstructionError::Raise();
711   }
712   if (!mySkface.IsNull() && Until.IsSame(mySkface)) {
713     Perform(Angle);
714     return;
715   }
716   myGluedF.Clear();
717   myPerfSelection = BRepFeat_NoSelection;
718   PerfSelectionValid();
719   mySFrom.Nullify();
720   ShapeFromValid();
721   mySUntil = Until;
722   Standard_Boolean Trf = TransformShapeFU(1);
723   ShapeUntilValid();
724
725   // Produce systematicallt an almost complete revolution
726   //  BRepSweep_Revol theRevol(myPbase,myAxis,2.*M_PI-10.*Precision::Angular());
727   LocOpe_Revol theRevol;
728   theRevol.Perform(myPbase, myAxis, Angle);
729   TopoDS_Shape VraiRevol = theRevol.Shape();
730
731   MajMap(myPbase,theRevol,myMap,myFShape,myLShape);
732
733   if(Trf) {
734     myGShape = VraiRevol;
735     GeneratedShapeValid();
736
737     TopoDS_Shape Base = theRevol.FirstShape();
738     exp.Init(Base, TopAbs_FACE);
739     TopoDS_Face theBase = TopoDS::Face(exp.Current());
740     exp.Next();
741     if(exp.More()) {
742       NotDone();
743       myStatusError = BRepFeat_InvFirstShape;
744       return;
745     }
746     GluedFacesValid();
747     //VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF);
748
749
750     theRevol.Curves(myCurves);
751     myBCurve = theRevol.BarycCurve();
752     GlobalPerform();
753   }
754   else {
755     TColGeom_SequenceOfCurve scur;
756     theRevol.Curves(myCurves);
757     myBCurve = theRevol.BarycCurve();    
758     scur.Clear();    
759     scur.Append(myBCurve);
760     LocOpe_CSIntersector ASI(mySUntil);
761     ASI.Perform(scur);
762     if (ASI.IsDone() && ASI.NbPoints(1) >=1) {
763       TopAbs_Orientation Or = ASI.Point(1,1).Orientation();
764       TopoDS_Face FUntil = ASI.Point(1,1).Face();
765       TopoDS_Shape Comp;
766       BRep_Builder B;
767       B.MakeCompound(TopoDS::Compound(Comp));
768       TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, Or);
769       if (!S.IsNull()) B.Add(Comp,S);
770       //modified by NIZNHY-PKV Thu Mar 21 18:20:14 2002 f
771       //BRepAlgo_Cut trP(VraiRevol,Comp);
772       BRepAlgoAPI_Cut trP(VraiRevol,Comp);
773       //modified by NIZNHY-PKV Thu Mar 21 18:20:19 2002 t
774       TopoDS_Shape Cutsh = trP.Shape();
775       TopExp_Explorer ex(Cutsh, TopAbs_SOLID);
776       for(; ex.More(); ex.Next()) {
777         TopExp_Explorer ex1(ex.Current(), TopAbs_FACE);
778         for(; ex1.More(); ex1.Next()) {
779           const TopoDS_Face& fac = TopoDS::Face(ex1.Current());
780           if(fac.IsSame(myPbase)) {
781             VraiRevol = ex.Current();
782             break;
783           }
784         }
785       }
786       if(myFuse == 1) {
787         //modified by NIZNHY-PKV Thu Mar 21 18:20:36 2002 f
788         //BRepAlgo_Fuse f(mySbase, VraiRevol);
789         //myShape = f.Shape();
790         //UpdateDescendants(f.Builder(), myShape, Standard_False);
791         BRepAlgoAPI_Fuse f(mySbase, VraiRevol);
792         myShape = f.Shape();
793         UpdateDescendants(f, myShape, Standard_False);
794         //modified by NIZNHY-PKV Thu Mar 21 18:20:40 2002 t
795         Done();
796       }
797       else if(myFuse == 0) {
798         //modified by NIZNHY-PKV Thu Mar 21 18:21:07 2002 f
799         //BRepAlgo_Cut c(mySbase, VraiRevol);
800         //myShape = c.Shape();
801         //UpdateDescendants(c.Builder(), myShape, Standard_False);
802         BRepAlgoAPI_Cut c(mySbase, VraiRevol);
803         myShape = c.Shape();
804         UpdateDescendants(c, myShape, Standard_False);
805         //modified by NIZNHY-PKV Thu Mar 21 18:21:26 2002 t
806         Done();
807       }
808       else {
809         myShape = VraiRevol;
810         Done(); 
811       }
812     }         
813   }
814 }
815 //=======================================================================
816 //function : Curves
817 //purpose  : circles parallel to the generating edge of revolution
818 //=======================================================================
819
820 void BRepFeat_MakeRevol::Curves(TColGeom_SequenceOfCurve& scur)
821 {
822   scur = myCurves;
823 }
824
825 //=======================================================================
826 //function : BarycCurve
827 //purpose  : pass through the center of mass of the primitive
828 //=======================================================================
829
830 Handle(Geom_Curve) BRepFeat_MakeRevol::BarycCurve()
831 {
832   return myBCurve;
833 }
834
835 //=======================================================================
836 //function : VerifGluedFaces
837 //purpose  : Check intersection Tool/theSkface = thePbase
838 //           if yes -> OK otherwise -> case without gluing
839 //=======================================================================
840
841 static void VerifGluedFaces(const TopoDS_Face& theSkface,
842   const TopoDS_Shape& thePbase,
843   Handle(Geom_Curve)& theBCurve,
844   TColGeom_SequenceOfCurve& theCurves,
845   LocOpe_Revol& theRevol,
846   TopTools_DataMapOfShapeShape& theMap)
847 {
848   Standard_Boolean GluedFaces = Standard_True;
849   TopoDS_Shape VraiRevol = theRevol.Shape();
850
851   TColGeom_SequenceOfCurve scur;
852   theRevol.Curves(theCurves);
853   theBCurve = theRevol.BarycCurve();    
854   scur.Clear();    
855   scur.Append(theBCurve);
856   LocOpe_CSIntersector ASI(theSkface);
857   ASI.Perform(scur);
858   if (ASI.IsDone() && ASI.NbPoints(1) >=1) {
859     TopAbs_Orientation Or = ASI.Point(1,1).Orientation();
860     TopoDS_Face FSk = ASI.Point(1,1).Face();
861     TopoDS_Shape Comp;
862     BRep_Builder B;
863     B.MakeCompound(TopoDS::Compound(Comp));
864     TopoDS_Solid S = BRepFeat::Tool(theSkface, FSk, Or);
865     if (!S.IsNull()) B.Add(Comp,S);
866     //modified by NIZNHY-PKV Thu Mar 21 18:21:54 2002 f
867     //BRepAlgo_Cut trP(VraiRevol,Comp);
868     BRepAlgoAPI_Cut trP(VraiRevol,Comp);
869     //modified by NIZNHY-PKV Thu Mar 21 18:21:58 2002 t
870     TopoDS_Shape Cutsh = trP.Shape();
871     TopExp_Explorer ex(Cutsh, TopAbs_SOLID);
872     for(; ex.More(); ex.Next()) {
873       TopExp_Explorer ex1(ex.Current(), TopAbs_FACE);
874       for(; ex1.More(); ex1.Next()) {
875         const TopoDS_Face& fac1 = TopoDS::Face(ex1.Current());
876         TopExp_Explorer ex2(thePbase, TopAbs_FACE);
877         for(; ex2.More(); ex2.Next()) {
878           const TopoDS_Face& fac2 = TopoDS::Face(ex2.Current());
879           if(fac1.IsSame(fac2)) break;
880         }
881         if (ex2.More()) break;
882       }
883       if (ex1.More()) continue;
884       GluedFaces = Standard_False;
885       break;
886     }
887     if (!GluedFaces) {
888 #ifdef OCCT_DEBUG
889       Standard_Boolean trc = BRepFeat_GettraceFEAT();
890       if (trc) cout << " Intersection Revol/skface : no gluing" << endl;
891 #endif
892       theMap.Clear();
893     }
894   }
895 }
896
897 //=======================================================================
898 //function : MajMap
899 //purpose  : management of descendants
900 //=======================================================================
901
902 static void MajMap(const TopoDS_Shape& theB,
903   const LocOpe_Revol& theP,
904   TopTools_DataMapOfShapeListOfShape& theMap, // myMap
905   TopoDS_Shape& theFShape,  // myFShape
906   TopoDS_Shape& theLShape) // myLShape
907 {
908   TopExp_Explorer exp(theP.FirstShape(),TopAbs_WIRE);
909   if (exp.More()) {
910     theFShape = exp.Current();
911     TopTools_ListOfShape thelist;
912     theMap.Bind(theFShape, thelist);
913     for (exp.Init(theP.FirstShape(),TopAbs_FACE);exp.More();exp.Next()) {
914       theMap(theFShape).Append(exp.Current());
915     }
916   }
917
918   exp.Init(theP.LastShape(),TopAbs_WIRE);
919   if (exp.More()) {
920     theLShape = exp.Current();
921     TopTools_ListOfShape thelist1;
922     theMap.Bind(theLShape, thelist1);
923     for (exp.Init(theP.LastShape(),TopAbs_FACE);exp.More();exp.Next()) {
924       theMap(theLShape).Append(exp.Current());
925     }
926   }
927
928   for (exp.Init(theB,TopAbs_EDGE); exp.More(); exp.Next()) {
929     if (!theMap.IsBound(exp.Current())) {
930       TopTools_ListOfShape thelist2;
931       theMap.Bind(exp.Current(), thelist2);
932       theMap(exp.Current()) = theP.Shapes(exp.Current());
933     }
934   }
935 }
936
937
938
939 //=======================================================================
940 //function : ToFuse
941 //purpose  : two faces samedomaine or not
942 //=======================================================================
943
944 Standard_Boolean ToFuse(const TopoDS_Face& F1,
945   const TopoDS_Face& F2)
946 {
947   if (F1.IsNull() || F2.IsNull()) {
948     return Standard_False;
949   }
950
951   Handle(Geom_Surface) S1,S2;
952   TopLoc_Location loc1, loc2;
953   Handle(Standard_Type) typS1,typS2;
954   const Standard_Real tollin = Precision::Confusion();
955   const Standard_Real tolang = Precision::Angular();
956
957   S1 = BRep_Tool::Surface(F1,loc1);
958   S2 = BRep_Tool::Surface(F2,loc2);
959
960   typS1 = S1->DynamicType();
961   typS2 = S2->DynamicType();
962
963   if (typS1 == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) {
964     S1 =  Handle(Geom_RectangularTrimmedSurface)::DownCast (S1)->BasisSurface();
965     typS1 = S1->DynamicType();
966   }
967
968   if (typS2 == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) {
969     S2 =  Handle(Geom_RectangularTrimmedSurface)::DownCast (S2)->BasisSurface();
970     typS2 = S2->DynamicType();
971   }
972
973   if (typS1 != typS2) {
974     return Standard_False;
975   }
976
977
978   Standard_Boolean ValRet = Standard_False;
979   if (typS1 == STANDARD_TYPE(Geom_Plane)) {
980     S1 = BRep_Tool::Surface(F1);  // to apply the location.
981     S2 = BRep_Tool::Surface(F2);
982     gp_Pln pl1( Handle(Geom_Plane)::DownCast (S1)->Pln());
983     gp_Pln pl2( Handle(Geom_Plane)::DownCast (S2)->Pln());
984
985     if (pl1.Position().IsCoplanar(pl2.Position(),tollin,tolang)) {
986       ValRet = Standard_True;
987     }
988   }
989
990   return ValRet;
991 }
992
993
994
995
996
997
998
999
1000