0024624: Lost word in license statement in source files
[occt.git] / src / BRepTest / BRepTest_TopologyCommands.cxx
1 // Created on: 1993-07-22
2 // Created by: Remi LEQUETTE
3 // Copyright (c) 1993-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 #ifdef HAVE_CONFIG_H
18 # include <config.h>
19 #endif
20 #include <BRepTest.hxx>
21 #include <DBRep.hxx>
22 #include <Draw_Interpretor.hxx>
23 #include <DrawTrSurf.hxx>
24 #include <Geom_Plane.hxx>
25 #include <Draw_Appli.hxx>
26 #include <BRep_Builder.hxx>
27
28 #include <BRepAlgo_Fuse.hxx>
29 #include <BRepAlgo_Common.hxx>
30 #include <BRepAlgo_Cut.hxx>
31 #include <BRepAlgo_Section.hxx>
32
33 #include <BRepFilletAPI_MakeFillet.hxx>
34 #include <BRepBuilderAPI_MakeVertex.hxx>
35 #include <BRepPrimAPI_MakeHalfSpace.hxx>
36 #include <BRepAlgo_FaceRestrictor.hxx>
37 #include <BRepExtrema_ExtPF.hxx>
38 #include <BRepLProp_SLProps.hxx>
39 #include <TopTools_ListIteratorOfListOfShape.hxx>
40 #include <TopoDS.hxx>
41 #include <TopoDS_Edge.hxx>
42 #include <TopoDS_Compound.hxx>
43 #include <TopoDS_Wire.hxx>
44 #include <TopExp_Explorer.hxx>
45 #include <TopOpeBRepBuild_HBuilder.hxx>
46 #include <TopOpeBRepDS_HDataStructure.hxx>
47 #include <gp.hxx>
48 #include <gp_Pln.hxx>
49 #include <TopTools_IndexedMapOfShape.hxx>
50 #include <TopExp.hxx>
51
52 #ifdef HAVE_STRINGS_H
53 # include <strings.h>
54 #endif
55
56 //=======================================================================
57 // topop
58 //=======================================================================
59
60 static Standard_Integer topop(Draw_Interpretor& , Standard_Integer n, const char** a)
61 {
62   if (n < 4) return 1;
63
64   TopoDS_Shape s1 = DBRep::Get(a[2]);
65   TopoDS_Shape s2 = DBRep::Get(a[3]);
66
67   if (s1.IsNull() || s2.IsNull()) return 1;
68
69   TopoDS_Shape res;
70
71   if (*a[0] == 'f')
72     res = BRepAlgo_Fuse(s1,s2);
73   else if (*(a[0]+1) == 'o')
74     res = BRepAlgo_Common(s1,s2);
75   else 
76     res = BRepAlgo_Cut(s1,s2);
77
78   DBRep::Set(a[1],res);
79
80   return 0;
81 }
82
83
84 //=======================================================================
85 // section
86 //=======================================================================
87
88 static Standard_Integer section(Draw_Interpretor& , Standard_Integer n, const char** a)
89 {
90
91   if (n < 4) return 1;
92
93   TopoDS_Shape s1 = DBRep::Get(a[2]);
94   TopoDS_Shape s2 = DBRep::Get(a[3]);
95   
96   if (s1.IsNull() || s2.IsNull())
97     return 1;
98
99   BRepAlgo_Section Sec(s1, s2, Standard_False);
100   TopoDS_Shape res;
101
102   for (int i=4; i < n; i++) {
103     if (!strcasecmp(a[i], "-2d"))
104     {
105       Sec.ComputePCurveOn1(Standard_True);
106       Sec.ComputePCurveOn2(Standard_True);
107     } 
108     else if (!strcasecmp(a[i], "-2d1")) 
109     {
110       Sec.ComputePCurveOn1(Standard_True);
111       Sec.ComputePCurveOn2(Standard_False);
112     } 
113     else if (!strcasecmp(a[i], "-2d2")) 
114     {
115       Sec.ComputePCurveOn1(Standard_False);
116       Sec.ComputePCurveOn2(Standard_True);
117     } 
118     else if (!strcasecmp(a[i], "-no2d"))
119     {
120       Sec.ComputePCurveOn1(Standard_False);
121       Sec.ComputePCurveOn2(Standard_False);
122     } 
123     else if (!strcasecmp(a[i], "-a")) 
124       Sec.Approximation(Standard_True);
125     else if (strcasecmp(a[i], "-p"))
126     {
127       cout << "Unknown option: " << a[i] << endl;
128       return 1;
129     }
130   }
131
132   res = Sec.Shape();
133   
134   DBRep::Set(a[1],res);
135   
136   return 0;
137 }
138
139 //=======================================================================
140 // psection
141 //=======================================================================
142
143 static Standard_Integer psection(Draw_Interpretor& , Standard_Integer n, const char** a)
144 {
145   if (n < 4) return 1;
146
147   TopoDS_Shape s = DBRep::Get(a[2]);
148   if (s.IsNull()) return 1;
149
150   Handle(Geom_Surface) ps = DrawTrSurf::GetSurface(a[3]);
151   if (ps.IsNull()) return 1;
152
153   Handle(Geom_Plane) pg = Handle(Geom_Plane)::DownCast(ps);
154   if (pg.IsNull()) return 1;
155
156   const gp_Pln& p = pg->Pln();
157
158   TopoDS_Shape res = BRepAlgo_Section(s,p);
159
160   DBRep::Set(a[1],res);
161
162   return 0;
163 }
164
165 static Standard_Integer halfspace(Draw_Interpretor& di,
166                                   Standard_Integer n, const char** a)
167 {
168   if (n < 6) return 1;
169
170   // Le point indiquant le cote "matiere".
171   gp_Pnt RefPnt = gp_Pnt(Draw::Atof(a[3]),Draw::Atof(a[4]),Draw::Atof(a[5]));
172
173   TopoDS_Shape Face = DBRep::Get(a[2],TopAbs_FACE);
174   if ( Face.IsNull()) {
175     TopoDS_Shape Shell  = DBRep::Get(a[2],TopAbs_SHELL);
176     if (Shell.IsNull()) {
177       //cout << a[2] << " must be a face or a shell" << endl;
178       di << a[2] << " must be a face or a shell" << "\n";
179       return 1;
180     }
181     else {
182       BRepPrimAPI_MakeHalfSpace Half(TopoDS::Shell(Shell),RefPnt);
183       if ( Half.IsDone()) {
184         DBRep::Set(a[1],Half.Solid());
185       }
186       else {
187         //cout << " HalfSpace NotDone" << endl;
188         di << " HalfSpace NotDone" << "\n";
189         return 1;
190       }
191     }
192   }
193   else {
194     BRepPrimAPI_MakeHalfSpace Half(TopoDS::Face(Face),RefPnt);
195     if ( Half.IsDone()) {
196       DBRep::Set(a[1],Half.Solid());
197     }
198     else {
199       //cout << " HalfSpace NotDone" << endl;
200       di << " HalfSpace NotDone" << "\n";
201       return 1;
202     }
203   }
204   return 0;
205 }
206
207 //=======================================================================
208 //function : buildfaces
209 //purpose  : 
210 //=======================================================================
211
212 static Standard_Integer buildfaces(Draw_Interpretor& , Standard_Integer narg, const char** a)
213 {
214   if (narg < 4) return 1;
215   
216   TopoDS_Shape InputShape(DBRep::Get( a[2] ,TopAbs_FACE));
217   TopoDS_Face F = TopoDS::Face(InputShape);
218 //  TopoDS_Face F = TopoDS::Face(DBRep::Get(a[2],TopAbs_FACE));
219   BRepAlgo_FaceRestrictor FR;
220   FR.Init(F);
221   
222   for (Standard_Integer i = 3 ; i < narg ; i++) {
223     TopoDS_Shape InputWire(DBRep::Get(a[i],TopAbs_WIRE));
224     TopoDS_Wire W = TopoDS::Wire(InputWire);
225 //    TopoDS_Wire W = TopoDS::Wire(DBRep::Get(a[i],TopAbs_WIRE));
226     FR.Add(W);
227   }
228   FR.Perform();
229   if (!FR.IsDone()) return 1;
230   
231   TopoDS_Compound Res;
232   BRep_Builder BB;
233   BB.MakeCompound(Res);
234
235   for (; FR.More(); FR.Next()) {
236     TopoDS_Face FF = FR.Current();
237     BB.Add(Res,FF);
238     DBRep::Set(a[1],Res);
239   }
240   return 0;
241 }
242
243 //=======================================================================
244 //function : TopologyCommands
245 //purpose  : 
246 //=======================================================================
247
248 void  BRepTest::TopologyCommands(Draw_Interpretor& theCommands)
249 {
250   static Standard_Boolean done = Standard_False;
251   if (done) return;
252   done = Standard_True;
253
254   DBRep::BasicCommands(theCommands);
255
256   const char* g = "TOPOLOGY Topological operation commands";
257   
258   theCommands.Add("fuse","fuse result s1 s2",__FILE__,topop,g);
259   theCommands.Add("common","common result s1 s2",__FILE__,topop,g);
260   theCommands.Add("cut","cut result part tool",__FILE__,topop,g);
261   theCommands.Add("section","section result s1 s2 [-no2d/-2d/-2d1/-2d2] [-p/-a]",__FILE__,section,g);
262   theCommands.Add("psection","psection result s plane",__FILE__,psection,g);
263   theCommands.Add("halfspace","halfspace result face/shell x y z",__FILE__,halfspace,g);
264   theCommands.Add("buildfaces","buildfaces result faceReference wire1 wire2 ...",__FILE__,buildfaces,g);
265 }