src/QABugs/QABugs_10.cxx

   1 // Created on: 2002-05-28
   2 // Created by: QA Admin
   3 // Copyright (c) 2002-2014 OPEN CASCADE SAS
   4 //
   5 // This file is part of Open CASCADE Technology software library.
   6 //
   7 // This library is free software; you can redistribute it and/or modify it under
   8 // the terms of the GNU Lesser General Public License version 2.1 as published
   9 // by the Free Software Foundation, with special exception defined in the file
  10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  11 // distribution for complete text of the license and disclaimer of any warranty.
  12 //
  13 // Alternatively, this file may be used under the terms of Open CASCADE
  14 // commercial license or contractual agreement.
  15
  16 #include <QABugs.hxx>
  17
  18 #include <Draw.hxx>
  19 #include <Draw_Interpretor.hxx>
  20 #include <DBRep.hxx>
  21 #include <DrawTrSurf.hxx>
  22 #include <AIS_InteractiveContext.hxx>
  23 #include <ViewerTest.hxx>
  24 #include <AIS_Shape.hxx>
  25 #include <TopoDS_Shape.hxx>
  26
  27 #include <BRepBuilderAPI_MakePolygon.hxx>
  28 #include <TopoDS_Face.hxx>
  29 #include <BRepBuilderAPI_MakeFace.hxx>
  30 #include <gp_Pnt.hxx>
  31 #include <gp_Dir.hxx>
  32 #include <gp_Ax1.hxx>
  33 #include <BRepPrimAPI_MakeRevol.hxx>
  34
  35 #include <BRepAlgoAPI_Fuse.hxx>
  36
  37 #include <GProp_GProps.hxx>
  38 #include <BRepGProp.hxx>
  39 #include <TopExp_Explorer.hxx>
  40 #include <TopoDS.hxx>
  41 #include <BRepMesh_IncrementalMesh.hxx>
  42 #include <BRep_Tool.hxx>
  43 #include <Poly_Triangulation.hxx>
  44 #include <TopExp.hxx>
  45 #include <ChFi3d_FilletShape.hxx>
  46 #include <BRepFilletAPI_MakeFillet.hxx>
  47 #include <Standard_ErrorHandler.hxx>
  48 #include <Geom_TrimmedCurve.hxx>
  49 #include <Geom_BSplineSurface.hxx>
  50
  51 #include <ShapeUpgrade_UnifySameDomain.hxx>
  52
  53 static Standard_Integer OCC426 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
  54 {
  55   if(argc != 8) {
  56     di << "Usage : " << argv[0] << " shape1 shape2 shape3 shape4 shape5 shape6 shape7\n";
  57     return 1;
  58   }
  59
  60   BRepBuilderAPI_MakePolygon W1;
  61   W1.Add(gp_Pnt(10, 0, 0));
  62   W1.Add(gp_Pnt(20, 0, 0));
  63   W1.Add(gp_Pnt(20, 0, 10));
  64   W1.Add(gp_Pnt(10, 0, 10));
  65   W1.Add(gp_Pnt(10, 0, 0));
  66
  67   Standard_Boolean OnlyPlane1 = Standard_False;
  68   TopoDS_Face F1 = BRepBuilderAPI_MakeFace(W1.Wire(), OnlyPlane1);
  69
  70   gp_Pnt P1(0, 0, 0);
  71   gp_Dir D1(0, 0, 30);
  72   gp_Ax1 A1(P1,D1);
  73   Standard_Real angle1 = 360 * (M_PI / 180.0);
  74   TopoDS_Shape rs1 = BRepPrimAPI_MakeRevol(F1, A1, angle1);
  75
  76   BRepBuilderAPI_MakePolygon W2;
  77   Standard_Real f1 =  7.0710678118654752440;
  78   Standard_Real f2 = 14.1421356237309504880;
  79   W2.Add(gp_Pnt(f1, f1, 10));
  80   W2.Add(gp_Pnt(f2, f2, 10));
  81   W2.Add(gp_Pnt(f2, f2, 20));
  82   W2.Add(gp_Pnt(f1, f1, 20));
  83   W2.Add(gp_Pnt(f1, f1, 10));
  84
  85   Standard_Boolean OnlyPlane2 = Standard_False;
  86   TopoDS_Face F2 = BRepBuilderAPI_MakeFace(W2.Wire(), OnlyPlane2);
  87
  88   gp_Pnt P2(0, 0, 0);
  89   gp_Dir D2(0, 0, 30);
  90   gp_Ax1 A2(P2,D2);
  91   Standard_Real angle2 = 270 * (M_PI / 180.0);
  92   TopoDS_Shape rs2 = BRepPrimAPI_MakeRevol(F2, A2, angle2);
  93
  94   BRepBuilderAPI_MakePolygon W3;
  95   W3.Add(gp_Pnt(10, 0, 20));
  96   W3.Add(gp_Pnt(20, 0, 20));
  97   W3.Add(gp_Pnt(20, 0, 30));
  98   W3.Add(gp_Pnt(10, 0, 30));
  99   W3.Add(gp_Pnt(10, 0, 20));
 100
 101   Standard_Boolean OnlyPlane3 = Standard_False;
 102   TopoDS_Face F3 = BRepBuilderAPI_MakeFace(W3.Wire(), OnlyPlane3);
 103
 104   gp_Pnt P3(0, 0, 0);
 105   gp_Dir D3(0, 0, 30);
 106   gp_Ax1 A3(P3,D3);
 107   Standard_Real angle3 = 360 * (M_PI / 180.0);
 108   TopoDS_Shape rs3 = BRepPrimAPI_MakeRevol(F3, A3, angle3);
 109
 110   di << "fuse32 = BRepAlgoAPI_Fuse(rs3, rs2)\n";
 111   di << "fuse321 = BRepAlgoAPI_Fuse(fuse32, rs1)\n";
 112   TopoDS_Shape fuse32 = BRepAlgoAPI_Fuse(rs3, rs2).Shape();
 113   TopoDS_Shape fuse321 = BRepAlgoAPI_Fuse(fuse32, rs1).Shape();
 114
 115   // unify the faces of the Fuse result
 116   ShapeUpgrade_UnifySameDomain anUnify(fuse321, Standard_True, Standard_True, Standard_True);
 117   anUnify.Build();
 118   const TopoDS_Shape& aFuseUnif = anUnify.Shape();
 119
 120   //Give the mass claculation of the shpae "aFuseUnif"
 121   GProp_GProps G;
 122   BRepGProp::VolumeProperties(aFuseUnif, G);
 123   di<<" \n";
 124   di<<"Mass: "<<G.Mass()<<"\n\n";
 125
 126   di << "Trianglating Faces .....\n";
 127   TopExp_Explorer ExpFace;
 128
 129   for (ExpFace.Init (aFuseUnif,TopAbs_FACE); ExpFace.More(); ExpFace.Next())
 130     {
 131       TopoDS_Face TopologicalFace = TopoDS::Face (ExpFace.Current());
 132       TopologicalFace.Orientation (TopAbs_FORWARD) ;
 133       BRepMesh_IncrementalMesh IM(TopologicalFace, 1);
 134       TopLoc_Location loc;
 135       Handle(Poly_Triangulation) facing = BRep_Tool::Triangulation(TopologicalFace, loc);
 136       if (facing.IsNull())
 137       {
 138         di << "Triangulation FAILED for this face\n";
 139         continue;
 140       }
 141       di << "No of Triangles = " << facing->NbTriangles() << "\n";
 142     }
 143   di<<"Triangulation of all Faces Completed. \n\n";
 144
 145   TopTools_IndexedDataMapOfShapeListOfShape edgemap;
 146   TopExp::MapShapesAndAncestors(aFuseUnif, TopAbs_EDGE, TopAbs_SOLID, edgemap);
 147   di << "No. of Edges: " << edgemap.Extent() << "\n";
 148   ChFi3d_FilletShape FShape = ChFi3d_Rational;
 149   BRepFilletAPI_MakeFillet blend(aFuseUnif,FShape);
 150   di << "Adding Edges ..... \n";
 151   for(int i = 1; i <= edgemap.Extent(); i++)
 152     {
 153       // std::cout << "Adding Edge : " << i << std::endl;
 154       TopoDS_Edge edg = TopoDS::Edge( edgemap.FindKey(i) );
 155       if(!edg.IsNull()) blend.Add(1, edg);
 156     }
 157   di << "All Edges added !  Now Building the Blend ... \n";
 158   di<<" \n";
 159   blend.Build();
 160
 161   //DBRep::Set ( argv[1], fuse321 );
 162   DBRep::Set ( argv[1], blend );
 163   DBRep::Set ( argv[2], rs1 );
 164   DBRep::Set ( argv[3], rs2 );
 165   DBRep::Set ( argv[4], rs3 );
 166   DBRep::Set ( argv[5], fuse32 );
 167   DBRep::Set ( argv[6], fuse321 );
 168   DBRep::Set ( argv[7], aFuseUnif );
 169
 170   return 0;
 171 }
 172
 173 #include <Geom_SurfaceOfRevolution.hxx>
 174 //=======================================================================
 175 //function : isPeriodic
 176 //purpose  :
 177 //=======================================================================
 178 static Standard_Integer isPeriodic(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 179 {
 180   try
 181     {
 182     OCC_CATCH_SIGNALS
 183       // 1. Verify amount of arguments of the command
 184       if (argc < 2) { di << "isperiodic FAULTY. Use : isperiodic surfaceOfRevolution"; return 0;}
 185       // 2. Retrieve surface
 186       Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
 187       if(aSurf.IsNull()) {di << "isperiodic FAULTY. argument of command is not a surface"; return 0;}
 188       Handle(Geom_SurfaceOfRevolution) aRevolSurf = Handle(Geom_SurfaceOfRevolution)::DownCast(aSurf);
 189       if(aRevolSurf.IsNull()) {di << "isperiodic FAULTY. argument of command is not a surface of revolution"; return 0;}
 190       // 3. Verify whether entry surface is u-periodic and v-periodic
 191       if(aRevolSurf->IsUPeriodic()) {di << "Surface is u-periodic \n";} else {di << "Surface is not u-periodic \n";}
 192       if(aRevolSurf->IsVPeriodic()) {di << "Surface is v-periodic \n";} else {di << "Surface is not v-periodic \n";}
 193     }
 194   catch (Standard_Failure const&) {di << "isperiodic Exception \n" ;return 0;}
 195
 196   return 0;
 197 }
 198
 199 #include <Precision.hxx>
 200 #include <Extrema_ExtPS.hxx>
 201 #include <GeomAdaptor_Surface.hxx>
 202 //=======================================================================
 203 //function : OCC486
 204 //purpose  :
 205 //=======================================================================
 206 static Standard_Integer OCC486(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 207 {
 208   try
 209     {
 210     OCC_CATCH_SIGNALS
 211       if (argc < 2) { di << "OCC486 FAULTY. Use : OCC486 surf x y z du dv"; return 1;}
 212
 213       Standard_Real du = 0;
 214       Standard_Real dv = 0;
 215
 216       Handle(Geom_Surface) GS;
 217       GS = DrawTrSurf::GetSurface(argv[1]);
 218       if (GS.IsNull()) { di << "OCC486 FAULTY. Null surface /n";return 1;}
 219       gp_Pnt P3D( Draw::Atof(argv[2]),Draw::Atof(argv[3]),Draw::Atof(argv[4]) );
 220
 221       Standard_Real Tol = Precision::PConfusion();
 222       Extrema_ExtPS myExtPS;
 223       if (argc > 5) du = Draw::Atof(argv[5]);
 224       if (argc > 6) dv = Draw::Atof(argv[6]);
 225
 226       Standard_Real uf, ul, vf, vl;
 227       GS->Bounds(uf, ul, vf, vl);
 228
 229       GeomAdaptor_Surface aSurf(GS);
 230       myExtPS.Initialize (aSurf, uf-du, ul+du, vf-dv, vl+dv, Tol, Tol );
 231       myExtPS.Perform ( P3D );
 232       Standard_Integer nPSurf = ( myExtPS.IsDone() ? myExtPS.NbExt() : 0 );
 233
 234       if ( nPSurf > 0 )
 235       {
 236         //Standard_Real distMin = myExtPS.Value ( 1 );
 237         Standard_Real distMin = myExtPS.SquareDistance ( 1 );
 238         Standard_Integer indMin=1;
 239         for (Standard_Integer sol = 2; sol <= nPSurf ; sol++)
 240         {
 241           //Standard_Real dist = myExtPS.Value(sol);
 242           Standard_Real dist = myExtPS.SquareDistance(sol);
 243           if ( distMin > dist )
 244           {
 245             distMin = dist;
 246             indMin = sol;
 247           }
 248         }
 249         distMin = sqrt(distMin);
 250         Standard_Real S, T;
 251         myExtPS.Point(indMin).Parameter ( S, T );
 252         gp_Pnt aCheckPnt = aSurf.Value( S, T );
 253         Standard_Real aCheckDist = P3D.Distance(aCheckPnt);
 254         di << "Solution is : U = "<< S << "\t V = "<< T << "\n";
 255         di << "Solution is : X = "<< aCheckPnt.X() << "\t Y = "<< aCheckPnt.Y() << "\t Z = "<< aCheckPnt.Z() << "\n";
 256         di << "ExtremaDistance = " << distMin  << "\n";
 257         di << "CheckDistance = " << aCheckDist << "\n";
 258
 259         if(fabs(distMin - aCheckDist) < Precision::Confusion()) return 0;
 260         else return 1;
 261       }
 262       else return 1;
 263     }
 264   catch (Standard_Failure const&) {di << "OCC486 Exception \n" ;return 1;}
 265 }
 266
 267 #include <GC_MakeArcOfCircle.hxx>
 268 #include <BRepBuilderAPI_MakeEdge.hxx>
 269 #include <BRepBuilderAPI_MakeWire.hxx>
 270 #include <TopoDS_Wire.hxx>
 271 #include <BRepBuilderAPI_MakeFace.hxx>
 272 #include <Geom_Plane.hxx>
 273 #include <gp_Pln.hxx>
 274 #include <BRepPrimAPI_MakePrism.hxx>
 275 #include <BRepOffsetAPI_DraftAngle.hxx>
 276 //=======================================================================
 277 //function : OCC712
 278 //purpose  :
 279 //=======================================================================
 280 static Standard_Integer OCC712 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv) {
 281   if (argc != 3) {
 282     di << "Usage : " << argv[0] << " draftAngle slabThick\n";
 283     return 1;
 284   }
 285   //NOTE: Case:1 - draftAngle = 15, slabThick = 30 --> Fails
 286   //      Case:2   draftAngle = 10, slabThick = 30 --> Ok
 287   //      Case:3   draftAngle = 10, slabThick = 40 --> Ok
 288   //
 289   //      --------------------------------------------------
 290   Standard_Real draftAngle = Draw::Atof(argv[1]);
 291   Standard_Real slabThick = Draw::Atof(argv[2]);
 292
 293   Standard_Real f1 = 75;
 294   Standard_Real f2 = 35;
 295
 296   gp_Pnt p1(-f2,  f2, 0);
 297   gp_Pnt p2(  0,  f1, 0);
 298   gp_Pnt p3( f2,  f2, 0);
 299   gp_Pnt p4( f1,   0, 0);
 300   gp_Pnt p5( f2, -f2, 0);
 301   gp_Pnt p6(  0, -f1, 0);
 302   gp_Pnt p7(-f2, -f2, 0);
 303   gp_Pnt p8(-f1,   0, 0);
 304
 305   GC_MakeArcOfCircle arc1(p1, p2, p3);
 306   GC_MakeArcOfCircle arc2(p3, p4, p5);
 307   GC_MakeArcOfCircle arc3(p5, p6, p7);
 308   GC_MakeArcOfCircle arc4(p7, p8, p1);
 309
 310   TopoDS_Edge e1 = BRepBuilderAPI_MakeEdge(arc1.Value());
 311   TopoDS_Edge e2 = BRepBuilderAPI_MakeEdge(arc2.Value());
 312   TopoDS_Edge e3 = BRepBuilderAPI_MakeEdge(arc3.Value());
 313   TopoDS_Edge e4 = BRepBuilderAPI_MakeEdge(arc4.Value());
 314
 315   BRepBuilderAPI_MakeWire MW;
 316   MW.Add(e1);
 317   MW.Add(e2);
 318   MW.Add(e3);
 319   MW.Add(e4);
 320
 321   if (!MW.IsDone())
 322     {
 323       di << "my Wire not done\n";
 324       return 1;
 325     }
 326   TopoDS_Wire W = MW.Wire();
 327
 328   TopoDS_Face F = BRepBuilderAPI_MakeFace(W);
 329   if ( F.IsNull())
 330     {
 331       di << " Error in Face creation \n";
 332       return 1;
 333     }
 334
 335   Handle(Geom_Surface) surf = BRep_Tool::Surface(F);
 336   Handle (Geom_Plane) P = Handle(Geom_Plane)::DownCast(surf);
 337   gp_Pln slabPln = P->Pln();
 338
 339   try
 340     {
 341     OCC_CATCH_SIGNALS
 342       gp_Dir slabDir(0, 0, 1);
 343       gp_Vec slabVect(slabDir);
 344       slabVect *= slabThick;
 345
 346       BRepPrimAPI_MakePrism slab(F, slabVect, Standard_True);
 347       if ( ! slab.IsDone() )
 348         {
 349           di << " Error in Slab creation \n";
 350           return 1;
 351         }
 352
 353       TopoDS_Shape slabShape = slab.Shape();
 354       if (fabs(draftAngle) > 0.01)
 355         {
 356           Standard_Real angle = draftAngle*(M_PI / 180.0);
 357           BRepOffsetAPI_DraftAngle draftSlab(slabShape);
 358
 359           TopoDS_Shape fShape = slab.FirstShape();
 360           TopoDS_Shape lShape = slab.LastShape();
 361
 362           TopExp_Explorer ex;
 363           for(ex.Init(slabShape, TopAbs_FACE); ex.More(); ex.Next())
 364             {
 365               TopoDS_Face aFace = TopoDS::Face(ex.Current());
 366               if(aFace.IsSame(fShape) || aFace.IsSame(lShape)) continue;
 367               draftSlab.Add(aFace, slabDir, angle, slabPln);
 368               if (!draftSlab.AddDone())
 369                 {
 370                   di << " Error in Add \n";
 371                   return 1;
 372                 }
 373             }
 374
 375           di << "All Faces added. Building... \n"; //std::cout.flush();
 376           draftSlab.Build();
 377           di << "Build done...\n"; //std::cout.flush();
 378           if (!draftSlab.IsDone())  //--------------> STEP:1
 379             {
 380               di << " Error in Build \n";
 381               return 1;
 382             }
 383           slabShape = draftSlab.Shape();
 384           DBRep::Set(argv[1], slabShape);
 385         }
 386     }
 387   catch ( Standard_Failure const& ) //--------------------> STEP:2
 388     {
 389       di << " Error in Draft Slab \n";
 390       return 1;
 391     }
 392   return 0;
 393 }
 394
 395 //=======================================================================
 396 //  performTriangulation
 397 //=======================================================================
 398
 399 Standard_Integer performTriangulation (TopoDS_Shape aShape, Draw_Interpretor& di)
 400 {
 401   int failed=0, total=0;
 402   TopExp_Explorer ExpFace;
 403   Handle(Poly_Triangulation) facing;
 404
 405   for (ExpFace.Init(aShape,TopAbs_FACE); ExpFace.More(); ExpFace.Next())
 406     {
 407       total++;
 408       TopoDS_Face TopologicalFace = TopoDS::Face (ExpFace.Current());
 409       TopologicalFace.Orientation (TopAbs_FORWARD) ;
 410       BRepMesh_IncrementalMesh IM(TopologicalFace, 1);
 411       TopLoc_Location loc;
 412       facing = BRep_Tool::Triangulation(TopologicalFace, loc);
 413       di << "Face " << total << " - ";
 414       if (facing.IsNull())
 415         {
 416           failed++;
 417           di << "******************** FAILED during Triangulation \n";
 418         }
 419       else
 420         {
 421           di << facing->NbTriangles() << " Triangles\n";
 422         }
 423     }
 424   di<<"Triangulation of all Faces Completed: \n\n";
 425   if (failed == 0) return 1;
 426   di<<"***************************************************\n";
 427   di<<"*******                                    ********\n";
 428   di<<"***** Triangulation FAILED for " << failed << " of " << total << " Faces ******\n";
 429   di<<"*******                                    ********\n";
 430   di<<"***************************************************\n";
 431   return 0;
 432 }
 433
 434 #include <BRepPrimAPI_MakeCylinder.hxx>
 435 #include <BRepPrimAPI_MakeCone.hxx>
 436 #include <BRepAlgoAPI_Cut.hxx>
 437 //=======================================================================
 438 //function : OCC822_1
 439 //purpose  :
 440 //=======================================================================
 441 static Standard_Integer OCC822_1 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv) {
 442
 443   if(argc != 4) {
 444     di << "Usage : " << argv[0] << " name1 name2 result\n";
 445     return 1;
 446   }
 447
 448   int index = 1;
 449
 450   gp_Pnt P1(0, 0, 0);
 451   gp_Dir D1(0, 0, 1);
 452   gp_Ax2 A1(P1,D1);
 453
 454   BRepPrimAPI_MakeCylinder cylMakerIn(A1, 40, 110);
 455   BRepPrimAPI_MakeCylinder cylMakerOut(A1, 50, 100);
 456   TopoDS_Shape cylIn = cylMakerIn.Shape();
 457   TopoDS_Shape cylOut = cylMakerOut.Shape();
 458
 459   gp_Pnt P2(0, 0, 0);
 460   gp_Dir D2(0, 0, -1);
 461   gp_Ax2 A2(P2,D2);
 462
 463   BRepPrimAPI_MakeCone conMakerIn(A2, 40, 60, 110);
 464   BRepPrimAPI_MakeCone conMakerOut(A2, 50, 70, 100);
 465   TopoDS_Shape conIn = conMakerIn.Shape();
 466   TopoDS_Shape conOut = conMakerOut.Shape();
 467
 468   di << "All primitives created.....  Creating Boolean\n";
 469
 470   try
 471   {
 472     OCC_CATCH_SIGNALS
 473
 474     di << "theIn = BRepAlgoAPI_Fuse(cylIn, conIn)\n";
 475     di << "theOut = BRepAlgoAPI_Fuse(cylOut, conOut)\n";
 476     di << "theRes = BRepAlgoAPI_Cut(theOut, theIn)\n";
 477     TopoDS_Shape theIn = BRepAlgoAPI_Fuse(cylIn, conIn).Shape();
 478     TopoDS_Shape theOut = BRepAlgoAPI_Fuse(cylOut, conOut).Shape();
 479     TopoDS_Shape theRes = BRepAlgoAPI_Cut(theOut, theIn).Shape();
 480
 481     if (index < argc) DBRep::Set(argv[index++], theIn);
 482     if (index < argc) DBRep::Set(argv[index++], theOut);
 483     if (index < argc) DBRep::Set(argv[index++], theRes);
 484     di << "Booleans Created !    Triangulating !\n";
 485
 486     performTriangulation(theRes, di);
 487   }
 488   catch ( Standard_Failure const& )
 489   {
 490     di << "*********************************************************\n";
 491     di << "*****                                              ******\n";
 492     di << "***** Standard_Failure : Exception in Shoe Function *****\n";
 493     di << "*****                                              ******\n";
 494     di << "*********************************************************\n";
 495     return 1;
 496   }
 497   return 0;
 498
 499 }
 500
 501 #include <BRepPrimAPI_MakeBox.hxx>
 502 #include <BRepPrimAPI_MakeSphere.hxx>
 503
 504 //=======================================================================
 505 //  OCC822_2
 506 //=======================================================================
 507
 508 static Standard_Integer OCC822_2 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 509 {
 510   if(argc != 4) {
 511     di << "Usage : " << argv[0] << " name1 name2 result\n";
 512     return 1;
 513   }
 514
 515   int index = 1;
 516
 517   gp_Dir xDir(1, 0, 0);
 518   gp_Dir zDir(0, 0, 1);
 519   gp_Pnt cen1(0, 0, 0);
 520   gp_Ax2 cor1(cen1, zDir, xDir);
 521   BRepPrimAPI_MakeBox boxMaker(cor1, 100, 100, 100);
 522   TopoDS_Shape box = boxMaker.Shape();
 523   if (index < argc) DBRep::Set(argv[index++], box);
 524
 525   BRepPrimAPI_MakeSphere sphereMaker(gp_Pnt(100.0, 50.0, 50.0), 25.0);
 526   TopoDS_Shape sph = sphereMaker.Shape();
 527   if (index < argc) DBRep::Set(argv[index++], sph);
 528
 529   di << "All primitives created.....  Creating Cut Objects\n";
 530
 531   try
 532   {
 533     OCC_CATCH_SIGNALS
 534
 535     di << "fuse = BRepAlgoAPI_Fuse(box, sph)\n";
 536     TopoDS_Shape fuse = BRepAlgoAPI_Fuse(box, sph).Shape();
 537
 538     if (index < argc) DBRep::Set(argv[index++], fuse);
 539     di << "Object Created !   Now Triangulating !";
 540
 541     performTriangulation(fuse, di);
 542   }
 543   catch ( Standard_Failure const& )
 544   {
 545     di << "*********************************************************\n";
 546     di << "*****                                              ******\n";
 547     di << "***** Standard_Failure : Exception in HSP Function ******\n";
 548     di << "*****                                              ******\n";
 549     di << "*********************************************************\n";
 550     return 1;
 551   }
 552
 553   return 0;
 554 }
 555
 556 //=======================================================================
 557 //  OCC823
 558 //=======================================================================
 559
 560 static Standard_Integer OCC823 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 561 {
 562   if(argc != 4) {
 563     di << "Usage : " << argv[0] << " name1 name2 result\n";
 564     return 1;
 565   }
 566
 567   int index = 1;
 568   Standard_Real size = 0.001;
 569
 570   gp_Pnt P1(40, 50, 0);
 571   gp_Dir D1(100, 0, 0);
 572   gp_Ax2 A1(P1,D1);
 573   BRepPrimAPI_MakeCylinder mkCyl1(A1, 20, 100);
 574   TopoDS_Shape cyl1 = mkCyl1.Shape();
 575   if (index < argc) DBRep::Set(argv[index++], cyl1);
 576
 577   gp_Pnt P2(100, 50, size);
 578   gp_Dir D2(0, size, 80);
 579   gp_Ax2 A2(P2,D2);
 580   BRepPrimAPI_MakeCylinder mkCyl2(A2, 20, 80);
 581   TopoDS_Shape cyl2 = mkCyl2.Shape();
 582   if (index < argc) DBRep::Set(argv[index++], cyl2);
 583
 584   di << "All primitives created.....  Creating Boolean\n";
 585
 586   try
 587   {
 588     OCC_CATCH_SIGNALS
 589
 590     di << "fuse = BRepAlgoAPI_Fuse(cyl2, cyl1)\n";
 591     TopoDS_Shape fuse = BRepAlgoAPI_Fuse(cyl2, cyl1).Shape();
 592
 593     if (index < argc) DBRep::Set(argv[index++], fuse);
 594     di << "Fuse Created !    Triangulating !\n";
 595
 596     performTriangulation(fuse, di);
 597   }
 598   catch (Standard_Failure const&)
 599   {
 600     di << "*********************************************************\n";
 601     di << "*****                                              ******\n";
 602     di << "***** Standard_Failure : Exception in TEE Function ******\n";
 603     di << "*****                                              ******\n";
 604     di << "*********************************************************\n";
 605     return 1;
 606   }
 607   return 0;
 608 }
 609
 610 //=======================================================================
 611 //  OCC824
 612 //=======================================================================
 613
 614 static Standard_Integer OCC824 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 615 {
 616   if(argc != 4) {
 617     di << "Usage : " << argv[0] << " name1 name2 result\n";
 618     return 1;
 619   }
 620
 621   int index = 1;
 622
 623   gp_Pnt P1(100, 0, 0);
 624   gp_Dir D1(-1, 0, 0);
 625   gp_Ax2 A1(P1,D1);
 626   BRepPrimAPI_MakeCylinder mkCyl(A1, 20, 100);
 627   TopoDS_Shape cyl = mkCyl.Shape();
 628   if (index < argc) DBRep::Set(argv[index++], cyl);
 629
 630   BRepPrimAPI_MakeSphere sphere(P1, 20.0);
 631   TopoDS_Shape sph = sphere.Shape();
 632   if (index < argc) DBRep::Set(argv[index++], sph);
 633
 634   di << "All primitives created.....  Creating Boolean\n";
 635
 636   try
 637   {
 638     OCC_CATCH_SIGNALS
 639
 640     di << "fuse = BRepAlgoAPI_Fuse(cyl, sph)\n";
 641     TopoDS_Shape fuse = BRepAlgoAPI_Fuse(cyl, sph).Shape();
 642
 643     di << "Fuse Created !    Triangulating !\n";
 644     if (index < argc) DBRep::Set(argv[index++], fuse);
 645
 646     performTriangulation(fuse, di);
 647   }
 648   catch (Standard_Failure const&)
 649   {
 650     di << "*********************************************************\n";
 651     di << "*****                                              ******\n";
 652     di << "***** Standard_Failure : Exception in YOU Function ******\n";
 653     di << "*****                                              ******\n";
 654     di << "*********************************************************\n";
 655     return 1;
 656   }
 657   return 0;
 658 }
 659
 660 #include <TColgp_Array2OfPnt.hxx>
 661 #include <GeomConvert.hxx>
 662 #include <Geom_BezierSurface.hxx>
 663 #include <BRepPrimAPI_MakeHalfSpace.hxx>
 664
 665 //=======================================================================
 666 //  OCC825
 667 //=======================================================================
 668
 669 static Standard_Integer OCC825 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 670 {
 671   if(argc != 6) {
 672     di << "Usage : " << argv[0] << " name1 name2 name3 result1 result2\n";
 673     return 1;
 674   }
 675
 676   int index = 1;
 677
 678   Standard_Real size = 50;
 679   TColgp_Array2OfPnt poles(1, 2, 1, 2);
 680
 681   poles(1, 1).SetCoord(-size, 0, -size);
 682   poles(1, 2).SetCoord(-size, 0,  size);
 683   poles(2, 1).SetCoord( size, 0, -size);
 684   poles(2, 2).SetCoord( size, 0,  size);
 685
 686   Handle(Geom_BezierSurface) BezSurf = new Geom_BezierSurface(poles);
 687   Handle(Geom_BSplineSurface) BSpSurf = GeomConvert::SurfaceToBSplineSurface(BezSurf);
 688   BRepBuilderAPI_MakeFace faceMaker(BSpSurf, Precision::Confusion());
 689   TopoDS_Face face = faceMaker.Face();
 690
 691   gp_Pnt pnt(0, size, 0);
 692   BRepPrimAPI_MakeHalfSpace *hSpace = new BRepPrimAPI_MakeHalfSpace(face,pnt);
 693   TopoDS_Shape hsp = hSpace->Solid();
 694   if (index < argc) DBRep::Set(argv[index++], hsp);
 695
 696   BRepPrimAPI_MakeSphere sphere1(gp_Pnt(0.0, 0.0, 0.0), 25.0);
 697   TopoDS_Shape sph1 = sphere1.Shape();
 698   if (index < argc) DBRep::Set(argv[index++], sph1);
 699
 700   BRepPrimAPI_MakeSphere sphere2(gp_Pnt(0.0, 0.00001, 0.0), 25.0);
 701   TopoDS_Shape sph2 = sphere2.Shape();
 702   if (index < argc) DBRep::Set(argv[index++], sph2);
 703
 704   di << "All primitives created.....  Creating Cut Objects\n";
 705
 706   try
 707   {
 708     OCC_CATCH_SIGNALS
 709
 710     di << "cut1 = BRepAlgoAPI_Cut(sph1, hsp)\n";
 711     TopoDS_Shape cut1 = BRepAlgoAPI_Cut(sph1, hsp).Shape();
 712
 713     if (index < argc) DBRep::Set(argv[index++], cut1);
 714     di << "CUT 1 Created !   " ;
 715
 716
 717     di << "cut2 = BRepAlgoAPI_Cut(sph2, hsp)\n";
 718     TopoDS_Shape cut2 = BRepAlgoAPI_Cut(sph2, hsp).Shape();
 719
 720     if (index < argc) DBRep::Set(argv[index++], cut2);
 721     di << "CUT 2 Created !\n\n";
 722
 723     GProp_GProps G;
 724     BRepGProp::VolumeProperties(cut1, G);
 725     di << "CUT 1 Mass = " << G.Mass() << "\n\n";
 726     BRepGProp::VolumeProperties(cut2, G);
 727     di << "CUT 2 Mass = " << G.Mass() << "\n\n";
 728
 729     di << "Trianglating Faces of CUT 1 .....\n";
 730     performTriangulation(cut1, di);
 731
 732     di << "Trianglating Faces of CUT 2 .....\n";
 733     performTriangulation(cut2, di);
 734   }
 735   catch (Standard_Failure const&)
 736   {
 737     di << "*********************************************************\n";
 738     di << "*****                                              ******\n";
 739     di << "***** Standard_Failure : Exception in HSP Function ******\n";
 740     di << "*****                                              ******\n";
 741     di << "*********************************************************\n";
 742     return 1;
 743   }
 744
 745   di << "*************************************************************\n";
 746   di << " CUT 1 and CUT 2 gives entirely different results during\n";
 747   di << " mass computation and face triangulation, eventhough the\n";
 748   di << " two spheres are located more or less at the same position.\n";
 749   di << "*************************************************************\n";
 750
 751   return 0;
 752 }
 753
 754 //=======================================================================
 755 //  OCC826
 756 //=======================================================================
 757
 758 static Standard_Integer OCC826 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 759 {
 760   if(argc != 4) {
 761     di << "Usage : " << argv[0] << " name1 name2 result\n";
 762     return 1;
 763   }
 764
 765   int index = 1;
 766
 767   Standard_Real x1 = 181.82808;
 768   Standard_Real x2 = 202.39390;
 769   Standard_Real y1 = 31.011970;
 770   Standard_Real y2 = 123.06856;
 771
 772   BRepBuilderAPI_MakePolygon W1;
 773   W1.Add(gp_Pnt(x1, y1, 0));
 774   W1.Add(gp_Pnt(x2, y1, 0));
 775   W1.Add(gp_Pnt(x2, y2, 0));
 776   W1.Add(gp_Pnt(x1, y2, 0));
 777   W1.Add(gp_Pnt(x1, y1, 0));
 778
 779   Standard_Boolean myFalse = Standard_False;
 780   TopoDS_Face F1 = BRepBuilderAPI_MakeFace(W1.Wire(), myFalse);
 781
 782   gp_Pnt P1(0, 0, 0);
 783   gp_Dir D1(0, 30, 0);
 784   gp_Ax1 A1(P1,D1);
 785   Standard_Real angle1 = 360 * (M_PI / 180.0);
 786   TopoDS_Shape rev = BRepPrimAPI_MakeRevol(F1, A1, angle1);
 787   if (index < argc) DBRep::Set(argv[index++], rev);
 788
 789   BRepPrimAPI_MakeSphere sphere(gp_Pnt(166.373, 77.0402, 96.0555), 23.218586);
 790   TopoDS_Shape sph = sphere.Shape();
 791   if (index < argc) DBRep::Set(argv[index++], sph);
 792
 793   di << "All primitives created.....  Creating Boolean\n";
 794
 795   try
 796   {
 797     OCC_CATCH_SIGNALS
 798
 799     di << "fuse = BRepAlgoAPI_Fuse(rev, sph)\n";
 800     TopoDS_Shape fuse = BRepAlgoAPI_Fuse(rev, sph).Shape();
 801
 802     if (index < argc) DBRep::Set(argv[index++], fuse);
 803     di << "Fuse Created !   Triangulating !\n";
 804     performTriangulation(fuse, di);
 805   }
 806   catch (Standard_Failure const&)
 807   {
 808     di << "*********************************************************\n";
 809     di << "*****                                              ******\n";
 810     di << "***** Standard_Failure : Exception in SPH Function ******\n";
 811     di << "*****                                              ******\n";
 812     di << "*********************************************************\n";
 813     return 1;
 814   }
 815   return 0;
 816 }
 817
 818 #include <BRepPrimAPI_MakeTorus.hxx>
 819 //=======================================================================
 820 //  OCC827
 821 //=======================================================================
 822
 823 static Standard_Integer OCC827 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 824 {
 825   if(argc != 6) {
 826     di << "Usage : " << argv[0] << " name1 name2 name3 result1 result2\n";
 827     return 1;
 828   }
 829
 830   int index = 1;
 831
 832   BRepBuilderAPI_MakePolygon W1;
 833   W1.Add(gp_Pnt(10, 0, 0));
 834   W1.Add(gp_Pnt(20, 0, 0));
 835   W1.Add(gp_Pnt(20, 0, 50));
 836   W1.Add(gp_Pnt(10, 0, 50));
 837   W1.Add(gp_Pnt(10, 0, 0));
 838
 839   Standard_Boolean myFalse = Standard_False;
 840   TopoDS_Face F1 = BRepBuilderAPI_MakeFace(W1.Wire(), myFalse);
 841
 842   gp_Pnt P1(0, 0, 0);
 843   gp_Dir D1(0, 0, 30);
 844   gp_Ax1 A1(P1,D1);
 845   Standard_Real angle1 = 360 * (M_PI / 180.0);
 846   TopoDS_Shape rev = BRepPrimAPI_MakeRevol(F1, A1, angle1);
 847   if (index < argc) DBRep::Set(argv[index++], rev);
 848
 849   gp_Pnt P2(0, 0, 50);
 850   gp_Dir D2(0, 0, 30);
 851   gp_Ax2 A2(P2,D2);
 852   Standard_Real majRad = 15;
 853   Standard_Real minRad = 5;
 854   BRepPrimAPI_MakeTorus Torus1(A2, majRad, minRad);
 855   TopoDS_Shape tor1 = Torus1.Shape();
 856   if (index < argc) DBRep::Set(argv[index++], tor1);
 857
 858   gp_Pnt P3(0, 0, 10);
 859   gp_Dir D3(0, 0, 30);
 860   gp_Ax2 A3(P3,D3);
 861   BRepPrimAPI_MakeTorus Torus2(A3, majRad, minRad);
 862   TopoDS_Shape tor2 = Torus2.Shape();
 863   if (index < argc) DBRep::Set(argv[index++], tor2);
 864
 865   di << "All primitives created.....  Creating Boolean\n";
 866
 867   try
 868   {
 869     OCC_CATCH_SIGNALS
 870
 871     di << "Fuse1 = BRepAlgoAPI_Fuse(tor1, rev)\n";
 872     TopoDS_Shape fuse1 = BRepAlgoAPI_Fuse(tor1, rev).Shape();
 873
 874     if (index < argc) DBRep::Set(argv[index++], fuse1);
 875     di << "Fuse1 Created !    Creating Fuse 2\n";
 876
 877     di << "Fuse2 = BRepAlgoAPI_Fuse(tor2, fuse1)\n";
 878     TopoDS_Shape fuse2 = BRepAlgoAPI_Fuse(tor2, fuse1).Shape();
 879
 880     if (index < argc) DBRep::Set(argv[index++], fuse2);
 881     di << "Fuse2 Created !    Triangulating !\n";
 882
 883     performTriangulation(fuse2, di);
 884   }
 885   catch (Standard_Failure const&)
 886   {
 887     di << "*********************************************************\n";
 888     di << "*****                                              ******\n";
 889     di << "***** Standard_Failure : Exception in REV Function ******\n";
 890     di << "*****                                              ******\n";
 891     di << "*********************************************************\n";
 892     return 1;
 893   }
 894   return 0;
 895 }
 896
 897 //=======================================================================
 898 //  performBlend
 899 //=======================================================================
 900
 901 int performBlend (TopoDS_Shape aShape, Standard_Real rad, TopoDS_Shape& bShape, Draw_Interpretor& di)
 902 {
 903   Standard_Integer status = 0;
 904   TopoDS_Shape newShape;
 905   TopTools_IndexedDataMapOfShapeListOfShape edgemap;
 906   TopExp::MapShapesAndAncestors(aShape,TopAbs_EDGE,TopAbs_SOLID,edgemap);
 907   di << "Blending All Edges: No. of Edges: " << edgemap.Extent() << "\n";
 908   ChFi3d_FilletShape FShape = ChFi3d_Rational;
 909   BRepFilletAPI_MakeFillet blend(aShape,FShape);
 910   for(int i = 1; i <= edgemap.Extent(); i++)
 911     {
 912       TopoDS_Edge edg = TopoDS::Edge( edgemap.FindKey(i) );
 913       if(!edg.IsNull()) blend.Add(rad, edg);
 914     }
 915
 916   try
 917     {
 918     OCC_CATCH_SIGNALS
 919       blend.Build();
 920       if(!blend.HasResult() || blend.Shape().IsNull()) {
 921         status = 1;
 922       }
 923     }
 924   catch ( Standard_Failure const& )
 925     {
 926       status = 1;
 927     }
 928   if(status) {
 929       di<<"*******************************************************\n";
 930       di<<"******                                          *******\n";
 931       di<<"****** Blending Failed (Radius = " << rad << ") *******\n";
 932       di<<"******                                          *******\n";
 933       di<<"*******************************************************\n";
 934       return 1;
 935   } else {
 936     di<<"Blending successfully performed on all Edges: \n\n";
 937   }
 938   bShape = blend.Shape();
 939   return 0;
 940 }
 941
 942 #include <GC_MakeSegment.hxx>
 943 //=======================================================================
 944 //  OCC828
 945 //=======================================================================
 946
 947 static Standard_Integer OCC828 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv)
 948 {
 949   if(argc != 2) {
 950     di << "Usage : " << argv[0] << " shape\n";
 951   }
 952   int index = 1;
 953
 954   Standard_Real slabThick = 111;
 955
 956   gp_Pnt p11(-27.598139, -7.0408573, 0.0);
 957   gp_Pnt p12(-28.483755, -17.487625, 0.0);
 958   gp_Pnt p13(-19.555504, -22.983587, 0.0);
 959   GC_MakeArcOfCircle arc1(p11, p12, p13);
 960
 961   gp_Pnt p21(12.125083, -22.983587, 0.0);
 962   gp_Pnt p22(21.1572, -17.27554, 0.0);
 963   gp_Pnt p23(19.878168, -6.6677585, 0.0);
 964   GC_MakeArcOfCircle arc2(p21, p22, p23);
 965
 966   gp_Pnt p31(3.265825, 13.724955, 0.0);
 967   gp_Pnt p32(-4.7233953, 17.406338, 0.0);
 968   gp_Pnt p33(-12.529893, 13.351856, 0.0);
 969   GC_MakeArcOfCircle arc3(p31, p32, p33);
 970
 971   GC_MakeSegment ln1(p13, p21);
 972   GC_MakeSegment ln2(p23, p31);
 973   GC_MakeSegment ln3(p33, p11);
 974
 975   TopoDS_Edge e1 = BRepBuilderAPI_MakeEdge(arc1.Value());
 976   TopoDS_Edge e2 = BRepBuilderAPI_MakeEdge(arc2.Value());
 977   TopoDS_Edge e3 = BRepBuilderAPI_MakeEdge(arc3.Value());
 978
 979   TopoDS_Edge e4 = BRepBuilderAPI_MakeEdge(ln1.Value());
 980   TopoDS_Edge e5 = BRepBuilderAPI_MakeEdge(ln2.Value());
 981   TopoDS_Edge e6 = BRepBuilderAPI_MakeEdge(ln3.Value());
 982
 983   BRepBuilderAPI_MakeWire MW;
 984   MW.Add(e1);
 985   MW.Add(e4);
 986   MW.Add(e2);
 987   MW.Add(e5);
 988   MW.Add(e3);
 989   MW.Add(e6);
 990
 991   if (!MW.IsDone())
 992     {
 993       di << "my Wire not done\n";
 994       return 1;
 995     }
 996
 997   TopoDS_Wire W = MW.Wire();
 998   TopoDS_Face F = BRepBuilderAPI_MakeFace(W);
 999   if ( F.IsNull())
1000     {
1001       di << " Error in Face creation \n";
1002       return 1;
1003     }
1004
1005   try
1006     {
1007     OCC_CATCH_SIGNALS
1008       gp_Dir slabDir(0, 0, 1);
1009       gp_Vec slabVect(slabDir);
1010       slabVect *= slabThick;
1011
1012       BRepPrimAPI_MakePrism slab(F, slabVect, Standard_True);
1013       if ( ! slab.IsDone() )
1014         {
1015           di << " Error in Slab creation \n";
1016           return 1;
1017         }
1018       if (index < argc) DBRep::Set(argv[index++], slab.Shape());
1019
1020 //       std::cout << "Slab Successfully Created !   Now Blending ..." << std::endl;
1021 //       TopoDS_Shape aShape;
1022 //       int ret = performBlend(slab.Shape(), radius, aShape);
1023 //       if (ret) return 1;
1024 //       if (index < argc) DBRep::Set(argv[index++], aShape);
1025
1026 //       std::cout << "Blending Successfully Done !   Now Triangulating ..." << std::endl;
1027 //       performTriangulation(aShape);
1028     }
1029   catch ( Standard_Failure const& )
1030     {
1031       di << " Error in Draft Slab \n";
1032       return 1;
1033     }
1034   return 0;
1035 }
1036
1037 void QABugs::Commands_10(Draw_Interpretor& theCommands) {
1038   const char *group = "QABugs";
1039
1040   theCommands.Add ("OCC426", "OCC426 shape1 shape2 shape3 shape4 shape5 shape6 shape7", __FILE__, OCC426, group);
1041
1042   theCommands.Add("isperiodic", "Use : isperiodic surfaceOfRevolution", __FILE__, isPeriodic, group);
1043   theCommands.Add("OCC486", "Use : OCC486 surf x y z du dv ", __FILE__, OCC486, group);
1044   theCommands.Add("OCC712", "OCC712 draftAngle slabThick", __FILE__, OCC712, group);
1045   theCommands.Add("OCC822_1", "OCC822_1 name1 name2 result", __FILE__,OCC822_1, group);
1046   theCommands.Add("OCC822_2", "OCC822_2 name1 name2 result", __FILE__,OCC822_2, group);
1047   theCommands.Add("OCC823", "OCC823 name1 name2 result", __FILE__,OCC823, group);
1048   theCommands.Add("OCC824", "OCC824 name1 name2 result", __FILE__,OCC824, group);
1049   theCommands.Add("OCC825", "OCC825 name1 name2 name3 name4 name5", __FILE__,OCC825, group);
1050   theCommands.Add("OCC826", "OCC826 name1 name2 result", __FILE__,OCC826, group);
1051   theCommands.Add("OCC827", "OCC827 name1 name2 name3 result1 result2", __FILE__,OCC827, group);
1052   theCommands.Add("OCC828", "OCC828 redius shape result ", __FILE__,OCC828, group);
1053
1054   return;
1055 }