0026585: Eliminate compile warnings obtained by building occt with vc14: 'type cast...
[occt.git] / src / BRepFeat / BRepFeat_MakeDPrism.cxx
1 // Created on: 1996-09-03
2 // Created by: Olga KOULECHOVA
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
18 #include <Bnd_Box.hxx>
19 #include <BRep_Builder.hxx>
20 #include <BRep_Tool.hxx>
21 #include <BRepAlgo.hxx>
22 #include <BRepAlgoAPI_Cut.hxx>
23 #include <BRepAlgoAPI_Fuse.hxx>
24 #include <BRepBndLib.hxx>
25 #include <BRepFeat.hxx>
26 #include <BRepFeat_MakeDPrism.hxx>
27 #include <BRepLib.hxx>
28 #include <BRepLib_MakeFace.hxx>
29 #include <BRepPrimAPI_MakeBox.hxx>
30 #include <BRepTools.hxx>
31 #include <Geom2d_Curve.hxx>
32 #include <Geom_Curve.hxx>
33 #include <Geom_Line.hxx>
34 #include <Geom_Plane.hxx>
35 #include <Geom_RectangularTrimmedSurface.hxx>
36 #include <Geom_Surface.hxx>
37 #include <gp_Pln.hxx>
38 #include <gp_Pnt.hxx>
39 #include <gp_Pnt2d.hxx>
40 #include <gp_Vec.hxx>
41 #include <LocOpe.hxx>
42 #include <LocOpe_CSIntersector.hxx>
43 #include <LocOpe_DPrism.hxx>
44 #include <LocOpe_PntFace.hxx>
45 #include <LocOpe_SequenceOfLin.hxx>
46 #include <Precision.hxx>
47 #include <Standard_ConstructionError.hxx>
48 #include <TColGeom_SequenceOfCurve.hxx>
49 #include <TColgp_SequenceOfPnt.hxx>
50 #include <TopAbs.hxx>
51 #include <TopExp.hxx>
52 #include <TopExp_Explorer.hxx>
53 #include <TopoDS.hxx>
54 #include <TopoDS_Compound.hxx>
55 #include <TopoDS_Edge.hxx>
56 #include <TopoDS_Face.hxx>
57 #include <TopoDS_Shape.hxx>
58 #include <TopoDS_Shell.hxx>
59 #include <TopoDS_Solid.hxx>
60 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
61 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
62 #include <TopTools_ListIteratorOfListOfShape.hxx>
63 #include <TopTools_MapIteratorOfMapOfShape.hxx>
64 #include <TopTools_MapOfShape.hxx>
65
66 //modified by NIZNHY-PKV Thu Mar 21 17:49:46 2002 f
67 //#include <BRepAlgo_Cut.hxx>
68 //#include <BRepAlgo_Fuse.hxx>
69 //modified by NIZNHY-PKV Thu Mar 21 17:50:04 2002 t
70 #ifdef OCCT_DEBUG
71 extern Standard_Boolean BRepFeat_GettraceFEAT();
72 #endif
73
74 static void MajMap(const TopoDS_Shape&,
75   const LocOpe_DPrism&,
76   TopTools_DataMapOfShapeListOfShape&, // myMap
77   TopoDS_Shape&,  // myFShape
78   TopoDS_Shape&); // myLShape
79
80 static Standard_Real HeightMax(const TopoDS_Shape& theSbase,
81   const TopoDS_Face& theSkface,
82   const TopoDS_Shape& theSFrom,
83   const TopoDS_Shape& theSUntil);
84
85 static Standard_Integer SensOfPrism(const Handle(Geom_Curve) C,
86   const TopoDS_Shape& Until);
87
88 static Handle(Geom_Curve) TestCurve(const TopoDS_Face&);
89
90
91 //=======================================================================
92 //function : Init
93 //purpose  : 
94 //=======================================================================
95
96 void BRepFeat_MakeDPrism::Init(const TopoDS_Shape&    Sbase,
97   const TopoDS_Face&     Pbase,
98   const TopoDS_Face&     Skface,
99   const Standard_Real    Angle,
100   const Standard_Integer Mode,
101   const Standard_Boolean Modify)
102
103 {
104 #ifdef OCCT_DEBUG
105   Standard_Boolean trc = BRepFeat_GettraceFEAT();
106   if (trc) cout << "BRepFeat_MakeDPrism::Init" << endl;
107 #endif
108   mySkface = Skface;
109   SketchFaceValid();
110   mySbase  = Sbase;
111   BasisShapeValid();
112   myPbase  = Pbase;
113   mySlface.Clear();
114   if(Mode == 0) {
115     myFuse   = Standard_False;
116     myJustFeat = Standard_False;
117   }
118   else if(Mode == 1) {
119     myFuse   = Standard_True;
120     myJustFeat = Standard_False;
121   }
122   else if(Mode == 2) {
123     myFuse   = Standard_True;
124     myJustFeat = Standard_True;
125   }
126   else {    
127   }
128   myModify = Modify;
129   myJustGluer = Standard_False;
130
131
132   //-------------- ifv
133   //mySkface.Nullify();
134   //-------------- ifv
135
136   myShape.Nullify();
137   myMap.Clear();
138   myFShape.Nullify();
139   myLShape.Nullify();
140   myTopEdges.Clear();
141   myLatEdges.Clear();
142   TopExp_Explorer exp;
143   for (exp.Init(mySbase,TopAbs_FACE);exp.More();exp.Next()) {
144     TopTools_ListOfShape thelist;
145     myMap.Bind(exp.Current(), thelist);
146     myMap(exp.Current()).Append(exp.Current());
147   }
148   myAngle = Angle;
149 #ifdef OCCT_DEBUG
150   if (trc) {
151     if (myJustFeat)  cout << " Just Feature" << endl;
152     if (myFuse)  cout << " Fuse" << endl;
153     if (!myFuse)  cout << " Cut" << endl;
154     if (!myModify) cout << " Modify = 0" << endl;
155     //    cout <<" Angle = " << myAngle << endl;
156   }
157 #endif 
158 }
159
160
161 //=======================================================================
162 //function : Add
163 //purpose  : add sliding faces and edges
164 //=======================================================================
165
166 void BRepFeat_MakeDPrism::Add(const TopoDS_Edge& E,
167   const TopoDS_Face& F)
168 {
169 #ifdef OCCT_DEBUG
170   Standard_Boolean trc = BRepFeat_GettraceFEAT();
171   if (trc) cout << "BRepFeat_MakeDPrism::Add(Edge,face)" << endl;
172 #endif
173   TopExp_Explorer exp;
174   for (exp.Init(mySbase,TopAbs_FACE);exp.More();exp.Next()) {
175     if (exp.Current().IsSame(F)) {
176       break;
177     }
178   }
179   if (!exp.More()) {
180     Standard_ConstructionError::Raise();
181   }
182
183   for (exp.Init(myPbase,TopAbs_EDGE);exp.More();exp.Next()) {
184     if (exp.Current().IsSame(E)) {
185       break;
186     }
187   }
188   if (!exp.More()) {
189     Standard_ConstructionError::Raise();
190   }
191
192   if (!mySlface.IsBound(F)) {
193     TopTools_ListOfShape thelist;
194     mySlface.Bind(F, thelist);
195   }
196   TopTools_ListIteratorOfListOfShape itl(mySlface(F));
197   for (; itl.More();itl.Next()) {
198     if (itl.Value().IsSame(E)) {
199       break;
200     }
201   }
202   if (!itl.More()) {
203     mySlface(F).Append(E);
204   }
205 }
206
207
208 //=======================================================================
209 //function : Perform
210 //purpose  : feature of Height
211 //=======================================================================
212
213 void BRepFeat_MakeDPrism::Perform(const Standard_Real Height)
214 {
215 #ifdef OCCT_DEBUG
216   Standard_Boolean trc = BRepFeat_GettraceFEAT();
217   if (trc) cout << "BRepFeat_MakeDPrism::Perform(Height)" << endl;
218 #endif
219   mySFrom.Nullify();
220   ShapeFromValid();
221   mySUntil.Nullify();
222   ShapeUntilValid();
223   myGluedF.Clear();
224   myPerfSelection = BRepFeat_NoSelection;
225   PerfSelectionValid();
226
227   Standard_Real theheight = Height/cos(myAngle);
228   //  myPbase.Orientation(TopAbs_FORWARD);
229
230   LocOpe_DPrism theDPrism(myPbase,theheight,myAngle);
231   TopoDS_Shape VraiDPrism = theDPrism.Shape();
232
233   MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
234
235   myGShape = VraiDPrism;
236   GeneratedShapeValid();  
237   TopoDS_Shape Base = theDPrism.FirstShape();
238   TopExp_Explorer exp(Base, TopAbs_FACE);
239   TopoDS_Face theBase = TopoDS::Face(exp.Current());
240   exp.Next();
241   if(exp.More()) {
242     NotDone();
243     myStatusError = BRepFeat_InvFirstShape;
244     return;
245   }
246
247   // management of gluing faces 
248
249   GluedFacesValid();
250
251   if(!myGluedF.IsEmpty()) {   // case gluing
252     myJustGluer = Standard_True;
253     theDPrism.Curves(myCurves);
254     myBCurve = theDPrism.BarycCurve();    
255     GlobalPerform();
256   }
257
258   // if there is no gluing -> call topological operations
259   if(!myJustGluer) {
260     if(myFuse == 1) {
261       //modified by NIZNHY-PKV Thu Mar 21 17:32:17 2002 f
262       //BRepAlgo_Fuse f(mySbase, myGShape);
263       //myShape = f.Shape();
264       //UpdateDescendants(f.Builder(), myShape, Standard_False);
265       BRepAlgoAPI_Fuse f(mySbase, myGShape);
266       myShape = f.Shape();
267       UpdateDescendants(f, myShape, Standard_False);
268       //modified by NIZNHY-PKV Thu Mar 21 17:32:23 2002 t
269       Done();
270     }
271     else if(myFuse == 0) {
272       //modified by NIZNHY-PKV Thu Mar 21 17:32:57 2002 f
273       //BRepAlgo_Cut c(mySbase, myGShape);
274       //myShape = c.Shape();
275       //UpdateDescendants(c.Builder(), myShape, Standard_False);
276       BRepAlgoAPI_Cut c(mySbase, myGShape);
277       myShape = c.Shape();
278       UpdateDescendants(c, myShape, Standard_False);
279       //modified by NIZNHY-PKV Thu Mar 21 17:33:00 2002 t
280       Done();
281     }
282     else {
283       myShape = myGShape;
284       Done();
285     }
286   }
287 }
288
289 //=======================================================================
290 //function : Perform
291 //purpose  : feature limited by the shape Until
292 //=======================================================================
293
294 void BRepFeat_MakeDPrism::Perform(const TopoDS_Shape& Until)
295 {
296 #ifdef OCCT_DEBUG
297   Standard_Boolean trc = BRepFeat_GettraceFEAT();
298   if (trc) cout << "BRepFeat_MakeDPrism::Perform(Until)" << endl;
299 #endif
300   if (Until.IsNull()) {
301     Standard_ConstructionError::Raise();
302   }
303   TopExp_Explorer exp(Until, TopAbs_FACE);
304   if (!exp.More()) {
305     Standard_ConstructionError::Raise();
306   }
307   //  myPbase.Orientation(TopAbs_FORWARD);
308
309   myGluedF.Clear();
310   myPerfSelection = BRepFeat_SelectionU;
311   PerfSelectionValid();
312   mySFrom.Nullify();
313   ShapeFromValid();
314   mySUntil = Until;
315   Standard_Boolean Trf = TransformShapeFU(1);  
316   ShapeUntilValid();
317   Handle(Geom_Curve) C = TestCurve(myPbase);  
318   Standard_Integer sens = SensOfPrism(C, mySUntil);
319
320   BRep_Builder bB;
321   Standard_Real Height = 
322     sens*HeightMax(mySbase, mySkface, mySFrom, mySUntil);
323   LocOpe_DPrism theDPrism(myPbase,Height,myAngle);
324   TopoDS_Shape VraiDPrism = theDPrism.Shape();
325   if(!Trf) {
326     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
327     myGShape = VraiDPrism;
328     GeneratedShapeValid();
329     TopoDS_Shape Base = theDPrism.FirstShape();
330     exp.Init(Base, TopAbs_FACE);
331     TopoDS_Face theBase = TopoDS::Face(exp.Current());
332     exp.Next();
333     if(exp.More()) {
334       NotDone();
335       myStatusError = BRepFeat_InvFirstShape;
336       return;
337     }
338
339     GluedFacesValid();
340     theDPrism.Curves(myCurves);
341     myBCurve = theDPrism.BarycCurve();
342     GlobalPerform();
343   }
344   else {
345     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
346     Handle(Geom_Curve) C1;
347     if(sens == -1) {
348       C1 = C->Reversed();
349     } 
350     else {
351       C1 = C;
352     }
353
354     TColGeom_SequenceOfCurve scur;
355     scur.Clear();    
356     scur.Append(C1);
357     LocOpe_CSIntersector ASI(mySUntil);
358     ASI.Perform(scur);
359     TopAbs_Orientation Or;
360     if (ASI.IsDone() && ASI.NbPoints(1) >=1) { 
361       if (myFuse == 1) {
362         Or = ASI.Point(1,1).Orientation();
363       }
364       else {
365         Or = ASI.Point(1,ASI.NbPoints(1)).Orientation();
366       }     
367       //      Standard_Real prm = ASI.Point(1,1).Parameter();
368       //      if(prm < 0) Or = TopAbs::Reverse(Or);
369       TopoDS_Face FUntil = ASI.Point(1,1).Face();
370       TopoDS_Shape Comp;
371       bB.MakeCompound(TopoDS::Compound(Comp));
372       TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, Or);
373       if (!S.IsNull()) bB.Add(Comp,S);
374
375       //modified by NIZNHY-PKV Thu Mar 21 17:33:27 2002
376       //BRepAlgo_Cut trP(VraiDPrism,Comp);
377       //UpdateDescendants(trP.Builder(),trP.Shape(), Standard_False);
378       BRepAlgoAPI_Cut trP(VraiDPrism,Comp);
379       UpdateDescendants(trP, trP.Shape(), Standard_False);
380       //modified by NIZNHY-PKV Thu Mar 21 17:33:32 2002 t
381
382       TopExp_Explorer ex(trP.Shape(), TopAbs_SOLID);
383       TopoDS_Shape Cutsh = ex.Current();
384       if(myFuse == 1) {
385         //modified by NIZNHY-PKV Thu Mar 21 17:34:09 2002 f
386         //BRepAlgo_Fuse f(mySbase, Cutsh);
387         //myShape = f.Shape();
388         //UpdateDescendants(f.Builder(), myShape, Standard_False);
389         BRepAlgoAPI_Fuse f(mySbase, Cutsh);
390         myShape = f.Shape();
391         UpdateDescendants(f, myShape, Standard_False);
392         //modified by NIZNHY-PKV Thu Mar 21 17:34:13 2002 t
393         Done();
394       }
395       else if(myFuse == 0) {
396         //modified by NIZNHY-PKV Thu Mar 21 17:34:56 2002 f
397         //BRepAlgo_Cut c(mySbase, Cutsh);
398         //myShape = c.Shape();
399         //UpdateDescendants(c.Builder(), myShape, Standard_False);
400         BRepAlgoAPI_Cut c(mySbase, Cutsh);
401         myShape = c.Shape();
402         UpdateDescendants(c, myShape, Standard_False);
403         //modified by NIZNHY-PKV Thu Mar 21 17:34:59 2002 t
404         Done();
405       }
406       else {
407         myShape = Cutsh;
408         Done(); 
409       }
410     }         
411   }
412   TopTools_ListIteratorOfListOfShape ited(myNewEdges);
413   for (; ited.More();ited.Next()) {
414     const TopoDS_Edge& ledg=TopoDS::Edge(ited.Value());
415     if (!BRepAlgo::IsValid(ledg)) {
416       bB.SameRange(ledg, Standard_False);
417       bB.SameParameter(ledg, Standard_False);
418       BRepLib::SameParameter(ledg, BRep_Tool::Tolerance(ledg));
419     }
420   }  
421 }
422
423
424 //=======================================================================
425 //function : Perform
426 //purpose  : feature limited by two shapes 
427 //=======================================================================
428
429 void BRepFeat_MakeDPrism::Perform(const TopoDS_Shape& From,
430   const TopoDS_Shape& Until)
431 {
432 #ifdef OCCT_DEBUG
433   Standard_Boolean trc = BRepFeat_GettraceFEAT();
434   if (trc) cout << "BRepFeat_MakeDPrism::Perform(From,Until)" << endl;
435 #endif
436   if (From.IsNull() || Until.IsNull()) {
437     Standard_ConstructionError::Raise();
438   }
439
440   if (!mySkface.IsNull()) {
441     if (From.IsSame(mySkface)) {
442       myJustGluer = Standard_True;
443       Perform(Until);
444       if (myJustGluer) return;
445     }
446     else if (Until.IsSame(mySkface)) {
447       myJustGluer = Standard_True;
448       Perform(From);
449       if (myJustGluer) return;
450     }
451   }
452   //  myPbase.Orientation(TopAbs_FORWARD);
453
454   myGluedF.Clear();
455   myPerfSelection = BRepFeat_SelectionFU;
456   PerfSelectionValid();
457
458   TopExp_Explorer exp(From, TopAbs_FACE);
459   if (!exp.More()) {
460     Standard_ConstructionError::Raise();
461   }
462   exp.Init(Until, TopAbs_FACE);
463   if (!exp.More()) {
464     Standard_ConstructionError::Raise();
465   }
466   mySFrom = From;
467   Standard_Boolean Trff = TransformShapeFU(0);
468   ShapeFromValid();
469   mySUntil = Until;
470   Standard_Boolean Trfu = TransformShapeFU(1);
471   ShapeUntilValid();  
472   if(Trfu != Trff) {
473     NotDone();
474     myStatusError = BRepFeat_IncTypes;
475     return;
476   }
477   Handle(Geom_Curve) C = TestCurve(myPbase);  
478   Standard_Integer sens;
479   if(From.IsSame(Until)) {
480     sens = 1;
481   }
482   else {
483     sens = SensOfPrism(C, mySUntil);
484   }
485
486   Standard_Real Height = 
487     sens*HeightMax(mySbase, myPbase, mySFrom, mySUntil);
488   LocOpe_DPrism theDPrism(myPbase, Height, Height, myAngle);
489   TopoDS_Shape VraiDPrism = theDPrism.Shape();
490
491
492   if(!Trff) {
493     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
494
495     // Make systematically dprism 
496     myGShape = VraiDPrism;
497     GeneratedShapeValid();
498
499     // management of gluing faces 
500     //    mySbase, myPbase, mySlface, theDPrism, myGluedF);
501     GluedFacesValid();
502     theDPrism.Curves(myCurves);
503     myBCurve = theDPrism.BarycCurve();
504
505     // topologic reconstruction 
506     GlobalPerform();
507   }
508   else {
509     // management of descendants
510     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
511     Handle(Geom_Curve) C1;
512     if(sens == -1) {
513       C1 = C->Reversed();
514     }
515     else {
516       C1 = C;
517     }
518     TColGeom_SequenceOfCurve scur;
519     scur.Clear();    
520     scur.Append(C1);
521     LocOpe_CSIntersector ASI1(mySUntil);
522     LocOpe_CSIntersector ASI2(mySFrom);
523     ASI1.Perform(scur);
524     ASI2.Perform(scur);
525     TopAbs_Orientation OrU, OrF;
526     TopoDS_Face FFrom, FUntil;
527     //direction of dprism
528     if (ASI1.IsDone() && ASI1.NbPoints(1) >=1) {
529       if (myFuse == 1) {
530         OrU = ASI1.Point(1,1).Orientation();
531       }
532       else {
533         OrU = ASI1.Point(1,ASI1.NbPoints(1)).Orientation();
534       }
535       //      Standard_Real prm = ASI1.Point(1,1).Parameter();
536       //      if(prm < 0) OrU = TopAbs::Reverse(OrU);
537       FUntil = ASI1.Point(1,1).Face();
538     }
539     else {
540       NotDone();
541       myStatusError = BRepFeat_NoIntersectU;
542       return;
543     }
544     if (ASI2.IsDone() && ASI2.NbPoints(1) >=1) {
545       OrF = ASI2.Point(1,1).Orientation();
546       //      Standard_Real prm = ASI2.Point(1,1).Parameter();
547       OrF = TopAbs::Reverse(OrF);
548       FFrom = ASI2.Point(1,1).Face();
549     }
550     else {
551       NotDone();
552       myStatusError = BRepFeat_NoIntersectF;
553       return;
554     }
555     TopoDS_Shape Comp;
556     BRep_Builder B;
557     B.MakeCompound(TopoDS::Compound(Comp));
558     TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, OrU);
559     if (!S.IsNull()) B.Add(Comp,S);
560     else {
561       NotDone();
562       myStatusError = BRepFeat_NullToolU;
563       return;
564     }
565     TopoDS_Solid SS = BRepFeat::Tool(mySFrom, FFrom, OrF);
566     if (!SS.IsNull()) B.Add(Comp,SS);
567     else {
568       NotDone();
569       myStatusError = BRepFeat_NullToolF;
570       return;
571     }
572
573     //modified by NIZNHY-PKV Thu Mar 21 17:35:48 2002 f
574     //BRepAlgo_Cut trP(VraiDPrism,Comp);
575     BRepAlgoAPI_Cut trP(VraiDPrism,Comp);
576     //modified by NIZNHY-PKV Thu Mar 21 17:35:59 2002 t
577
578     if(myFuse == 1) {
579       //modified by NIZNHY-PKV Thu Mar 21 17:36:06 2002 f
580       //BRepAlgo_Fuse f(mySbase, trP.Shape());
581       //myShape = f.Shape();
582       //UpdateDescendants(f.Builder(), myShape, Standard_False);
583       BRepAlgoAPI_Fuse f(mySbase, trP.Shape());
584       myShape = f.Shape();
585       UpdateDescendants(f, myShape, Standard_False);
586       //modified by NIZNHY-PKV Thu Mar 21 17:36:33 2002 t
587       Done();
588     }
589     //
590     else if(myFuse == 0) {
591       //modified by NIZNHY-PKV Thu Mar 21 17:37:16 2002 f
592       //BRepAlgo_Cut c(mySbase, trP.Shape());
593       //myShape = c.Shape();
594       //UpdateDescendants(c.Builder(), myShape, Standard_False);
595       BRepAlgoAPI_Cut c(mySbase, trP.Shape());
596       myShape = c.Shape();
597       UpdateDescendants(c, myShape, Standard_False);
598       //modified by NIZNHY-PKV Thu Mar 21 17:37:49 2002 t
599       Done();
600     }
601     else {
602       myShape = trP.Shape();
603       Done();   
604     }
605   }
606 }
607
608
609 //=======================================================================
610 //function : Perform
611 //purpose  : feature semi-infinie
612 //=======================================================================
613
614 void BRepFeat_MakeDPrism::PerformUntilEnd()
615 {  
616 #ifdef OCCT_DEBUG
617   Standard_Boolean trc = BRepFeat_GettraceFEAT();
618   if (trc) cout << "BRepFeat_MakeDPrism::PerformUntilEnd()" << endl;
619 #endif 
620   myPerfSelection = BRepFeat_SelectionSh;
621   PerfSelectionValid();
622   myGluedF.Clear();
623   mySUntil.Nullify();
624   ShapeUntilValid();
625   mySFrom.Nullify();
626   ShapeFromValid();
627   Standard_Real Height = HeightMax(mySbase, mySkface, mySFrom, mySUntil);
628   //  myPbase.Orientation(TopAbs_FORWARD);
629
630   LocOpe_DPrism theDPrism(myPbase, Height, myAngle);  
631   TopoDS_Shape VraiDPrism = theDPrism.Shape();
632
633   MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
634
635   myGShape = VraiDPrism;
636   GeneratedShapeValid();
637   GluedFacesValid();
638   theDPrism.Curves(myCurves);
639   myBCurve = theDPrism.BarycCurve();
640   GlobalPerform();
641 }
642
643 //=======================================================================
644 //function : PerformFromEnd
645 //purpose  : feature semiinfinite limited by the shape Until from the other side
646 //=======================================================================
647
648 void BRepFeat_MakeDPrism::PerformFromEnd(const TopoDS_Shape& Until)
649 {
650 #ifdef OCCT_DEBUG
651   Standard_Boolean trc = BRepFeat_GettraceFEAT();
652   if (trc) cout << "BRepFeat_MakeDPrism::PerformFromEnd(From,Until)" << endl;
653 #endif
654   if (Until.IsNull()) {
655     Standard_ConstructionError::Raise();
656   }
657   if (!mySkface.IsNull() && Until.IsSame(mySkface)) {
658     PerformUntilEnd();
659     return;
660   }
661   TopExp_Explorer exp;
662   exp.Init(Until, TopAbs_FACE);
663   if (!exp.More()) {
664     Standard_ConstructionError::Raise();
665   }
666   //  myPbase.Orientation(TopAbs_FORWARD);
667   myPerfSelection = BRepFeat_SelectionShU;
668   PerfSelectionValid();
669   mySFrom.Nullify();
670   ShapeFromValid();
671   mySUntil = Until;
672   Standard_Boolean Trf = TransformShapeFU(1);
673   ShapeUntilValid();
674   Handle(Geom_Curve) C = TestCurve(myPbase);  
675   Standard_Integer sens = SensOfPrism(C, mySUntil);
676   Standard_Real Height = 
677     sens*HeightMax(mySbase, mySkface, mySFrom, mySUntil);
678
679   LocOpe_DPrism theDPrism(myPbase, Height, Height, myAngle);
680   TopoDS_Shape VraiDPrism = theDPrism.Shape();
681   if(VraiDPrism.IsNull()) {
682     NotDone();
683     myStatusError = BRepFeat_NullRealTool;
684     return;
685   }
686
687   if(!Trf) {   // case finite face
688     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);    
689     myGShape = VraiDPrism;
690     GeneratedShapeValid();
691     myGluedF.Clear();
692     GluedFacesValid();
693     theDPrism.Curves(myCurves);
694     myBCurve = theDPrism.BarycCurve();
695     GlobalPerform();
696   }
697   else {   // case support
698     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);    
699     Handle(Geom_Curve) C2;
700     if(sens == -1) {
701       C2 = C->Reversed();
702     }
703     else {
704       C2 = C;
705     }
706     TColGeom_SequenceOfCurve scur;
707     scur.Clear();    
708     scur.Append(C2);
709     LocOpe_CSIntersector ASI1(mySUntil);
710     LocOpe_CSIntersector ASI2(mySbase);
711     ASI1.Perform(scur);
712     ASI2.Perform(scur);
713     TopAbs_Orientation OrU = TopAbs_FORWARD, OrF = TopAbs_FORWARD;
714     TopoDS_Face FUntil, FFrom;
715     if (ASI1.IsDone() && ASI1.NbPoints(1) >=1) {
716       OrU = ASI1.Point(1,1).Orientation();
717       Standard_Real prm = ASI1.Point(1,1).Parameter();
718       if(prm < 0) OrU = TopAbs::Reverse(OrU);
719       FUntil = ASI1.Point(1,1).Face();
720     }
721
722     if (ASI2.IsDone() && ASI2.NbPoints(1) >=1) {
723       Standard_Integer jj = ASI2.NbPoints(1);
724       Standard_Real prm = ASI2.Point(1,1).Parameter();
725       FFrom = ASI2.Point(1, 1).Face();
726       OrF = ASI2.Point(1,1).Orientation();
727       OrF = TopAbs::Reverse(OrF);         
728       for(Standard_Integer iii = 1; iii <= jj; iii++) {
729         if(ASI2.Point(1,iii).Parameter() < prm) {
730           prm = ASI2.Point(1,iii).Parameter();
731           FFrom = ASI2.Point(1, iii).Face(); 
732           OrF = ASI2.Point(1,iii).Orientation();
733           OrF = TopAbs::Reverse(OrF);     
734         }
735       }
736       Handle(Geom_Surface) S = BRep_Tool::Surface(FFrom);
737       if (S->DynamicType() == 
738         STANDARD_TYPE(Geom_RectangularTrimmedSurface)) {
739           S = Handle(Geom_RectangularTrimmedSurface)::
740             DownCast(S)->BasisSurface();
741       }
742       BRepLib_MakeFace fac(S, Precision::Confusion());
743       mySFrom = fac.Face();  
744       Trf = TransformShapeFU(0);
745       //      FFrom = TopoDS::Face(mySFrom);      
746     }
747
748     TopoDS_Shape Comp;
749     BRep_Builder B;
750     B.MakeCompound(TopoDS::Compound(Comp));
751     TopoDS_Solid Sol = BRepFeat::Tool(mySUntil, FUntil, OrU);
752     if (!Sol.IsNull()) B.Add(Comp,Sol);
753     else {
754       NotDone();
755       myStatusError = BRepFeat_NullToolU;
756       return;
757     }
758
759     TopoDS_Solid Sol1 = BRepFeat::Tool(mySFrom, FFrom, OrF);
760     if (!Sol1.IsNull()) B.Add(Comp,Sol1);
761     else {
762       NotDone();
763       myStatusError = BRepFeat_NullToolF;
764       return;
765     }
766
767     //modified by NIZNHY-PKV Thu Mar 21 17:38:33 2002 f
768     //BRepAlgo_Cut trP(VraiDPrism,Comp);
769     BRepAlgoAPI_Cut trP(VraiDPrism,Comp);
770     //modified by NIZNHY-PKV Thu Mar 21 17:38:37 2002 t
771     if(myFuse == 1) {
772       //modified by NIZNHY-PKV Thu Mar 21 17:40:03 2002 f
773       //BRepAlgo_Fuse f(mySbase, trP.Shape());
774       //myShape = f.Shape();
775       //UpdateDescendants(f.Builder(), myShape, Standard_False);
776       BRepAlgoAPI_Fuse f(mySbase, trP.Shape());
777       myShape = f.Shape();
778       UpdateDescendants(f, myShape, Standard_False);
779       //modified by NIZNHY-PKV Thu Mar 21 17:40:07 2002 t
780       Done();
781     }
782     else if(myFuse == 0) {
783       //modified by NIZNHY-PKV Thu Mar 21 17:40:33 2002 t
784       //BRepAlgo_Cut c(mySbase, trP.Shape());
785       //myShape = c.Shape();
786       //UpdateDescendants(c.Builder(), myShape, Standard_False);
787       BRepAlgoAPI_Cut c(mySbase, trP.Shape());
788       myShape = c.Shape();
789       UpdateDescendants(c, myShape, Standard_False);
790       //modified by NIZNHY-PKV Thu Mar 21 17:40:44 2002 t
791       Done();
792     }
793     else {
794       myShape = trP.Shape();
795       Done();   
796     }
797   }
798 }
799
800
801
802 //=======================================================================
803 //function : PerformThruAll
804 //purpose  : feature throughout the entire initial shape 
805 //=======================================================================
806
807 void BRepFeat_MakeDPrism::PerformThruAll()
808 {
809 #ifdef OCCT_DEBUG
810   Standard_Boolean trc = BRepFeat_GettraceFEAT();
811   if (trc) cout << "BRepFeat_MakeDPrism::PerformThruAll()" << endl;
812 #endif
813   mySUntil.Nullify();
814   ShapeUntilValid();
815   mySFrom.Nullify();
816   ShapeFromValid();
817   if(myFuse == 0) {
818     myPerfSelection = BRepFeat_NoSelection;
819   }
820   else {
821     myPerfSelection = BRepFeat_SelectionSh;
822   }
823
824   PerfSelectionValid();
825   //  myPbase.Orientation(TopAbs_FORWARD);
826   myGluedF.Clear();
827   GluedFacesValid();
828
829   Standard_Real Height = HeightMax(mySbase, mySkface, mySFrom, mySUntil);
830   LocOpe_DPrism theDPrism(myPbase, Height, Height, myAngle);
831   TopoDS_Shape VraiDPrism = theDPrism.Shape();
832   MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
833
834   myGShape = VraiDPrism;
835   GeneratedShapeValid();
836
837   if(myFuse == 0) {
838     //modified by NIZNHY-PKV Thu Mar 21 17:46:16 2002 f
839     //BRepAlgo_Cut c(mySbase, myGShape);
840     BRepAlgoAPI_Cut c(mySbase, myGShape);
841     //modified by NIZNHY-PKV Thu Mar 21 17:46:26 2002 t
842     if (c.IsDone()) {
843       myShape = c.Shape();
844       //modified by NIZNHY-PKV Thu Mar 21 17:46:39 2002 f
845       //UpdateDescendants(c.Builder(), myShape, Standard_False);
846       UpdateDescendants(c, myShape, Standard_False);
847       //modified by NIZNHY-PKV Thu Mar 21 17:46:43 2002 t
848       Done();
849     }
850   }
851   else {  
852     theDPrism.Curves(myCurves);
853     myBCurve = theDPrism.BarycCurve();
854     GlobalPerform();
855   }
856 }
857
858
859 //=======================================================================
860 //function : PerformUntilHeight
861 //purpose  : feature until the shape is of the given height
862 //=======================================================================
863
864 void BRepFeat_MakeDPrism::PerformUntilHeight(const TopoDS_Shape& Until,
865   const Standard_Real Height)
866 {
867 #ifdef OCCT_DEBUG
868   Standard_Boolean trc = BRepFeat_GettraceFEAT();
869   if (trc) cout << "BRepFeat_MakeDPrism::PerformUntilHeight(Until,Height)" << endl;
870 #endif
871   if (Until.IsNull()) {
872     Perform(Height);
873   }
874   if(Height == 0) {
875     Perform(Until);
876   }
877   TopExp_Explorer exp(Until, TopAbs_FACE);
878   if (!exp.More()) {
879     Standard_ConstructionError::Raise();
880   }
881   //  myPbase.Orientation(TopAbs_FORWARD);
882   myGluedF.Clear();
883   myPerfSelection = BRepFeat_NoSelection;
884   PerfSelectionValid();
885   mySFrom.Nullify();
886   ShapeFromValid();
887   mySUntil = Until;
888   Standard_Boolean Trf = TransformShapeFU(1);
889   ShapeUntilValid();
890   Handle(Geom_Curve) C = TestCurve(myPbase);  
891   Standard_Integer sens = SensOfPrism(C, mySUntil);
892
893   LocOpe_DPrism theDPrism(myPbase,sens*Height,myAngle);
894   TopoDS_Shape VraiDPrism = theDPrism.Shape();
895
896   if(!Trf) {    // case face finished
897     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
898     myGShape = VraiDPrism;
899     GeneratedShapeValid();
900     TopoDS_Shape Base = theDPrism.FirstShape();
901     exp.Init(Base, TopAbs_FACE);
902     TopoDS_Face theBase = TopoDS::Face(exp.Current());
903     exp.Next();
904     if(exp.More()) {
905       NotDone();
906       myStatusError = BRepFeat_InvFirstShape;
907       return;
908     }
909
910     GluedFacesValid();
911     theDPrism.Curves(myCurves);
912     myBCurve = theDPrism.BarycCurve();
913     GlobalPerform();
914   }
915   else {      // case support
916     MajMap(myPbase,theDPrism,myMap,myFShape,myLShape);
917     Handle(Geom_Curve) C1;
918     if(sens == -1) {
919       C1 = C->Reversed();      
920     }
921     else {
922       C1 = C;
923     }
924     TColGeom_SequenceOfCurve scur;
925     scur.Clear();    
926     scur.Append(C1);
927     LocOpe_CSIntersector ASI(mySUntil);
928     ASI.Perform(scur);
929     TopAbs_Orientation Or;
930     if (ASI.IsDone() && ASI.NbPoints(1) >=1) {
931       if (myFuse == 1) {
932         Or = ASI.Point(1,1).Orientation();
933       }
934       else {
935         Or = ASI.Point(1,ASI.NbPoints(1)).Orientation();
936       }
937       //      Standard_Real prm = ASI.Point(1,1).Parameter();
938       //      if(prm < 0) Or = TopAbs::Reverse(Or);
939       TopoDS_Face FUntil = ASI.Point(1,1).Face();
940       TopoDS_Shape Comp;
941       BRep_Builder B;
942       B.MakeCompound(TopoDS::Compound(Comp));
943       TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, Or);
944       if (!S.IsNull()) B.Add(Comp,S);
945       //modified by NIZNHY-PKV Thu Mar 21 17:47:14 2002 f
946       //BRepAlgo_Cut trP(VraiDPrism,Comp);
947       BRepAlgoAPI_Cut trP(VraiDPrism,Comp);
948       //modified by NIZNHY-PKV Thu Mar 21 17:47:21 2002 t
949       if(myFuse == 1) {
950         //modified by NIZNHY-PKV Thu Mar 21 17:47:42 2002 f
951         //BRepAlgo_Fuse f(mySbase, trP.Shape());
952         //myShape = f.Shape();
953         //UpdateDescendants(f.Builder(), myShape, Standard_False);
954         BRepAlgoAPI_Fuse f(mySbase, trP.Shape());
955         myShape = f.Shape();
956         UpdateDescendants(f, myShape, Standard_False);
957         //modified by NIZNHY-PKV Thu Mar 21 17:47:49 2002 t
958         Done();
959       }
960       else if(myFuse == 0) {
961         //modified by NIZNHY-PKV Thu Mar 21 17:48:18 2002 f
962         //BRepAlgo_Cut c(mySbase, trP.Shape());
963         //myShape = c.Shape();
964         //UpdateDescendants(c.Builder(), myShape, Standard_False);
965         BRepAlgoAPI_Cut c(mySbase, trP.Shape());
966         myShape = c.Shape();
967         UpdateDescendants(c, myShape, Standard_False);
968         //modified by NIZNHY-PKV Thu Mar 21 17:48:26 2002 t
969         Done();
970       }
971       else {
972         myShape = trP.Shape();
973         Done(); 
974       }
975     }             
976   }
977 }
978
979
980 //=======================================================================
981 //function : Curves
982 //purpose  : curves parallel to the axis of the prism
983 //=======================================================================
984
985 void BRepFeat_MakeDPrism::Curves(TColGeom_SequenceOfCurve& scur)
986 {
987   scur = myCurves;
988 }
989
990 //============================================================================
991 // function : BRepFeat_BossEgdes
992 // purpose: Determination of TopEdges and LatEdges.
993 //          sig = 1 -> TopEdges = FirstShape of the DPrism
994 //          sig = 2 -> TOpEdges = LastShape of the DPrism
995 //============================================================================
996 void BRepFeat_MakeDPrism::BossEdges (const Standard_Integer signature) 
997 {
998 #ifdef OCCT_DEBUG
999   Standard_Boolean trc = BRepFeat_GettraceFEAT();
1000   if (trc) cout << "BRepFeat_MakeDPrism::BossEdges (integer)" << endl;
1001 #endif
1002   TopTools_ListOfShape theLastShape;
1003   theLastShape.Clear();
1004   if (signature == 1 || signature == -1) {
1005     theLastShape = FirstShape();
1006   }
1007   else if (signature == 2 || signature == -2) {
1008     theLastShape = LastShape();
1009   }
1010   else {
1011     return;
1012   }
1013
1014   // Edges Top
1015   TopTools_ListIteratorOfListOfShape itLS;
1016   TopExp_Explorer ExpE;
1017   for (itLS.Initialize(theLastShape);itLS.More();itLS.Next()) {
1018     const TopoDS_Face& FF = TopoDS::Face(itLS.Value());
1019     for (ExpE.Init(FF,TopAbs_EDGE);ExpE.More();ExpE.Next()) {
1020       const TopoDS_Edge& EE = TopoDS::Edge(ExpE.Current());
1021       myTopEdges.Append(EE);  
1022     }
1023   }
1024
1025   // Edges Bottom    
1026   if (signature < 0) {
1027     // Attention check if TgtEdges is important
1028     myLatEdges = NewEdges();
1029   }
1030   else if (signature > 0) {
1031     if ( !myShape.IsNull() ) {
1032       TopTools_MapOfShape MapE;
1033       Standard_Boolean Found;
1034
1035       TopExp_Explorer ExpF;
1036       for (ExpF.Init(myShape,TopAbs_FACE);ExpF.More();ExpF.Next()) {
1037         Found = Standard_False;
1038         const TopoDS_Face& FF = TopoDS::Face(ExpF.Current());   
1039         for (itLS.Initialize(theLastShape);itLS.More();itLS.Next()) {
1040           const TopoDS_Face& TopFace = TopoDS::Face(itLS.Value());
1041           if (!FF.IsSame(TopFace)) {
1042             for (ExpE.Init(FF,TopAbs_EDGE);ExpE.More() && !Found ;ExpE.Next()) {
1043               const TopoDS_Edge& E1 = TopoDS::Edge(ExpE.Current());
1044               TopoDS_Vertex V1,V2;
1045               TopExp::Vertices (E1,V1,V2);
1046               TopTools_ListIteratorOfListOfShape it(myTopEdges);
1047               for (;it.More() && !Found ; it.Next()) {
1048                 TopoDS_Edge E2 = TopoDS::Edge(it.Value());
1049                 TopoDS_Vertex VT1,VT2;
1050                 TopExp::Vertices (E2,VT1,VT2);
1051
1052                 if (V1.IsSame(VT1) || V1.IsSame(VT2) || V2.IsSame(VT1) || V2.IsSame(VT2)) {
1053                   Found = Standard_True;
1054                   TopExp_Explorer ExpE2;
1055                   for (ExpE2.Init(FF,TopAbs_EDGE);ExpE2.More();ExpE2.Next()) {
1056                     const TopoDS_Edge& E3 = TopoDS::Edge(ExpE2.Current());
1057                     if (MapE.Contains(E3)) { 
1058                       MapE.Remove(E3); }
1059                     else  {
1060                       MapE.Add(E3);}
1061                   }
1062                 }
1063               }
1064             }
1065           }
1066         }
1067       }
1068
1069       TopTools_ListIteratorOfListOfShape it(myTopEdges);
1070       for (;it.More() ; it.Next()) {
1071         if (MapE.Contains(it.Value())) {MapE.Remove(it.Value()); }
1072       }
1073
1074       TopTools_MapIteratorOfMapOfShape itMap;
1075       for (itMap.Initialize(MapE);itMap.More();itMap.Next()) {
1076         if (!BRep_Tool::Degenerated(TopoDS::Edge(itMap.Key())))
1077           myLatEdges.Append(itMap.Key());
1078       }
1079     }
1080   }
1081 }
1082
1083
1084
1085 //============================================================================
1086 // function : BRepFeat_TopEgdes
1087 // Purpose: Returns the list of TopoDS Edges of the top of the boss           
1088 //============================================================================
1089 const TopTools_ListOfShape& BRepFeat_MakeDPrism::TopEdges () 
1090 {    
1091   return myTopEdges;
1092 }
1093
1094 //============================================================================
1095 // function : BRepFeat_LatEgdes
1096 // Purpose: Returns the list of TopoDS Edges of the top of the boss           
1097 //============================================================================
1098 const TopTools_ListOfShape& BRepFeat_MakeDPrism::LatEdges () 
1099 {
1100   return myLatEdges;
1101 }
1102
1103 //=======================================================================
1104 //function : BarycCurve
1105 //purpose  : passe par le centre de masses de la primitive
1106 //=======================================================================
1107
1108 Handle(Geom_Curve) BRepFeat_MakeDPrism::BarycCurve()
1109 {
1110   return myBCurve;
1111 }
1112
1113
1114 //=======================================================================
1115 //function : HeightMax
1116 //purpose  : Calculate the height of the prism following the parameters of the bounding box 
1117 //=======================================================================
1118
1119 static Standard_Real HeightMax(const TopoDS_Shape& theSbase,  // shape initial
1120   const TopoDS_Face&  theSkface, // face de sketch
1121   const TopoDS_Shape& theSFrom,  // shape from
1122   const TopoDS_Shape& theSUntil) // shape until
1123 {
1124   Bnd_Box Box;
1125   BRepBndLib::Add(theSbase,Box);
1126   BRepBndLib::Add(theSkface,Box);
1127   if(!theSFrom.IsNull()) {
1128     BRepBndLib::Add(theSFrom,Box);
1129   }
1130   if(!theSUntil.IsNull()) {
1131     BRepBndLib::Add(theSUntil,Box);
1132   }
1133   Standard_Real c[6];
1134
1135   Box.Get(c[0],c[2],c[4],c[1],c[3],c[5]);
1136   //  Standard_Real parmin=c[0], parmax = c[0];
1137   //  for(Standard_Integer i = 0 ; i < 6; i++) {
1138   //    if(c[i] > parmax) parmax = c[i];
1139   //    if(c[i] < parmin ) parmin = c[i];    
1140   //  }
1141   //  Standard_Real Height = abs(2.*(parmax - parmin));
1142   //  return(2.*Height);
1143   //#ifndef OCCT_DEBUG
1144   Standard_Real par = Max(  Max( fabs(c[1] - c[0]), fabs(c[3] - c[2]) ), fabs(c[5] - c[4]) );
1145   //#else
1146   //  Standard_Real par = Max(  Max( abs(c[1] - c[0]), abs(c[3] - c[2]) ), abs(c[5] - c[4]) );
1147   //#endif
1148 #ifdef OCCT_DEBUG
1149   cout << "Height = > " <<  par  << endl;
1150 #endif
1151   return par;
1152 }
1153
1154
1155
1156
1157
1158 //=======================================================================
1159 //function : SensOfPrism
1160 //purpose  : determine the direction of prism generation
1161 //=======================================================================
1162 Standard_Integer SensOfPrism(const Handle(Geom_Curve) C,
1163   const TopoDS_Shape& Until)
1164 {
1165   LocOpe_CSIntersector ASI1(Until);
1166   TColGeom_SequenceOfCurve scur;
1167   scur.Append(C);
1168   ASI1.Perform(scur);
1169   Standard_Integer sens = 1;
1170   if(ASI1.IsDone() && ASI1.NbPoints(1) >= 1) {
1171     Standard_Integer nb = ASI1.NbPoints(1);
1172     Standard_Real prm1 = ASI1.Point(1, 1).Parameter();
1173     Standard_Real prm2 = ASI1.Point(1, nb).Parameter();
1174     if(prm1 < 0. && prm2 < 0.) {
1175       sens = -1;
1176     }
1177   }
1178   else if(BRepFeat::ParametricBarycenter(Until,C) < 0) {
1179     sens = -1;
1180   }
1181   else {}
1182   return sens;
1183 }
1184
1185 //=======================================================================
1186 //function : MajMap
1187 //purpose  : 
1188 //=======================================================================
1189
1190 static void MajMap(const TopoDS_Shape& theB,
1191   const LocOpe_DPrism& theP,
1192   TopTools_DataMapOfShapeListOfShape& theMap, // myMap
1193   TopoDS_Shape& theFShape,  // myFShape
1194   TopoDS_Shape& theLShape) // myLShape
1195 {
1196   TopExp_Explorer exp;
1197   if(!theP.FirstShape().IsNull()) {    
1198     exp.Init(theP.FirstShape(),TopAbs_WIRE);
1199     if (exp.More()) {
1200       theFShape = exp.Current();
1201       TopTools_ListOfShape thelist;
1202       theMap.Bind(theFShape, thelist);
1203       for (exp.Init(theP.FirstShape(),TopAbs_FACE);exp.More();exp.Next()) {
1204         theMap(theFShape).Append(exp.Current());
1205       }
1206     }
1207   }
1208
1209   if(!theP.LastShape().IsNull()) {    
1210     exp.Init(theP.LastShape(),TopAbs_WIRE);
1211     if (exp.More()) {
1212       theLShape = exp.Current();
1213       TopTools_ListOfShape thelist1;
1214       theMap.Bind(theLShape, thelist1);
1215       for (exp.Init(theP.LastShape(),TopAbs_FACE);exp.More();exp.Next()) {
1216         theMap(theLShape).Append(exp.Current());
1217       }
1218     }
1219   }
1220
1221   for (exp.Init(theB,TopAbs_EDGE); exp.More(); exp.Next()) {
1222     if (!theMap.IsBound(exp.Current())) {
1223       const TopoDS_Edge& edg = TopoDS::Edge(exp.Current());
1224       TopTools_ListOfShape thelist2;
1225       theMap.Bind(edg, thelist2);
1226       theMap(edg) = theP.Shapes(edg);
1227     }
1228   }
1229 }
1230
1231
1232 //=======================================================================
1233 //function : MajMap
1234 //purpose  : 
1235 //=======================================================================
1236
1237 static Handle(Geom_Curve) TestCurve(const TopoDS_Face& Base)
1238 {
1239   gp_Pnt bar(0., 0., 0.);
1240   TColgp_SequenceOfPnt spt;
1241   LocOpe::SampleEdges(Base,spt);
1242   for (Standard_Integer jj=1;jj<=spt.Length(); jj++) {
1243     const gp_Pnt& pvt = spt(jj);
1244     bar.ChangeCoord() += pvt.XYZ();
1245   }
1246   bar.ChangeCoord().Divide(spt.Length());
1247   Handle(Geom_Surface) s = BRep_Tool::Surface(Base);
1248   if (s->DynamicType() == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) {
1249     s = Handle(Geom_RectangularTrimmedSurface)::
1250       DownCast(s)->BasisSurface();
1251   }
1252   Handle (Geom_Plane)   P = Handle(Geom_Plane)::DownCast(s);
1253   if(P.IsNull())  {
1254     Handle(Geom_Curve) toto;
1255     return toto;
1256   }
1257   gp_Pln pp = P->Pln();
1258   gp_Dir Normale(pp.Position().XDirection()^pp.Position().YDirection());
1259   gp_Ax1 theAx(bar, Normale);
1260   Handle(Geom_Line) theLin = new Geom_Line(theAx);
1261   return theLin;
1262 }
1263