0026150: BRepOffsetAPI_ThruSections doesn't implement history method Generated().
[occt.git] / src / BRepTest / BRepTest_SweepCommands.cxx
index 4c1203b..2cd7a1c 100644 (file)
@@ -43,6 +43,7 @@
 #include <TColgp_Array1OfPnt2d.hxx>
 
 static BRepOffsetAPI_MakePipeShell* Sweep= 0;
+static BRepOffsetAPI_ThruSections* Generator = 0;
 
 #include <stdio.h>
 #include <Geom_Curve.hxx>
@@ -130,11 +131,11 @@ static Standard_Integer pipe(Draw_Interpretor& di,
 {
   if (n == 1)
   {
-    di << "pipe result Wire_spine Profile [Mode [Approx]]" << "\n";
-    di << "Mode = 0 - CorrectedFrenet," << "\n";
-    di << "     = 1 - Frenet," << "\n";
-    di << "     = 2 - DiscreteTrihedron" << "\n";
-    di << "Approx - force C1-approximation if result is C0" << "\n";
+    di << "pipe result Wire_spine Profile [Mode [Approx]]\n";
+    di << "Mode = 0 - CorrectedFrenet,\n";
+    di << "     = 1 - Frenet,\n";
+    di << "     = 2 - DiscreteTrihedron\n";
+    di << "Approx - force C1-approximation if result is C0\n";
     return 0;
   }
   
@@ -200,6 +201,11 @@ static Standard_Integer geompipe(Draw_Interpretor& ,
     rotate = (Draw::Atoi(a[k++])==1);
   GeomFill_Pipe aPipe(ProfileCurve,aAdaptCurve,cStart,ByACR,rotate);
   aPipe.Perform(Standard_True);
+  if (!aPipe.IsDone())
+  {
+    cout << "GeomFill_Pipe cannot make a surface" << endl;
+    return 1;
+  }
   Handle(Geom_Surface) Sur=aPipe.Surface();
   TopoDS_Face F;
   if(!Sur.IsNull())
@@ -221,11 +227,11 @@ Standard_Integer evolved(Draw_Interpretor& di, Standard_Integer n, const char**
     //cout << "        is given in the referencial axis. " << endl;
     //cout << " 2) evolved result base profil o : "<< endl;
     //cout << "        This position is automatically computed." << endl;
-    di << " 1) evolved result base profil : "<< "\n";
-    di << "        The relative position of the profil on the base" << "\n";
-    di << "        is given in the referencial axis. " << "\n";
-    di << " 2) evolved result base profil o : "<< "\n";
-    di << "        This position is automatically computed." << "\n";
+    di << " 1) evolved result base profil : \n";
+    di << "        The relative position of the profil on the base\n";
+    di << "        is given in the referencial axis. \n";
+    di << " 2) evolved result base profil o : \n";
+    di << "        This position is automatically computed.\n";
     return 0;
   }
 
@@ -318,19 +324,19 @@ Standard_Integer gener(Draw_Interpretor&, Standard_Integer n, const char** a)
 
   TopoDS_Shape Shape; 
 
-  BRepFill_Generator Generator;
+  BRepFill_Generator aGenerator;
   
   for ( Standard_Integer i = 2; i<= n-1 ; i++) {
     Shape = DBRep::Get(a[i],TopAbs_WIRE);
     if ( Shape.IsNull()) 
       return 1;
 
-    Generator.AddWire(TopoDS::Wire(Shape));
+    aGenerator.AddWire(TopoDS::Wire(Shape));
   }
 
-  Generator.Perform();
+  aGenerator.Perform();
 
-  TopoDS_Shell Shell = Generator.Shell();
+  TopoDS_Shell Shell = aGenerator.Shell();
   
   DBRep::Set(a[1], Shell);
 
@@ -363,7 +369,12 @@ Standard_Integer thrusections(Draw_Interpretor&, Standard_Integer n, const char*
   Standard_Boolean issolid = ( Draw::Atoi(a[index]) == 1 );
   Standard_Boolean isruled = ( Draw::Atoi(a[index+1]) == 1 );
 
-  BRepOffsetAPI_ThruSections Generator(issolid,isruled);
+  if (Generator != 0)
+  {
+    delete Generator; 
+    Generator = 0;
+  }
+  Generator = new BRepOffsetAPI_ThruSections(issolid,isruled);
   
   Standard_Integer NbEdges = 0;
   Standard_Boolean IsFirstWire = Standard_False;
@@ -372,7 +383,7 @@ Standard_Integer thrusections(Draw_Interpretor&, Standard_Integer n, const char*
     Shape = DBRep::Get(a[i], TopAbs_WIRE);
     if (!Shape.IsNull())
       {
-       Generator.AddWire(TopoDS::Wire(Shape));
+       Generator->AddWire(TopoDS::Wire(Shape));
        if (!IsFirstWire)
          IsFirstWire = Standard_True;
        else
@@ -383,7 +394,7 @@ Standard_Integer thrusections(Draw_Interpretor&, Standard_Integer n, const char*
        Shape = DBRep::Get(a[i], TopAbs_VERTEX);
        IsWire = Standard_False;
        if (!Shape.IsNull())
-         Generator.AddVertex(TopoDS::Vertex(Shape));
+         Generator->AddVertex(TopoDS::Vertex(Shape));
        else
          return 1;
       }
@@ -402,12 +413,12 @@ Standard_Integer thrusections(Draw_Interpretor&, Standard_Integer n, const char*
   }
 
   check = (check || !samenumber);
-  Generator.CheckCompatibility(check);
+  Generator->CheckCompatibility(check);
 
-  Generator.Build();
+  Generator->Build();
 
-  if (Generator.IsDone()) {
-    TopoDS_Shape Shell = Generator.Shape();
+  if (Generator->IsDone()) {
+    TopoDS_Shape Shell = Generator->Shape();
     DBRep::Set(a[index-1], Shell);
   }
   else {
@@ -416,6 +427,55 @@ Standard_Integer thrusections(Draw_Interpretor&, Standard_Integer n, const char*
   return 0;
 }
 
+//============================================================================
+//function : genthrus
+//purpose  : returns generated shape for subshape of a section of thrusections
+//           Thrusections must be done previously
+//============================================================================
+static Standard_Integer genthrus(Draw_Interpretor& di,
+                                 Standard_Integer n, const char** a)
+{
+  if (n != 3)
+  {
+    di << "genthrus: it is called after thrusections command\n";
+    di << "returns:\n";
+    di << "- chain of generated faces for sub-edge of a profile;\n";
+    di << "- chain of generated edges for sub-vertex of a profile;\n";
+    di << "- bunch of chains of generated edges for start or end vertex if it is degenerated section.\n";
+    di << "Usage: genthrus res subshape_of_profile, thrusections must be done\n";
+    return 1;
+  }
+  if (Generator == 0)
+  {
+    di << "You have forgotten the <<thrusections>> command !\n";
+    return 1;
+  }
+  if (!Generator->IsDone())
+  {
+    di << "Thrusections is not done\n";
+    return 1;
+  }
+  TopoDS_Shape aShape = DBRep::Get(a[2]);
+  if (aShape.IsNull())
+  {
+    cout<<"Null subshape"<<endl;
+    return 1;
+  }
+  const TopTools_ListOfShape& Edges = Generator->Generated(aShape);
+  TopoDS_Compound aCompound;
+  BRep_Builder BB;
+  BB.MakeCompound(aCompound);
+  TopTools_ListIteratorOfListOfShape iter(Edges);
+  for (; iter.More(); iter.Next())
+  {
+    const TopoDS_Shape& anEdge = iter.Value();
+    BB.Add(aCompound, anEdge);
+  }
+
+  DBRep::Set(a[1], aCompound);
+  return 0;
+}
+
 //=======================================================================
 //  mksweep
 //=======================================================================
@@ -450,22 +510,22 @@ static Standard_Integer setsweep(Draw_Interpretor& di,
     //cout << "   -CN dx dy dz : BiNormal is given by dx dy dz" << endl;
     //cout << "   -FX Tx Ty TZ [Nx Ny Nz] : Tangent and Normal are fixed" <<endl;
     //cout << "   -G guide  0|1(ACR|Plan)  0|1(contact|no contact) : with guide"<<endl;
-    di << "setsweep options [arg1 [arg2 [...]]] : options are :" << "\n";
-    di << "   -FR : Tangent and Normal are given by Frenet trihedron" <<"\n";
-    di << "   -CF : Tangente is given by Frenet," << "\n";
-    di << "         the Normal is computed to minimize the torsion " << "\n";
-    di << "   -DT : discrete trihedron" << "\n";
-    di << "   -DX Surf : Tangent and Normal are given by Darboux trihedron," <<"\n";     
-    di << "       Surf have to be a shell or a face" <<"\n";
-    di << "   -CN dx dy dz : BiNormal is given by dx dy dz" << "\n";
-    di << "   -FX Tx Ty TZ [Nx Ny Nz] : Tangent and Normal are fixed" <<"\n";
-    di << "   -G guide  0|1(Plan|ACR)  0|1|2(no contact|contact|contact on border) : with guide"<<"\n";
+    di << "setsweep options [arg1 [arg2 [...]]] : options are :\n";
+    di << "   -FR : Tangent and Normal are given by Frenet trihedron\n";
+    di << "   -CF : Tangente is given by Frenet,\n";
+    di << "         the Normal is computed to minimize the torsion \n";
+    di << "   -DT : discrete trihedron\n";
+    di << "   -DX Surf : Tangent and Normal are given by Darboux trihedron,\n";     
+    di << "       Surf have to be a shell or a face\n";
+    di << "   -CN dx dy dz : BiNormal is given by dx dy dz\n";
+    di << "   -FX Tx Ty TZ [Nx Ny Nz] : Tangent and Normal are fixed\n";
+    di << "   -G guide  0|1(Plan|ACR)  0|1|2(no contact|contact|contact on border) : with guide\n";
     return 0;
   }
 
    if (Sweep ==0) {
      //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
-     di << "You have forgotten the <<mksweep>> command  !"<< "\n";
+     di << "You have forgotten the <<mksweep>> command  !\n";
      return 1;
    }
   if (!strcmp(a[1],"-FR")) {
@@ -480,14 +540,14 @@ static Standard_Integer setsweep(Draw_Interpretor& di,
   else if (!strcmp(a[1],"-DX")) {
     if (n!=3) {
       //cout << "bad arguments !" << endl;
-      di << "bad arguments !" << "\n";
+      di << "bad arguments !\n";
       return 1;
     }
     TopoDS_Shape Surf;
     Surf = DBRep::Get(a[2],TopAbs_SHAPE);
     if (Surf.IsNull()) {
        //cout << a[2] <<"is not a shape !" << endl;
-       di << a[2] <<"is not a shape !" << "\n";
+       di << a[2] <<"is not a shape !\n";
       return 1;
     }
     Sweep->SetMode(Surf);
@@ -495,7 +555,7 @@ static Standard_Integer setsweep(Draw_Interpretor& di,
   else if (!strcmp(a[1],"-CN")) {
     if (n!=5) {
       //cout << "bad arguments !" << endl;
-      di << "bad arguments !" << "\n";
+      di << "bad arguments !\n";
       return 1;
     }
     gp_Dir D(Draw::Atof(a[2]), Draw::Atof(a[3]), Draw::Atof(a[4]));
@@ -504,7 +564,7 @@ static Standard_Integer setsweep(Draw_Interpretor& di,
   else if (!strcmp(a[1],"-FX")) {
     if ((n!=5)&&(n!=8)) {
       //cout << "bad arguments !" << endl;
-      di << "bad arguments !" << "\n";
+      di << "bad arguments !\n";
       return 1;
     }
     gp_Dir D(Draw::Atof(a[2]), Draw::Atof(a[3]), Draw::Atof(a[4]));
@@ -523,13 +583,13 @@ static Standard_Integer setsweep(Draw_Interpretor& di,
      if (n != 5)
        {
         //cout << "bad arguments !" << endl;
-        di << "bad arguments !" << "\n";
+        di << "bad arguments !\n";
         return 1; 
        }
      else
        {  
          TopoDS_Shape Guide = DBRep::Get(a[2],TopAbs_WIRE);
-          Standard_Integer CurvilinearEquivalence = Draw::Atoi(a[3]);
+          Standard_Boolean CurvilinearEquivalence = Draw::Atoi(a[3]) != 0;
           Standard_Integer KeepContact = Draw::Atoi(a[4]);
           Sweep->SetMode(TopoDS::Wire(Guide),
                          CurvilinearEquivalence,
@@ -539,7 +599,7 @@ static Standard_Integer setsweep(Draw_Interpretor& di,
  
   else {
     //cout << "The option "<< a[1] << " is unknown !" << endl;
-    di << "The option "<< a[1] << " is unknown !" << "\n";
+    di << "The option "<< a[1] << " is unknown !\n";
     return 1;
   }
   return 0;
@@ -558,17 +618,17 @@ static Standard_Integer addsweep(Draw_Interpretor& di,
     //cout << "        with the spine" <<endl;
     //cout << "   -R : the wire have to be rotated to assume orthogonality"<<endl;
     //cout << "        with the spine's tangent" << endl;
-    di << "addsweep wire/vertex [Vertex] [-T] [-R] [u0 v0 u1 v1 [...[uN vN]]] : options are :" << "\n";
-    di << "   -T : the wire/vertex have to be translated to assume contact"<< "\n";
-    di << "        with the spine" <<"\n";
-    di << "   -R : the wire have to be rotated to assume orthogonality"<<"\n";
-    di << "        with the spine's tangent" << "\n";
+    di << "addsweep wire/vertex [Vertex] [-T] [-R] [u0 v0 u1 v1 [...[uN vN]]] : options are :\n";
+    di << "   -T : the wire/vertex have to be translated to assume contact\n";
+    di << "        with the spine\n";
+    di << "   -R : the wire have to be rotated to assume orthogonality\n";
+    di << "        with the spine's tangent\n";
     return 0;
   }
 
   if (Sweep ==0) {
     //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
-    di << "You have forgotten the <<mksweep>> command  !"<< "\n";
+    di << "You have forgotten the <<mksweep>> command  !\n";
     return 1;
   }
 
@@ -577,11 +637,12 @@ static Standard_Integer addsweep(Draw_Interpretor& di,
   Handle(Law_Interpol) thelaw;
 
   Section = DBRep::Get(a[1], TopAbs_SHAPE);
-  if (Section.ShapeType() != TopAbs_WIRE &&
-      Section.ShapeType() != TopAbs_VERTEX)
+  if (Section.IsNull() ||
+      (Section.ShapeType() != TopAbs_WIRE &&
+       Section.ShapeType() != TopAbs_VERTEX))
     {
       //cout << a[1] <<"is not a wire and is not a vertex!" << endl;
-      di << a[1] <<"is not a wire and is not a vertex!" << "\n";
+      di << a[1] <<" is not a wire and is not a vertex!\n";
       return 1;
     }
 
@@ -659,7 +720,7 @@ static Standard_Integer deletesweep(Draw_Interpretor& di,
 //  Section = TopoDS::Wire(DBRep::Get(a[1],TopAbs_SHAPE));
   if (Section.IsNull()) {
     //cout << a[1] <<"is not a wire !" << endl;
-    di << a[1] <<"is not a wire !" << "\n";
+    di << a[1] <<"is not a wire !\n";
     return 1;
   }  
 
@@ -683,27 +744,27 @@ static Standard_Integer buildsweep(Draw_Interpretor& di,
     //cout << "   -R : Discontinuities are treated like Round Corner" << endl;
     //cout << "        Treatement is Intersect and Fill" << endl;
     //cout << "   -S : To build a Solid" << endl;
-    di << "build sweep result [-M/-C/-R] [-S] [tol] : options are" << "\n";
-    di << "   -M : Discontinuities are treated by Modfication of"<< "\n"; 
-    di << "        the sweeping mode : it is the default" <<"\n";
-    di << "   -C : Discontinuities are treated like Right Corner" << "\n";
-    di << "        Treatement is Extent && Intersect" << "\n";
-    di << "   -R : Discontinuities are treated like Round Corner" << "\n";
-    di << "        Treatement is Intersect and Fill" << "\n";
-    di << "   -S : To build a Solid" << "\n";
+    di << "build sweep result [-M/-C/-R] [-S] [tol] : options are\n";
+    di << "   -M : Discontinuities are treated by Modfication of\n"; 
+    di << "        the sweeping mode : it is the default\n";
+    di << "   -C : Discontinuities are treated like Right Corner\n";
+    di << "        Treatement is Extent && Intersect\n";
+    di << "   -R : Discontinuities are treated like Round Corner\n";
+    di << "        Treatement is Intersect and Fill\n";
+    di << "   -S : To build a Solid\n";
     return 0;
   }
 
   Standard_Boolean mksolid = Standard_False;
   if (Sweep ==0) {
     //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
-    di << "You have forgotten the <<mksweep>> command  !"<< "\n";
+    di << "You have forgotten the <<mksweep>> command  !\n";
     return 1;
   }
 
   if (!Sweep->IsReady()) {
     //cout << "You have forgotten the <<addsweep>> command  !"<< endl;
-    di << "You have forgotten the <<addsweep>> command  !"<< "\n";
+    di << "You have forgotten the <<addsweep>> command  !\n";
     return 1;
   }
 
@@ -730,15 +791,15 @@ static Standard_Integer buildsweep(Draw_Interpretor& di,
   Sweep->Build();
   if (!Sweep->IsDone()) {
     //cout << "Buildsweep : Not Done" << endl;
-    di << "Buildsweep : Not Done" << "\n";
+    di << "Buildsweep : Not Done\n";
     BRepBuilderAPI_PipeError Stat = Sweep->GetStatus(); 
     if (Stat == BRepBuilderAPI_PlaneNotIntersectGuide) {
       //cout << "Buildsweep : One Plane not intersect the guide" << endl;
-      di << "Buildsweep : One Plane not intersect the guide" << "\n";
+      di << "Buildsweep : One Plane not intersect the guide\n";
     }
     if (Stat == BRepBuilderAPI_ImpossibleContact) {
       //cout << "BuildSweep : One section can not be in contact with the guide" << endl;
-      di << "BuildSweep : One section can not be in contact with the guide" << "\n";
+      di << "BuildSweep : One section can not be in contact with the guide\n";
     }
     return 1;
   }
@@ -747,7 +808,7 @@ static Standard_Integer buildsweep(Draw_Interpretor& di,
       Standard_Boolean B;
       B = Sweep->MakeSolid();
       //if (!B) cout << " BuildSweep : It is impossible to make a solid !" << endl;
-      if (!B) di << " BuildSweep : It is impossible to make a solid !" << "\n";
+      if (!B) di << " BuildSweep : It is impossible to make a solid !\n";
     }
     result = Sweep->Shape();
     DBRep::Set(a[1],result);
@@ -756,6 +817,45 @@ static Standard_Integer buildsweep(Draw_Interpretor& di,
   return 0;
 }
 
+//=======================================================================
+//function : gensweep
+//purpose  : returns generated shape for subshape of a section of sweep
+//           Sweep must be done previously
+//=======================================================================
+static Standard_Integer gensweep(Draw_Interpretor&,
+                                 Standard_Integer n, const char** a)
+{
+  if (n != 3)
+  {
+    cout<<"Usage: gensweep res subshape_of_profile, sweep must be done"<<endl;
+    return 1;
+  }
+  if (!Sweep->IsDone())
+  {
+    cout<<"Sweep is not done"<<endl;
+    return 1;
+  }
+  TopoDS_Shape aShape = DBRep::Get(a[2]);
+  if (aShape.IsNull())
+  {
+    cout<<"Null subshape"<<endl;
+    return 1;
+  }
+  TopTools_ListOfShape Shells = Sweep->Generated(aShape);
+  TopoDS_Compound aCompound;
+  BRep_Builder BB;
+  BB.MakeCompound(aCompound);
+  TopTools_ListIteratorOfListOfShape itsh(Shells);
+  for (; itsh.More(); itsh.Next())
+  {
+    const TopoDS_Shape& aShell = itsh.Value();
+    BB.Add(aCompound, aShell);
+  }
+
+  DBRep::Set(a[1], aCompound);
+  return 0;
+}
+
 //=======================================================================
 //  simulsweep
 //=======================================================================
@@ -766,13 +866,13 @@ static Standard_Integer simulsweep(Draw_Interpretor& di,
   
   if (Sweep ==0) {
     //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
-    di << "You have forgotten the <<mksweep>> command  !"<< "\n";
+    di << "You have forgotten the <<mksweep>> command  !\n";
     return 1;
   }
   
   if (!Sweep->IsReady()) {
     //cout << "You have forgotten the <<addsweep>> command  !"<< endl;
-    di << "You have forgotten the <<addsweep>> command  !"<< "\n";
+    di << "You have forgotten the <<addsweep>> command  !\n";
     return 1;
   }
   
@@ -875,6 +975,8 @@ void  BRepTest::SweepCommands(Draw_Interpretor& theCommands)
   theCommands.Add("thrusections", "thrusections [-N] result issolid isruled shape1 shape2 [..shape..], the option -N means no check on wires, shapes must be wires or vertices (only first or last)",
                  __FILE__,thrusections,g);
 
+  theCommands.Add("genthrus", "genthrus res subshape_of_profile",
+                 __FILE__,genthrus,g);
   
   theCommands.Add("mksweep", "mksweep wire",
                  __FILE__,mksweep,g);
@@ -890,9 +992,12 @@ void  BRepTest::SweepCommands(Draw_Interpretor& theCommands)
                  "deletesweep wire, To delete a section",
                  __FILE__,deletesweep,g);
 
-  theCommands.Add("buildsweep", "builsweep [r] [option] [Tol] , no args to get help"
+ theCommands.Add("buildsweep", "builsweep [r] [option] [Tol] , no args to get help",
                  __FILE__,buildsweep,g);
 
+ theCommands.Add("gensweep", "gensweep res subshape_of_profile",
+                 __FILE__,gensweep,g);
+
   theCommands.Add("simulsweep", "simulsweep r [n] [option]"
                  __FILE__,simulsweep,g);
   theCommands.Add("geompipe", "geompipe r spineedge profileedge radius [byACR [byrotate]]"