src/QABugs/QABugs_19.cxx

   1 // Created on: 2002-05-21
   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_Interpretor.hxx>
  19 #include <DBRep.hxx>
  20 #include <DrawTrSurf.hxx>
  21 #include <ViewerTest.hxx>
  22 #include <TopoDS_Shape.hxx>
  23 #include <AIS_InteractiveContext.hxx>
  24 #include <AIS_TexturedShape.hxx>
  25 #include <Image_PixMap.hxx>
  26 #include <Image_Color.hxx>
  27
  28 #include <gp_Pnt2d.hxx>
  29 #include <gp_Ax1.hxx>
  30 #include <GCE2d_MakeSegment.hxx>
  31 #include <Geom2d_TrimmedCurve.hxx>
  32 #include <DrawTrSurf.hxx>
  33
  34 #include <Precision.hxx>
  35
  36 #include <PCollection_HAsciiString.hxx>
  37
  38 #include <cstdio>
  39 #include <cmath>
  40 #include <iostream>
  41 #include <OSD_PerfMeter.hxx>
  42 #include <OSD_Timer.hxx>
  43 #include <BRepPrimAPI_MakeBox.hxx>
  44 #include <BRepPrimAPI_MakeSphere.hxx>
  45 #include <BRepAlgo_Cut.hxx>
  46 #include <NCollection_Map.hxx>
  47 #include <NCollection_Handle.hxx>
  48 #include <TCollection_HAsciiString.hxx>
  49 #include <GeomFill_Trihedron.hxx>
  50 #include <BRepOffsetAPI_MakePipe.hxx>
  51
  52 #include <Standard_Version.hxx>
  53
  54 #define QCOMPARE(val1, val2) \
  55   di << "Checking " #val1 " == " #val2 << \
  56         ((val1) == (val2) ? ": OK\n" : ": Error\n")
  57
  58 static Standard_Integer OCC230 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
  59 {
  60   if ( argc != 4) {
  61     di << "ERROR OCC230: Usage : " << argv[0] << " TrimmedCurve Pnt2d Pnt2d" << "\n";
  62     return 1;
  63   }
  64
  65   gp_Pnt2d P1, P2;
  66   if ( !DrawTrSurf::GetPoint2d(argv[2],P1)) {
  67     di << "ERROR OCC230: " << argv[2] << " is not Pnt2d" << "\n";
  68     return 1;
  69   }
  70   if ( !DrawTrSurf::GetPoint2d(argv[3],P2)) {
  71     di << "ERROR OCC230: " << argv[3] << " is not Pnt2d" << "\n";
  72     return 1;
  73   }
  74
  75   GCE2d_MakeSegment MakeSegment(P1,P2);
  76   Handle(Geom2d_TrimmedCurve) TrimmedCurve = MakeSegment.Value();
  77   DrawTrSurf::Set(argv[1], TrimmedCurve);
  78   return 0;
  79 }
  80
  81 static Standard_Integer OCC142 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
  82 {
  83   for(Standard_Integer i= 0;i <= 20;i++){
  84     Handle(PCollection_HAsciiString) pstr = new PCollection_HAsciiString("TEST");
  85     pstr->Clear();
  86   }
  87   di << "OCC142: OK" << "\n";
  88   return 0;
  89 }
  90
  91 static Standard_Integer OCC23361 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
  92 {
  93   gp_Pnt p(0, 0, 2);
  94
  95   gp_Trsf t1, t2;
  96   t1.SetRotation(gp_Ax1(p, gp_Dir(0, 1, 0)), -0.49328285294022267);
  97   t2.SetRotation(gp_Ax1(p, gp_Dir(0, 0, 1)), 0.87538474718473880);
  98
  99   gp_Trsf tComp = t2 * t1;
 100
 101   gp_Pnt p1(10, 3, 4);
 102   gp_Pnt p2 = p1.Transformed(tComp);
 103   gp_Pnt p3 = p1.Transformed(t1);
 104   p3.Transform(t2);
 105
 106   // points must be equal
 107   if ( ! p2.IsEqual(p3, Precision::Confusion()) )
 108     di << "ERROR OCC23361: equivalent transformations does not produce equal points" << "\n";
 109   else
 110     di << "OCC23361: OK" << "\n";
 111
 112   return 0;
 113 }
 114
 115 static Standard_Integer OCC23237 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
 116 {
 117   OSD_PerfMeter aPM("TestMeter",0);
 118   OSD_Timer aTM;
 119
 120   // run some operation in cycle for about 2 seconds to have good values of times to compare
 121   int count = 0;
 122   printf("OSD_PerfMeter test.\nRunning Boolean operation on solids in loop.\n");
 123   for (; aTM.ElapsedTime() < 2.; count++)
 124   {
 125     aPM.Start();
 126     aTM.Start();
 127
 128     // do some operation that will take considerable time compared with time of starting / stopping timers
 129     BRepPrimAPI_MakeBox aBox (10., 10., 10.);
 130     BRepPrimAPI_MakeSphere aSphere (10.);
 131     BRepAlgo_Cut aCutter (aBox.Shape(), aSphere.Shape());
 132
 133     aTM.Stop();
 134     aPM.Stop();
 135   }
 136
 137   int aNbEnters = 0;
 138   Standard_Real aPerfMeter_CPUtime = 0., aTimer_ElapsedTime = aTM.ElapsedTime();
 139
 140   perf_get_meter("TestMeter", &aNbEnters, &aPerfMeter_CPUtime);
 141
 142   Standard_Real aTimeDiff = (fabs(aTimer_ElapsedTime - aPerfMeter_CPUtime) / aTimer_ElapsedTime);
 143
 144   printf("\nMeasurement results (%d cycles):\n", count);
 145   printf("\nOSD_PerfMeter CPU time: %lf\nOSD_Timer elapsed time: %lf\n", aPerfMeter_CPUtime, aTimer_ElapsedTime);
 146   printf("Time delta is: %.3lf %%\n", aTimeDiff * 100);
 147
 148   if (aTimeDiff > 0.2)
 149     di << "OCC23237: Error: too much difference between CPU and elapsed times";
 150   else if (aNbEnters != count)
 151     di << "OCC23237: Error: counter reported by PerfMeter (" << aNbEnters << ") does not correspond to actual number of cycles";
 152   else
 153     di << "OCC23237: OK";
 154
 155   return 0;
 156 }
 157
 158 #ifdef HAVE_TBB
 159
 160 #include <Standard_Atomic.hxx>
 161 #include <tbb/blocked_range.h>
 162 #include <tbb/parallel_for.h>
 163
 164 class IncrementerDecrementer
 165 {
 166 public:
 167     IncrementerDecrementer (Standard_Integer* theVal, Standard_Boolean thePositive) : myVal (theVal), myPositive (thePositive)
 168     {}
 169     void operator() (const tbb::blocked_range<size_t>& r) const
 170     {
 171         if (myPositive)
 172             for (size_t i = r.begin(); i != r.end(); ++i)
 173                 Standard_Atomic_Increment (myVal);
 174         else
 175             for (size_t i = r.begin(); i != r.end(); ++i)
 176                 Standard_Atomic_Decrement (myVal);
 177     }
 178 private:
 179     Standard_Integer*   myVal;
 180     Standard_Boolean   myPositive;
 181 };
 182 #endif
 183
 184 #ifdef HAVE_TBB
 185 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
 186 {
 187   int aSum = 0;
 188
 189   //check returned value
 190   QCOMPARE (Standard_Atomic_Decrement (&aSum), -1);
 191   QCOMPARE (Standard_Atomic_Increment (&aSum), 0);
 192   QCOMPARE (Standard_Atomic_Increment (&aSum), 1);
 193   QCOMPARE (Standard_Atomic_Increment (&aSum), 2);
 194 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 0);
 195 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 1);
 196
 197   //check atomicity
 198   aSum = 0;
 199   const int N = 1 << 24; //big enough to ensure concurrency
 200
 201   //increment
 202   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, true));
 203   QCOMPARE (aSum, N);
 204
 205   //decrement
 206   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
 207   QCOMPARE (aSum, 0);
 208
 209   return 0;
 210 }
 211
 212 #else /* HAVE_TBB */
 213
 214 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
 215 {
 216   di << "Test skipped: command " << argv[0] << " requires TBB library\n";
 217   return 0;
 218 }
 219
 220 #endif /* HAVE_TBB */
 221
 222 #include <TDocStd_Application.hxx>
 223 #include <XCAFApp_Application.hxx>
 224 #include <TDocStd_Document.hxx>
 225 #include <XCAFDoc_ShapeTool.hxx>
 226 #include <XCAFDoc_DocumentTool.hxx>
 227 #include <TDF_Label.hxx>
 228 #include <TDataStd_Name.hxx>
 229
 230 static Standard_Integer OCC23595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
 231 {
 232   const Handle(TDocStd_Application)& anApp = XCAFApp_Application::GetApplication();
 233   Handle(TDocStd_Document) aDoc;
 234   anApp->NewDocument ("XmlXCAF", aDoc);
 235   QCOMPARE (!aDoc.IsNull(), Standard_True);
 236
 237   Handle(XCAFDoc_ShapeTool) aShTool = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main());
 238
 239   //check default value
 240   Standard_Boolean aValue = XCAFDoc_ShapeTool::AutoNaming();
 241   QCOMPARE (aValue, Standard_True);
 242
 243   //true
 244   XCAFDoc_ShapeTool::SetAutoNaming (Standard_True);
 245   TopoDS_Shape aShape = BRepPrimAPI_MakeBox (100., 200., 300.).Shape();
 246   TDF_Label aLabel = aShTool->AddShape (aShape);
 247   Handle(TDataStd_Name) anAttr;
 248   QCOMPARE (aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
 249
 250   //false
 251   XCAFDoc_ShapeTool::SetAutoNaming (Standard_False);
 252   aShape = BRepPrimAPI_MakeBox (300., 200., 100.).Shape();
 253   aLabel = aShTool->AddShape (aShape);
 254   QCOMPARE (!aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
 255
 256   //restore
 257   XCAFDoc_ShapeTool::SetAutoNaming (aValue);
 258
 259   return 0;
 260 }
 261
 262 #include <ExprIntrp_GenExp.hxx>
 263 Standard_Integer OCC22611 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 264 {
 265
 266   if (argc != 3) {
 267     di << "Usage : " << argv[0] << " string nb\n";
 268     return 1;
 269   }
 270
 271   TCollection_AsciiString aToken = argv[1];
 272   Standard_Integer aNb = atoi(argv[2]);
 273
 274   Handle(ExprIntrp_GenExp) aGen = ExprIntrp_GenExp::Create();
 275   for (Standard_Integer i=0; i < aNb; i++)
 276   {
 277     aGen->Process(aToken);
 278     Handle(Expr_GeneralExpression) aExpr = aGen->Expression();
 279   }
 280
 281   return 0;
 282 }
 283
 284 Standard_Integer OCC22595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
 285 {
 286   gp_Mat M0;
 287   di << "M0 = "
 288   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}"
 289   << " {" << M0(2,1) << "} {" << M0(2,2) << "} {" << M0(2,3) <<"}"
 290   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}";
 291   return 0;
 292 }
 293
 294 #include <TopoDS_Face.hxx>
 295 #include <TopoDS_Face.hxx>
 296 #include <TopoDS.hxx>
 297 #include <BRepBuilderAPI_Transform.hxx>
 298 #include <BRepExtrema_DistShapeShape.hxx>
 299 #include <BRepTools.hxx>
 300
 301 static Standard_Boolean OCC23774Test(const TopoDS_Face& grossPlateFace, const TopoDS_Shape& originalWire, Draw_Interpretor& di)
 302 {
 303   BRepExtrema_DistShapeShape distShapeShape(grossPlateFace,originalWire,Extrema_ExtFlag_MIN);
 304   if(!distShapeShape.IsDone()) {
 305     di <<"Distance ShapeShape is Not Done\n";
 306     return Standard_False;
 307   }
 308
 309   if(distShapeShape.Value() > 0.01) {
 310     di << "Wrong Dist = " <<distShapeShape.Value() << "\n";
 311     return Standard_False;
 312   } else
 313     di << "Dist0 = " <<distShapeShape.Value() <<"\n";
 314
 315   //////////////////////////////////////////////////////////////////////////
 316   /// First Flip Y
 317   const gp_Pnt2d axis1P1(1474.8199035519228,1249.9995745636970);
 318   const gp_Pnt2d axis1P2(1474.8199035519228,1250.9995745636970);
 319
 320   gp_Vec2d mirrorVector1(axis1P1,axis1P2);
 321
 322   gp_Trsf2d mirror1;
 323   mirror1.SetMirror(gp_Ax2d(axis1P1,mirrorVector1));
 324
 325   BRepBuilderAPI_Transform transformer1(mirror1);
 326   transformer1.Perform(originalWire);
 327   if(!transformer1.IsDone()) {
 328     di << "Not Done1 " << "\n";
 329     return Standard_False;
 330   }
 331   const TopoDS_Shape& step1ModifiedShape = transformer1.ModifiedShape(originalWire);
 332
 333   BRepExtrema_DistShapeShape distShapeShape1(grossPlateFace,step1ModifiedShape,Extrema_ExtFlag_MIN);
 334   if(!distShapeShape1.IsDone())
 335     return Standard_False;
 336   if(distShapeShape1.Value() > 0.01) {
 337     di << "Dist = " <<distShapeShape1.Value() <<"\n";
 338     return Standard_False;
 339   } else
 340     di << "Dist1 = " <<distShapeShape1.Value() <<"\n";
 341
 342   //////////////////////////////////////////////////////////////////////////
 343   /// Second flip Y
 344   transformer1.Perform(step1ModifiedShape);
 345   if(!transformer1.IsDone()) {
 346     di << "Not Done1 \n";
 347     return Standard_False;
 348   }
 349   const TopoDS_Shape& step2ModifiedShape = transformer1.ModifiedShape(step1ModifiedShape);
 350
 351   //This is identity matrix for values but for type is gp_Rotation ?!
 352   gp_Trsf2d mirror11 = mirror1;
 353   mirror11.PreMultiply(mirror1);
 354
 355   BRepExtrema_DistShapeShape distShapeShape2(grossPlateFace,step2ModifiedShape);//,Extrema_ExtFlag_MIN);
 356   if(!distShapeShape2.IsDone())
 357     return Standard_False;
 358
 359   //This last test case give error (the value is 1008.8822038689706)
 360   if(distShapeShape2.Value() > 0.01) {
 361     di  << "Wrong Dist2 = " <<distShapeShape2.Value() <<"\n";
 362     Standard_Integer N = distShapeShape2.NbSolution();
 363     di << "Nb = " <<N <<"\n";
 364     for (Standard_Integer i=1;i <= N;i++)
 365         di <<"Sol(" <<i<<") = " <<distShapeShape2.PointOnShape1(i).Distance(distShapeShape2.PointOnShape2(i)) <<"\n";
 366     return Standard_False;
 367   }
 368   di << "Distance2 = " <<distShapeShape2.Value() <<"\n";
 369
 370   return Standard_True;
 371 }
 372 static Standard_Integer OCC23774(Draw_Interpretor& di, Standard_Integer n, const char** a)
 373 {
 374
 375   if (n != 3) {
 376         di <<"OCC23774: invalid number of input parameters\n";
 377         return 1;
 378   }
 379
 380   const char *ns1 = (a[1]), *ns2 = (a[2]);
 381   TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2));
 382   if (S1.IsNull() || S2.IsNull()) {
 383         di <<"OCC23774: Null input shapes\n";
 384         return 1;
 385   }
 386   const TopoDS_Face& aFace  = TopoDS::Face(S1);
 387   if(!OCC23774Test(aFace, S2, di))
 388         di << "Something is wrong\n";
 389
 390  return 0;
 391 }
 392
 393 #include <GeomConvert_ApproxSurface.hxx>
 394 #include <Geom_BSplineSurface.hxx>
 395 #include <Draw.hxx>
 396 #include <OSD_Thread.hxx>
 397 static void GeomConvertTest (Draw_Interpretor& di, Standard_Integer theTargetNbUPoles, Standard_CString theFileName)
 398 {
 399         Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(theFileName);
 400         GeomConvert_ApproxSurface aGAS (aSurf, 1e-4, GeomAbs_C1, GeomAbs_C1, 9, 9, 100, 1);
 401         if (!aGAS.IsDone()) {
 402                 di << "ApproxSurface is not done!" << "\n";
 403                 return;
 404         }
 405         const Handle(Geom_BSplineSurface)& aBSurf = aGAS.Surface();
 406         if (aBSurf.IsNull()) {
 407                 di << "BSplineSurface is not created!" << "\n";
 408                 return;
 409         }
 410         di << "Number of UPoles:" << aBSurf->NbUPoles() << "\n";
 411         QCOMPARE (aBSurf->NbUPoles(), theTargetNbUPoles);
 412 }
 413
 414 struct aData {
 415         Draw_Interpretor* di;
 416         Standard_Integer nbupoles;
 417         Standard_CString filename;
 418 };
 419
 420 Standard_EXPORT Standard_Address convert(Standard_Address data)
 421 {
 422         aData* info = (aData*) data;
 423         GeomConvertTest(*(info->di),info->nbupoles,info->filename);
 424         return NULL;
 425 }
 426
 427 static Standard_Integer OCC23952sweep (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 428 {
 429         if (argc != 3) {
 430                 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
 431                 return 1;
 432         }
 433         struct aData aStorage;
 434         aStorage.di = &di;
 435         aStorage.nbupoles = Draw::Atoi(argv[1]);
 436         aStorage.filename = argv[2];
 437
 438         OSD_Thread aThread1(convert);
 439         aThread1.Run(&aStorage);
 440         GeomConvertTest(di,aStorage.nbupoles,aStorage.filename);
 441         cout << "result of thread: " << aThread1.Wait() << endl;
 442
 443         return 0;
 444 }
 445
 446 #include <GeomInt_IntSS.hxx>
 447 static void GeomIntSSTest (Draw_Interpretor& di, Standard_Integer theNbSol, Standard_CString theFileName1, Standard_CString theFileName2)
 448 {
 449         Handle(Geom_Surface) aSurf1 = DrawTrSurf::GetSurface(theFileName1);
 450         Handle(Geom_Surface) aSurf2 = DrawTrSurf::GetSurface(theFileName2);
 451         GeomInt_IntSS anInter;
 452         anInter.Perform(aSurf1, aSurf2, Precision::Confusion(), Standard_True);
 453         if (!anInter.IsDone()) {
 454                 di << "An intersection is not done!" << "\n";
 455                 return;
 456         }
 457
 458         di << "Number of Lines:" << anInter.NbLines() << "\n";
 459         QCOMPARE (anInter.NbLines(), theNbSol);
 460 }
 461
 462 struct aNewData {
 463         Draw_Interpretor* di;
 464         Standard_Integer nbsol;
 465         Standard_CString filename1;
 466         Standard_CString filename2;
 467 };
 468 Standard_EXPORT Standard_Address convert_inter(Standard_Address data)
 469 {
 470         aNewData* info = (aNewData*) data;
 471         GeomIntSSTest(*(info->di),info->nbsol,info->filename1,info->filename2);
 472         return NULL;
 473 }
 474
 475 static Standard_Integer OCC23952intersect (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 476 {
 477         if (argc != 4) {
 478                 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
 479                 return 1;
 480         }
 481         struct aNewData aStorage;
 482         aStorage.di = &di;
 483         aStorage.nbsol = Draw::Atoi(argv[1]);
 484         aStorage.filename1 = argv[2];
 485         aStorage.filename2 = argv[3];
 486
 487         OSD_Thread aThread1(convert_inter);
 488         aThread1.Run(&aStorage);
 489         GeomIntSSTest(di,aStorage.nbsol,aStorage.filename1,aStorage.filename2);
 490         cout << "result of thread: " << aThread1.Wait() << endl;
 491
 492         return 0;
 493 }
 494
 495 #include <Geom_SurfaceOfRevolution.hxx>
 496 static Standard_Integer OCC23683 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
 497 {
 498   if (argc < 2) {
 499     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
 500     return 1;
 501   }
 502
 503   Standard_Integer ucontinuity = 1;
 504   Standard_Integer vcontinuity = 1;
 505   Standard_Boolean iscnu = false;
 506   Standard_Boolean iscnv = false;
 507
 508   Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
 509
 510   QCOMPARE (aSurf->IsCNu (ucontinuity), iscnu);
 511   QCOMPARE (aSurf->IsCNv (vcontinuity), iscnv);
 512
 513   return 0;
 514 }
 515
 516 #include <gp_Ax1.hxx>
 517 #include <gp_Ax22d.hxx>
 518 #include <Geom_Plane.hxx>
 519 #include <Geom2d_Circle.hxx>
 520 #include <Geom2d_TrimmedCurve.hxx>
 521 #include <BRepBuilderAPI_MakeEdge.hxx>
 522 #include <BRepPrimAPI_MakeRevol.hxx>
 523 #include <Geom2d_OffsetCurve.hxx>
 524
 525 static int test_offset(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
 526 {
 527   // Check the command arguments
 528   if ( argc != 1 )
 529   {
 530     di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
 531     di << "Usage: type help " << argv[0] << "\n";
 532     return 1; // TCL_ERROR
 533   }
 534
 535   gp_Ax1 RotoAx( gp::Origin(), gp::DZ() );
 536   gp_Ax22d Ax2( gp::Origin2d(), gp::DY2d(), gp::DX2d() );
 537   Handle(Geom_Surface) Plane = new Geom_Plane( gp::YOZ() );
 538
 539   di << "<<<< Preparing sample surface of revolution based on trimmed curve >>>>" << "\n";
 540   di << "-----------------------------------------------------------------------" << "\n";
 541
 542   Handle(Geom2d_Circle) C2d1 = new Geom2d_Circle(Ax2, 1.0);
 543   Handle(Geom2d_TrimmedCurve) C2d1Trimmed = new Geom2d_TrimmedCurve(C2d1, 0.0, M_PI/2.0);
 544   TopoDS_Edge E1 = BRepBuilderAPI_MakeEdge(C2d1Trimmed, Plane);
 545
 546   DBRep::Set("e1", E1);
 547
 548   BRepPrimAPI_MakeRevol aRevolBuilder1(E1, RotoAx);
 549   TopoDS_Face F1 = TopoDS::Face( aRevolBuilder1.Shape() );
 550
 551   DBRep::Set("f1", F1);
 552
 553   di << "Result: f1" << "\n";
 554
 555   di << "<<<< Preparing sample surface of revolution based on offset curve  >>>>" << "\n";
 556   di << "-----------------------------------------------------------------------" << "\n";
 557
 558   Handle(Geom2d_OffsetCurve) C2d2Offset = new Geom2d_OffsetCurve(C2d1Trimmed, -0.5);
 559   TopoDS_Edge E2 = BRepBuilderAPI_MakeEdge(C2d2Offset, Plane);
 560
 561   DBRep::Set("e2", E2);
 562
 563   BRepPrimAPI_MakeRevol aRevolBuilder2(E2, RotoAx);
 564   TopoDS_Face F2 = TopoDS::Face( aRevolBuilder2.Shape() );
 565
 566   DBRep::Set("f2", F2);
 567
 568   di << "Result: f2" << "\n";
 569
 570   return 0;
 571 }
 572
 573 #include <Geom_Curve.hxx>
 574 #include <Geom_Surface.hxx>
 575 #include <Precision.hxx>
 576 #include <ShapeConstruct_ProjectCurveOnSurface.hxx>
 577 //=======================================================================
 578 //function : OCC24008
 579 //purpose  :
 580 //=======================================================================
 581 static Standard_Integer OCC24008 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 582 {
 583   if (argc != 3) {
 584     di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
 585     return 1;
 586   }
 587   Handle(Geom_Curve) aCurve = DrawTrSurf::GetCurve(argv[1]);
 588   Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[2]);
 589   if (aCurve.IsNull()) {
 590     di << "Curve was not read" << "\n";
 591         return 1;
 592   }
 593   if (aSurf.IsNull()) {
 594         di << "Surface was not read" << "\n";
 595         return 1;
 596   }
 597   ShapeConstruct_ProjectCurveOnSurface aProj;
 598   aProj.Init (aSurf, Precision::Confusion());
 599   try {
 600     Handle(Geom2d_Curve) aPCurve;
 601     aProj.Perform (aCurve, aCurve->FirstParameter(), aCurve->LastParameter(), aPCurve);
 602     if (aPCurve.IsNull()) {
 603           di << "PCurve was not created" << "\n";
 604           return 1;
 605     }
 606   } catch (...) {
 607     di << "Exception was caught" << "\n";
 608   }
 609   return 0;
 610 }
 611
 612 #include <GeomAdaptor_Surface.hxx>
 613 #include <Draw.hxx>
 614 //=======================================================================
 615 //function : OCC23945
 616 //purpose  :
 617 //=======================================================================
 618
 619 static Standard_Integer OCC23945 (Draw_Interpretor& /*di*/,Standard_Integer n, const char** a)
 620 {
 621   if (n < 5) return 1;
 622
 623   Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(a[1]);
 624   if (aS.IsNull()) return 1;
 625
 626   GeomAdaptor_Surface GS(aS);
 627
 628   Standard_Real U = Draw::Atof(a[2]);
 629   Standard_Real V = Draw::Atof(a[3]);
 630
 631   Standard_Boolean DrawPoint = ( n%3 == 2);
 632   if ( DrawPoint) n--;
 633
 634   gp_Pnt P;
 635   if (n >= 13) {
 636     gp_Vec DU,DV;
 637     if (n >= 22) {
 638       gp_Vec D2U,D2V,D2UV;
 639       GS.D2(U,V,P,DU,DV,D2U,D2V,D2UV);
 640       Draw::Set(a[13],D2U.X());
 641       Draw::Set(a[14],D2U.Y());
 642       Draw::Set(a[15],D2U.Z());
 643       Draw::Set(a[16],D2V.X());
 644       Draw::Set(a[17],D2V.Y());
 645       Draw::Set(a[18],D2V.Z());
 646       Draw::Set(a[19],D2UV.X());
 647       Draw::Set(a[20],D2UV.Y());
 648       Draw::Set(a[21],D2UV.Z());
 649     }
 650     else
 651       GS.D1(U,V,P,DU,DV);
 652
 653     Draw::Set(a[7],DU.X());
 654     Draw::Set(a[8],DU.Y());
 655     Draw::Set(a[9],DU.Z());
 656     Draw::Set(a[10],DV.X());
 657     Draw::Set(a[11],DV.Y());
 658     Draw::Set(a[12],DV.Z());
 659   }
 660   else
 661     GS.D0(U,V,P);
 662
 663   if ( n > 6) {
 664     Draw::Set(a[4],P.X());
 665     Draw::Set(a[5],P.Y());
 666     Draw::Set(a[6],P.Z());
 667   }
 668   if ( DrawPoint) {
 669     DrawTrSurf::Set(a[n],P);
 670   }
 671
 672   return 0;
 673 }
 674
 675 #include <Voxel_BoolDS.hxx>
 676 #include <Voxel_FastConverter.hxx>
 677 #include <Voxel_BooleanOperation.hxx>
 678 static Standard_Integer OCC24019 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
 679 {
 680   if ( argc != 2 ) {
 681     di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
 682         return 1;
 683   }
 684
 685   TCollection_AsciiString aFileName = argv[1];
 686   TopoDS_Shape aShape;
 687   BRep_Builder aBuilder;
 688
 689   if (!BRepTools::Read(aShape, aFileName.ToCString(), aBuilder)) {
 690     di << "Error: Could not read a shape!" << "\n";
 691         return 1;
 692   }
 693
 694   TopoDS_Solid aShape1 = BRepPrimAPI_MakeSphere(gp_Pnt(20,25,35), 7);
 695
 696   Standard_Real deflection = 0.005;
 697   Standard_Integer nbThreads = 1;
 698   Standard_Integer nbx = 200, nby = 200, nbz = 200;
 699   Voxel_BoolDS theVoxels(0,0,0, 50, 50, 50, nbx, nby, nbz);
 700   Voxel_BoolDS theVoxels1(0,0,0, 50, 50, 50, nbx, nby, nbz);
 701
 702   Standard_Integer progress = 0;
 703   Voxel_FastConverter fcp(aShape, theVoxels, deflection, nbx, nby, nbz, nbThreads);
 704   fcp.ConvertUsingSAT(progress, 1);
 705   fcp.FillInVolume(1);
 706
 707   Voxel_FastConverter fcp1(aShape1, theVoxels1, deflection, nbx, nby, nbz, nbThreads);
 708   fcp1.ConvertUsingSAT(progress, 1);
 709   fcp1.FillInVolume(1);
 710
 711   Voxel_BooleanOperation op;
 712   Standard_Boolean result = op.Cut(theVoxels1, theVoxels);
 713   if ( result != 1 ) {
 714     di << "Error: invalid boolean operation" << "\n";
 715   } else {
 716         di << "OK: boolean operation is ok" << "\n";
 717   }
 718
 719   return 0;
 720 }
 721
 722 //=======================================================================
 723 //function : OCC11758
 724 //purpose  :
 725 //=======================================================================
 726 static Standard_Integer OCC11758 (Draw_Interpretor& di, Standard_Integer n, const char**)
 727 {
 728   if (n != 1) return 1;
 729
 730   const char* theStr = "0123456789";
 731   Standard_Integer i, j;
 732   for ( i = 0; i < 5; ++i ) {
 733     // TCollection_AsciiString(const Standard_CString astring)
 734     TCollection_AsciiString a(theStr+i);
 735     // IsEqual (const Standard_CString other)const
 736     //assert( a == theStr+i );
 737     QCOMPARE ( a , theStr+i );
 738
 739     //TCollection_AsciiString(const Standard_CString astring,const Standard_Integer aLen )
 740     TCollection_AsciiString b(theStr+i, 3);
 741     //assert( b.Length() == 3 );
 742     //assert( strncmp( b.ToCString(), theStr+i, 3 ) == 0 );
 743     //assert( strlen( b.ToCString() ) == 3 );
 744     QCOMPARE ( b.Length() , 3 );
 745     QCOMPARE ( strncmp( b.ToCString() , theStr+i, 3 ) , 0 );
 746     QCOMPARE ( b.Length() , 3 );
 747
 748     //TCollection_AsciiString(const Standard_Integer aValue)
 749     TCollection_AsciiString c(i);
 750     //assert( c.IsIntegerValue() );
 751     //assert( c.IntegerValue() == i );
 752     QCOMPARE ( c.IsIntegerValue() , Standard_True );
 753     QCOMPARE ( c.IntegerValue() , i );
 754
 755     //TCollection_AsciiString(const Standard_Real aValue)
 756     TCollection_AsciiString d( 0.1*i );
 757     //assert( d.IsRealValue() );
 758     //assert( TCollection_AsciiString(3.3) == "3.3");
 759     QCOMPARE ( d.IsRealValue() , Standard_True );
 760     QCOMPARE ( TCollection_AsciiString(3.3) , "3.3" );
 761
 762     //TCollection_AsciiString(const TCollection_AsciiString& astring)
 763     TCollection_AsciiString e(d);
 764     //assert( e == d );
 765     //assert( e.Length() == d.Length() );
 766     //assert( strcmp( e.ToCString(), d.ToCString() ) == 0 );
 767     QCOMPARE ( e ,d  );
 768     QCOMPARE ( e.Length() , d.Length() );
 769     QCOMPARE ( strcmp( e.ToCString(), d.ToCString() ) , 0 );
 770
 771     // TCollection_AsciiString(const TCollection_AsciiString& astring ,
 772     //                         const Standard_Character other )
 773     TCollection_AsciiString f(e,'\a');
 774     //assert( f.Length() == e.Length() + 1 );
 775     //assert( strncmp( f.ToCString(), e.ToCString(), e.Length() ) == 0 );
 776     //assert( f.Value( f.Length() ) == '\a');
 777     QCOMPARE ( f.Length() , e.Length() + 1 );
 778     QCOMPARE ( strncmp( f.ToCString(), e.ToCString(), e.Length() ) , 0 );
 779     QCOMPARE ( f.Value( f.Length() ) , '\a' );
 780
 781     // TCollection_AsciiString(const TCollection_AsciiString& astring ,
 782     //                         const Standard_CString other )
 783     TCollection_AsciiString g(f, theStr);
 784     //assert( g.Length() == f.Length() + strlen( theStr ));
 785     //assert( strncmp( g.ToCString(), f.ToCString(), f.Length() ) == 0 );
 786     //assert( g.Search( theStr ) == f.Length() + 1 );
 787     QCOMPARE ( g.Length() , f.Length() + (Standard_Integer)strlen( theStr ) );
 788     QCOMPARE ( strncmp( g.ToCString(), f.ToCString(), f.Length() ) , 0 );
 789     QCOMPARE ( g.Search( theStr ) , f.Length() + 1 );
 790
 791     // TCollection_AsciiString(const TCollection_AsciiString& astring ,
 792     //                         const TCollection_AsciiString& other )
 793     TCollection_AsciiString h(d,a);
 794     //assert( h.Length() == d.Length() + a.Length() );
 795     //assert( strncmp( h.ToCString(), d.ToCString(), d.Length() ) == 0 );
 796     //assert( strncmp( h.ToCString() + d.Length(), a.ToCString(), a.Length() ) == 0 );
 797     QCOMPARE ( h.Length() , d.Length() + a.Length() );
 798     QCOMPARE ( strncmp( h.ToCString(), d.ToCString(), d.Length() ) , 0 );
 799     QCOMPARE ( strncmp( h.ToCString() + d.Length(), a.ToCString(), a.Length() ) , 0 );
 800
 801     // AssignCat(const Standard_CString other)
 802     c.AssignCat( a.ToCString() );
 803     //assert( c.Length() == 1 + a.Length() );
 804     //assert( c.Search( a ) == 2 );
 805     QCOMPARE ( c.Length() , 1 + a.Length() );
 806     QCOMPARE ( c.Search( a ) , 2 );
 807
 808     // AssignCat(const TCollection_AsciiString& other)
 809     Standard_Integer dl = d.Length();
 810     d.AssignCat( a );
 811     //assert( d.Length() == dl + a.Length() );
 812     //assert( d.Search( a ) == dl + 1 );
 813     QCOMPARE ( d.Length() , dl + a.Length() );
 814     QCOMPARE ( d.Search( a ) , dl + 1 );
 815
 816     // Capitalize()
 817     TCollection_AsciiString capitalize("aBC");
 818     capitalize.Capitalize();
 819     //assert( capitalize == "Abc" );
 820     QCOMPARE ( capitalize , "Abc" );
 821
 822     // Copy(const Standard_CString fromwhere)
 823     d = theStr+i;
 824     //assert( d == theStr+i );
 825     QCOMPARE ( d , theStr+i );
 826
 827     // Copy(const TCollection_AsciiString& fromwhere)
 828     d = h;
 829     // IsEqual (const TCollection_AsciiString& other)const
 830     //assert( d == h );
 831     QCOMPARE ( d , h );
 832
 833     // Insert(const Standard_Integer where, const Standard_CString what)
 834     dl = d.Length();
 835     d.Insert( 2, theStr );
 836     //assert( d.Length() == dl + strlen( theStr ));
 837     //assert( strncmp( d.ToCString() + 1, theStr, strlen( theStr )) == 0 );
 838     QCOMPARE ( d.Length() , dl + (Standard_Integer)strlen( theStr ) );
 839     QCOMPARE ( strncmp( d.ToCString() + 1, theStr, strlen( theStr )) , 0 );
 840
 841     //Insert(const Standard_Integer where,const Standard_Character what)
 842     d = theStr;
 843     d.Insert( i+1, 'i' );
 844     //assert( d.Length() == strlen( theStr ) + 1 );
 845     //assert( d.Value( i+1 ) == 'i');
 846     //assert( strcmp( d.ToCString() + i + 1, theStr+i ) == 0 );
 847     QCOMPARE ( d.Length() , (Standard_Integer)strlen( theStr ) + 1 );
 848     QCOMPARE ( d.Value( i+1 ) , 'i' );
 849     QCOMPARE ( strcmp( d.ToCString() + i + 1, theStr+i ) , 0 );
 850
 851     //Insert(const Standard_Integer where,const TCollection_AsciiString& what)
 852     d = theStr;
 853     d.Insert( i+1, TCollection_AsciiString( "i" ));
 854     //assert( d.Length() == strlen( theStr ) + 1 );
 855     //assert( d.Value( i+1 ) == 'i');
 856     //assert( strcmp( d.ToCString() + i + 1, theStr+i ) == 0 );
 857     QCOMPARE ( d.Length() , (Standard_Integer)strlen( theStr ) + 1 );
 858     QCOMPARE ( d.Value( i+1 ) , 'i' );
 859     QCOMPARE ( strcmp( d.ToCString() + i + 1, theStr+i ) , 0 );
 860
 861     // IsDifferent (const Standard_CString other)const
 862     //assert( d.IsDifferent( theStr ));
 863     //assert( d.IsDifferent( "theStr" ));
 864     //assert( d.IsDifferent( "" ));
 865     //assert( !d.IsDifferent( d.ToCString() ));
 866     QCOMPARE ( d.IsDifferent( theStr ) , Standard_True );
 867     QCOMPARE ( d.IsDifferent( "theStr" ) , Standard_True );
 868     QCOMPARE ( d.IsDifferent( "" ) , Standard_True );
 869     QCOMPARE ( !d.IsDifferent( d.ToCString() ) , Standard_True );
 870
 871     // IsDifferent (const TCollection_AsciiString& other)const
 872     //assert( d.IsDifferent( TCollection_AsciiString() ));
 873     //assert( d.IsDifferent( a ));
 874     //assert( d.IsDifferent( h ));
 875     //assert( !d.IsDifferent( d ));
 876     QCOMPARE ( d.IsDifferent( TCollection_AsciiString() ) , Standard_True );
 877     QCOMPARE ( d.IsDifferent( a ) , Standard_True );
 878     QCOMPARE ( d.IsDifferent( h ) , Standard_True );
 879     QCOMPARE ( !d.IsDifferent( d ) , Standard_True );
 880
 881     // IsLess (const Standard_CString other)const
 882     //assert( TCollection_AsciiString ("0"). IsLess("1"));
 883     //assert( TCollection_AsciiString ("0"). IsLess("00"));
 884     //assert( TCollection_AsciiString ("").  IsLess("0"));
 885     //assert( !TCollection_AsciiString("1"). IsLess("0"));
 886     //assert( !TCollection_AsciiString("00").IsLess("0"));
 887     //assert( !TCollection_AsciiString("0"). IsLess(""));
 888     //assert( TCollection_AsciiString (theStr+i).IsLess(theStr+i+1));
 889     QCOMPARE ( TCollection_AsciiString ("0"). IsLess("1") , Standard_True );
 890     QCOMPARE ( TCollection_AsciiString ("0"). IsLess("00") , Standard_True );
 891     QCOMPARE ( TCollection_AsciiString ("").  IsLess("0") , Standard_True );
 892     QCOMPARE ( !TCollection_AsciiString("1"). IsLess("0"), Standard_True );
 893     QCOMPARE ( !TCollection_AsciiString("00").IsLess("0") , Standard_True );
 894     QCOMPARE ( !TCollection_AsciiString("0"). IsLess("") , Standard_True );
 895     QCOMPARE ( TCollection_AsciiString (theStr+i).IsLess(theStr+i+1) , Standard_True );
 896
 897     // IsLess (const TCollection_AsciiString& other)const
 898     //assert( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("1" )));
 899     //assert( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("00")));
 900     //assert( TCollection_AsciiString ("").  IsLess(TCollection_AsciiString("0" )));
 901     //assert( !TCollection_AsciiString("1"). IsLess(TCollection_AsciiString("0" )));
 902     //assert( !TCollection_AsciiString("00").IsLess(TCollection_AsciiString("0" )));
 903     //assert( !TCollection_AsciiString("0"). IsLess(TCollection_AsciiString(""  )));
 904     //assert( TCollection_AsciiString (theStr+i).IsLess(TCollection_AsciiString(theStr+i+1)));
 905     QCOMPARE ( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("1" )) , Standard_True );
 906     QCOMPARE ( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("00")) , Standard_True );
 907     QCOMPARE ( TCollection_AsciiString ("").  IsLess(TCollection_AsciiString("0" )) , Standard_True );
 908     QCOMPARE ( !TCollection_AsciiString("1"). IsLess(TCollection_AsciiString("0" )) , Standard_True );
 909     QCOMPARE ( !TCollection_AsciiString("00").IsLess(TCollection_AsciiString("0" )) , Standard_True );
 910     QCOMPARE ( !TCollection_AsciiString("0"). IsLess(TCollection_AsciiString(""  )) , Standard_True );
 911     QCOMPARE ( TCollection_AsciiString (theStr+i).IsLess(TCollection_AsciiString(theStr+i+1)) , Standard_True );
 912
 913     // IsGreater (const Standard_CString other)const
 914     //assert( !TCollection_AsciiString("0"). IsGreater("1"));
 915     //assert( !TCollection_AsciiString("0"). IsGreater("00"));
 916     //assert( !TCollection_AsciiString("").  IsGreater("0"));
 917     //assert( TCollection_AsciiString ("1"). IsGreater("0"));
 918     //assert( TCollection_AsciiString ("00").IsGreater("0"));
 919     //assert( TCollection_AsciiString ("0"). IsGreater(""));
 920     //assert( TCollection_AsciiString (theStr+i+1).IsGreater(theStr+i));
 921     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater("1") , Standard_True );
 922     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater("00") , Standard_True );
 923     QCOMPARE ( !TCollection_AsciiString("").  IsGreater("0") , Standard_True );
 924     QCOMPARE ( TCollection_AsciiString ("1"). IsGreater("0") , Standard_True );
 925     QCOMPARE ( TCollection_AsciiString ("00").IsGreater("0") , Standard_True );
 926     QCOMPARE ( TCollection_AsciiString ("0"). IsGreater("") , Standard_True );
 927     QCOMPARE ( TCollection_AsciiString (theStr+i+1).IsGreater(theStr+i) , Standard_True );
 928
 929     // IsGreater (const TCollection_AsciiString& other)const
 930     //assert( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("1" )));
 931     //assert( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("00")));
 932     //assert( !TCollection_AsciiString("").  IsGreater(TCollection_AsciiString("0" )));
 933     //assert( TCollection_AsciiString ("1"). IsGreater(TCollection_AsciiString("0" )));
 934     //assert( TCollection_AsciiString ("00").IsGreater(TCollection_AsciiString("0" )));
 935     //assert( TCollection_AsciiString ("0"). IsGreater(TCollection_AsciiString(""  )));
 936     //assert( TCollection_AsciiString (theStr+i+1).IsGreater(TCollection_AsciiString(theStr+i)));
 937     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("1" )) , Standard_True );
 938     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("00")) , Standard_True );
 939     QCOMPARE ( !TCollection_AsciiString("").  IsGreater(TCollection_AsciiString("0" )) , Standard_True );
 940     QCOMPARE ( TCollection_AsciiString ("1"). IsGreater(TCollection_AsciiString("0" )) , Standard_True );
 941     QCOMPARE ( TCollection_AsciiString ("00").IsGreater(TCollection_AsciiString("0" )) , Standard_True );
 942     QCOMPARE ( TCollection_AsciiString ("0"). IsGreater(TCollection_AsciiString(""  )) , Standard_True );
 943     QCOMPARE ( TCollection_AsciiString (theStr+i+1).IsGreater(TCollection_AsciiString(theStr+i)) , Standard_True );
 944
 945     // void Read(Standard_IStream& astream)
 946     std::istringstream is( theStr );
 947     e.Read( is );
 948     //assert( e == theStr );
 949     QCOMPARE ( e , theStr );
 950
 951     // Standard_Integer SearchFromEnd (const Standard_CString what)const
 952     //assert( e.SearchFromEnd( theStr + i ) == i + 1 );
 953     QCOMPARE ( e.SearchFromEnd( theStr + i ) , i + 1 );
 954
 955     // SetValue(const Standard_Integer where, const Standard_CString what)
 956     e.SetValue( i+1, "what");
 957     //assert( e.Search( "what" ) == i+1 );
 958     //assert( e.Length() == strlen( theStr ));
 959     QCOMPARE ( e.Search( "what" ) , i+1 );
 960     QCOMPARE ( e.Length() , (Standard_Integer)strlen( theStr ) );
 961
 962     // TCollection_AsciiString Split (const Standard_Integer where)
 963     e = theStr;
 964     d = e.Split( i+1 );
 965     //assert( d.Length() + e.Length() == strlen( theStr ));
 966     QCOMPARE ( d.Length() + e.Length() , (Standard_Integer)strlen( theStr ) );
 967
 968     // TCollection_AsciiString SubString (const Standard_Integer FromIndex,
 969     //                                    const Standard_Integer ToIndex) const
 970     e = theStr;
 971     d = e.SubString( (unsigned int)i+1, (unsigned int)i+3 );
 972     //assert( d.Length() == 3 );
 973     //assert( d.Value(1) == theStr[ i ]);
 974     QCOMPARE ( d.Length() , 3 );
 975     QCOMPARE ( d.Value(1) , theStr[ i ] );
 976
 977     // TCollection_AsciiString Token (const Standard_CString separators,
 978     //                                const Standard_Integer whichone) const
 979     e = " ";
 980     for ( j = 0; j < i; ++j ) {
 981       e += TCollection_AsciiString( theStr[j] ) + " ";
 982       //assert( e.Token(" ", j+1 ) == TCollection_AsciiString( theStr+j, 1 ));
 983       QCOMPARE ( e.Token(" ", j+1 ) , TCollection_AsciiString( theStr+j, 1 ) );
 984     }
 985   }
 986   for ( i = 0; i < 5; ++i )
 987   {
 988     // TCollection_ExtendedString (const Standard_CString astring,
 989     //                             const Standard_Boolean isMultiByte)
 990     const TCollection_ExtendedString a( theStr+i );
 991     //assert( TCollection_AsciiString( a ) == theStr+i );
 992     QCOMPARE ( TCollection_AsciiString( a ) , theStr+i );
 993
 994     //TCollection_ExtendedString (const Standard_ExtString astring)
 995     const TCollection_ExtendedString b( a.ToExtString() );
 996     //assert( a == b );
 997     QCOMPARE ( a , b );
 998
 999     // TCollection_ExtendedString (const Standard_Integer      length,
1000     //                             const Standard_ExtCharacter filler )
1001     const TCollection_ExtendedString c( i, 1 );
1002     //assert( c.Length() == i );
1003     QCOMPARE ( c.Length() , i );
1004     if ( c.Length() > 0 ) {
1005       //assert( c.Value( i ) == 1 );
1006       QCOMPARE ( c.Value( i ) , 1 );
1007     }
1008
1009     // TCollection_ExtendedString (const Standard_Integer aValue)
1010     TCollection_ExtendedString d( i );
1011     const TCollection_AsciiString da( d );
1012     //assert( da.IsIntegerValue() );
1013     //assert( da.IntegerValue() == i );
1014     QCOMPARE ( da.IsIntegerValue() , Standard_True );
1015     QCOMPARE (  da.IntegerValue(), i );
1016
1017     // TCollection_ExtendedString (const Standard_Real aValue)
1018     const TCollection_ExtendedString e( 0.1 * i );
1019     const TCollection_AsciiString ea( e );
1020     //assert( ea.IsRealValue() );
1021     //assert( Abs( ea.RealValue() - 0.1 * i ) < 1e-10 );
1022     QCOMPARE ( ea.IsRealValue() , Standard_True );
1023     QCOMPARE ( Abs( ea.RealValue() - 0.1 * i ) < 1e-10 , Standard_True );
1024
1025     // TCollection_ExtendedString (const TCollection_ExtendedString& astring)
1026     const TCollection_ExtendedString f(e);
1027     //assert( f.Length() == e.Length());
1028     //assert( f == e );
1029     QCOMPARE ( f.Length() , e.Length() );
1030     QCOMPARE ( f , e );
1031
1032     // TCollection_ExtendedString (const TCollection_AsciiString& astring)
1033     const TCollection_ExtendedString g( ea );
1034     //assert( g.Length() == ea.Length() );
1035     //assert( TCollection_AsciiString( g ) == ea );
1036     QCOMPARE ( g.Length() , ea.Length() );
1037     QCOMPARE ( TCollection_AsciiString( g ) , ea );
1038
1039     // AssignCat (const TCollection_ExtendedString& other)
1040     const TCollection_ExtendedString sep(",");
1041     d.AssignCat( sep );
1042     d.AssignCat( g );
1043     //assert( d.Length() == 2 + g.Length() );
1044     //assert( d.Token( sep.ToExtString(), 1 ) == TCollection_ExtendedString( i ));
1045     //assert( d.Token( sep.ToExtString(), 2 ) == g );
1046     QCOMPARE ( d.Length() , 2 + g.Length() );
1047     QCOMPARE ( d.Token( sep.ToExtString(), 1 ) , TCollection_ExtendedString( i ) );
1048     QCOMPARE ( d.Token( sep.ToExtString(), 2 ) , g );
1049
1050     // TCollection_ExtendedString Cat (const TCollection_ExtendedString& other) const
1051     const TCollection_ExtendedString cat = a.Cat( sep );
1052     //assert( cat.Length() == a.Length() + sep.Length() );
1053     //assert( cat.Search( a ) == 1 );
1054     //assert( cat.Search( sep ) == a.Length() + 1 );
1055     QCOMPARE ( cat.Length() , a.Length() + sep.Length() );
1056     QCOMPARE ( cat.Search( a ) , 1 );
1057     QCOMPARE ( cat.Search( sep ) , a.Length() + 1 );
1058
1059     // Copy (const TCollection_ExtendedString& fromwhere)
1060     d = cat;
1061     //assert( d.Length() == cat.Length() );
1062     //assert( d == cat );
1063     QCOMPARE ( d.Length() , cat.Length() );
1064     QCOMPARE ( d , cat );
1065
1066     // IsEqual (const Standard_ExtString other) const
1067     //assert( d.IsEqual( d.ToExtString() ));
1068     QCOMPARE ( d.IsEqual( d.ToExtString() ) , Standard_True );
1069
1070     // IsDifferent (const Standard_ExtString other ) const
1071     //assert( d.IsDifferent( a.ToExtString() ));
1072     QCOMPARE ( d.IsDifferent( a.ToExtString() ) , Standard_True );
1073
1074     // IsDifferent (const TCollection_ExtendedString& other) const
1075     //assert( d.IsDifferent( a ));
1076     QCOMPARE ( d.IsDifferent( a ) , Standard_True );
1077
1078     // IsLess (const Standard_ExtString other) const
1079     const TCollection_ExtendedString l0("0"), l1("1"), l00("00"), l, ls(theStr+i), ls1(theStr+i+1);
1080     //assert( l0. IsLess( l1.ToExtString() ));
1081     //assert( l0. IsLess( l00.ToExtString() ));
1082     //assert( l.  IsLess( l0.ToExtString() ));
1083     //assert( ! l1. IsLess( l0.ToExtString() ));
1084     //assert( ! l00.IsLess( l0.ToExtString() ));
1085     //assert( ! l0. IsLess( l.ToExtString() ));
1086     //assert( ls.IsLess( ls1.ToExtString() ));
1087     QCOMPARE ( l0. IsLess( l1.ToExtString() ) , Standard_True );
1088     QCOMPARE ( l0. IsLess( l00.ToExtString() ) , Standard_True );
1089     QCOMPARE ( l.  IsLess( l0.ToExtString() ) , Standard_True );
1090     QCOMPARE ( ! l1. IsLess( l0.ToExtString() ) , Standard_True );
1091     QCOMPARE ( ! l00.IsLess( l0.ToExtString() ) , Standard_True );
1092     QCOMPARE ( ! l0. IsLess( l.ToExtString() ) , Standard_True );
1093     QCOMPARE ( ls.IsLess( ls1.ToExtString() ) , Standard_True );
1094
1095     // IsLess (const TCollection_ExtendedString& other) const
1096     //assert( l0. IsLess( l1 ));
1097     //assert( l0. IsLess( l00 ));
1098     //assert( l.  IsLess( l0 ));
1099     //assert( ! l1. IsLess( l0 ));
1100     //assert( ! l00.IsLess( l0 ));
1101     //assert( ! l0. IsLess( l ));
1102     //assert( ls.IsLess( ls1 ));
1103     QCOMPARE ( l0. IsLess( l1 ) , Standard_True );
1104     QCOMPARE ( l0. IsLess( l00 ) , Standard_True );
1105     QCOMPARE ( l.  IsLess( l0 ) , Standard_True );
1106     QCOMPARE ( ! l1. IsLess( l0 ) , Standard_True );
1107     QCOMPARE ( ! l00.IsLess( l0 ) , Standard_True );
1108     QCOMPARE ( ! l0. IsLess( l ) , Standard_True );
1109     QCOMPARE ( ls.IsLess( ls1 ) , Standard_True );
1110
1111     // IsGreater (const Standard_ExtString other) const
1112     //assert( ! l0.IsGreater( l1.ToExtString() ));
1113     //assert( ! l0.IsGreater( l00.ToExtString() ));
1114     //assert( ! l. IsGreater( l0.ToExtString() ));
1115     //assert(  l1. IsGreater( l0.ToExtString() ));
1116     //assert(  l00.IsGreater( l0.ToExtString() ));
1117     //assert(  l0. IsGreater( l.ToExtString() ));
1118     //assert(  ls1.IsGreater( ls.ToExtString() ));
1119     QCOMPARE ( ! l0.IsGreater( l1.ToExtString() ) , Standard_True );
1120     QCOMPARE ( ! l0.IsGreater( l00.ToExtString() ) , Standard_True );
1121     QCOMPARE ( ! l. IsGreater( l0.ToExtString() ) , Standard_True );
1122     QCOMPARE ( l1. IsGreater( l0.ToExtString() ) , Standard_True );
1123     QCOMPARE ( l00.IsGreater( l0.ToExtString() ) , Standard_True );
1124     QCOMPARE ( l0. IsGreater( l.ToExtString() ) , Standard_True );
1125     QCOMPARE ( ls1.IsGreater( ls.ToExtString() ) ,Standard_True  );
1126
1127     // IsGreater (const TCollection_ExtendedString& other) const
1128     //assert( ! l0.IsGreater( l1));
1129     //assert( ! l0.IsGreater( l00));
1130     //assert( ! l. IsGreater( l0));
1131     //assert(  l1. IsGreater( l0));
1132     //assert(  l00.IsGreater( l0));
1133     //assert(  l0. IsGreater( l));
1134     //assert(  ls1.IsGreater( ls));
1135     QCOMPARE ( ! l0.IsGreater( l1) , Standard_True );
1136     QCOMPARE ( ! l0.IsGreater( l00) , Standard_True );
1137     QCOMPARE ( ! l. IsGreater( l0) , Standard_True );
1138     QCOMPARE ( l1. IsGreater( l0) , Standard_True );
1139     QCOMPARE ( l00.IsGreater( l0) , Standard_True );
1140     QCOMPARE ( l0. IsGreater( l) , Standard_True );
1141     QCOMPARE ( ls1.IsGreater( ls) , Standard_True );
1142
1143     // ==========================
1144     //TCollection_HAsciiString::
1145     // ==========================
1146
1147     // IsDifferent(const Handle(TCollection_HAsciiString)& S)
1148     Handle(TCollection_HAsciiString) ha1 = new TCollection_HAsciiString( theStr+i );
1149     Handle(TCollection_HAsciiString) ha2 = new TCollection_HAsciiString( theStr+i+1 );
1150     //assert( ha1->IsDifferent( ha2 ));
1151     //assert( !ha1->IsDifferent( ha1 ));
1152     QCOMPARE ( ha1->IsDifferent( ha2 ) , Standard_True );
1153     QCOMPARE ( !ha1->IsDifferent( ha1 ) , Standard_True );
1154
1155     // IsSameString (const Handle(TCollection_HAsciiString)& S)
1156     //assert( !ha1->IsSameString( ha2 ));
1157     //assert( ha1->IsSameString( ha1 ));
1158     QCOMPARE ( !ha1->IsSameString( ha2 ) , Standard_True );
1159     QCOMPARE ( ha1->IsSameString( ha1 ) , Standard_True );
1160
1161     // IsSameState (const Handle(TCollection_HAsciiString)& other) const
1162     //assert( !ha1->IsSameState( ha2 ));
1163     //assert( ha1->IsSameState( ha1 ));
1164     QCOMPARE ( !ha1->IsSameState( ha2 ) , Standard_True );
1165     QCOMPARE ( ha1->IsSameState( ha1 ) , Standard_True );
1166
1167     // IsSameString (const Handle(TCollection_HAsciiString)& S ,
1168     //               const Standard_Boolean CaseSensitive) const
1169     //assert( !ha1->IsSameString( ha2, true ));
1170     //assert( ha1->IsSameString( ha1, true ));
1171     //assert( !ha1->IsSameString( ha2, false ));
1172     //assert( ha1->IsSameString( ha1, false ));
1173     QCOMPARE ( !ha1->IsSameString( ha2, Standard_True ) , Standard_True );
1174     QCOMPARE ( ha1->IsSameString( ha1, Standard_True ) , Standard_True );
1175     QCOMPARE ( !ha1->IsSameString( ha2, Standard_False ) , Standard_True );
1176     QCOMPARE ( ha1->IsSameString( ha1, Standard_False ) , Standard_True );
1177
1178     ha1->SetValue( 1, "AbC0000000");
1179     ha2->SetValue( 1, "aBc0000000");
1180     //assert( !ha1->IsSameString( ha2, true ));
1181     //assert( ha1->IsSameString( ha2, false ));
1182     QCOMPARE ( !ha1->IsSameString( ha2, Standard_True ) , Standard_True );
1183     QCOMPARE (  ha1->IsSameString( ha2, Standard_False ), Standard_True );
1184   }
1185   return 0;
1186 }
1187
1188 #include <Geom_CylindricalSurface.hxx>
1189 #include <IntTools_FaceFace.hxx>
1190 #include <IntTools_Curve.hxx>
1191 #include <IntTools_PntOn2Faces.hxx>
1192
1193 static Standard_Integer OCC24005 (Draw_Interpretor& theDI, Standard_Integer theNArg, const char** theArgv)
1194 {
1195   if(theNArg < 2)
1196   {
1197     theDI << "Wrong a number of arguments!\n";
1198     return 1;
1199   }
1200
1201   Handle_Geom_Plane plane(new Geom_Plane(
1202                                   gp_Ax3( gp_Pnt(-72.948737453424499, 754.30437716359393, 259.52151854671678),
1203                                   gp_Dir(6.2471473085930200e-007, -0.99999999999980493, 0.00000000000000000),
1204                                   gp_Dir(0.99999999999980493, 6.2471473085930200e-007, 0.00000000000000000))));
1205   Handle(Geom_CylindricalSurface) cylinder(
1206                   new Geom_CylindricalSurface(
1207                                   gp_Ax3(gp_Pnt(-6.4812490053250649, 753.39408794522092, 279.16400974257465),
1208                                   gp_Dir(1.0000000000000000, 0.0, 0.00000000000000000),
1209                                   gp_Dir(0.0, 1.0000000000000000, 0.00000000000000000)),
1210                                                                                           19.712534607908712));
1211
1212   DrawTrSurf::Set("pln", plane);
1213   theDI << "pln\n";
1214   DrawTrSurf::Set("cyl", cylinder);
1215   theDI << "cyl\n";
1216
1217   BRep_Builder builder;
1218   TopoDS_Face face1, face2;
1219   builder.MakeFace(face1, plane, Precision::Confusion());
1220   builder.MakeFace(face2, cylinder, Precision::Confusion());
1221   IntTools_FaceFace anInters;
1222   anInters.SetParameters(false, true, true, Precision::Confusion());
1223   anInters.Perform(face1, face2);
1224
1225   if (!anInters.IsDone())
1226   {
1227     theDI<<"No intersections found!"<<"\n";
1228
1229     return 1;
1230   }
1231
1232   //Handle(Geom_Curve) aResult;
1233   //gp_Pnt             aPoint;
1234
1235   const IntTools_SequenceOfCurves& aCvsX=anInters.Lines();
1236   const IntTools_SequenceOfPntOn2Faces& aPntsX=anInters.Points();
1237
1238   char buf[1024];
1239   Standard_Integer aNbCurves, aNbPoints;
1240
1241   aNbCurves=aCvsX.Length();
1242   aNbPoints=aPntsX.Length();
1243
1244   if (aNbCurves >= 2)
1245   {
1246     for (Standard_Integer i=1; i<=aNbCurves; ++i)
1247     {
1248       Sprintf(buf, "%s_%d",theArgv[1],i);
1249       theDI << buf << " ";
1250
1251       const IntTools_Curve& aIC = aCvsX(i);
1252       const Handle(Geom_Curve)& aC3D= aIC.Curve();
1253       DrawTrSurf::Set(buf,aC3D);
1254     }
1255   }
1256   else if (aNbCurves == 1)
1257   {
1258     const IntTools_Curve& aIC = aCvsX(1);
1259     const Handle(Geom_Curve)& aC3D= aIC.Curve();
1260     Sprintf(buf, "%s",theArgv[1]);
1261     theDI << buf << " ";
1262     DrawTrSurf::Set(buf,aC3D);
1263   }
1264
1265   for (Standard_Integer i = 1; i<=aNbPoints; ++i)
1266   {
1267     const IntTools_PntOn2Faces& aPi=aPntsX(i);
1268     const gp_Pnt& aP=aPi.P1().Pnt();
1269
1270     Sprintf(buf,"%s_p_%d",theArgv[1],i);
1271     theDI << buf << " ";
1272     DrawTrSurf::Set(buf, aP);
1273   }
1274
1275   return 0;
1276 }
1277
1278 #include <BRepAlgo_NormalProjection.hxx>
1279 static Standard_Integer OCC24012 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1280 {
1281         if (argc != 3) {
1282                 di << "Usage : " << argv[0] << " should be 2 arguments (face and edge)";
1283                 return 1;
1284         }
1285
1286         Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
1287         if(myAISContext.IsNull()) {
1288                 di << "use 'vinit' command before " << argv[0] << "\n";
1289                 return 1;
1290         }
1291
1292         TopoDS_Face m_Face1 = TopoDS::Face(DBRep::Get(argv[1]));
1293         TopoDS_Edge m_Edge = TopoDS::Edge(DBRep::Get(argv[2]));
1294
1295         BRepAlgo_NormalProjection anormpro(m_Face1);
1296     anormpro.Add(m_Edge);
1297     anormpro.SetDefaultParams();
1298
1299     //anormpro.Compute3d();
1300     //anormpro.SetLimit();
1301
1302     anormpro.Build();
1303
1304     if (anormpro.IsDone())
1305     {
1306         TopoDS_Shape rshape = anormpro.Projection();
1307                 Handle(AIS_InteractiveObject) myShape = new AIS_Shape (rshape);
1308                 myAISContext->SetColor(myShape, Quantity_Color(Quantity_NOC_YELLOW));
1309                 myAISContext->Display(myShape, Standard_True);
1310     }
1311
1312         return 0;
1313 }
1314
1315 #include <Voxel_FastConverter.hxx>
1316 static Standard_Integer OCC24051 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1317 {
1318         if (argc != 1) {
1319                 di << "Usage : " << argv[0] << " should be one argument (command name only)";
1320                 return 1;
1321         }
1322
1323         TopoDS_Shape shape = BRepPrimAPI_MakeBox(gp_Pnt(5, 10, 10), 10, 20, 30).Shape();
1324         Standard_Integer progress = 0;
1325         Standard_Real deflection = 0.005;
1326         Standard_Integer nbx = 200, nby = 200, nbz = 200;
1327         Voxel_BoolDS theVoxels(-50,-50,-30, 100, 100, 100, nbx, nby, nbz);
1328         Voxel_BoolDS theVoxels1(-50,-50,-30, 100, 100, 100, nbx, nby, nbz);
1329         Standard_Integer nbThreads = 5;
1330         Voxel_FastConverter fcp(shape, theVoxels, deflection, nbx, nby, nbz, nbThreads, Standard_True);
1331
1332         #ifdef WNT
1333         #pragma omp parallel for
1334         for(int i = 0; i < nbThreads; i++)
1335                         fcp.ConvertUsingSAT(progress, i+1);
1336         #endif
1337
1338         fcp.ConvertUsingSAT(progress);
1339
1340         return 0;
1341 }
1342
1343 #include <BRepFeat_SplitShape.hxx>
1344 #include <ShapeAnalysis_ShapeContents.hxx>
1345 #include <BRepAlgo.hxx>
1346 static Standard_Integer OCC24086 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1347 {
1348         if (argc != 3) {
1349                 di << "Usage : " << argv[0] << " should be 2 arguments (face and wire)";
1350                 return 1;
1351         }
1352
1353         Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
1354         if(myAISContext.IsNull()) {
1355                 di << "use 'vinit' command before " << argv[0] << "\n";
1356                 return 1;
1357         }
1358
1359         TopoDS_Shape result;
1360         TopoDS_Face face = TopoDS::Face(DBRep::Get(argv[1]));
1361         TopoDS_Wire wire = TopoDS::Wire(DBRep::Get(argv[2]));
1362
1363         BRepFeat_SplitShape asplit(face);
1364         asplit.Add(wire, face);
1365         asplit.Build();
1366     result = asplit.Shape();
1367     ShapeAnalysis_ShapeContents ana;
1368     ana.Perform(result);
1369     ana.NbFaces();
1370
1371         if (!(BRepAlgo::IsValid(result))) {
1372                 di << "Result was checked and it is INVALID" << "\n";
1373         } else {
1374                 di << "Result was checked and it is VALID" << "\n";
1375         }
1376
1377         Handle(AIS_InteractiveObject) myShape = new AIS_Shape (result);
1378         myAISContext->Display(myShape, Standard_True);
1379
1380         return 0;
1381 }
1382
1383 #include <Extrema_FuncExtPS.hxx>
1384 #include <math_FunctionSetRoot.hxx>
1385 #include <math_Vector.hxx>
1386 #include <BRepBuilderAPI_MakeVertex.hxx>
1387 static Standard_Integer OCC24137 (Draw_Interpretor& theDI, Standard_Integer theNArg, const char** theArgv)
1388 {
1389   Standard_Integer anArgIter = 1;
1390   if (theNArg < 5)
1391     {
1392       theDI <<"Usage: " << theArgv[0] << " face vertex U V [N]"<<"\n";
1393       return 1;
1394     }
1395
1396   // get target shape
1397   Standard_CString aFaceName = theArgv[anArgIter++];
1398   Standard_CString aVertName = theArgv[anArgIter++];
1399   const TopoDS_Shape aShapeF = DBRep::Get (aFaceName);
1400   const TopoDS_Shape aShapeV = DBRep::Get (aVertName);
1401   const Standard_Real aUFrom = Atof (theArgv[anArgIter++]);
1402   const Standard_Real aVFrom = Atof (theArgv[anArgIter++]);
1403   const Standard_Integer aNbIts = (anArgIter < theNArg) ? atol (theArgv[anArgIter++]) : 100;
1404   if (aShapeF.IsNull() || aShapeF.ShapeType() != TopAbs_FACE)
1405     {
1406       std::cout << "Error: " << aFaceName << " shape is null / not a face" << std::endl;
1407       return 1;
1408     }
1409   if (aShapeV.IsNull() || aShapeV.ShapeType() != TopAbs_VERTEX)
1410     {
1411       std::cout << "Error: " << aVertName << " shape is null / not a vertex" << std::endl;
1412       return 1;
1413     }
1414   const TopoDS_Face   aFace = TopoDS::Face   (aShapeF);
1415   const TopoDS_Vertex aVert = TopoDS::Vertex (aShapeV);
1416   GeomAdaptor_Surface aSurf (BRep_Tool::Surface (aFace));
1417
1418   gp_Pnt aPnt = BRep_Tool::Pnt (aVert), aRes;
1419
1420   Extrema_FuncExtPS    anExtFunc;
1421   math_FunctionSetRoot aRoot (anExtFunc, aNbIts);
1422
1423   math_Vector aTolUV (1, 2), aUVinf  (1, 2), aUVsup  (1, 2), aFromUV (1, 2);
1424   aTolUV (1) =  Precision::Confusion(); aTolUV (2) =  Precision::Confusion();
1425   aUVinf (1) = -Precision::Infinite();  aUVinf (2) = -Precision::Infinite();
1426   aUVsup (1) =  Precision::Infinite();  aUVsup (2) =  Precision::Infinite();
1427   aFromUV(1) =  aUFrom; aFromUV(2) = aVFrom;
1428
1429   anExtFunc.Initialize (aSurf);
1430   anExtFunc.SetPoint (aPnt);
1431   aRoot.SetTolerance (aTolUV);
1432   aRoot.Perform (anExtFunc, aFromUV, aUVinf, aUVsup);
1433   if (!aRoot.IsDone())
1434     {
1435       std::cerr << "No results!\n";
1436       return 1;
1437     }
1438
1439   theDI << aRoot.Root()(1) << " " << aRoot.Root()(2) << "\n";
1440
1441   aSurf.D0 (aRoot.Root()(1), aRoot.Root()(2), aRes);
1442   DBRep::Set ("result", BRepBuilderAPI_MakeVertex (aRes));
1443   return 0;
1444 }
1445
1446 //! Check boolean operations on NCollection_Map
1447 static Standard_Integer OCC24271 (Draw_Interpretor& di,
1448                                   Standard_Integer  /*theArgNb*/,
1449                                   const char**      /*theArgVec*/)
1450 {
1451   // input data
1452   const Standard_Integer aLeftLower  = 1;
1453   const Standard_Integer aLeftUpper  = 10;
1454   const Standard_Integer aRightLower = 5;
1455   const Standard_Integer aRightUpper = 15;
1456
1457   // define arguments
1458   NCollection_Map<Standard_Integer> aMapLeft;
1459   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1460   {
1461     aMapLeft.Add (aKeyIter);
1462   }
1463
1464   NCollection_Map<Standard_Integer> aMapRight;
1465   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aRightUpper; ++aKeyIter)
1466   {
1467     aMapRight.Add (aKeyIter);
1468   }
1469
1470   QCOMPARE (aMapLeft .Contains (aMapRight), Standard_False);
1471   QCOMPARE (aMapRight.Contains (aMapLeft),  Standard_False);
1472
1473   // validate Union operation
1474   NCollection_Map<Standard_Integer> aMapUnion;
1475   aMapUnion.Union (aMapLeft, aMapRight);
1476   QCOMPARE (aMapUnion.Extent(), aRightUpper - aLeftLower + 1);
1477   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter <= aRightUpper; ++aKeyIter)
1478   {
1479     QCOMPARE (aMapUnion.Contains (aKeyIter), Standard_True);
1480   }
1481
1482   // validate Intersection operation
1483   NCollection_Map<Standard_Integer> aMapSect;
1484   aMapSect.Intersection (aMapLeft, aMapRight);
1485   QCOMPARE (aMapSect.Extent(), aLeftUpper - aRightLower + 1);
1486   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1487   {
1488     QCOMPARE (aMapSect.Contains (aKeyIter), Standard_True);
1489   }
1490   QCOMPARE (aMapLeft .Contains (aMapSect), Standard_True);
1491   QCOMPARE (aMapRight.Contains (aMapSect), Standard_True);
1492
1493   // validate Substruction operation
1494   NCollection_Map<Standard_Integer> aMapSubsLR;
1495   aMapSubsLR.Subtraction (aMapLeft, aMapRight);
1496   QCOMPARE (aMapSubsLR.Extent(), aRightLower - aLeftLower);
1497   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter < aRightLower; ++aKeyIter)
1498   {
1499     QCOMPARE (aMapSubsLR.Contains (aKeyIter), Standard_True);
1500   }
1501
1502   NCollection_Map<Standard_Integer> aMapSubsRL;
1503   aMapSubsRL.Subtraction (aMapRight, aMapLeft);
1504   QCOMPARE (aMapSubsRL.Extent(), aRightUpper - aLeftUpper);
1505   for (Standard_Integer aKeyIter = aLeftUpper + 1; aKeyIter < aRightUpper; ++aKeyIter)
1506   {
1507     QCOMPARE (aMapSubsRL.Contains (aKeyIter), Standard_True);
1508   }
1509
1510   // validate Difference operation
1511   NCollection_Map<Standard_Integer> aMapDiff;
1512   aMapDiff.Difference (aMapLeft, aMapRight);
1513   QCOMPARE (aMapDiff.Extent(), aRightLower - aLeftLower + aRightUpper - aLeftUpper);
1514   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter < aRightLower; ++aKeyIter)
1515   {
1516     QCOMPARE (aMapDiff.Contains (aKeyIter), Standard_True);
1517   }
1518   for (Standard_Integer aKeyIter = aLeftUpper + 1; aKeyIter < aRightUpper; ++aKeyIter)
1519   {
1520     QCOMPARE (aMapDiff.Contains (aKeyIter), Standard_True);
1521   }
1522
1523   // validate Exchange operation
1524   NCollection_Map<Standard_Integer> aMapSwap;
1525   aMapSwap.Exchange (aMapSect);
1526   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1527   {
1528     QCOMPARE (aMapSwap.Contains (aKeyIter), Standard_True);
1529   }
1530   QCOMPARE (aMapSect.IsEmpty(), Standard_True);
1531   aMapSwap.Add (34);
1532   aMapSect.Add (43);
1533
1534   NCollection_Map<Standard_Integer> aMapCopy (aMapSwap);
1535   QCOMPARE (aMapCopy.IsEqual (aMapSwap), Standard_True);
1536   aMapCopy.Remove (34);
1537   aMapCopy.Add    (43);
1538   QCOMPARE (aMapCopy.IsEqual (aMapSwap), Standard_False);
1539
1540   return 0;
1541 }
1542
1543 #include <ShapeFix_EdgeProjAux.hxx>
1544 static Standard_Integer OCC24370 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
1545 {
1546   if (argc < 5) {
1547     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
1548     return 1;
1549   }
1550
1551   TopoDS_Shape aSh = DBRep::Get(argv[1]);
1552   if (aSh.IsNull()) {
1553     di << argv[0] << " Error: Null input edge\n";
1554     return 1;
1555   }
1556   const TopoDS_Edge& anEdge = TopoDS::Edge (aSh);
1557
1558   Handle(Geom2d_Curve) aC = DrawTrSurf::GetCurve2d(argv[2]);
1559   if (aC.IsNull()) {
1560     di << argv[0] << " Error: Null input curve\n";
1561     return 1;
1562   }
1563
1564   Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(argv[3]);
1565   if (aS.IsNull()) {
1566     di << argv[0] << " Error: Null input surface\n";
1567     return 1;
1568   }
1569
1570   Standard_Real prec = Draw::Atof(argv[4]);
1571
1572   //prepare data
1573   TopoDS_Face aFace;
1574   BRep_Builder aB;
1575   aB.MakeFace (aFace, aS, Precision::Confusion());
1576   aB.UpdateEdge (anEdge, aC, aFace, Precision::Confusion());
1577   aB.Range (anEdge, aFace, aC->FirstParameter(), aC->LastParameter());
1578
1579   //call algorithm
1580   ShapeFix_EdgeProjAux aProj (aFace, anEdge);
1581   aProj.Compute (prec);
1582
1583   Standard_Boolean isfirstdone = aProj.IsFirstDone();
1584   Standard_Boolean islastdone = aProj.IsLastDone();
1585
1586   Standard_Real first = 0.;
1587   Standard_Real last = 0.;
1588   Standard_Integer isfirstdoneInteger = 0;
1589   Standard_Integer islastdoneInteger = 0;
1590
1591
1592   if (isfirstdone) {
1593     first = aProj.FirstParam();
1594     isfirstdoneInteger = 1;
1595   }
1596
1597   if (islastdone) {
1598     last= aProj.LastParam();
1599     islastdoneInteger = 1;
1600   }
1601
1602   di << isfirstdoneInteger << " "<< islastdoneInteger << " "<< first << " "<< last << " \n";
1603
1604   return 0;
1605 }
1606
1607 template<typename T, typename HT>
1608 static void DoIsNull(Draw_Interpretor& di)
1609 {
1610   HT aHandle;
1611   //    QVERIFY (aHandle.IsNull());
1612   QCOMPARE (aHandle.IsNull(), Standard_True);
1613   const T* p = aHandle.Access();
1614 #if OCC_VERSION_HEX > 0x060700
1615   //QVERIFY (!p);
1616   //QVERIFY (p == 0);
1617   QCOMPARE (!p, Standard_True);
1618   QCOMPARE (p == 0, Standard_True);
1619 #endif
1620
1621   aHandle = new T;
1622   //QVERIFY (!aHandle.IsNull());
1623   QCOMPARE (!aHandle.IsNull(), Standard_True);
1624   p = aHandle.Access();
1625   //QVERIFY (p);
1626   //QVERIFY (p != 0);
1627   QCOMPARE (p != NULL, Standard_True);
1628   QCOMPARE (p != 0, Standard_True);
1629 }
1630
1631 //=======================================================================
1632 //function : OCC24533
1633 //purpose  :
1634 //=======================================================================
1635 static Standard_Integer OCC24533 (Draw_Interpretor& di, Standard_Integer n, const char**)
1636 {
1637   if (n != 1) return 1;
1638
1639   DoIsNull<Standard_Transient, Handle_Standard_Transient>(di);
1640   DoIsNull<Standard_Persistent, Handle_Standard_Persistent>(di);
1641
1642   return 0;
1643 }
1644
1645 // Dummy class to test interface for compilation issues
1646 class QABugs_HandleClass : public Standard_Transient
1647 {
1648 public:
1649   Standard_Integer HandleProc (Draw_Interpretor& , Standard_Integer  , const char** theArgVec)
1650   {
1651     std::cerr << "QABugs_HandleClass[" << this << "] " << theArgVec[0] << "\n";
1652     return 0;
1653   }
1654   DEFINE_STANDARD_RTTI(QABugs_HandleClass) // Type definition
1655 };
1656 DEFINE_STANDARD_HANDLE    (QABugs_HandleClass, Standard_Transient)
1657 IMPLEMENT_STANDARD_HANDLE (QABugs_HandleClass, Standard_Transient)
1658 IMPLEMENT_STANDARD_RTTIEXT(QABugs_HandleClass, Standard_Transient)
1659
1660 // Dummy class to test interface for compilation issues
1661 struct QABugs_NHandleClass
1662 {
1663   Standard_Integer NHandleProc (Draw_Interpretor& , Standard_Integer  , const char** theArgVec)
1664   {
1665     std::cerr << "QABugs_NHandleClass[" << this << "] " << "" << theArgVec[0] << "\n";
1666     return 0;
1667   }
1668 };
1669
1670 //=======================================================================
1671 //function : OCC24622
1672 //purpose  : The command tests sourcing Image_PixMap to AIS_TexturedShape
1673 //=======================================================================
1674 static Standard_Integer OCC24622 (Draw_Interpretor& /*theDi*/, Standard_Integer theArgNb, const char** theArgVec)
1675 {
1676   if (theArgNb != 2)
1677   {
1678     std::cout << "Usage : " << theArgVec[0] << " texture={1D|2D}";
1679     return 1;
1680   }
1681
1682   const Handle(AIS_InteractiveContext)& anAISContext = ViewerTest::GetAISContext();
1683   if (anAISContext.IsNull())
1684   {
1685     std::cout << "Please initialize view with \"vinit\".\n";
1686     return 1;
1687   }
1688
1689   Handle(Image_PixMap) anImage = new Image_PixMap();
1690
1691   static const Image_ColorRGB aBitmap[8] =
1692   {
1693     {{255,   0, 0}}, {{0,  148, 255}}, {{ 0, 148, 255}}, {{255,  94, 0}},
1694     {{255, 121, 0}}, {{76, 255,   0}}, {{76, 255,   0}}, {{255, 202, 0}}
1695   };
1696
1697   TCollection_AsciiString aTextureTypeArg (theArgVec[1]);
1698   aTextureTypeArg.UpperCase();
1699   if (aTextureTypeArg == "1D")
1700   {
1701     anImage->InitWrapper (Image_PixMap::ImgRGB, (Standard_Byte*)aBitmap, 8, 1);
1702   }
1703   else if (aTextureTypeArg == "2D")
1704   {
1705     anImage->InitTrash (Image_PixMap::ImgRGB, 8, 8);
1706     Image_PixMapData<Image_ColorRGB>& anImageData = anImage->EditData<Image_ColorRGB>();
1707     for (Standard_Integer aRow = 0; aRow < 8; ++aRow)
1708     {
1709       for (Standard_Integer aCol = 0; aCol < 8; ++aCol)
1710       {
1711         anImageData.ChangeValue (aRow, aCol) = aBitmap[aRow];
1712       }
1713     }
1714   }
1715   else
1716   {
1717     std::cout << "Please specify type of texture to test {1D|2D}.\n";
1718     return 1;
1719   }
1720
1721   TopoDS_Shape aBlankShape = BRepPrimAPI_MakeBox (10.0, 10.0, 10.0).Shape();
1722
1723   Handle(AIS_TexturedShape) aTexturedShape = new AIS_TexturedShape (aBlankShape);
1724   aTexturedShape->SetTexturePixMap (anImage);
1725   anAISContext->Display (aTexturedShape, 3, 0);
1726
1727   return 0;
1728 }
1729
1730 //=======================================================================
1731 //function : OCC24667
1732 //purpose  :
1733 //=======================================================================
1734 static Standard_Integer OCC24667 (Draw_Interpretor& di, Standard_Integer n, const char** a)
1735 {
1736   if (n == 1)
1737   {
1738     di << "OCC24667 result Wire_spine Profile [Mode [Approx]]" << "\n";
1739     di << "Mode = 0 - CorrectedFrenet," << "\n";
1740     di << "     = 1 - Frenet," << "\n";
1741     di << "     = 2 - DiscreteTrihedron" << "\n";
1742     di << "Approx - force C1-approximation if result is C0" << "\n";
1743     return 0;
1744   }
1745
1746   if (n > 1 && n < 4) return 1;
1747
1748   TopoDS_Shape Spine = DBRep::Get(a[2],TopAbs_WIRE);
1749   if ( Spine.IsNull()) return 1;
1750
1751   TopoDS_Shape Profile = DBRep::Get(a[3]);
1752   if ( Profile.IsNull()) return 1;
1753
1754   GeomFill_Trihedron Mode = GeomFill_IsCorrectedFrenet;
1755   if (n >= 5)
1756   {
1757     Standard_Integer iMode = atoi(a[4]);
1758     if (iMode == 1)
1759       Mode = GeomFill_IsFrenet;
1760     else if (iMode == 2)
1761       Mode = GeomFill_IsDiscreteTrihedron;
1762   }
1763
1764   Standard_Boolean ForceApproxC1 = Standard_False;
1765   if (n >= 6)
1766     ForceApproxC1 = Standard_True;
1767
1768   BRepOffsetAPI_MakePipe aPipe(TopoDS::Wire(Spine),
1769                                           Profile,
1770                                           Mode,
1771                                           ForceApproxC1);
1772
1773   TopoDS_Shape S = aPipe.Shape();
1774   TopoDS_Shape aSF = aPipe.FirstShape();
1775   TopoDS_Shape aSL = aPipe.LastShape();
1776
1777   DBRep::Set(a[1],S);
1778
1779   TCollection_AsciiString aStrF(a[1], "_f");
1780   TCollection_AsciiString aStrL(a[1], "_l");
1781
1782   DBRep::Set(aStrF.ToCString(), aSF);
1783   DBRep::Set(aStrL.ToCString(), aSL);
1784
1785   return 0;
1786 }
1787
1788 void QABugs::Commands_19(Draw_Interpretor& theCommands) {
1789   const char *group = "QABugs";
1790
1791   Handle(QABugs_HandleClass) aClassPtr = new QABugs_HandleClass();
1792   theCommands.Add ("OCC24202_1", "Test Handle-based procedure",
1793                    __FILE__, aClassPtr, &QABugs_HandleClass::HandleProc, group);
1794   NCollection_Handle<QABugs_NHandleClass> aNClassPtr = new QABugs_NHandleClass();
1795   theCommands.Add ("OCC24202_2", "Test NCollection_Handle-based procedure",
1796                    __FILE__, aNClassPtr, &QABugs_NHandleClass::NHandleProc, group);
1797
1798   theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
1799   theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
1800   theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
1801   theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group);
1802   theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
1803   theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
1804   theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
1805   theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
1806   theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
1807   theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
1808   theCommands.Add ("OCC23952sweep", "OCC23952sweep nbupoles shape", __FILE__, OCC23952sweep, group);
1809   theCommands.Add ("OCC23952intersect", "OCC23952intersect nbsol shape1 shape2", __FILE__, OCC23952intersect, group);
1810   theCommands.Add ("test_offset", "test_offset", __FILE__, test_offset, group);
1811   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);
1812   theCommands.Add ("OCC24008", "OCC24008 curve surface", __FILE__, OCC24008, group);
1813   theCommands.Add ("OCC24019", "OCC24019 aShape", __FILE__, OCC24019, group);
1814   theCommands.Add ("OCC11758", "OCC11758", __FILE__, OCC11758, group);
1815   theCommands.Add ("OCC24005", "OCC24005 result", __FILE__, OCC24005, group);
1816   theCommands.Add ("OCC24137", "OCC24137 face vertex U V [N]", __FILE__, OCC24137, group);
1817   theCommands.Add ("OCC24271", "Boolean operations on NCollection_Map", __FILE__, OCC24271, group);
1818   theCommands.Add ("OCC24370", "OCC24370 edge pcurve surface prec", __FILE__, OCC24370, group);
1819   theCommands.Add ("OCC24533", "OCC24533", __FILE__, OCC24533, group);
1820   theCommands.Add ("OCC24012", "OCC24012 face edge", __FILE__, OCC24012, group);
1821   theCommands.Add ("OCC24051", "OCC24051", __FILE__, OCC24051, group);
1822   theCommands.Add ("OCC24086", "OCC24086 face wire", __FILE__, OCC24086, group);
1823   theCommands.Add ("OCC24622", "OCC24622 texture={1D|2D}\n Tests sourcing of 1D/2D pixmaps for AIS_TexturedShape", __FILE__, OCC24622, group);
1824   theCommands.Add ("OCC24667", "OCC24667 result Wire_spine Profile [Mode [Approx]], no args to get help", __FILE__, OCC24667, group);
1825   return;
1826 }