src/BRepTest/BRepTest_BasicCommands.cxx

   1 // Created on: 1994-12-13
   2 // Created by: Jacques GOUSSARD
   3 // Copyright (c) 1994-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 #include <Standard_Stream.hxx>
  23 #include <Standard_Macro.hxx>
  24
  25 #include <BRepTest.hxx>
  26
  27 #include <DBRep.hxx>
  28 #include <Draw_Appli.hxx>
  29 #include <Draw_Interpretor.hxx>
  30 #include <Draw_Box.hxx>
  31
  32 #include <BRepBuilderAPI.hxx>
  33 #include <BRepBuilderAPI_FindPlane.hxx>
  34 #include <BRepBuilderAPI_Copy.hxx>
  35 #include <BRepBuilderAPI_Transform.hxx>
  36 #include <BRepBuilderAPI_GTransform.hxx>
  37 #include <BRepBuilderAPI_NurbsConvert.hxx>
  38 #include <gp_Ax2.hxx>
  39 #include <gp_Mat.hxx>
  40 #include <gp_GTrsf.hxx>
  41 #include <BRepOffsetAPI_NormalProjection.hxx>
  42 #include <BRepLib.hxx>
  43 #include <BRep_Builder.hxx>
  44 #include <BRepBndLib.hxx>
  45 #include <Bnd_Box.hxx>
  46 #include <TopExp_Explorer.hxx>
  47 #include <TopoDS.hxx>
  48 #include <BRepTools_WireExplorer.hxx>
  49
  50 #include <GCPnts_QuasiUniformAbscissa.hxx>
  51 #include <Geom2dAdaptor_Curve.hxx>
  52 #include <GeomAdaptor_Curve.hxx>
  53 #include <ProjLib_ComputeApproxOnPolarSurface.hxx>
  54 #include <DrawTrSurf.hxx>
  55
  56
  57 #include <Draw_Segment3D.hxx>
  58 #include <Draw_Marker3D.hxx>
  59 #include <Draw_MarkerShape.hxx>
  60
  61 #include <stdio.h>
  62
  63 Standard_IMPORT Draw_Viewer dout;
  64
  65
  66
  67 //=======================================================================
  68 // addpcurve
  69 //=======================================================================
  70
  71 static Standard_Integer addpcurve(Draw_Interpretor& , Standard_Integer n, const char** a)
  72 {
  73   if(n < 4) return 1;
  74   TopoDS_Shape E = DBRep::Get(a[1]);
  75   if (E.IsNull()) return 1;
  76   Handle(Geom2d_Curve) PC = DrawTrSurf::GetCurve2d(a[2]);
  77   TopoDS_Shape F = DBRep::Get(a[3]);
  78   Standard_Real tol = 1.e-7;
  79   if (n > 4) {
  80     tol = atof(a[4]);
  81   }
  82   BRep_Builder BB;
  83   BB.UpdateEdge(TopoDS::Edge(E), PC, TopoDS::Face(F),tol);
  84   DBRep::Set(a[1], E);
  85   return 0;
  86 }
  87
  88
  89 //=======================================================================
  90 // transform
  91 //=======================================================================
  92
  93 static Standard_Integer transform(Draw_Interpretor& di,Standard_Integer n,const char** a)
  94 {
  95   if (n <= 1) return 1;
  96
  97   gp_Trsf T;
  98   Standard_Integer last = n;
  99
 100   if (!strcmp(a[0],"reset")) {
 101   }
 102   else if (!strcmp(a[0],"tmove")) {
 103     if (n < 3) return 1;
 104     TopoDS_Shape SL = DBRep::Get(a[n-1]);
 105     if (SL.IsNull()) return 0;
 106     T = SL.Location().Transformation();
 107     last = n-1;
 108   }
 109   else if (!strcmp(a[0],"ttranslate")) {
 110     if (n < 5) return 1;
 111     T.SetTranslation(gp_Vec(atof(a[n-3]),atof(a[n-2]),atof(a[n-1])));
 112     last = n-3;
 113   }
 114   else if (!strcmp(a[0],"trotate")) {
 115     if (n < 9) return 1;
 116     T.SetRotation(gp_Ax1(gp_Pnt(atof(a[n-7]),atof(a[n-6]),atof(a[n-5])),
 117                          gp_Vec(atof(a[n-4]),atof(a[n-3]),atof(a[n-2]))),
 118                   atof(a[n-1])* (M_PI / 180.0));
 119     last = n-7;
 120   }
 121   else if (!strcmp(a[0],"tmirror")) {
 122     if (n < 8) return 1;
 123     T.SetMirror(gp_Ax2(gp_Pnt(atof(a[n-6]),atof(a[n-5]),atof(a[n-4])),
 124                        gp_Vec(atof(a[n-3]),atof(a[n-2]),atof(a[n-1]))));
 125     last = n-6;
 126   }
 127   else if (!strcmp(a[0],"tscale")) {
 128     if (n < 6) return 1;
 129     T.SetScale(gp_Pnt(atof(a[n-4]),atof(a[n-3]),atof(a[n-2])),atof(a[n-1]));
 130     last = n-4;
 131
 132   }
 133
 134   if (T.Form() == gp_Identity) {
 135     TopLoc_Location L;
 136     for (Standard_Integer i = 1; i < last; i++) {
 137       TopoDS_Shape S = DBRep::Get(a[i]);
 138       if (S.IsNull()) {
 139         //cout << a[i] << " is not a valid shape" << endl;
 140         di << a[i] << " is not a valid shape" << "\n";
 141       }
 142       else {
 143         DBRep::Set(a[i],S.Located(L));
 144       }
 145     }
 146   }
 147   else {
 148     BRepBuilderAPI_Transform trf(T);
 149     for (Standard_Integer i = 1; i < last; i++) {
 150       TopoDS_Shape S = DBRep::Get(a[i]);
 151       if (S.IsNull()) {
 152         //cout << a[i] << " is not a valid shape" << endl;
 153         di << a[i] << " is not a valid shape" << "\n";
 154       }
 155       else {
 156         trf.Perform(S);
 157         if (trf.IsDone()){
 158           DBRep::Set(a[i],trf.Shape());
 159         }
 160         else {
 161           return 1;
 162         }
 163       }
 164     }
 165   }
 166   return 0;
 167 }
 168
 169
 170 ///=======================================================================
 171 // gtransform
 172 //=======================================================================
 173
 174 static Standard_Integer deform(Draw_Interpretor& di,Standard_Integer n,const char** a)
 175 {
 176   if (n <= 1) return 1;
 177
 178   Standard_Integer last = n;
 179
 180   gp_Trsf T;
 181   gp_GTrsf GT(T);
 182
 183 //  gp_Mat rot(atof(a[last-3]),0,0,0,atof(a[last-2]),0,0,0,atof(a[last-1]));
 184   gp_Mat rot(atof(a[3]),0,0,0,atof(a[4]),0,0,0,atof(a[5]));
 185   GT.SetVectorialPart(rot);
 186   last -= 3;
 187   BRepBuilderAPI_GTransform gtrf(GT);
 188   BRepBuilderAPI_NurbsConvert nbscv;
 189   //  for (Standard_Integer i = 1; i < last; i++) {
 190   //    TopoDS_Shape S = DBRep::Get(a[i]);
 191   TopoDS_Shape S = DBRep::Get(a[2]);
 192   if (S.IsNull()) {
 193     //cout << a[2] << " is not a valid shape" << endl;
 194     di << a[2] << " is not a valid shape" << "\n";
 195   }
 196   else {
 197     gtrf.Perform(S);
 198     if (gtrf.IsDone()){
 199       DBRep::Set(a[1],gtrf.Shape());
 200     }
 201     else {
 202       return 1;
 203     }
 204   }
 205
 206   return 0;
 207 }
 208
 209 //=======================================================================
 210 // tcopy
 211 //=======================================================================
 212
 213 static Standard_Integer tcopy(Draw_Interpretor& di,Standard_Integer n,const char** a)
 214 {
 215   Standard_Boolean copyGeom = Standard_True;
 216   Standard_Integer iFirst = 1; // index of first shape argument
 217
 218   if (n > 1 && a[1][0] == '-' && a[1][1] == 'n' )
 219   {
 220     copyGeom = Standard_False;
 221     iFirst = 2;
 222   }
 223
 224   if (n < 3 || (n - iFirst) % 2) {
 225     cout << "Use: " << a[0] << " [-n(ogeom)] shape1 copy1 [shape2 copy2 [...]]" << endl;
 226     cout << "Option -n forbids copying of geometry (it will be shared)" << endl;
 227     return 1;
 228   }
 229
 230   BRepBuilderAPI_Copy cop;
 231   Standard_Integer nbPairs = (n - iFirst) / 2;
 232   for (Standard_Integer i=0; i < nbPairs; i++) {
 233     cop.Perform(DBRep::Get(a[i+iFirst]), copyGeom);
 234     DBRep::Set(a[i+iFirst+1],cop.Shape());
 235     di << a[i+iFirst+1] << " ";
 236   }
 237   return 0;
 238 }
 239
 240
 241 //=======================================================================
 242 // NurbsConvert
 243 //=======================================================================
 244
 245 static Standard_Integer nurbsconvert(Draw_Interpretor& di,Standard_Integer n,const char** a)
 246 {
 247   if (n < 3) return 1;
 248   if ((n-1)%2 != 0) return 1;
 249   BRepBuilderAPI_NurbsConvert nbscv;
 250   for (Standard_Integer i=0; i<(n-1)/2; i++) {
 251     TopoDS_Shape S = DBRep::Get(a[2*i+2]);
 252     if (S.IsNull()) {
 253       //cout << a[2*i+2] << " is not a valid shape" << endl;
 254       di << a[2*i+2] << " is not a valid shape" << "\n";
 255     }
 256     else {
 257       nbscv.Perform(S);
 258       if (nbscv.IsDone()){
 259         DBRep::Set(a[2*i+1],nbscv.Shape());
 260       }
 261       else {
 262         return 1;
 263       }
 264     }
 265   }
 266
 267   return 0;
 268
 269 }
 270
 271 //=======================================================================
 272 // make a 3D edge curve
 273 //=======================================================================
 274
 275 static Standard_Integer mkedgecurve (Draw_Interpretor& ,Standard_Integer n,const char** a)
 276 {
 277
 278   Standard_Boolean CurveDone ;
 279
 280   if (n < 3) return 1;
 281   Standard_Real Tolerance = atof(a[2]) ;
 282
 283   TopoDS_Shape S = DBRep::Get(a[1]);
 284
 285   if (S.IsNull()) return 1;
 286
 287    CurveDone =
 288    BRepLib::BuildCurves3d(S,
 289                           Tolerance) ;
 290    return 0 ;
 291 }
 292
 293 //=======================================================================
 294 // sameparameter
 295 //=======================================================================
 296
 297 static Standard_Integer sameparameter(Draw_Interpretor& ,Standard_Integer n,const char** a)
 298 {
 299   if (n < 2) return 1;
 300   Standard_Real tol = 1.e-7;
 301   TopoDS_Shape S = DBRep::Get(a[1]);
 302   if (S.IsNull()) return 1;
 303   Standard_Boolean force  = !strcmp(a[0],"fsameparameter");
 304   if (n == 3) tol = atof(a[2]);
 305
 306   BRepLib::SameParameter(S,tol,force);
 307
 308   DBRep::Set(a[1],S);
 309   return 0;
 310 }
 311 //=======================================================================
 312 //function : updatetol
 313 //purpose  :
 314 //=======================================================================
 315 static Standard_Integer updatetol(Draw_Interpretor& ,Standard_Integer n,const char** a)
 316 {
 317   if (n < 2) return 1;
 318
 319   TopoDS_Shape S = DBRep::Get(a[1]);
 320   if (S.IsNull()) return 1;
 321
 322   if (n==2) BRepLib::UpdateTolerances(S);
 323   else BRepLib::UpdateTolerances(S,Standard_True);
 324   DBRep::Set(a[1],S);
 325   return 0;
 326
 327 }
 328
 329 //=======================================================================
 330 //function : OrienSolid
 331 //purpose  :
 332 //=======================================================================
 333 static Standard_Integer orientsolid(Draw_Interpretor& ,Standard_Integer n,const char** a)
 334 {
 335   if (n < 2) return 1;
 336
 337   TopoDS_Shape S = DBRep::Get(a[1]);
 338   if (S.IsNull()) return 1;
 339   if (S.ShapeType()!=TopAbs_SOLID) return 1;
 340
 341   BRepLib::OrientClosedSolid(TopoDS::Solid(S));
 342
 343   DBRep::Set(a[1],S);
 344   return 0;
 345
 346 }
 347
 348 //=======================================================================
 349 //function : boundingstr
 350 //purpose  :
 351 //=======================================================================
 352 static Standard_Integer boundingstr(Draw_Interpretor& di,Standard_Integer n,const char** a)
 353 {
 354   if (n < 2) return 1;
 355   TopoDS_Shape S = DBRep::Get(a[1]);
 356   if (S.IsNull()) return 1;
 357   Bnd_Box B;
 358   BRepBndLib::Add(S,B);
 359   Standard_Real axmin,aymin,azmin,axmax,aymax,azmax;
 360   B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
 361   di << axmin<<" "<< aymin<<" "<< azmin<<" "<< axmax<<" "<< aymax<<" "<< azmax;
 362   if (n >= 8) {
 363     Draw::Set(a[2],axmin) ;
 364     Draw::Set(a[3],aymin) ;
 365     Draw::Set(a[4],azmin) ;
 366     Draw::Set(a[5],axmax) ;
 367     Draw::Set(a[6],aymax) ;
 368     Draw::Set(a[7],azmax) ;
 369   }
 370   return 0;
 371 }
 372
 373 //=======================================================================
 374 //function : bounding
 375 //purpose  :
 376 //=======================================================================
 377 static Standard_Integer bounding(Draw_Interpretor& di,Standard_Integer n,const char** a)
 378 {
 379   if (n < 2) return 1;
 380   Standard_Real axmin,aymin,azmin,axmax,aymax,azmax;
 381   Bnd_Box B; Handle(Draw_Box) DB;
 382
 383   if (n == 2) {
 384     TopoDS_Shape S = DBRep::Get(a[1]);
 385     if (S.IsNull()) return 1;
 386     BRepBndLib::Add(S,B);
 387     B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
 388     DB = new Draw_Box(gp_Pnt(axmin,aymin,azmin),gp_Pnt(axmax,aymax,azmax),Draw_orange);
 389     dout<<DB;
 390     di << axmin<<" "<< aymin<<" "<< azmin<<" "<< axmax<<" "<< aymax<<" "<< azmax;
 391   }
 392   else if (n == 7) {
 393     axmin=atof(a[1]);
 394     aymin=atof(a[2]);
 395     azmin=atof(a[3]);
 396     axmax=atof(a[4]);
 397     aymax=atof(a[5]);
 398     azmax=atof(a[6]);
 399     DB = new Draw_Box(gp_Pnt(axmin,aymin,azmin),gp_Pnt(axmax,aymax,azmax),Draw_orange);
 400     dout<<DB;
 401   }
 402   return 0;
 403 }
 404 //=======================================================================
 405 //function : findplane
 406 //purpose  :
 407 //=======================================================================
 408 static Standard_Integer findplane(Draw_Interpretor& di,Standard_Integer n,const char** a)
 409 {
 410   if (n < 3) return 1;
 411   TopoDS_Shape S = DBRep::Get(a[1]);
 412   if (S.IsNull()) return 1;
 413   Standard_Real tolerance = 1.0e-5 ;
 414   BRepBuilderAPI_FindPlane a_plane_finder(S,
 415                                    tolerance) ;
 416   if (a_plane_finder.Found()) {
 417     //cout << " a plane is found "   ;
 418     di << " a plane is found \n";
 419     DrawTrSurf::Set(a[2],a_plane_finder.Plane()) ;
 420   }
 421   return 0 ;
 422 }
 423 //=======================================================================
 424 //function : precision
 425 //purpose  :
 426 //=======================================================================
 427
 428 static Standard_Integer precision(Draw_Interpretor& di,Standard_Integer n,const char** a)
 429 {
 430   n--;
 431
 432   if ( n == 0) {
 433     //cout << " Current Precision = " << BRepBuilderAPI::Precision() << endl;
 434     di << " Current Precision = " << BRepBuilderAPI::Precision() << "\n";
 435   }
 436   else {
 437     BRepBuilderAPI::Precision(atof(a[1]));
 438   }
 439   return 0;
 440 }
 441
 442
 443 //=======================================================================
 444 //function : reperage shape (Int lin Shape) + pointe double click   + maxtol
 445 //purpose  :
 446 //=======================================================================
 447 #include <IntCurvesFace_ShapeIntersector.hxx>
 448 #include <gp_Lin.hxx>
 449
 450 static Standard_Integer reperageshape(Draw_Interpretor& di, Standard_Integer narg , const char** a)
 451 {
 452   Standard_Integer details=0;
 453   if(narg<2) return 1;
 454   if(narg==3) details=1;
 455   const char *id1 = a[1];
 456   TopoDS_Shape TheShape1 = DBRep::Get(id1);
 457
 458   //cout << "Pick positions with button "<<endl;
 459   di << "Pick positions with button "<< "\n";
 460   Standard_Integer id,X,Y,b;
 461   gp_Trsf T;
 462   gp_Pnt P1,P2;
 463   dout.Select(id,X,Y,b);
 464
 465   dout.GetTrsf(id,T);
 466   T.Invert();
 467   Standard_Real z = dout.Zoom(id);
 468   P2.SetCoord((Standard_Real)X /z,(Standard_Real)Y /z, 0.0);
 469   P2.Transform(T);
 470   P1.SetCoord((Standard_Real)X /z,(Standard_Real)Y /z,-1.0);
 471   P1.Transform(T);
 472
 473
 474   gp_Ax1 Axe(P1,gp_Vec(P1,P2));
 475   IntCurvesFace_ShapeIntersector Inter;
 476   Inter.Load(TheShape1,1e-7);
 477
 478   Inter.Perform(Axe,-RealLast(),RealLast());
 479
 480   //cout<<"\n --> ";
 481   di <<"\n --> ";
 482   if(Inter.NbPnt()) {
 483     for(Standard_Integer i=1; i<=Inter.NbPnt(); i++) {
 484       Standard_Integer numface=1;
 485       TopExp_Explorer ExF;
 486       for(ExF.Init(TheShape1,TopAbs_FACE);
 487           ExF.More();
 488           ExF.Next(),numface++) {
 489         TopoDS_Face Face=TopoDS::Face(ExF.Current());
 490         if(Face.IsEqual(Inter.Face(i))) {
 491           //cout<<" "<<a[1]<<"_"<<numface;
 492           di<<" "<<a[1]<<"_"<<numface;
 493           continue;
 494         }
 495       }
 496       const gp_Pnt& P = Inter.Pnt(i);
 497       Standard_Real PMin = Inter.WParameter(i);
 498       if(details) {
 499         //cout<<" w:"<<PMin<<endl;
 500         di<<" w:"<<PMin<< "\n";
 501       }
 502       if(Inter.Transition(i) == IntCurveSurface_In) {
 503         if(Inter.State(i) == TopAbs_IN) {
 504           Handle(Draw_Marker3D) p = new Draw_Marker3D(P, Draw_Square, Draw_rouge,2);
 505           dout << p;   dout.Flush();
 506         }
 507         else if(Inter.State(i) == TopAbs_ON) {
 508           Handle(Draw_Marker3D) p = new Draw_Marker3D(P, Draw_Square, Draw_vert,2);
 509           dout << p;   dout.Flush();
 510         }
 511       }
 512       else {
 513         if(Inter.Transition(i) == IntCurveSurface_Out) {
 514           if(Inter.State(i) == TopAbs_IN) {
 515             Handle(Draw_Marker3D) p = new Draw_Marker3D(P, Draw_X, Draw_rouge,2);
 516             dout << p;   dout.Flush();
 517           }
 518           else if(Inter.State(i) == TopAbs_ON) {
 519             Handle(Draw_Marker3D) p = new Draw_Marker3D(P, Draw_X, Draw_vert,2);
 520             dout << p;   dout.Flush();
 521           }
 522         }
 523       }
 524     }
 525   }
 526   //cout<<endl;
 527   di << "\n";
 528   return(0);
 529 }
 530
 531
 532 static Standard_Integer maxtolerance(Draw_Interpretor& theCommands,
 533                                      Standard_Integer n, const char** a) {
 534   if(n<2) return(1);
 535   TopoDS_Shape TheShape = DBRep::Get(a[1]);
 536   if(TheShape.IsNull()) return(1);
 537
 538   Standard_Real T,TMF,TME,TMV,TmF,TmE,TmV;
 539   Standard_Integer nbF,nbE,nbV;
 540   TMF=TME=TMV=-RealLast();
 541   TmF=TmE=TmV=RealLast();
 542
 543   TopTools_MapOfShape mapS;
 544   mapS.Clear();
 545
 546   for(TopExp_Explorer ex(TheShape,TopAbs_FACE);
 547       ex.More();
 548       ex.Next()) {
 549     TopoDS_Face Face=TopoDS::Face(ex.Current());
 550     T=BRep_Tool::Tolerance(Face);
 551     if(T>TMF) TMF=T;
 552     if(T<TmF) TmF=T;
 553     mapS.Add(Face);
 554   }
 555
 556   nbF = mapS.Extent();
 557   mapS.Clear();
 558
 559   for(TopExp_Explorer ex(TheShape,TopAbs_EDGE);
 560       ex.More();
 561       ex.Next()) {
 562     TopoDS_Edge Edge=TopoDS::Edge(ex.Current());
 563     T=BRep_Tool::Tolerance(Edge);
 564     if(T>TME) TME=T;
 565     if(T<TmE) TmE=T;
 566     mapS.Add(Edge);
 567   }
 568
 569   nbE = mapS.Extent();
 570   mapS.Clear();
 571
 572   for(TopExp_Explorer ex(TheShape,TopAbs_VERTEX);
 573       ex.More();
 574       ex.Next()) {
 575     TopoDS_Vertex Vertex=TopoDS::Vertex(ex.Current());
 576     T=BRep_Tool::Tolerance(Vertex);
 577     if(T>TMV) TMV=T;
 578     if(T<TmV) TmV=T;
 579     mapS.Add(Vertex);
 580   }
 581
 582   nbV = mapS.Extent();
 583
 584 //#ifndef WNT
 585   Standard_SStream sss;
 586   sss << "\n## Tolerances on the shape " << a[1] << "  (nbFaces:" << nbF
 587       << "  nbEdges:" << nbE << " nbVtx:" << nbV << ")\n" ;
 588   sss.precision(5);
 589   sss.setf(ios::scientific);
 590   if(TmF<=TMF) sss << "\n    Face   : Min " << setw(8) << TmF <<"    Max  " << setw(8) << TMF << " \n ";
 591   if(TmE<=TME) sss << "\n    Edge   : Min " << setw(8) << TmE <<"    Max  " << setw(8) << TME << " \n ";
 592   if(TmV<=TMV) sss << "\n    Vertex : Min " << setw(8) << TmV <<"    Max  " << setw(8) << TMV << " \n ";
 593   theCommands << sss;
 594   //#endif*/
 595   return 0;
 596 }
 597
 598
 599 static Standard_Integer vecdc(Draw_Interpretor& di,Standard_Integer ,const char** ) {
 600   //cout << "Pick positions with button "<<endl;
 601   di << "Pick positions with button "<< "\n";
 602
 603   Standard_Integer id,X,Y,b;
 604   gp_Trsf T;
 605   gp_Pnt P1,P2,PP1,PP2;
 606
 607   //-----------------------------------------------------------
 608   dout.Select(id,X,Y,b);    dout.GetTrsf(id,T);
 609   T.Invert();
 610   Standard_Real z = dout.Zoom(id);
 611   P1.SetCoord((Standard_Real)X /z,(Standard_Real)Y /z,0.0);
 612   P1.Transform(T);
 613
 614   dout.Select(id,X,Y,b);  dout.GetTrsf(id,T);
 615   T.Invert();  z = dout.Zoom(id);
 616
 617   P2.SetCoord((Standard_Real)X /z,(Standard_Real)Y /z,0.0);
 618   P2.Transform(T);
 619   Standard_Real xa,ya,za;
 620   if(Abs(P1.X())>Abs(P2.X())) xa = P1.X(); else xa = P2.X();
 621   if(Abs(P1.Y())>Abs(P2.Y())) ya = P1.Y(); else ya = P2.Y();
 622   if(Abs(P1.Z())>Abs(P2.Z())) za = P1.Z(); else za = P2.Z();
 623   P1.SetCoord(xa,ya,za);
 624   Handle(Draw_Marker3D) D0 = new Draw_Marker3D(gp_Pnt(P1.X(),
 625                                                       P1.Y(),
 626                                                       P1.Z()),
 627                                                Draw_Square,Draw_blanc,1);
 628
 629   dout << D0;
 630   dout.Flush();
 631   //-----------------------------------------------------------
 632   dout.Select(id,X,Y,b);
 633   dout.GetTrsf(id,T);
 634   T.Invert();
 635   z = dout.Zoom(id);
 636   PP1.SetCoord((Standard_Real)X /z,(Standard_Real)Y /z,0.0);
 637   PP1.Transform(T);
 638   dout.Select(id,X,Y,b);
 639   dout.GetTrsf(id,T);
 640   T.Invert();
 641   z = dout.Zoom(id);
 642   PP2.SetCoord((Standard_Real)X /z,(Standard_Real)Y /z,0.0);
 643   PP2.Transform(T);
 644   if(Abs(PP1.X())>Abs(PP2.X())) xa = PP1.X(); else xa = PP2.X();
 645   if(Abs(PP1.Y())>Abs(PP2.Y())) ya = PP1.Y(); else ya = PP2.Y();
 646   if(Abs(PP1.Z())>Abs(PP2.Z())) za = PP1.Z(); else za = PP2.Z();
 647   PP1.SetCoord(xa,ya,za);
 648   Handle(Draw_Segment3D) d = new Draw_Segment3D(P1,PP1,Draw_blanc);
 649   dout << d;
 650   dout.Flush();
 651   //cout<<"\nttran   "<<PP1.X()-P1.X()<<" "<<PP1.Y()-P1.Y()<<" "<<PP1.Z()-P1.Z()<<endl;
 652   di <<"\nttran   "<<PP1.X()-P1.X()<<" "<<PP1.Y()-P1.Y()<<" "<<PP1.Z()-P1.Z()<< "\n";
 653
 654   static Standard_Integer nboxvecdp=0;
 655   //cout<<"\nbox  b"<<++nboxvecdp<<" "<<Min(P1.X(),PP1.X())<<" "<<Min(P1.Y(),PP1.Y())<<" "<<Min(PP1.Z(),P1.Z());
 656   //cout<<"  "<<Abs(PP1.X()-P1.X())<<" "<<Abs(PP1.Y()-P1.Y())<<" "<<Abs(PP1.Z()-P1.Z())<<endl;
 657
 658   //cout<<"\nDistance :"<<sqrt( (PP1.X()-P1.X())*(PP1.X()-P1.X())
 659         //                   +(PP1.Y()-P1.Y())*(PP1.Y()-P1.Y())
 660         //                   +(PP1.Z()-P1.Z())*(PP1.Z()-P1.Z()))<<endl;
 661
 662   di <<"\nbox  b"<<++nboxvecdp<<" "<<Min(P1.X(),PP1.X())<<" "<<Min(P1.Y(),PP1.Y())<<" "<<Min(PP1.Z(),P1.Z());
 663   di <<"  "<<Abs(PP1.X()-P1.X())<<" "<<Abs(PP1.Y()-P1.Y())<<" "<<Abs(PP1.Z()-P1.Z())<< "\n";
 664
 665   di <<"\nDistance :"<<sqrt( (PP1.X()-P1.X())*(PP1.X()-P1.X())
 666                              +(PP1.Y()-P1.Y())*(PP1.Y()-P1.Y())
 667                              +(PP1.Z()-P1.Z())*(PP1.Z()-P1.Z()))<< "\n";
 668   return(0);
 669 }
 670 //=======================================================================
 671 // nproject
 672 //=======================================================================
 673
 674 #include <TopTools_SequenceOfShape.hxx>
 675  static Standard_Integer nproject(Draw_Interpretor& di, Standard_Integer n, const char** a)
 676 {
 677   if ( n < 4) return 1;
 678   TopoDS_Shape InpShape;
 679   Standard_Integer arg = 2, i;
 680   TopTools_SequenceOfShape Args;
 681
 682   Standard_Real Tol = 1.e-4;
 683   Standard_Real Tol2d;
 684   Standard_Real MaxDistance = 1.e-3;
 685   GeomAbs_Shape Continuity = GeomAbs_C2;
 686   Standard_Integer MaxDeg = 14;
 687   Standard_Integer MaxSeg = 16;
 688
 689   while((n > arg) && !(InpShape = DBRep::Get(a[arg])).IsNull()){
 690     Args.Append(InpShape);
 691     arg++;
 692   }
 693   if(Args.Length() < 2) return 1;
 694
 695   BRepOffsetAPI_NormalProjection OrtProj(Args.Last());
 696
 697   for(i = 1; i < Args.Length(); i++)
 698     OrtProj.Add(Args(i));
 699
 700   if(n > arg)
 701     if (!strcmp(a[arg],"-g")) {
 702       OrtProj.SetLimit(Standard_False);
 703       arg++;
 704     }
 705
 706   if(n > arg)
 707     if (!strcmp(a[arg],"-d")) {
 708       arg++;
 709       if(n > arg)
 710         MaxDistance = atof(a[arg++]);
 711       OrtProj.SetMaxDistance(MaxDistance);
 712     }
 713   if(n > arg) {
 714     Tol = Max(atof(a[arg++]),1.e-10);
 715   }
 716
 717   if(n > arg) {
 718     if (atoi(a[arg]) == 0) Continuity = GeomAbs_C0;
 719     else if (atoi(a[arg]) == 1) Continuity = GeomAbs_C1;
 720     arg++;
 721   }
 722
 723
 724   if(n > arg) {
 725     MaxDeg = atoi(a[arg++]);
 726     if (MaxDeg<1 || MaxDeg>14) MaxDeg = 14;
 727   }
 728
 729   if(n > arg) MaxSeg = atoi(a[arg]);
 730
 731   Tol2d = Pow(Tol, 2./3);
 732
 733   OrtProj.SetParams(Tol, Tol2d, Continuity, MaxDeg, MaxSeg);
 734   OrtProj.Build();
 735   TopTools_ListOfShape Wire;
 736   Standard_Boolean IsWire=OrtProj.BuildWire(Wire);
 737   if (IsWire) {
 738     //cout << " BuildWire OK " << endl;
 739     di << " BuildWire OK " << "\n";
 740   }
 741   DBRep::Set(a[1], OrtProj.Shape());
 742   return 0;
 743 }
 744
 745 //==========================================================================
 746 //function : wexplo
 747 //           exploration of a wire
 748 //==========================================================================
 749 static Standard_Integer wexplo (Draw_Interpretor&,
 750                                 Standard_Integer argc, const char** argv)
 751 {
 752   char name[100];
 753   if (argc < 2) return 1;
 754
 755   TopoDS_Shape C1 = DBRep::Get (argv[1],TopAbs_WIRE);
 756   TopoDS_Shape C2 ;
 757
 758   if (argc > 2)  C2 = DBRep::Get (argv[2],TopAbs_FACE);
 759
 760   if (C1.IsNull()) return 1;
 761
 762   BRepTools_WireExplorer we;
 763   if (C2.IsNull()) we.Init(TopoDS::Wire(C1));
 764   else             we.Init(TopoDS::Wire(C1),TopoDS::Face(C2));
 765
 766   Standard_Integer k = 1;
 767   while (we.More()) {
 768     TopoDS_Edge E = we.Current();
 769     sprintf(name,"WEDGE_%d",k);
 770           DBRep::Set(name,E);
 771     we.Next();
 772     k++;
 773   }
 774
 775   return 0;
 776 }
 777
 778 static Standard_Integer scalexyz(Draw_Interpretor& di, Standard_Integer n, const char** a)
 779 {
 780   if (n < 6) return 1;
 781
 782   TopoDS_Shape aShapeBase = DBRep::Get(a[2]);
 783   if (aShapeBase.IsNull()) return 1;
 784
 785   Standard_Real aFactorX = atof(a[3]);
 786   Standard_Real aFactorY = atof(a[4]);
 787   Standard_Real aFactorZ = atof(a[5]);
 788
 789   gp_GTrsf aGTrsf;
 790   gp_Mat rot (aFactorX, 0, 0,
 791               0, aFactorY, 0,
 792               0, 0, aFactorZ);
 793   aGTrsf.SetVectorialPart(rot);
 794   BRepBuilderAPI_GTransform aBRepGTrsf (aShapeBase, aGTrsf, Standard_False);
 795   if (!aBRepGTrsf.IsDone())
 796     Standard_ConstructionError::Raise("Scaling not done");
 797   TopoDS_Shape Result = aBRepGTrsf.Shape();
 798
 799   DBRep::Set(a[1], Result);
 800   return 0;
 801 }
 802
 803 void  BRepTest::BasicCommands(Draw_Interpretor& theCommands)
 804 {
 805   static Standard_Boolean done = Standard_False;
 806   if (done) return;
 807   done = Standard_True;
 808
 809   DBRep::BasicCommands(theCommands);
 810
 811   const char* g = "TOPOLOGY Basic shape commands";
 812
 813   theCommands.Add("addpcurve",
 814                   "addpcurve edge 2dcurve face [tol (default 1.e-7)]",
 815                   __FILE__,
 816                   addpcurve,g);
 817
 818   theCommands.Add("reset",
 819                   "reset name1 name2..., remove location",
 820                   __FILE__,
 821                   transform,g);
 822
 823   theCommands.Add("tmove",
 824                   "tmove name1 name2 ...  name(location)",
 825                   __FILE__,
 826                   transform,g);
 827
 828   theCommands.Add("ttranslate",
 829                   "tmove name1 name2 ...  dx dy dz",
 830                   __FILE__,
 831                   transform,g);
 832
 833   theCommands.Add("trotate","trotate name1 name2 ... x y z dx dy dz angle",
 834                   __FILE__,
 835                   transform,g);
 836
 837   theCommands.Add("tmirror",
 838                   "tmirror name x y z dx dy dz",
 839                   __FILE__,
 840                   transform,g);
 841
 842   theCommands.Add("tscale",
 843                   "tscale name  x y z scale",
 844                   __FILE__,
 845                   transform,g);
 846
 847   theCommands.Add("tcopy",
 848     "tcopy [-n(ogeom)] name1 result1 [name2 result2 ...]",
 849                   __FILE__,
 850                   tcopy,g);
 851
 852   theCommands.Add("precision",
 853                   "precision [preci]",
 854                   __FILE__,
 855                   precision,g);
 856
 857   theCommands.Add("mkedgecurve",
 858                   "mkedgecurve name tolerance",
 859                   __FILE__,
 860                   mkedgecurve,g);
 861
 862   theCommands.Add("fsameparameter",
 863                   "fsameparameter shapename [tol (default 1.e-7)], \nforce sameparameter on all edges of the shape",
 864                   __FILE__,
 865                   sameparameter,g);
 866
 867   theCommands.Add("sameparameter",
 868                   "sameparameter shapename [tol (default 1.e-7)]",
 869                   __FILE__,
 870                   sameparameter,g);
 871
 872   theCommands.Add("updatetolerance",
 873                   "updatetolerance myShape [param] \n  if [param] is absent - not verify of face tolerance, else - perform it",
 874                   __FILE__,
 875                   updatetol,g);
 876
 877   theCommands.Add("solidorientation",
 878                   "orientsolid myClosedSolid",
 879                   __FILE__,
 880                   orientsolid,g);
 881
 882   theCommands.Add("bounding",
 883                   "bounding shape [ xmin ymin zmin xmax ymax zmax] ; draw bounds",
 884                   __FILE__,
 885                   bounding,g);
 886
 887   theCommands.Add("boundingstr",
 888                   "boundingstr shape [ xmin ymin zmin xmax ymax zmax] ; print bounding box",
 889                   __FILE__,
 890                   boundingstr,g);
 891
 892   theCommands.Add("nurbsconvert",
 893                   "nurbsconvert result name [result name]",
 894                   __FILE__,
 895                   nurbsconvert,g);
 896
 897   theCommands.Add("deform",
 898                   "deform newname name CoeffX CoeffY CoeffZ",
 899                   __FILE__,
 900                   deform,g);
 901
 902   theCommands.Add("findplane",
 903                   "findplane name planename ",
 904                   __FILE__,
 905                   findplane,g) ;
 906
 907   theCommands.Add("maxtolerance",
 908                   "maxtolerance shape ",
 909                   __FILE__,
 910                   maxtolerance,g) ;
 911
 912   theCommands.Add("reperageshape",
 913                   "reperage shape -> list of shape (result of interstion shape , line)",
 914                   __FILE__,
 915                   reperageshape,g) ;
 916
 917   theCommands.Add("vecdc",
 918                   "vecdc + Pointe double click ",
 919                   __FILE__,
 920                   vecdc,g) ;
 921
 922   theCommands.Add("nproject","nproject pj e1 e2 e3 ... surf -g -d [dmax] [Tol [continuity [maxdeg [maxseg]]]",
 923                   __FILE__,
 924                   nproject,g);
 925
 926   theCommands.Add("wexplo","wexplo wire [face] create WEDGE_i",
 927                   __FILE__,
 928                   wexplo,g);
 929
 930   theCommands.Add("scalexyz",
 931                   "scalexyz res shape factor_x factor_y factor_z",
 932                   __FILE__,
 933                   scalexyz, g);
 934 }