0030346: Modeling Algorithms - BRepPrimAPI_MakeRevol throws "BRepSweep_Translation...
[occt.git] / src / BRepTest / BRepTest_OtherCommands.cxx
1 // Created on: 1995-04-13
2 // Created by: Jean Yves LEBEY
3 // Copyright (c) 1995-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 <BRepTest.hxx>
18
19 #include <Draw.hxx>
20 #include <Draw_Interpretor.hxx>
21 #include <DBRep.hxx>
22 #include <DrawTrSurf.hxx>
23
24 #include <string.h>
25 #include <stdio.h>
26
27 #include <Standard_ErrorHandler.hxx>
28 #include <Precision.hxx>
29 #include <TCollection_AsciiString.hxx>
30 #include <gp_Lin.hxx>
31 #include <gp_Pnt.hxx>
32 #include <gp_Dir.hxx>
33 #include <gp_Ax1.hxx>
34
35 #include <ElCLib.hxx>
36
37 #include <TColgp_SequenceOfPnt.hxx>
38
39 #include <GeomAbs_JoinType.hxx>
40 #include <Geom_Line.hxx>
41
42 #include <IntCurvesFace_Intersector.hxx>
43
44 #include <TopAbs.hxx>
45 #include <TopAbs_Orientation.hxx>
46
47 #include <TopoDS.hxx>
48 #include <TopoDS_Shape.hxx>
49 #include <TopoDS_Compound.hxx>
50 #include <TopoDS_CompSolid.hxx>
51 #include <TopoDS_Solid.hxx>
52 #include <TopoDS_Shell.hxx>
53 #include <TopoDS_Face.hxx>
54 #include <TopoDS_Wire.hxx>
55 #include <TopoDS_Edge.hxx>
56 #include <TopoDS_Vertex.hxx>
57 #include <TopoDS_Iterator.hxx>
58
59 #include <BRep_Builder.hxx>
60 #include <BRep_Tool.hxx>
61
62 #include <TopExp_Explorer.hxx>
63
64 #include <TopTools_ListOfShape.hxx>
65 #include <TopTools_ListIteratorOfListOfShape.hxx>
66 #include <TopTools_MapOfShape.hxx>
67
68 #include <LocOpe_CSIntersector.hxx>
69 #include <LocOpe_SequenceOfLin.hxx>
70 #include <LocOpe_PntFace.hxx>
71 #include <BRepFeat_MakeDPrism.hxx>
72
73 #include <BRepTools.hxx>
74 #include <BRepIntCurveSurface_Inter.hxx>
75 #include <BRepOffset_MakeOffset.hxx>
76 #include <BRepClass3d_SolidClassifier.hxx>
77 #include <GeomAdaptor_Curve.hxx>
78
79 static 
80   void SampleEdges (const TopoDS_Shape&   theShape, 
81                     TColgp_SequenceOfPnt& theSeq);
82 static 
83   TopoDS_Face NextFaceForPrism (const TopoDS_Shape& shape, 
84                                 const TopoDS_Shape& basis,
85                                 const gp_Ax1&       ax1);
86 static
87   void PrintState (Draw_Interpretor& aDI, 
88                    const TopAbs_State& aState);
89 //
90 static Standard_Integer emptyshape(Draw_Interpretor&, Standard_Integer, const char** );
91 static Standard_Integer subshape  (Draw_Interpretor&, Standard_Integer, const char** );
92 static Standard_Integer brepintcs (Draw_Interpretor&, Standard_Integer, const char** );
93 static Standard_Integer MakeBoss  (Draw_Interpretor&, Standard_Integer, const char** );
94 static Standard_Integer MakeShell (Draw_Interpretor&, Standard_Integer, const char** );
95 static Standard_Integer xbounds   (Draw_Interpretor&, Standard_Integer, const char** );
96 static Standard_Integer xclassify (Draw_Interpretor&, Standard_Integer, const char** );
97
98 //=======================================================================
99 //function : OtherCommands
100 //purpose  : 
101 //=======================================================================
102 void  BRepTest::OtherCommands(Draw_Interpretor& theCommands)
103 {
104   static Standard_Boolean done = Standard_False;
105   if (done) return;
106   done = Standard_True;
107
108   const char* g = "TOPOLOGY other commands";
109
110   theCommands.Add("shape",
111                   "shape name V/E/W/F/Sh/So/CS/C; make a empty shape",__FILE__,emptyshape,g);
112
113   theCommands.Add("subshape",
114                   "subshape name V/E/W/F/Sh/So/CS/C index; get subsshape <index> of given type"
115                   ,__FILE__,subshape,g);
116
117   theCommands.Add("BRepIntCS",
118     "Calcul d'intersection entre face et curve : BRepIntCS curve1 [curve2 ...] shape [res] [tol]"
119                   ,__FILE__,brepintcs,g);
120
121   theCommands.Add("makeboss",  "create a boss on the shape myS", __FILE__, MakeBoss, g);
122   theCommands.Add("mksh", "create a shell on Shape", __FILE__, MakeShell, g);
123   theCommands.Add("xbounds",  "xbounds face", __FILE__, xbounds, g);
124   theCommands.Add("xclassify",  "use xclassify Solid [Tolerance=1.e-7]", __FILE__, xclassify, g);
125   
126
127 }
128 //=======================================================================
129 //function : emptyshape
130 //purpose  : shape : shape name V/E/W/F/SH/SO/CS/C
131 //=======================================================================
132 Standard_Integer emptyshape(Draw_Interpretor& , Standard_Integer n, const char** a)
133 {
134   if (n <= 1) return 1;
135
136   BRep_Builder B;
137   TopoDS_Shape S;
138
139   if (n == 3) {
140     TCollection_AsciiString as(a[2]); as.LowerCase();
141     const char* a2 = as.ToCString();
142
143     if      ( ! strcmp(a2,"c") ) {
144       TopoDS_Compound SS; B.MakeCompound(SS); S = SS;
145     }
146     else if ( ! strcmp(a2,"cs") ) {
147       TopoDS_CompSolid SS; B.MakeCompSolid(SS); S = SS;
148     }
149     else if ( ! strcmp(a2,"so") ) {
150       TopoDS_Solid SS; B.MakeSolid(SS); S = SS;
151     }
152     else if ( ! strcmp(a2,"sh") ) {
153       TopoDS_Shell SS; B.MakeShell(SS); S = SS;
154     }
155     else if ( ! strcmp(a2,"f") ) {
156       TopoDS_Face SS; B.MakeFace(SS); S = SS;
157     }
158     else if ( ! strcmp(a2,"w") ) {
159       TopoDS_Wire SS; B.MakeWire(SS); S = SS;
160     }
161     else if ( ! strcmp(a2,"e") ) {
162       TopoDS_Edge SS; B.MakeEdge(SS); S = SS;
163     }
164     else if ( ! strcmp(a2,"v") ) {
165       TopoDS_Vertex SS; B.MakeVertex(SS); S = SS;
166     }
167     else {
168       return 1;
169     }
170   }
171
172   const char *shapename = a[1];
173   DBRep::Set(shapename,S);
174   return 0;
175 }
176
177 //=======================================================================
178 //function : subshape
179 //purpose  : 
180 //=======================================================================
181 Standard_Integer subshape(Draw_Interpretor& di, Standard_Integer n, const char** a)
182 {
183   if (n <= 2) return 1;
184
185   
186   TopoDS_Shape S = DBRep::Get(a[1]);
187   if (S.IsNull()) return 0;
188   char newname[1024];
189   strcpy(newname,a[1]);
190   char* p = newname;
191   while (*p != '\0') p++;
192   *p = '_';
193   p++;
194   Standard_Integer i = 0;
195   if (n == 3) {
196     Standard_Integer isub = Draw::Atoi(a[2]);
197     TopoDS_Iterator itr(S);
198     while (itr.More()) {
199       i++;
200       if ( i == isub ) {
201         Sprintf(p,"%d",i);
202         DBRep::Set(newname,itr.Value());
203         di.AppendElement(newname);
204         break;
205       }
206       itr.Next();
207     }
208   }
209   else {
210     // explode a type
211     TopAbs_ShapeEnum typ;
212     switch (a[2][0]) {
213       
214     case 'C' :
215     case 'c' :
216       typ = TopAbs_COMPSOLID;
217       break;
218       
219     case 'S' :
220     case 's' :
221       if ((a[2][1] == 'O')||(a[2][1] == 'o')) 
222         typ = TopAbs_SOLID;
223       else if ((a[2][1] == 'H')||(a[2][1] == 'h')) 
224         typ = TopAbs_SHELL;
225       else
226         return 1;
227       break;
228       
229     case 'F' :
230     case 'f' :
231       typ = TopAbs_FACE;
232       break;
233       
234     case 'W' :
235     case 'w' :
236       typ = TopAbs_WIRE;
237       break;
238       
239     case 'E' :
240     case 'e' :
241       typ = TopAbs_EDGE;
242       break;
243       
244     case 'V' :
245     case 'v' :
246       typ = TopAbs_VERTEX;
247       break;
248       
249       default :
250         return 1;
251     }
252     
253     Standard_Integer isub = Draw::Atoi(a[3]);
254     TopTools_MapOfShape M;
255     M.Add(S);
256     TopExp_Explorer ex(S,typ);
257     while (ex.More()) {
258       if (M.Add(ex.Current())) {
259         i++;
260         if ( i == isub ) {
261           Sprintf(p,"%d",i);
262           DBRep::Set(newname,ex.Current());
263           di.AppendElement(newname);
264           break;
265         }
266       }
267       ex.Next();
268     }
269   }
270   return 0;
271 }
272 //=======================================================================
273 //function : brepintcs
274 //purpose  : 
275 //=======================================================================
276 Standard_Integer brepintcs(Draw_Interpretor& di, Standard_Integer n, const char** a)
277 {
278   if (n <= 2) 
279   {
280     cout<<"Invalid input arguments. Should be: curve1 [curve2 ...] shape [result] [tol]"<<endl;
281     return 1;
282   }
283   Standard_Integer indshape = 2;
284   TopoDS_Shape S;
285   for( ; indshape <= n-1 ; indshape++)
286   {
287     S = DBRep::Get(a[indshape]);
288     if(!S.IsNull())
289       break;
290   }
291   if (S.IsNull()) 
292   {
293     cout<<"Invalid input shape"<<endl;
294     return 1;
295   }
296
297   BRepIntCurveSurface_Inter theAlg;
298   double tol=1e-6;
299   if( indshape < n-1)
300   {
301     Standard_Real preci = atof(a[n-1]);
302     if(preci >= Precision::Confusion())
303       tol = preci;
304   }
305   int nbpi=0;
306   gp_Pnt curp;
307   TopoDS_Compound aComp;
308   BRep_Builder aB;
309   aB.MakeCompound(aComp);
310   if (indshape == 2) {
311     Handle(Geom_Curve) C= DrawTrSurf::GetCurve(a[1]);
312     if (C.IsNull()) return 2;
313     GeomAdaptor_Curve acur(C);
314     theAlg.Init(S, acur, tol);
315
316     for (; theAlg.More(); theAlg.Next()) {
317       curp=theAlg.Pnt();
318       TopoDS_Vertex aV;
319
320       aB.MakeVertex(aV, curp, 0);
321       aB.Add(aComp, aV);
322       nbpi++;
323       di<<"Point "<<nbpi<<" : "<<curp.X()<<" "<<curp.Y()<<" "<<curp.Z()<<"\n";
324       char name[64];
325       char* temp = name; // pour portage WNT
326       Sprintf(temp, "%s_%d", "brics", nbpi); 
327       DrawTrSurf::Set(temp, curp);
328     }
329   }
330   else {
331     theAlg.Load(S,tol );
332     for (Standard_Integer il = 1; il<indshape ; il++) 
333     {
334       Handle(Geom_Curve) hl= DrawTrSurf::GetCurve(a[il]);
335       if (!hl.IsNull()) {
336         theAlg.Init(hl);
337         for (; theAlg.More(); theAlg.Next()) {
338           curp=theAlg.Pnt();
339           nbpi++;
340           TopoDS_Vertex aV;
341           aB.MakeVertex(aV, curp, 0);
342           aB.Add(aComp, aV);
343           di<<"Point "<<nbpi<<" : "<<curp.X()<<" "<<curp.Y()<<" "<<curp.Z()<<"\n";
344           char name[64];
345           char* temp = name; // pour portage WNT
346           Sprintf(temp, "%s_%d", "brics", nbpi); 
347           DrawTrSurf::Set(temp, curp);
348         }
349       }
350     }
351   }
352   if(!nbpi)
353     di<<"Points of intersections are not found\n";
354   if(indshape < n-1)
355     DBRep::Set(a[n-1], aComp);
356   //POP pour NT
357   return 0;
358 }
359 //=======================================================================
360 //function : MakeBoss
361 //purpose  : 
362 //=======================================================================
363 Standard_Integer MakeBoss(Draw_Interpretor& , Standard_Integer , const char** a)
364 {
365   TopoDS_Shape myS = DBRep::Get( a[2] );
366
367   TopoDS_Shape myBasis = DBRep::Get( a[3] ) ;
368
369   Standard_Real ang = -0.05235987901687622;
370
371
372   TopoDS_Face basis = TopoDS::Face(myBasis);
373
374   BRepFeat_MakeDPrism DPRISM(myS, basis, basis, ang, 1, Standard_True);
375
376   TopoDS_Shape myFaceOnShape;
377   gp_Pnt Pnt(0.0, 0.0, 50.0);
378   gp_Dir Dir(-0.0, -0.0, -1.0);
379   gp_Ax1 ax(Pnt, Dir);
380   
381   myFaceOnShape = NextFaceForPrism(myS, myBasis, ax);
382
383   DPRISM.Perform (myFaceOnShape);
384   DPRISM.Build();
385
386   if( DPRISM.IsDone() )  DBRep::Set( a[1], DPRISM.Shape() );
387
388   return 0;
389 }
390 //=======================================================================
391 //function : MakeShell
392 //purpose  : 
393 //=======================================================================
394 Standard_Integer MakeShell(Draw_Interpretor& , Standard_Integer , const char** a)
395 {
396
397   TopoDS_Shape aShape = DBRep::Get( a[1] ); 
398   TopTools_ListOfShape Lst;
399   TopExp_Explorer Exp(aShape, TopAbs_FACE);
400   TopoDS_Shape InputShape(DBRep::Get( a[2] ));
401   TopoDS_Face F = TopoDS::Face(InputShape);
402 //  TopoDS_Face F = TopoDS::Face(DBRep::Get( a[2] ));
403   
404   Standard_Real Off = -Draw::Atof( a[3] );
405
406   BRepOffset_MakeOffset Offset;
407
408   Offset.Initialize( aShape, Off,  1.0e-3, BRepOffset_Skin, 
409                                Standard_True , Standard_False , GeomAbs_Arc );
410   Offset.AddFace( F );
411   Offset.MakeThickSolid();
412
413   if( Offset.IsDone() ) {
414   //    SaveShape::Save(Offset.Shape(), "ss");
415     DBRep::Set( a[1], Offset.Shape() );
416   }
417   return 0;
418 }
419 //=======================================================================
420 //function : xbounds
421 //purpose  : 
422 //=======================================================================
423 Standard_Integer xbounds(Draw_Interpretor& di, Standard_Integer n, const char** a)
424 {
425   if (n<2) {
426     di << "Usage : " << a[0] << " face\n";
427     return 0;
428   }
429   //
430  
431   Standard_Real aUMin, aUMax, aVMin, aVMax;
432   TopoDS_Shape aS;
433   TopoDS_Face aF;
434   //
435   aS=DBRep::Get(a[1]);
436   if (aS.IsNull()) {
437     di << " null shapes is not allowed here\n";
438     return 0;
439   }
440   if (aS.ShapeType()!=TopAbs_FACE) {
441     di << " shape" << a[1] <<" must be a face\n";
442     return 0;
443   }
444   //
445   aF=*((TopoDS_Face*)&aS);
446   //
447   BRepTools::UVBounds(aF, aUMin, aUMax, aVMin, aVMax);
448   //
449   TCollection_AsciiString aStr;
450   TCollection_AsciiString sUMin(aUMin);
451   TCollection_AsciiString sUMax(aUMax);
452   TCollection_AsciiString sVMin(aVMin);
453   TCollection_AsciiString sVMax(aVMax);
454   //
455   aStr=aStr+sUMin + "\n";
456   aStr=aStr+sUMax + "\n";
457   aStr=aStr+sVMin + "\n";
458   aStr=aStr+sVMax + "\n";
459   di <<aStr.ToCString();
460   //
461   return 0;
462 }
463 //=======================================================================
464 //function : xclassify
465 //purpose  : 
466 //=======================================================================
467 Standard_Integer xclassify (Draw_Interpretor& aDI, Standard_Integer n, const char** a)
468 {
469   if (n < 2) {
470     aDI<<" use xclassify Solid [Tolerance=1.e-7]\n";
471     return 1;
472   }
473   
474   TopoDS_Shape aS = DBRep::Get(a[1]);
475   if (aS.IsNull()) {
476     aDI<<" Null Shape is not allowed here\n";
477     return 0;
478   }
479   
480   if (aS.ShapeType()!=TopAbs_SOLID) {
481     aDI<< " Shape type must be SOLID\n";
482     return 0;
483   }
484   //
485   Standard_Real aTol=1.e-7;
486   TopAbs_State aState = TopAbs_UNKNOWN;
487   //
488   aTol=1.e-7; 
489   if (n==3) {
490     aTol=Draw::Atof(a[2]);
491   }
492   //
493   BRepClass3d_SolidClassifier aSC(aS);
494   aSC.PerformInfinitePoint(aTol);
495   
496   aState = aSC.State();
497   PrintState(aDI, aState);
498   //
499   return 0;
500 }
501 //=======================================================================
502 //function : PrintState
503 //purpose  : 
504 //=======================================================================
505 void PrintState (Draw_Interpretor& aDI, 
506                  const TopAbs_State& aState)
507 {
508   aDI<<"state is: ";
509   switch (aState) {
510   case TopAbs_IN:
511     aDI<<"IN\n"; 
512     break;
513   case TopAbs_OUT:
514     aDI<<"OUT\n";               
515     break;
516   case TopAbs_ON:        
517     aDI<<"ON\n";        
518     break;
519   case TopAbs_UNKNOWN:
520   default:
521     aDI<<"UNKNOWN\n";           
522     break;
523   }
524 }
525 //=======================================================================
526 //function : NextFaceForPrism
527 //purpose  : Search a face from <shape> which intersects with a line of
528 //           direction <ax1> and location a point of <basis>.
529 //=======================================================================
530 TopoDS_Face NextFaceForPrism (const TopoDS_Shape& shape, 
531                               const TopoDS_Shape& basis,
532                               const gp_Ax1&       ax1)
533 {
534   TopoDS_Face nextFace;
535
536   TColgp_SequenceOfPnt seqPnts;
537   SampleEdges(basis, seqPnts);
538   
539   for (Standard_Integer i=1; i<=seqPnts.Length(); i++) {
540     const gp_Pnt& pt = seqPnts(i);
541     // find a axis through a face
542     gp_Dir dir = ax1.Direction();
543     gp_Ax1 ax1b(pt, dir);
544     
545     LocOpe_CSIntersector ASI(shape);
546     LocOpe_SequenceOfLin slin;
547     slin.Append(ax1b);
548     ASI.Perform(slin);
549     
550     if (ASI.IsDone()) {
551       Standard_Integer no=1, IndFrom, IndTo;
552       TopAbs_Orientation theOr;
553       Standard_Real min = 1.e-04, Tol = -Precision::Confusion();
554       if (ASI.LocalizeAfter (no, min, Tol, theOr, IndFrom, IndTo))  {
555         nextFace = ASI.Point(no, IndFrom).Face();  
556         break;
557       }
558     }
559   }
560   
561   return nextFace;
562 }
563
564
565 //=======================================================================
566 //function : SampleEdges
567 //purpose  : Sampling of <theShape>.
568 //design   : Collect the vertices and points on the edges
569 //=======================================================================
570 void SampleEdges (const TopoDS_Shape&   theShape, TColgp_SequenceOfPnt& theSeq)
571 {
572
573   theSeq.Clear();
574
575    
576   TopTools_MapOfShape theMap;
577   TopExp_Explorer exp;
578   
579   // Adds all vertices/pnt
580   for (exp.Init(theShape,TopAbs_VERTEX); exp.More(); exp.Next()) {
581     if (theMap.Add(exp.Current())) {
582       theSeq.Append (BRep_Tool::Pnt(TopoDS::Vertex(exp.Current())));
583     }
584   }
585
586   // Computes points on edge, but does not take the extremities into account
587   Standard_Integer NECHANT = 5;
588   Handle(Geom_Curve) C;
589   Standard_Real f,l,prm;
590   for (exp.Init (theShape,TopAbs_EDGE); exp.More(); exp.Next()) {
591     const TopoDS_Edge& edg = TopoDS::Edge(exp.Current());
592     if (theMap.Add(edg)) {
593       if (!BRep_Tool::Degenerated(edg)) {
594         C = BRep_Tool::Curve(edg,f,l);
595         for (Standard_Integer i=1; i < NECHANT; i++) {
596           prm = ((NECHANT-i)*f+i*l)/NECHANT;
597           theSeq.Append (C->Value(prm));
598         }
599       }
600     }
601   }
602 }
603