src/QABugs/QABugs_19.cxx

   1 // Created on: 2002-05-21
   2 // Created by: QA Admin
   3 // Copyright (c) 2002-2012 OPEN CASCADE SAS
   4 //
   5 // The content of this file is subject to the Open CASCADE Technology Public
   6 // License Version 6.5 (the "License"). You may not use the content of this file
   7 // except in compliance with the License. Please obtain a copy of the License
   8 // at http://www.opencascade.org and read it completely before using this file.
   9 //
  10 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
  11 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
  12 //
  13 // The Original Code and all software distributed under the License is
  14 // distributed on an "AS IS" basis, without warranty of any kind, and the
  15 // Initial Developer hereby disclaims all such warranties, including without
  16 // limitation, any warranties of merchantability, fitness for a particular
  17 // purpose or non-infringement. Please see the License for the specific terms
  18 // and conditions governing the rights and limitations under the License.
  19
  20
  21 #include <QABugs.hxx>
  22
  23 #include <Draw_Interpretor.hxx>
  24 #include <DBRep.hxx>
  25 #include <DrawTrSurf.hxx>
  26 #include <AIS_InteractiveContext.hxx>
  27 #include <ViewerTest.hxx>
  28 #include <AIS_Shape.hxx>
  29 #include <TopoDS_Shape.hxx>
  30
  31 #include <gp_Pnt2d.hxx>
  32 #include <gp_Ax1.hxx>
  33 #include <GCE2d_MakeSegment.hxx>
  34 #include <Geom2d_TrimmedCurve.hxx>
  35 #include <DrawTrSurf.hxx>
  36
  37 #include <Precision.hxx>
  38
  39 #include <PCollection_HAsciiString.hxx>
  40
  41 #include <cstdio>
  42 #include <cmath>
  43 #include <iostream>
  44 #include <OSD_PerfMeter.hxx>
  45 #include <OSD_Timer.hxx>
  46 #include <BRepPrimAPI_MakeBox.hxx>
  47 #include <BRepPrimAPI_MakeSphere.hxx>
  48 #include <BRepAlgo_Cut.hxx>
  49
  50
  51 static Standard_Integer OCC230 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
  52 {
  53   if ( argc != 4) {
  54     di << "ERROR OCC230: Usage : " << argv[0] << " TrimmedCurve Pnt2d Pnt2d" << "\n";
  55     return 1;
  56   }
  57
  58   gp_Pnt2d P1, P2;
  59   if ( !DrawTrSurf::GetPoint2d(argv[2],P1)) {
  60     di << "ERROR OCC230: " << argv[2] << " is not Pnt2d" << "\n";
  61     return 1;
  62   }
  63   if ( !DrawTrSurf::GetPoint2d(argv[3],P2)) {
  64     di << "ERROR OCC230: " << argv[3] << " is not Pnt2d" << "\n";
  65     return 1;
  66   }
  67
  68   GCE2d_MakeSegment MakeSegment(P1,P2);
  69   Handle(Geom2d_TrimmedCurve) TrimmedCurve = MakeSegment.Value();
  70   DrawTrSurf::Set(argv[1], TrimmedCurve);
  71   return 0;
  72 }
  73
  74 static Standard_Integer OCC142 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
  75 {
  76   for(Standard_Integer i= 0;i <= 20;i++){
  77     Handle(PCollection_HAsciiString) pstr = new PCollection_HAsciiString("TEST");
  78     pstr->Clear();
  79   }
  80   di << "OCC142: OK" << "\n";
  81   return 0;
  82 }
  83
  84 static Standard_Integer OCC23361 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
  85 {
  86   gp_Pnt p(0, 0, 2);
  87
  88   gp_Trsf t1, t2;
  89   t1.SetRotation(gp_Ax1(p, gp_Dir(0, 1, 0)), -0.49328285294022267);
  90   t2.SetRotation(gp_Ax1(p, gp_Dir(0, 0, 1)), 0.87538474718473880);
  91
  92   gp_Trsf tComp = t2 * t1;
  93
  94   gp_Pnt p1(10, 3, 4);
  95   gp_Pnt p2 = p1.Transformed(tComp);
  96   gp_Pnt p3 = p1.Transformed(t1);
  97   p3.Transform(t2);
  98
  99   // points must be equal
 100   if ( ! p2.IsEqual(p3, Precision::Confusion()) )
 101     di << "ERROR OCC23361: equivalent transformations does not produce equal points" << "\n";
 102   else
 103     di << "OCC23361: OK" << "\n";
 104
 105   return 0;
 106 }
 107
 108 static Standard_Integer OCC23237 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
 109 {
 110   OSD_PerfMeter aPM("TestMeter",0);
 111   OSD_Timer aTM;
 112
 113   // run some operation in cycle for about 2 seconds to have good values of times to compare
 114   int count = 0;
 115   printf("OSD_PerfMeter test.\nRunning Boolean operation on solids in loop.\n");
 116   for (; aTM.ElapsedTime() < 2.; count++)
 117   {
 118     aPM.Start();
 119     aTM.Start();
 120
 121     // do some operation that will take considerable time compared with time or starting / stopping timers
 122     BRepPrimAPI_MakeBox aBox (10., 10., 10.);
 123     BRepPrimAPI_MakeSphere aSphere (10.);
 124     BRepAlgo_Cut (aBox.Shape(), aSphere.Shape());
 125
 126     aTM.Stop();
 127     aPM.Stop();
 128   }
 129
 130   int aNbEnters = 0;
 131   Standard_Real aPerfMeter_CPUtime = 0., aTimer_ElapsedTime = aTM.ElapsedTime();
 132
 133   perf_get_meter("TestMeter", &aNbEnters, &aPerfMeter_CPUtime);
 134
 135   Standard_Real aTimeDiff = (fabs(aTimer_ElapsedTime - aPerfMeter_CPUtime) / aTimer_ElapsedTime);
 136
 137   printf("\nMeasurement results (%d cycles):\n", count);
 138   printf("\nOSD_PerfMeter CPU time: %lf\nOSD_Timer elapsed time: %lf\n", aPerfMeter_CPUtime, aTimer_ElapsedTime);
 139   printf("Time delta is: %.3lf %%\n", aTimeDiff * 100);
 140
 141   if (aTimeDiff > 0.2)
 142     di << "OCC23237: Error: too much difference between CPU and elapsed times";
 143   else if (aNbEnters != count)
 144     di << "OCC23237: Error: counter reported by PerfMeter (" << aNbEnters << ") does not correspond to actual number of cycles";
 145   else
 146     di << "OCC23237: OK";
 147
 148   return 0;
 149 }
 150
 151 #ifdef HAVE_TBB
 152
 153 #include <Standard_Atomic.hxx>
 154 #include <tbb/blocked_range.h>
 155 #include <tbb/parallel_for.h>
 156
 157 class IncrementerDecrementer
 158 {
 159 public:
 160     IncrementerDecrementer (Standard_Integer* theVal, Standard_Boolean thePositive) : myVal (theVal), myPositive (thePositive)
 161     {}
 162     void operator() (const tbb::blocked_range<size_t>& r) const
 163     {
 164         if (myPositive)
 165             for (size_t i = r.begin(); i != r.end(); ++i)
 166                 Standard_Atomic_Increment (myVal);
 167         else
 168             for (size_t i = r.begin(); i != r.end(); ++i)
 169                 Standard_Atomic_Decrement (myVal);
 170     }
 171 private:
 172     Standard_Integer*   myVal;
 173     Standard_Boolean   myPositive;
 174 };
 175 #endif
 176
 177 #define QCOMPARE(val1, val2) \
 178   di << "Checking " #val1 " == " #val2 << \
 179         ((val1) == (val2) ? ": OK\n" : ": Error\n")
 180
 181 #ifdef HAVE_TBB
 182 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
 183 {
 184   int aSum = 0;
 185
 186   //check returned value
 187   QCOMPARE (Standard_Atomic_Decrement (&aSum), -1);
 188   QCOMPARE (Standard_Atomic_Increment (&aSum), 0);
 189   QCOMPARE (Standard_Atomic_Increment (&aSum), 1);
 190   QCOMPARE (Standard_Atomic_Increment (&aSum), 2);
 191 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 0);
 192 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 1);
 193
 194   //check atomicity
 195   aSum = 0;
 196   const int N = 1 << 24; //big enough to ensure concurrency
 197
 198   //increment
 199   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, true));
 200   QCOMPARE (aSum, N);
 201
 202   //decrement
 203   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
 204   QCOMPARE (aSum, 0);
 205
 206   return 0;
 207 }
 208
 209 #else /* HAVE_TBB */
 210
 211 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
 212 {
 213   di << "Test skipped: command " << argv[0] << " requires TBB library\n";
 214   return 0;
 215 }
 216
 217 #endif /* HAVE_TBB */
 218
 219 #include <TDocStd_Application.hxx>
 220 #include <XCAFApp_Application.hxx>
 221 #include <TDocStd_Document.hxx>
 222 #include <XCAFDoc_ShapeTool.hxx>
 223 #include <XCAFDoc_DocumentTool.hxx>
 224 #include <TDF_Label.hxx>
 225 #include <TDataStd_Name.hxx>
 226
 227 static Standard_Integer OCC23595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
 228 {
 229   const Handle(TDocStd_Application)& anApp = XCAFApp_Application::GetApplication();
 230   Handle(TDocStd_Document) aDoc;
 231   anApp->NewDocument ("XmlXCAF", aDoc);
 232   QCOMPARE (!aDoc.IsNull(), Standard_True);
 233
 234   Handle(XCAFDoc_ShapeTool) aShTool = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main());
 235
 236   //check default value
 237   Standard_Boolean aValue = XCAFDoc_ShapeTool::AutoNaming();
 238   QCOMPARE (aValue, Standard_True);
 239
 240   //true
 241   XCAFDoc_ShapeTool::SetAutoNaming (Standard_True);
 242   TopoDS_Shape aShape = BRepPrimAPI_MakeBox (100., 200., 300.).Shape();
 243   TDF_Label aLabel = aShTool->AddShape (aShape);
 244   Handle(TDataStd_Name) anAttr;
 245   QCOMPARE (aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
 246
 247   //false
 248   XCAFDoc_ShapeTool::SetAutoNaming (Standard_False);
 249   aShape = BRepPrimAPI_MakeBox (300., 200., 100.).Shape();
 250   aLabel = aShTool->AddShape (aShape);
 251   QCOMPARE (!aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
 252
 253   //restore
 254   XCAFDoc_ShapeTool::SetAutoNaming (aValue);
 255
 256   return 0;
 257 }
 258
 259 #include <ExprIntrp_GenExp.hxx>
 260 Standard_Integer OCC22611 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 261 {
 262
 263   if (argc != 3) {
 264     di << "Usage : " << argv[0] << " string nb\n";
 265     return 1;
 266   }
 267
 268   TCollection_AsciiString aToken = argv[1];
 269   Standard_Integer aNb = atoi(argv[2]);
 270
 271   Handle(ExprIntrp_GenExp) aGen = ExprIntrp_GenExp::Create();
 272   for (Standard_Integer i=0; i < aNb; i++)
 273   {
 274     aGen->Process(aToken);
 275     Handle(Expr_GeneralExpression) aExpr = aGen->Expression();
 276   }
 277
 278   return 0;
 279 }
 280
 281 Standard_Integer OCC22595 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 282 {
 283   gp_Mat M0;
 284   di << "M0 = "
 285   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}"
 286   << " {" << M0(2,1) << "} {" << M0(2,2) << "} {" << M0(2,3) <<"}"
 287   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}";
 288   return 0;
 289 }
 290
 291 #include <TopoDS_Face.hxx>
 292 #include <TopoDS_Face.hxx>
 293 #include <TopoDS.hxx>
 294 #include <BRepBuilderAPI_Transform.hxx>
 295 #include <BRepExtrema_DistShapeShape.hxx>
 296 #include <BRepTools.hxx>
 297
 298 Standard_Boolean static OCC23774Test(const TopoDS_Face& grossPlateFace, const TopoDS_Shape& originalWire, Draw_Interpretor& di)
 299 {
 300   BRepExtrema_DistShapeShape distShapeShape(grossPlateFace,originalWire,Extrema_ExtFlag_MIN);
 301   if(!distShapeShape.IsDone()) {
 302     di <<"Distance ShapeShape is Not Done\n";
 303     return Standard_False;
 304   }
 305
 306   if(distShapeShape.Value() > 0.01) {
 307     di << "Wrong Dist = " <<distShapeShape.Value() << "\n";
 308     return Standard_False;
 309   } else
 310     di << "Dist0 = " <<distShapeShape.Value() <<"\n";
 311
 312   //////////////////////////////////////////////////////////////////////////
 313   /// First Flip Y
 314   const gp_Pnt2d axis1P1(1474.8199035519228,1249.9995745636970);
 315   const gp_Pnt2d axis1P2(1474.8199035519228,1250.9995745636970);
 316
 317   gp_Vec2d mirrorVector1(axis1P1,axis1P2);
 318
 319   gp_Trsf2d mirror1;
 320   mirror1.SetMirror(gp_Ax2d(axis1P1,mirrorVector1));
 321
 322   BRepBuilderAPI_Transform transformer1(mirror1);
 323   transformer1.Perform(originalWire);
 324   if(!transformer1.IsDone()) {
 325     di << "Not Done1 " << "\n";
 326     return Standard_False;
 327   }
 328   const TopoDS_Shape& step1ModifiedShape = transformer1.ModifiedShape(originalWire);
 329
 330   BRepExtrema_DistShapeShape distShapeShape1(grossPlateFace,step1ModifiedShape,Extrema_ExtFlag_MIN);
 331   if(!distShapeShape1.IsDone())
 332     return Standard_False;
 333   if(distShapeShape1.Value() > 0.01) {
 334     di << "Dist = " <<distShapeShape1.Value() <<"\n";
 335     return Standard_False;
 336   } else
 337     di << "Dist1 = " <<distShapeShape1.Value() <<"\n";
 338
 339   //////////////////////////////////////////////////////////////////////////
 340   /// Second flip Y
 341   transformer1.Perform(step1ModifiedShape);
 342   if(!transformer1.IsDone()) {
 343     di << "Not Done1 \n";
 344     return Standard_False;
 345   }
 346   const TopoDS_Shape& step2ModifiedShape = transformer1.ModifiedShape(step1ModifiedShape);
 347
 348   //This is identity matrix for values but for type is gp_Rotation ?!
 349   gp_Trsf2d mirror11 = mirror1;
 350   mirror11.PreMultiply(mirror1);
 351
 352   BRepExtrema_DistShapeShape distShapeShape2(grossPlateFace,step2ModifiedShape);//,Extrema_ExtFlag_MIN);
 353   if(!distShapeShape2.IsDone())
 354     return Standard_False;
 355
 356   //This last test case give error (the value is 1008.8822038689706)
 357   if(distShapeShape2.Value() > 0.01) {
 358     di  << "Wrong Dist2 = " <<distShapeShape2.Value() <<"\n";
 359     Standard_Integer N = distShapeShape2.NbSolution();
 360     di << "Nb = " <<N <<"\n";
 361     for (Standard_Integer i=1;i <= N;i++)
 362         di <<"Sol(" <<i<<") = " <<distShapeShape2.PointOnShape1(i).Distance(distShapeShape2.PointOnShape2(i)) <<"\n";
 363     return Standard_False;
 364   }
 365   di << "Distance2 = " <<distShapeShape2.Value() <<"\n";
 366
 367   return Standard_True;
 368 }
 369 static Standard_Integer OCC23774(Draw_Interpretor& di, Standard_Integer n, const char** a)
 370 {
 371
 372   if (n != 3) {
 373         di <<"OCC23774: invalid number of input parameters\n";
 374         return 1;
 375   }
 376
 377   const char *ns1 = (a[1]), *ns2 = (a[2]);
 378   TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2));
 379   if (S1.IsNull() || S2.IsNull()) {
 380         di <<"OCC23774: Null input shapes\n";
 381         return 1;
 382   }
 383   const TopoDS_Face& aFace  = TopoDS::Face(S1);
 384   if(!OCC23774Test(aFace, S2, di))
 385         di << "Something is wrong\n";
 386
 387  return 0;
 388 }
 389
 390 #include <GeomConvert_ApproxSurface.hxx>
 391 #include <Geom_BSplineSurface.hxx>
 392 #include <Draw.hxx>
 393 #include <OSD_Thread.hxx>
 394 static void GeomConvertTest (Draw_Interpretor& di, Standard_Integer theTargetNbUPoles, Standard_CString theFileName)
 395 {
 396         Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(theFileName);
 397         GeomConvert_ApproxSurface aGAS (aSurf, 1e-4, GeomAbs_C1, GeomAbs_C1, 9, 9, 100, 1);
 398         if (!aGAS.IsDone()) {
 399                 di << "ApproxSurface is not done!" << "\n";
 400                 return;
 401         }
 402         const Handle(Geom_BSplineSurface)& aBSurf = aGAS.Surface();
 403         if (aBSurf.IsNull()) {
 404                 di << "BSplineSurface is not created!" << "\n";
 405                 return;
 406         }
 407         di << "Number of UPoles:" << aBSurf->NbUPoles() << "\n";
 408         QCOMPARE (aBSurf->NbUPoles(), theTargetNbUPoles);
 409 }
 410
 411 struct aData {
 412         Draw_Interpretor* di;
 413         Standard_Integer nbupoles;
 414         Standard_CString filename;
 415 };
 416
 417 Standard_EXPORT Standard_Address convert(Standard_Address data)
 418 {
 419         aData* info = (aData*) data;
 420         GeomConvertTest(*(info->di),info->nbupoles,info->filename);
 421         return NULL;
 422 }
 423
 424 static Standard_Integer OCC23952sweep (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 425 {
 426         if (argc != 3) {
 427                 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
 428                 return 1;
 429         }
 430         struct aData aStorage;
 431         aStorage.di = &di;
 432         aStorage.nbupoles = Draw::Atoi(argv[1]);
 433         aStorage.filename = argv[2];
 434
 435         OSD_Thread aThread1(convert);
 436         aThread1.Run(&aStorage);
 437         GeomConvertTest(di,aStorage.nbupoles,aStorage.filename);
 438         cout << "result of thread: " << aThread1.Wait() << endl;
 439
 440         return 0;
 441 }
 442
 443 #include <GeomInt_IntSS.hxx>
 444 static void GeomIntSSTest (Draw_Interpretor& di, Standard_Integer theNbSol, Standard_CString theFileName1, Standard_CString theFileName2)
 445 {
 446         Handle(Geom_Surface) aSurf1 = DrawTrSurf::GetSurface(theFileName1);
 447         Handle(Geom_Surface) aSurf2 = DrawTrSurf::GetSurface(theFileName2);
 448         GeomInt_IntSS anInter;
 449         anInter.Perform(aSurf1, aSurf2, Precision::Confusion(), Standard_True);
 450         if (!anInter.IsDone()) {
 451                 di << "An intersection is not done!" << "\n";
 452                 return;
 453         }
 454
 455         di << "Number of Lines:" << anInter.NbLines() << "\n";
 456         QCOMPARE (anInter.NbLines(), theNbSol);
 457 }
 458
 459 struct aNewData {
 460         Draw_Interpretor* di;
 461         Standard_Integer nbsol;
 462         Standard_CString filename1;
 463         Standard_CString filename2;
 464 };
 465 Standard_EXPORT Standard_Address convert_inter(Standard_Address data)
 466 {
 467         aNewData* info = (aNewData*) data;
 468         GeomIntSSTest(*(info->di),info->nbsol,info->filename1,info->filename2);
 469         return NULL;
 470 }
 471
 472 static Standard_Integer OCC23952intersect (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 473 {
 474         if (argc != 4) {
 475                 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
 476                 return 1;
 477         }
 478         struct aNewData aStorage;
 479         aStorage.di = &di;
 480         aStorage.nbsol = Draw::Atoi(argv[1]);
 481         aStorage.filename1 = argv[2];
 482         aStorage.filename2 = argv[3];
 483
 484         OSD_Thread aThread1(convert_inter);
 485         aThread1.Run(&aStorage);
 486         GeomIntSSTest(di,aStorage.nbsol,aStorage.filename1,aStorage.filename2);
 487         cout << "result of thread: " << aThread1.Wait() << endl;
 488
 489         return 0;
 490 }
 491
 492 #include <Geom_SurfaceOfRevolution.hxx>
 493 static Standard_Integer OCC23683 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
 494 {
 495   if (argc < 2) {
 496     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
 497     return 1;
 498   }
 499
 500   Standard_Integer ucontinuity = 1;
 501   Standard_Integer vcontinuity = 1;
 502   Standard_Boolean iscnu = false;
 503   Standard_Boolean iscnv = false;
 504
 505   Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
 506
 507   QCOMPARE (aSurf->IsCNu (ucontinuity), iscnu);
 508   QCOMPARE (aSurf->IsCNv (vcontinuity), iscnv);
 509
 510   return 0;
 511 }
 512
 513 #include <gp_Ax1.hxx>
 514 #include <gp_Ax22d.hxx>
 515 #include <Geom_Plane.hxx>
 516 #include <Geom2d_Circle.hxx>
 517 #include <Geom2d_TrimmedCurve.hxx>
 518 #include <BRepBuilderAPI_MakeEdge.hxx>
 519 #include <BRepPrimAPI_MakeRevol.hxx>
 520 #include <Geom2d_OffsetCurve.hxx>
 521
 522 static int test_offset(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
 523 {
 524   // Check the command arguments
 525   if ( argc != 1 )
 526   {
 527     di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
 528     di << "Usage: type help " << argv[0] << "\n";
 529     return 1; // TCL_ERROR
 530   }
 531
 532   gp_Ax1 RotoAx( gp::Origin(), gp::DZ() );
 533   gp_Ax22d Ax2( gp::Origin2d(), gp::DY2d(), gp::DX2d() );
 534   Handle(Geom_Surface) Plane = new Geom_Plane( gp::YOZ() );
 535
 536   di << "<<<< Preparing sample surface of revolution based on trimmed curve >>>>" << "\n";
 537   di << "-----------------------------------------------------------------------" << "\n";
 538
 539   Handle(Geom2d_Circle) C2d1 = new Geom2d_Circle(Ax2, 1.0);
 540   Handle(Geom2d_TrimmedCurve) C2d1Trimmed = new Geom2d_TrimmedCurve(C2d1, 0.0, M_PI/2.0);
 541   TopoDS_Edge E1 = BRepBuilderAPI_MakeEdge(C2d1Trimmed, Plane);
 542
 543   DBRep::Set("e1", E1);
 544
 545   BRepPrimAPI_MakeRevol aRevolBuilder1(E1, RotoAx);
 546   TopoDS_Face F1 = TopoDS::Face( aRevolBuilder1.Shape() );
 547
 548   DBRep::Set("f1", F1);
 549
 550   di << "Result: f1" << "\n";
 551
 552   di << "<<<< Preparing sample surface of revolution based on offset curve  >>>>" << "\n";
 553   di << "-----------------------------------------------------------------------" << "\n";
 554
 555   Handle(Geom2d_OffsetCurve) C2d2Offset = new Geom2d_OffsetCurve(C2d1Trimmed, -0.5);
 556   TopoDS_Edge E2 = BRepBuilderAPI_MakeEdge(C2d2Offset, Plane);
 557
 558   DBRep::Set("e2", E2);
 559
 560   BRepPrimAPI_MakeRevol aRevolBuilder2(E2, RotoAx);
 561   TopoDS_Face F2 = TopoDS::Face( aRevolBuilder2.Shape() );
 562
 563   DBRep::Set("f2", F2);
 564
 565   di << "Result: f2" << "\n";
 566
 567   return 0;
 568 }
 569
 570 #include <GeomAdaptor_Surface.hxx>
 571 #include <Draw.hxx>
 572 //=======================================================================
 573 //function : OCC23945
 574 //purpose  :
 575 //=======================================================================
 576
 577 static Standard_Integer OCC23945 (Draw_Interpretor& di,Standard_Integer n, const char** a)
 578 {
 579   if (n < 5) return 1;
 580
 581   Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(a[1]);
 582   if (aS.IsNull()) return 1;
 583
 584   GeomAdaptor_Surface GS(aS);
 585
 586   Standard_Real U = Draw::Atof(a[2]);
 587   Standard_Real V = Draw::Atof(a[3]);
 588
 589   Standard_Boolean DrawPoint = ( n%3 == 2);
 590   if ( DrawPoint) n--;
 591
 592   gp_Pnt P;
 593   if (n >= 13) {
 594     gp_Vec DU,DV;
 595     if (n >= 22) {
 596       gp_Vec D2U,D2V,D2UV;
 597       GS.D2(U,V,P,DU,DV,D2U,D2V,D2UV);
 598       Draw::Set(a[13],D2U.X());
 599       Draw::Set(a[14],D2U.Y());
 600       Draw::Set(a[15],D2U.Z());
 601       Draw::Set(a[16],D2V.X());
 602       Draw::Set(a[17],D2V.Y());
 603       Draw::Set(a[18],D2V.Z());
 604       Draw::Set(a[19],D2UV.X());
 605       Draw::Set(a[20],D2UV.Y());
 606       Draw::Set(a[21],D2UV.Z());
 607     }
 608     else
 609       GS.D1(U,V,P,DU,DV);
 610
 611     Draw::Set(a[7],DU.X());
 612     Draw::Set(a[8],DU.Y());
 613     Draw::Set(a[9],DU.Z());
 614     Draw::Set(a[10],DV.X());
 615     Draw::Set(a[11],DV.Y());
 616     Draw::Set(a[12],DV.Z());
 617   }
 618   else
 619     GS.D0(U,V,P);
 620
 621   if ( n > 6) {
 622     Draw::Set(a[4],P.X());
 623     Draw::Set(a[5],P.Y());
 624     Draw::Set(a[6],P.Z());
 625   }
 626   if ( DrawPoint) {
 627     DrawTrSurf::Set(a[n],P);
 628   }
 629
 630   return 0;
 631 }
 632
 633 #include <Voxel_BoolDS.hxx>
 634 #include <Voxel_FastConverter.hxx>
 635 #include <Voxel_BooleanOperation.hxx>
 636 static Standard_Integer OCC24019 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
 637 {
 638   if ( argc != 2 ) {
 639     di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
 640         return 1;
 641   }
 642
 643   TCollection_AsciiString aFileName = argv[1];
 644   TopoDS_Shape aShape;
 645   BRep_Builder aBuilder;
 646
 647   if (!BRepTools::Read(aShape, aFileName.ToCString(), aBuilder)) {
 648     di << "Error: Could not read a shape!" << "\n";
 649         return 1;
 650   }
 651
 652   TopoDS_Solid aShape1 = BRepPrimAPI_MakeSphere(gp_Pnt(20,25,35), 7);
 653
 654   Standard_Real deflection = 0.005;
 655   Standard_Integer nbThreads = 1;
 656   Standard_Integer nbx = 200, nby = 200, nbz = 200;
 657   Voxel_BoolDS theVoxels(0,0,0, 50, 50, 50, nbx, nby, nbz);
 658   Voxel_BoolDS theVoxels1(0,0,0, 50, 50, 50, nbx, nby, nbz);
 659
 660   Standard_Integer progress = 0;
 661   Voxel_FastConverter fcp(aShape, theVoxels, deflection, nbx, nby, nbz, nbThreads);
 662   fcp.ConvertUsingSAT(progress, 1);
 663   fcp.FillInVolume(1);
 664
 665   Voxel_FastConverter fcp1(aShape1, theVoxels1, deflection, nbx, nby, nbz, nbThreads);
 666   fcp1.ConvertUsingSAT(progress, 1);
 667   fcp1.FillInVolume(1);
 668
 669   Voxel_BooleanOperation op;
 670   Standard_Boolean result = op.Cut(theVoxels1, theVoxels);
 671   if ( result != 1 ) {
 672     di << "Error: invalid boolean operation" << "\n";
 673   } else {
 674         di << "OK: boolean operation is ok" << "\n";
 675   }
 676
 677   return 0;
 678 }
 679
 680 void QABugs::Commands_19(Draw_Interpretor& theCommands) {
 681   const char *group = "QABugs";
 682
 683   theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
 684   theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
 685   theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
 686   theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group);
 687   theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
 688   theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
 689   theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
 690   theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
 691   theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
 692   theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
 693   theCommands.Add ("OCC23952sweep", "OCC23952sweep nbupoles shape", __FILE__, OCC23952sweep, group);
 694   theCommands.Add ("OCC23952intersect", "OCC23952intersect nbsol shape1 shape2", __FILE__, OCC23952intersect, group);
 695   theCommands.Add ("test_offset", "test_offset", __FILE__, test_offset, group);
 696   theCommands.Add("OCC23945", "OCC23945 surfname U V X Y Z [DUX DUY DUZ DVX DVY DVZ [D2UX D2UY D2UZ D2VX D2VY D2VZ D2UVX D2UVY D2UVZ]]", __FILE__, OCC23945,group);
 697   theCommands.Add ("OCC24019", "OCC24019 aShape", __FILE__, OCC24019, group);
 698
 699   return;
 700 }
 701
 702