src/BRepTest/BRepTest_SweepCommands.cxx

   1 // Created on: 1993-07-22
   2 // Created by: Remi LEQUETTE
   3 // Copyright (c) 1993-1999 Matra Datavision
   4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
   5 //
   6 // The content of this file is subject to the Open CASCADE Technology Public
   7 // License Version 6.5 (the "License"). You may not use the content of this file
   8 // except in compliance with the License. Please obtain a copy of the License
   9 // at http://www.opencascade.org and read it completely before using this file.
  10 //
  11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
  12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
  13 //
  14 // The Original Code and all software distributed under the License is
  15 // distributed on an "AS IS" basis, without warranty of any kind, and the
  16 // Initial Developer hereby disclaims all such warranties, including without
  17 // limitation, any warranties of merchantability, fitness for a particular
  18 // purpose or non-infringement. Please see the License for the specific terms
  19 // and conditions governing the rights and limitations under the License.
  20
  21
  22
  23
  24 #include <BRepTest.hxx>
  25 #include <DBRep.hxx>
  26 #include <Draw_Interpretor.hxx>
  27 #include <Draw_Appli.hxx>
  28
  29 #include <BRepFill.hxx>
  30 #include <BRepBuilderAPI_PipeError.hxx>
  31 #include <BRepFill_Generator.hxx>
  32 #include <BRepPrimAPI_MakePrism.hxx>
  33 #include <BRepPrimAPI_MakeRevol.hxx>
  34 #include <BRepOffsetAPI_MakePipe.hxx>
  35 #include <BRepOffsetAPI_MakeEvolved.hxx>
  36 #include <BRepOffsetAPI_ThruSections.hxx>
  37 #include <BRepOffsetAPI_MakePipeShell.hxx>
  38 #include <BRepOffsetAPI_MiddlePath.hxx>
  39
  40 #include <BRepLib_MakeWire.hxx>
  41 #include <TopoDS.hxx>
  42 #include <TopTools_ListIteratorOfListOfShape.hxx>
  43 #include <TopExp_Explorer.hxx>
  44
  45 #include <Precision.hxx>
  46 #include <Law_Interpol.hxx>
  47 #include <gp_Ax1.hxx>
  48 #include <gp_Ax2.hxx>
  49 #include <gp_Pnt2d.hxx>
  50 #include <TColgp_Array1OfPnt2d.hxx>
  51
  52 static BRepOffsetAPI_MakePipeShell* Sweep= 0;
  53
  54 //#ifdef WNT
  55 #include <stdio.h>
  56 #include <Geom_Curve.hxx>
  57 #include <GeomAdaptor_HCurve.hxx>
  58 #include <GeomFill_Pipe.hxx>
  59 #include <Geom_Surface.hxx>
  60 #include <BRepBuilderAPI_MakeFace.hxx>
  61 #include <BRep_Tool.hxx>
  62 #include <gp_Pnt.hxx>
  63 //#endi#include <gp_Vec.hxx>
  64 #include <gp_Vec.hxx>
  65 #include <Geom_Circle.hxx>
  66 #include <gp_Ax2.hxx>
  67
  68 //=======================================================================
  69 // prism
  70 //=======================================================================
  71
  72 static Standard_Integer prism(Draw_Interpretor& , Standard_Integer n, const char** a)
  73 {
  74   if (n < 6) return 1;
  75
  76   TopoDS_Shape base = DBRep::Get(a[2]);
  77   if (base.IsNull()) return 1;
  78
  79   gp_Vec V(atof(a[3]),atof(a[4]),atof(a[5]));
  80
  81   Standard_Boolean copy = Standard_False;
  82   Standard_Boolean inf  = Standard_False;
  83   Standard_Boolean sinf = Standard_False;
  84
  85   if (n > 6) {
  86     copy = (*a[6] == 'c') || (*a[6] == 'C');
  87     inf  = (*a[6] == 'i') || (*a[6] == 'I');
  88     sinf = (*a[6] == 's') || (*a[6] == 'S');
  89   }
  90
  91   TopoDS_Shape res;
  92
  93   if (inf || sinf)
  94     res = BRepPrimAPI_MakePrism(base,gp_Dir(V),inf);
  95   else
  96     res = BRepPrimAPI_MakePrism(base,V,copy);
  97
  98  DBRep::Set(a[1],res);
  99
 100   return 0;
 101 }
 102
 103
 104 //=======================================================================
 105 // revol
 106 //=======================================================================
 107
 108 static Standard_Integer revol(Draw_Interpretor& ,
 109                               Standard_Integer n, const char** a)
 110 {
 111   if (n < 10) return 1;
 112
 113   TopoDS_Shape base = DBRep::Get(a[2]);
 114   if (base.IsNull()) return 1;
 115
 116   gp_Pnt P(atof(a[3]),atof(a[4]),atof(a[5]));
 117   gp_Dir D(atof(a[6]),atof(a[7]),atof(a[8]));
 118   gp_Ax1 A(P,D);
 119
 120   Standard_Real angle = atof(a[9]) * (M_PI / 180.0);
 121
 122   Standard_Boolean copy = n > 10;
 123
 124   TopoDS_Shape res = BRepPrimAPI_MakeRevol(base,A,angle,copy);
 125
 126   DBRep::Set(a[1],res);
 127
 128   return 0;
 129 }
 130
 131
 132 //=======================================================================
 133 // pipe
 134 //=======================================================================
 135
 136 static Standard_Integer pipe(Draw_Interpretor& ,
 137                              Standard_Integer n, const char** a)
 138 {
 139   if ( n < 4) return 1;
 140
 141   TopoDS_Shape Spine = DBRep::Get(a[2],TopAbs_WIRE);
 142   if ( Spine.IsNull()) return 1;
 143
 144   TopoDS_Shape Profile = DBRep::Get(a[3]);
 145   if ( Profile.IsNull()) return 1;
 146
 147   TopoDS_Shape S = BRepOffsetAPI_MakePipe(TopoDS::Wire(Spine),Profile);
 148
 149   DBRep::Set(a[1],S);
 150
 151   return 0;
 152 }
 153 //=======================================================================
 154
 155 static Standard_Integer geompipe(Draw_Interpretor& ,
 156                              Standard_Integer n, const char** a)
 157 {
 158   TopoDS_Shape Spine = DBRep::Get(a[2],TopAbs_EDGE);
 159   if ( Spine.IsNull()) return 1;
 160   if ( n < 5) return 1;
 161   TopoDS_Shape Profile = DBRep::Get(a[3],TopAbs_EDGE);
 162   if ( Profile.IsNull()) return 1;
 163   Standard_Real aSpFirst,aSpLast,aPrFirst,aPrLast;
 164   Handle(Geom_Curve) SpineCurve = BRep_Tool::Curve(TopoDS::Edge(Spine),aSpFirst,aSpLast);
 165   Handle(Geom_Curve) ProfileCurve = BRep_Tool::Curve(TopoDS::Edge(Profile),aPrFirst,aPrLast);
 166   Handle(GeomAdaptor_HCurve) aAdaptCurve = new GeomAdaptor_HCurve(SpineCurve,aSpFirst,aSpLast);
 167   Standard_Boolean ByACR = Standard_False;
 168   Standard_Boolean rotate = Standard_False;
 169   Standard_Real Radius = atof(a[4]);
 170   gp_Pnt ctr;
 171   gp_Vec norm;
 172   ProfileCurve->D1(aSpFirst,ctr,norm);
 173   gp_Vec xAxisStart(ctr,SpineCurve->Value(aSpFirst));
 174   gp_Ax2 aAx2Start(ctr,norm,xAxisStart);
 175   Handle(Geom_Circle) cStart=new Geom_Circle(aAx2Start,Radius);
 176   Standard_Integer k =5;
 177   if(n > k)
 178     ByACR = (atoi(a[k++]) ==1);
 179   if(n > k)
 180     rotate = (atoi(a[k++])==1);
 181   GeomFill_Pipe aPipe(ProfileCurve,aAdaptCurve,cStart,ByACR,rotate);
 182   aPipe.Perform(Standard_True);
 183   Handle(Geom_Surface) Sur=aPipe.Surface();
 184   TopoDS_Face F;
 185   if(!Sur.IsNull())
 186     F = BRepBuilderAPI_MakeFace(Sur, Precision::Confusion());
 187   DBRep::Set(a[1],F);
 188   return 0;
 189 }
 190
 191 //=======================================================================
 192 //function : evolved
 193 //purpose  :
 194 //=======================================================================
 195
 196 Standard_Integer evolved(Draw_Interpretor& di, Standard_Integer n, const char** a)
 197 {
 198   if ( n == 1) {
 199     //cout << " 1) evolved result base profil : "<< endl;
 200     //cout << "        The relative position of the profil on the base" << endl;
 201     //cout << "        is given in the referencial axis. " << endl;
 202     //cout << " 2) evolved result base profil o : "<< endl;
 203     //cout << "        This position is automatically computed." << endl;
 204     di << " 1) evolved result base profil : "<< "\n";
 205     di << "        The relative position of the profil on the base" << "\n";
 206     di << "        is given in the referencial axis. " << "\n";
 207     di << " 2) evolved result base profil o : "<< "\n";
 208     di << "        This position is automatically computed." << "\n";
 209     return 0;
 210   }
 211
 212   if ( n < 4 ) return 1;
 213   Standard_Boolean IsAFace = Standard_False;
 214   Standard_Boolean Solid   = (!strcmp(a[0],"evolvedsolid"));
 215
 216
 217
 218   TopoDS_Shape Base = DBRep::Get(a[2],TopAbs_WIRE,Standard_False);
 219   if ( Base.IsNull()) {
 220     Base = DBRep::Get(a[2],TopAbs_FACE,Standard_False);
 221     IsAFace = Standard_True;
 222   }
 223   if ( Base.IsNull()) return 1;
 224
 225   TopoDS_Shape InpuTShape(DBRep::Get(a[3],TopAbs_WIRE,Standard_False));
 226   TopoDS_Wire Prof = TopoDS::Wire(InpuTShape);
 227 //  TopoDS_Wire Prof =
 228 //    TopoDS::Wire(DBRep::Get(a[3],TopAbs_WIRE,Standard_False));
 229   if ( Prof.IsNull()) return 1;
 230
 231   if (IsAFace) {
 232     TopoDS_Shape Volevo
 233       = BRepOffsetAPI_MakeEvolved(TopoDS::Face(Base),Prof,GeomAbs_Arc,n == 4,Solid);
 234     DBRep::Set(a[1],Volevo);
 235   }
 236   else {
 237     TopoDS_Shape Volevo
 238       = BRepOffsetAPI_MakeEvolved(TopoDS::Wire(Base),Prof,GeomAbs_Arc,n == 4,Solid);
 239     DBRep::Set(a[1],Volevo);
 240   }
 241
 242   return 0;
 243 }
 244
 245
 246 //=======================================================================
 247 //function : pruled
 248 //purpose  :
 249 //=======================================================================
 250
 251 static Standard_Integer pruled(Draw_Interpretor& ,
 252                                Standard_Integer n, const char** a)
 253 {
 254   if ( n != 4) return 1;
 255
 256   Standard_Boolean YaWIRE = Standard_False;
 257   TopoDS_Shape S1 = DBRep::Get(a[2],TopAbs_EDGE);
 258   if ( S1.IsNull()) {
 259     S1 = DBRep::Get(a[2],TopAbs_WIRE);
 260     if (S1.IsNull()) return 1;
 261     YaWIRE = Standard_True;
 262   }
 263
 264   TopoDS_Shape S2 = DBRep::Get(a[3],TopAbs_EDGE);
 265   if ( S2.IsNull()) {
 266     S2 = DBRep::Get(a[3],TopAbs_WIRE);
 267     if ( S2.IsNull()) return 1;
 268     if (!YaWIRE) {
 269       S1 = BRepLib_MakeWire(TopoDS::Edge(S1));
 270       YaWIRE = Standard_True;
 271     }
 272   }
 273   else if ( YaWIRE) {
 274     S2 = BRepLib_MakeWire(TopoDS::Edge(S2));
 275   }
 276
 277   TopoDS_Shape Result;
 278   if ( YaWIRE) {
 279     Result = BRepFill::Shell(TopoDS::Wire(S1),TopoDS::Wire(S2));
 280   }
 281   else {
 282     Result = BRepFill::Face(TopoDS::Edge(S1),TopoDS::Edge(S2));
 283   }
 284
 285   DBRep::Set(a[1],Result);
 286   return 0;
 287 }
 288
 289
 290 //=======================================================================
 291 //function : gener
 292 //purpose  : Create a surface between generating wires
 293 //=======================================================================
 294
 295 Standard_Integer gener(Draw_Interpretor&, Standard_Integer n, const char** a)
 296 {
 297   if ( n < 4) return 1;
 298
 299   TopoDS_Shape Shape;
 300
 301   BRepFill_Generator Generator;
 302
 303   for ( Standard_Integer i = 2; i<= n-1 ; i++) {
 304     Shape = DBRep::Get(a[i],TopAbs_WIRE);
 305     if ( Shape.IsNull())
 306       return 1;
 307
 308     Generator.AddWire(TopoDS::Wire(Shape));
 309   }
 310
 311   Generator.Perform();
 312
 313   TopoDS_Shell Shell = Generator.Shell();
 314
 315   DBRep::Set(a[1], Shell);
 316
 317
 318   return 0;
 319 }
 320
 321
 322 //=======================================================================
 323 //function : thrusections
 324 //purpose  :
 325 //=======================================================================
 326
 327 Standard_Integer thrusections(Draw_Interpretor&, Standard_Integer n, const char** a)
 328 {
 329   if (n<6) return 1;
 330
 331   Standard_Boolean check = Standard_True;
 332   Standard_Boolean samenumber = Standard_True;
 333   Standard_Integer index = 2;
 334     // Lecture option
 335   if (!strcmp(a[1],"-N")) {
 336     if (n<7) return 1;
 337     check = Standard_False;
 338     index++;
 339   }
 340
 341   TopoDS_Shape Shape;
 342
 343   Standard_Boolean issolid = ( atoi(a[index]) == 1 );
 344   Standard_Boolean isruled = ( atoi(a[index+1]) == 1 );
 345
 346   BRepOffsetAPI_ThruSections Generator(issolid,isruled);
 347
 348   Standard_Integer NbEdges = 0;
 349   Standard_Boolean IsFirstWire = Standard_False;
 350   for ( Standard_Integer i = index+2; i<= n-1 ; i++) {
 351     Standard_Boolean IsWire = Standard_True;
 352     Shape = DBRep::Get(a[i], TopAbs_WIRE);
 353     if (!Shape.IsNull())
 354       {
 355         Generator.AddWire(TopoDS::Wire(Shape));
 356         if (!IsFirstWire)
 357           IsFirstWire = Standard_True;
 358         else
 359           IsFirstWire = Standard_False;
 360       }
 361     else
 362       {
 363         Shape = DBRep::Get(a[i], TopAbs_VERTEX);
 364         IsWire = Standard_False;
 365         if (!Shape.IsNull())
 366           Generator.AddVertex(TopoDS::Vertex(Shape));
 367         else
 368           return 1;
 369       }
 370
 371     Standard_Integer cpt = 0;
 372     TopExp_Explorer PE;
 373     for (PE.Init(Shape, TopAbs_EDGE); PE.More(); PE.Next()) {
 374       cpt++;
 375     }
 376     if (IsFirstWire)
 377       NbEdges = cpt;
 378     else
 379       if (IsWire && cpt != NbEdges)
 380         samenumber = Standard_False;
 381
 382   }
 383
 384   check = (check || !samenumber);
 385   Generator.CheckCompatibility(check);
 386
 387   Generator.Build();
 388
 389   TopoDS_Shape Shell = Generator.Shape();
 390
 391   DBRep::Set(a[index-1], Shell);
 392
 393   return 0;
 394 }
 395
 396 //=======================================================================
 397 //  mksweep
 398 //=======================================================================
 399 static Standard_Integer mksweep(Draw_Interpretor& ,
 400                              Standard_Integer n, const char** a)
 401 {
 402   if ( n != 2) return 1;
 403   TopoDS_Shape Spine = DBRep::Get(a[1],TopAbs_WIRE);
 404   if ( Spine.IsNull()) return 1;
 405   if (Sweep !=0)  {
 406     delete Sweep;
 407     Sweep = 0;
 408   }
 409   Sweep = new BRepOffsetAPI_MakePipeShell(TopoDS::Wire(Spine));
 410   return 0;
 411 }
 412
 413 //=======================================================================
 414 //  setsweep
 415 //=======================================================================
 416 static Standard_Integer setsweep(Draw_Interpretor& di,
 417                                  Standard_Integer n, const char** a)
 418 {
 419   if ( n == 1) {
 420     //cout << "setsweep options [arg1 [arg2 [...]]] : options are :" << endl;
 421     //cout << "   -FR : Tangent and Normal are given by Frenet trihedron" <<endl;
 422     //cout << "   -CF : Tangente is given by Frenet," << endl;
 423     //cout << "         the Normal is computed to minimize the torsion " << endl;
 424     //cout << "   -DX Surf : Tangent and Normal are given by Darboux trihedron,"
 425     //  <<endl;
 426     //cout << "       Surf have to be a shell or a face" <<endl;
 427     //cout << "   -CN dx dy dz : BiNormal is given by dx dy dz" << endl;
 428     //cout << "   -FX Tx Ty TZ [Nx Ny Nz] : Tangent and Normal are fixed" <<endl;
 429     //cout << "   -G guide  0|1(ACR|Plan)  0|1(contact|no contact) : with guide"<<endl;
 430     di << "setsweep options [arg1 [arg2 [...]]] : options are :" << "\n";
 431     di << "   -FR : Tangent and Normal are given by Frenet trihedron" <<"\n";
 432     di << "   -CF : Tangente is given by Frenet," << "\n";
 433     di << "         the Normal is computed to minimize the torsion " << "\n";
 434     di << "   -DX Surf : Tangent and Normal are given by Darboux trihedron," <<"\n";
 435     di << "       Surf have to be a shell or a face" <<"\n";
 436     di << "   -CN dx dy dz : BiNormal is given by dx dy dz" << "\n";
 437     di << "   -FX Tx Ty TZ [Nx Ny Nz] : Tangent and Normal are fixed" <<"\n";
 438     di << "   -G guide  0|1(ACR|Plan)  0|1(contact|no contact) : with guide"<<"\n";
 439     return 0;
 440   }
 441
 442    if (Sweep ==0) {
 443      //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
 444      di << "You have forgotten the <<mksweep>> command  !"<< "\n";
 445      return 1;
 446    }
 447   if (!strcmp(a[1],"-FR")) {
 448     Sweep->SetMode(Standard_True);
 449   }
 450   else if (!strcmp(a[1],"-CF")) {
 451     Sweep->SetMode(Standard_False);
 452   }
 453   else if (!strcmp(a[1],"-DX")) {
 454     if (n!=3) {
 455       //cout << "bad arguments !" << endl;
 456       di << "bad arguments !" << "\n";
 457       return 1;
 458     }
 459     TopoDS_Shape Surf;
 460     Surf = DBRep::Get(a[2],TopAbs_SHAPE);
 461     if (Surf.IsNull()) {
 462        //cout << a[2] <<"is not a shape !" << endl;
 463        di << a[2] <<"is not a shape !" << "\n";
 464       return 1;
 465     }
 466     Sweep->SetMode(Surf);
 467   }
 468   else if (!strcmp(a[1],"-CN")) {
 469     if (n!=5) {
 470       //cout << "bad arguments !" << endl;
 471       di << "bad arguments !" << "\n";
 472       return 1;
 473     }
 474     gp_Dir D(atof(a[2]), atof(a[3]), atof(a[4]));
 475     Sweep->SetMode(D);;
 476   }
 477   else if (!strcmp(a[1],"-FX")) {
 478     if ((n!=5)&&(n!=8)) {
 479       //cout << "bad arguments !" << endl;
 480       di << "bad arguments !" << "\n";
 481       return 1;
 482     }
 483     gp_Dir D(atof(a[2]), atof(a[3]), atof(a[4]));
 484     if (n==8) {
 485       gp_Dir DN(atof(a[5]), atof(a[6]), atof(a[7]));
 486       gp_Ax2 Axe(gp_Pnt(0., 0., 0.), D, DN);
 487       Sweep->SetMode(Axe);
 488     }
 489     else {
 490       gp_Ax2 Axe(gp_Pnt(0., 0., 0.), D);
 491       Sweep->SetMode(Axe);
 492     }
 493   }
 494   else if (!strcmp(a[1],"-G"))  // contour guide
 495     {
 496      if (n != 5)
 497        {
 498          //cout << "bad arguments !" << endl;
 499          di << "bad arguments !" << "\n";
 500          return 1;
 501        }
 502      else
 503         {
 504           TopoDS_Shape Guide = DBRep::Get(a[2],TopAbs_WIRE);
 505           Sweep->SetMode(TopoDS::Wire(Guide), atoi(a[3]), atoi(a[4]));
 506         }
 507     }
 508
 509   else {
 510     //cout << "The option "<< a[1] << " is unknown !" << endl;
 511     di << "The option "<< a[1] << " is unknown !" << "\n";
 512     return 1;
 513   }
 514   return 0;
 515 }
 516
 517
 518 //=======================================================================
 519 //  addsweep
 520 //=======================================================================
 521 static Standard_Integer addsweep(Draw_Interpretor& di,
 522                              Standard_Integer n, const char** a)
 523 {
 524   if ( n == 1) {
 525     //cout << "addsweep wire/vertex [Vertex] [-T] [-R] [u0 v0 u1 v1 [...[uN vN]]] : options are :" << endl;
 526     //cout << "   -T : the wire/vertex have to be translated to assume contact"<< endl;
 527     //cout << "        with the spine" <<endl;
 528     //cout << "   -R : the wire have to be rotated to assume orthogonality"<<endl;
 529     //cout << "        with the spine's tangent" << endl;
 530     di << "addsweep wire/vertex [Vertex] [-T] [-R] [u0 v0 u1 v1 [...[uN vN]]] : options are :" << "\n";
 531     di << "   -T : the wire/vertex have to be translated to assume contact"<< "\n";
 532     di << "        with the spine" <<"\n";
 533     di << "   -R : the wire have to be rotated to assume orthogonality"<<"\n";
 534     di << "        with the spine's tangent" << "\n";
 535     return 0;
 536   }
 537
 538   if (Sweep ==0) {
 539     //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
 540     di << "You have forgotten the <<mksweep>> command  !"<< "\n";
 541     return 1;
 542   }
 543
 544   TopoDS_Shape  Section;
 545   TopoDS_Vertex Vertex;
 546   Handle(Law_Interpol) thelaw;
 547
 548   Section = DBRep::Get(a[1], TopAbs_SHAPE);
 549   if (Section.ShapeType() != TopAbs_WIRE &&
 550       Section.ShapeType() != TopAbs_VERTEX)
 551     {
 552       //cout << a[1] <<"is not a wire and is not a vertex!" << endl;
 553       di << a[1] <<"is not a wire and is not a vertex!" << "\n";
 554       return 1;
 555     }
 556
 557   Standard_Boolean HasVertex=Standard_False,
 558                    isT=Standard_False,
 559                    isR=Standard_False;
 560
 561   if (n > 2) {
 562     Standard_Integer cur = 2;
 563     // Reading of Vertex
 564     TopoDS_Shape InputVertex(DBRep::Get(a[cur],TopAbs_VERTEX));
 565     Vertex = TopoDS::Vertex(InputVertex);
 566 //    Vertex = TopoDS::Vertex(DBRep::Get(a[cur],TopAbs_VERTEX));
 567     if (!Vertex.IsNull()) {
 568       cur++;
 569       HasVertex = Standard_True;
 570     }
 571
 572     // Reading of the translation option
 573     if ((n>cur) && !strcmp(a[cur],"-T")) {
 574       cur++;
 575       isT = Standard_True;
 576     }
 577
 578     // Reading of the rotation option
 579     if ((n>cur) && !strcmp(a[cur],"-R")) {
 580       cur++;
 581       isR = Standard_True;
 582     }
 583
 584     // law ?
 585     if (n>cur) {
 586       Standard_Integer nbreal = n-cur;
 587       if ( (nbreal < 4) || (nbreal % 2 != 0) ) {
 588         //cout << "bad arguments ! :" <<a[cur] << endl;
 589         di << "bad arguments ! :" <<a[cur] << "\n";
 590       } else { //law of interpolation
 591         Standard_Integer ii, L= nbreal/2;
 592         TColgp_Array1OfPnt2d ParAndRad(1, L);
 593         for (ii=1; ii<=L; ii++, cur+=2) {
 594            ParAndRad(ii).SetX(atof(a[cur]));
 595            ParAndRad(ii).SetY(atof(a[cur+1]));
 596          }
 597         thelaw = new (Law_Interpol) ();
 598         thelaw->Set(ParAndRad,
 599                    Abs(ParAndRad(1).Y() - ParAndRad(L).Y()) < Precision::Confusion());
 600       }
 601     }
 602   }
 603
 604   if (thelaw.IsNull()) {
 605     if (HasVertex) Sweep->Add(Section, Vertex, isT, isR);
 606     else           Sweep->Add(Section, isT, isR);
 607   }
 608   else {
 609     if (HasVertex) Sweep->SetLaw(Section, thelaw, Vertex, isT, isR);
 610     else           Sweep->SetLaw(Section, thelaw, isT, isR);
 611   }
 612
 613   return 0;
 614 }
 615
 616 //=======================================================================
 617 //  deletesweep
 618 //=======================================================================
 619 static Standard_Integer deletesweep(Draw_Interpretor& di,
 620                                     Standard_Integer n, const char** a)
 621 {
 622   if ( n != 2) {
 623     return 1;
 624   }
 625   TopoDS_Wire Section;
 626   TopoDS_Shape InputShape(DBRep::Get(a[1],TopAbs_SHAPE));
 627   Section = TopoDS::Wire(InputShape);
 628 //  Section = TopoDS::Wire(DBRep::Get(a[1],TopAbs_SHAPE));
 629   if (Section.IsNull()) {
 630     //cout << a[1] <<"is not a wire !" << endl;
 631     di << a[1] <<"is not a wire !" << "\n";
 632     return 1;
 633   }
 634
 635   Sweep->Delete(Section);
 636
 637   return 0;
 638 }
 639
 640 //=======================================================================
 641 //  buildsweep
 642 //=======================================================================
 643 static Standard_Integer buildsweep(Draw_Interpretor& di,
 644                                    Standard_Integer n, const char** a)
 645 {
 646   if ( n == 1) {
 647     //cout << "build sweep result [-M/-C/-R] [-S] [tol] : options are" << endl;
 648     //cout << "   -M : Discontinuities are treated by Modfication of"<< endl;
 649     //cout << "        the sweeping mode : it is the default" <<endl;
 650     //cout << "   -C : Discontinuities are treated like Right Corner" << endl;
 651     //cout << "        Treatement is Extent && Intersect" << endl;
 652     //cout << "   -R : Discontinuities are treated like Round Corner" << endl;
 653     //cout << "        Treatement is Intersect and Fill" << endl;
 654     //cout << "   -S : To build a Solid" << endl;
 655     di << "build sweep result [-M/-C/-R] [-S] [tol] : options are" << "\n";
 656     di << "   -M : Discontinuities are treated by Modfication of"<< "\n";
 657     di << "        the sweeping mode : it is the default" <<"\n";
 658     di << "   -C : Discontinuities are treated like Right Corner" << "\n";
 659     di << "        Treatement is Extent && Intersect" << "\n";
 660     di << "   -R : Discontinuities are treated like Round Corner" << "\n";
 661     di << "        Treatement is Intersect and Fill" << "\n";
 662     di << "   -S : To build a Solid" << "\n";
 663     return 0;
 664   }
 665
 666   Standard_Boolean mksolid = Standard_False;
 667   if (Sweep ==0) {
 668     //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
 669     di << "You have forgotten the <<mksweep>> command  !"<< "\n";
 670     return 1;
 671   }
 672
 673   if (!Sweep->IsReady()) {
 674     //cout << "You have forgotten the <<addsweep>> command  !"<< endl;
 675     di << "You have forgotten the <<addsweep>> command  !"<< "\n";
 676     return 1;
 677   }
 678
 679   TopoDS_Shape result;
 680   Standard_Integer cur=2;
 681   if (n>cur) {
 682     BRepBuilderAPI_TransitionMode Transition = BRepBuilderAPI_Transformed;
 683
 684     // Reading Transition
 685     if (!strcmp(a[cur],"-C")) {
 686       Transition = BRepBuilderAPI_RightCorner;
 687       cur++;
 688     }
 689     else if (!strcmp(a[cur],"-R")) {
 690       Transition = BRepBuilderAPI_RoundCorner;
 691       cur++;
 692     }
 693     Sweep->SetTransitionMode(Transition);
 694   }
 695   // Reading solid ?
 696   if ((n>cur) && (!strcmp(a[cur],"-S")) ) mksolid = Standard_True;
 697
 698   // Calcul le resultat
 699   Sweep->Build();
 700   if (!Sweep->IsDone()) {
 701     //cout << "Buildsweep : Not Done" << endl;
 702     di << "Buildsweep : Not Done" << "\n";
 703     BRepBuilderAPI_PipeError Stat = Sweep->GetStatus();
 704     if (Stat == BRepBuilderAPI_PlaneNotIntersectGuide) {
 705       //cout << "Buildsweep : One Plane not intersect the guide" << endl;
 706       di << "Buildsweep : One Plane not intersect the guide" << "\n";
 707     }
 708     if (Stat == BRepBuilderAPI_ImpossibleContact) {
 709       //cout << "BuildSweep : One section can not be in contact with the guide" << endl;
 710       di << "BuildSweep : One section can not be in contact with the guide" << "\n";
 711     }
 712     return 1;
 713   }
 714   else {
 715     if (mksolid) {
 716       Standard_Boolean B;
 717       B = Sweep->MakeSolid();
 718       //if (!B) cout << " BuildSweep : It is impossible to make a solid !" << endl;
 719       if (!B) di << " BuildSweep : It is impossible to make a solid !" << "\n";
 720     }
 721     result = Sweep->Shape();
 722     DBRep::Set(a[1],result);
 723   }
 724
 725   return 0;
 726 }
 727
 728 //=======================================================================
 729 //  simulsweep
 730 //=======================================================================
 731 static Standard_Integer simulsweep(Draw_Interpretor& di,
 732                                    Standard_Integer n, const char** a)
 733 {
 734   if ( (n!=3) && (n!=4) ) return 1;
 735
 736   if (Sweep ==0) {
 737     //cout << "You have forgotten the <<mksweep>> command  !"<< endl;
 738     di << "You have forgotten the <<mksweep>> command  !"<< "\n";
 739     return 1;
 740   }
 741
 742   if (!Sweep->IsReady()) {
 743     //cout << "You have forgotten the <<addsweep>> command  !"<< endl;
 744     di << "You have forgotten the <<addsweep>> command  !"<< "\n";
 745     return 1;
 746   }
 747
 748   char name[100];
 749   TopTools_ListOfShape List;
 750   TopTools_ListIteratorOfListOfShape it;
 751   Standard_Integer N, ii;
 752   N = atoi(a[2]);
 753
 754   if (n>3) {
 755     BRepBuilderAPI_TransitionMode Transition = BRepBuilderAPI_Transformed;
 756     // Lecture Transition
 757     if (!strcmp(a[3],"-C")) {
 758       Transition = BRepBuilderAPI_RightCorner;
 759     }
 760     else if (!strcmp(a[3],"-R")) {
 761       Transition = BRepBuilderAPI_RoundCorner;
 762     }
 763     Sweep->SetTransitionMode(Transition);
 764   }
 765
 766   // Calculate the result
 767   Sweep->Simulate(N, List);
 768   for (ii=1, it.Initialize(List); it.More(); it.Next(), ii++) {
 769     sprintf(name,"%s_%d",a[1],ii);
 770     DBRep::Set(name, it.Value());
 771   }
 772
 773   return 0;
 774 }
 775
 776 //=======================================================================
 777 //  middlepath
 778 //=======================================================================
 779 static Standard_Integer middlepath(Draw_Interpretor& di,
 780                                    Standard_Integer n, const char** a)
 781 {
 782   if (n < 5) return 1;
 783
 784   TopoDS_Shape aShape = DBRep::Get(a[2]);
 785   if (aShape.IsNull()) return 1;
 786
 787   TopoDS_Shape StartShape = DBRep::Get(a[3]);
 788   if (StartShape.IsNull()) return 1;
 789
 790   TopoDS_Shape EndShape   = DBRep::Get(a[4]);
 791   if (EndShape.IsNull()) return 1;
 792
 793   BRepOffsetAPI_MiddlePath Builder(aShape, StartShape, EndShape);
 794   Builder.Build();
 795
 796   TopoDS_Shape Result = Builder.Shape();
 797   DBRep::Set(a[1], Result);
 798
 799   return 0;
 800 }
 801
 802 //=======================================================================
 803 //function : SweepCommands
 804 //purpose  :
 805 //=======================================================================
 806
 807 void  BRepTest::SweepCommands(Draw_Interpretor& theCommands)
 808 {
 809   static Standard_Boolean done = Standard_False;
 810   if (done) return;
 811   done = Standard_True;
 812
 813   DBRep::BasicCommands(theCommands);
 814
 815   const char* g = "Sweep commands";
 816
 817   theCommands.Add("prism",
 818                   "prism result base dx dy dz [Copy | Inf | Seminf]",
 819                   __FILE__,prism,g);
 820
 821   theCommands.Add("revol",
 822                   "revol result base px py pz dx dy dz angle [Copy]",
 823                   __FILE__,revol,g);
 824
 825   theCommands.Add("pipe",
 826                   "pipe result Wire_spine Profile",
 827                   __FILE__,pipe,g);
 828
 829   theCommands.Add("evolved",
 830                   "evolved , no args to get help",
 831                   __FILE__,evolved,g);
 832
 833   theCommands.Add("evolvedsolid",
 834                   "evolved , no args to get help",
 835                   __FILE__,evolved,g);
 836
 837   theCommands.Add("pruled",
 838                   "pruled result Edge1/Wire1 Edge2/Wire2",
 839                   __FILE__,pruled,g);
 840
 841   theCommands.Add("gener", "gener result wire1 wire2 [..wire..]",
 842                   __FILE__,gener,g);
 843
 844   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)",
 845                   __FILE__,thrusections,g);
 846
 847
 848   theCommands.Add("mksweep", "mksweep wire",
 849                   __FILE__,mksweep,g);
 850
 851   theCommands.Add("setsweep", "setsweep  no args to get help",
 852                   __FILE__,setsweep,g);
 853
 854   theCommands.Add("addsweep",
 855                   "addsweep wire [vertex] [-M ] [-C] [auxiilaryshape]:no args to get help",
 856                   __FILE__,addsweep,g);
 857
 858  theCommands.Add("deletesweep",
 859                   "deletesweep wire, To delete a section",
 860                   __FILE__,deletesweep,g);
 861
 862   theCommands.Add("buildsweep", "builsweep [r] [option] [Tol] , no args to get help"
 863                   __FILE__,buildsweep,g);
 864
 865   theCommands.Add("simulsweep", "simulsweep r [n] [option]"
 866                   __FILE__,simulsweep,g);
 867   theCommands.Add("geompipe", "geompipe r spineedge profileedge radius [byACR [byrotate]]"
 868                   __FILE__,geompipe,g);
 869
 870   theCommands.Add("middlepath", "middlepath res shape startshape endshape",
 871                   __FILE__,middlepath,g);
 872 }
 873