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