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 <Geom_Circle.hxx>
1384 #include <GeomAdaptor_Curve.hxx>
1385 #include <Extrema_ExtPC.hxx>
1386 #include <gp_Cylinder.hxx>
1387 #include <ElSLib.hxx>
1388 static Standard_Integer OCC24945 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1389 {
1390   if (argc != 1) {
1391     di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
1392     return 1;
1393   }
1394
1395   gp_Pnt aP3D( -1725.97, 843.257, -4.22741e-013 );
1396   gp_Ax2 aAxis( gp_Pnt( 0, 843.257, 0 ), gp_Dir( 0, -1, 0 ), gp::DX() );
1397   Handle(Geom_Circle) aCircle = new Geom_Circle( aAxis, 1725.9708621929999 );
1398   GeomAdaptor_Curve aC3D( aCircle );
1399
1400   Extrema_ExtPC aExtPC( aP3D, aC3D );
1401   //Standard_Real aParam = (aExtPC.Point(1)).Parameter();
1402   gp_Pnt aProj = (aExtPC.Point(1)).Value();
1403   di << "Projected point: X = " << aProj.X() << "; Y = " << aProj.Y() << "; Z = " << aProj.Z() << "\n";
1404
1405   // Result of deviation
1406   gp_Ax2 aCylAxis( gp_Pnt( 0, 2103.87, 0 ), -gp::DY(), -gp::DX() );
1407   gp_Cylinder aCylinder( aCylAxis, 1890. );
1408
1409   Standard_Real aU = 0., aV = 0.;
1410   ElSLib::Parameters( aCylinder, aProj, aU, aV );
1411   di << "Parameters on cylinder: U = " << aU << "; V = " << aV << "\n";
1412
1413   return 0;
1414 }
1415
1416 #include <Extrema_FuncExtPS.hxx>
1417 #include <math_FunctionSetRoot.hxx>
1418 #include <math_Vector.hxx>
1419 #include <BRepBuilderAPI_MakeVertex.hxx>
1420 static Standard_Integer OCC24137 (Draw_Interpretor& theDI, Standard_Integer theNArg, const char** theArgv)
1421 {
1422   Standard_Integer anArgIter = 1;
1423   if (theNArg < 5)
1424     {
1425       theDI <<"Usage: " << theArgv[0] << " face vertex U V [N]"<<"\n";
1426       return 1;
1427     }
1428
1429   // get target shape
1430   Standard_CString aFaceName = theArgv[anArgIter++];
1431   Standard_CString aVertName = theArgv[anArgIter++];
1432   const TopoDS_Shape aShapeF = DBRep::Get (aFaceName);
1433   const TopoDS_Shape aShapeV = DBRep::Get (aVertName);
1434   const Standard_Real aUFrom = Atof (theArgv[anArgIter++]);
1435   const Standard_Real aVFrom = Atof (theArgv[anArgIter++]);
1436   const Standard_Integer aNbIts = (anArgIter < theNArg) ? atol (theArgv[anArgIter++]) : 100;
1437   if (aShapeF.IsNull() || aShapeF.ShapeType() != TopAbs_FACE)
1438     {
1439       std::cout << "Error: " << aFaceName << " shape is null / not a face" << std::endl;
1440       return 1;
1441     }
1442   if (aShapeV.IsNull() || aShapeV.ShapeType() != TopAbs_VERTEX)
1443     {
1444       std::cout << "Error: " << aVertName << " shape is null / not a vertex" << std::endl;
1445       return 1;
1446     }
1447   const TopoDS_Face   aFace = TopoDS::Face   (aShapeF);
1448   const TopoDS_Vertex aVert = TopoDS::Vertex (aShapeV);
1449   GeomAdaptor_Surface aSurf (BRep_Tool::Surface (aFace));
1450
1451   gp_Pnt aPnt = BRep_Tool::Pnt (aVert), aRes;
1452
1453   Extrema_FuncExtPS    anExtFunc;
1454   math_FunctionSetRoot aRoot (anExtFunc, aNbIts);
1455
1456   math_Vector aTolUV (1, 2), aUVinf  (1, 2), aUVsup  (1, 2), aFromUV (1, 2);
1457   aTolUV (1) =  Precision::Confusion(); aTolUV (2) =  Precision::Confusion();
1458   aUVinf (1) = -Precision::Infinite();  aUVinf (2) = -Precision::Infinite();
1459   aUVsup (1) =  Precision::Infinite();  aUVsup (2) =  Precision::Infinite();
1460   aFromUV(1) =  aUFrom; aFromUV(2) = aVFrom;
1461
1462   anExtFunc.Initialize (aSurf);
1463   anExtFunc.SetPoint (aPnt);
1464   aRoot.SetTolerance (aTolUV);
1465   aRoot.Perform (anExtFunc, aFromUV, aUVinf, aUVsup);
1466   if (!aRoot.IsDone())
1467     {
1468       std::cerr << "No results!\n";
1469       return 1;
1470     }
1471
1472   theDI << aRoot.Root()(1) << " " << aRoot.Root()(2) << "\n";
1473
1474   aSurf.D0 (aRoot.Root()(1), aRoot.Root()(2), aRes);
1475   DBRep::Set ("result", BRepBuilderAPI_MakeVertex (aRes));
1476   return 0;
1477 }
1478
1479 //! Check boolean operations on NCollection_Map
1480 static Standard_Integer OCC24271 (Draw_Interpretor& di,
1481                                   Standard_Integer  /*theArgNb*/,
1482                                   const char**      /*theArgVec*/)
1483 {
1484   // input data
1485   const Standard_Integer aLeftLower  = 1;
1486   const Standard_Integer aLeftUpper  = 10;
1487   const Standard_Integer aRightLower = 5;
1488   const Standard_Integer aRightUpper = 15;
1489
1490   // define arguments
1491   NCollection_Map<Standard_Integer> aMapLeft;
1492   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1493   {
1494     aMapLeft.Add (aKeyIter);
1495   }
1496
1497   NCollection_Map<Standard_Integer> aMapRight;
1498   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aRightUpper; ++aKeyIter)
1499   {
1500     aMapRight.Add (aKeyIter);
1501   }
1502
1503   QCOMPARE (aMapLeft .Contains (aMapRight), Standard_False);
1504   QCOMPARE (aMapRight.Contains (aMapLeft),  Standard_False);
1505
1506   // validate Union operation
1507   NCollection_Map<Standard_Integer> aMapUnion;
1508   aMapUnion.Union (aMapLeft, aMapRight);
1509   QCOMPARE (aMapUnion.Extent(), aRightUpper - aLeftLower + 1);
1510   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter <= aRightUpper; ++aKeyIter)
1511   {
1512     QCOMPARE (aMapUnion.Contains (aKeyIter), Standard_True);
1513   }
1514
1515   // validate Intersection operation
1516   NCollection_Map<Standard_Integer> aMapSect;
1517   aMapSect.Intersection (aMapLeft, aMapRight);
1518   QCOMPARE (aMapSect.Extent(), aLeftUpper - aRightLower + 1);
1519   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1520   {
1521     QCOMPARE (aMapSect.Contains (aKeyIter), Standard_True);
1522   }
1523   QCOMPARE (aMapLeft .Contains (aMapSect), Standard_True);
1524   QCOMPARE (aMapRight.Contains (aMapSect), Standard_True);
1525
1526   // validate Substruction operation
1527   NCollection_Map<Standard_Integer> aMapSubsLR;
1528   aMapSubsLR.Subtraction (aMapLeft, aMapRight);
1529   QCOMPARE (aMapSubsLR.Extent(), aRightLower - aLeftLower);
1530   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter < aRightLower; ++aKeyIter)
1531   {
1532     QCOMPARE (aMapSubsLR.Contains (aKeyIter), Standard_True);
1533   }
1534
1535   NCollection_Map<Standard_Integer> aMapSubsRL;
1536   aMapSubsRL.Subtraction (aMapRight, aMapLeft);
1537   QCOMPARE (aMapSubsRL.Extent(), aRightUpper - aLeftUpper);
1538   for (Standard_Integer aKeyIter = aLeftUpper + 1; aKeyIter < aRightUpper; ++aKeyIter)
1539   {
1540     QCOMPARE (aMapSubsRL.Contains (aKeyIter), Standard_True);
1541   }
1542
1543   // validate Difference operation
1544   NCollection_Map<Standard_Integer> aMapDiff;
1545   aMapDiff.Difference (aMapLeft, aMapRight);
1546   QCOMPARE (aMapDiff.Extent(), aRightLower - aLeftLower + aRightUpper - aLeftUpper);
1547   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter < aRightLower; ++aKeyIter)
1548   {
1549     QCOMPARE (aMapDiff.Contains (aKeyIter), Standard_True);
1550   }
1551   for (Standard_Integer aKeyIter = aLeftUpper + 1; aKeyIter < aRightUpper; ++aKeyIter)
1552   {
1553     QCOMPARE (aMapDiff.Contains (aKeyIter), Standard_True);
1554   }
1555
1556   // validate Exchange operation
1557   NCollection_Map<Standard_Integer> aMapSwap;
1558   aMapSwap.Exchange (aMapSect);
1559   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1560   {
1561     QCOMPARE (aMapSwap.Contains (aKeyIter), Standard_True);
1562   }
1563   QCOMPARE (aMapSect.IsEmpty(), Standard_True);
1564   aMapSwap.Add (34);
1565   aMapSect.Add (43);
1566
1567   NCollection_Map<Standard_Integer> aMapCopy (aMapSwap);
1568   QCOMPARE (aMapCopy.IsEqual (aMapSwap), Standard_True);
1569   aMapCopy.Remove (34);
1570   aMapCopy.Add    (43);
1571   QCOMPARE (aMapCopy.IsEqual (aMapSwap), Standard_False);
1572
1573   return 0;
1574 }
1575
1576 #define QVERIFY(val1) \
1577   di << "Checking " #val1 " == Standard_True" << \
1578         ((val1) == Standard_True ? ": OK\n" : ": Error\n")
1579
1580 #include <GeomInt_IntSS.hxx>
1581 #include <Geom_ConicalSurface.hxx>
1582 #include <Standard_ErrorHandler.hxx>
1583 //=======================================================================
1584 //function : OCC23972
1585 //purpose  :
1586 //=======================================================================
1587 static void DoGeomIntSSTest (const Handle(Geom_Surface)& theSurf1,
1588                              const Handle(Geom_Surface)& theSurf2,
1589                              const Standard_Integer theNbSol,
1590                              Draw_Interpretor& di)
1591 {
1592   try {
1593     OCC_CATCH_SIGNALS
1594          GeomInt_IntSS anInter;
1595          anInter.Perform (theSurf1, theSurf2, Precision::Confusion(), Standard_True);
1596          QVERIFY (anInter.IsDone());
1597          QCOMPARE (anInter.NbLines(), theNbSol);
1598   } catch (...) {
1599     QVERIFY (Standard_False);
1600   }
1601 }
1602
1603 namespace {
1604   static Handle(Geom_ConicalSurface) CreateCone (const gp_Pnt& theLoc,
1605                                                  const gp_Dir& theDir,
1606                                                  const gp_Dir& theXDir,
1607                                                  const Standard_Real theRad,
1608                                                  const Standard_Real theSin,
1609                                                  const Standard_Real theCos)
1610   {
1611     const Standard_Real anA = atan (theSin / theCos);
1612     gp_Ax3 anAxis (theLoc, theDir, theXDir);
1613     Handle(Geom_ConicalSurface) aSurf = new Geom_ConicalSurface (anAxis, anA, theRad);
1614     return aSurf;
1615   }
1616 }
1617
1618 static Standard_Integer OCC23972 (Draw_Interpretor& di,Standard_Integer n, const char**)
1619 {
1620   if (n != 1) return 1;
1621
1622   //process specific cones, cannot read them from files because due to rounding the original error
1623   //in math_FunctionRoots gets hidden
1624   Handle(Geom_Surface) aS1 = CreateCone (
1625                                          gp_Pnt (123.694345356663, 789.9, 68.15),
1626                                          gp_Dir (-1, 3.48029791472957e-016, -8.41302743359754e-017),
1627                                          gp_Dir (-3.48029791472957e-016, -1, -3.17572289932207e-016),
1628                                          3.28206830417112,
1629                                          0.780868809443031,
1630                                          0.624695047554424);
1631   Handle(Geom_Surface) aS2 = CreateCone (
1632                                          gp_Pnt (123.694345356663, 784.9, 68.15),
1633                                          gp_Dir (-1, -2.5209507537117e-016, -1.49772808948866e-016),
1634                                          gp_Dir (1.49772808948866e-016, 3.17572289932207e-016, -1),
1635                                          3.28206830417112,
1636                                          0.780868809443031,
1637                                          0.624695047554424);
1638
1639   DoGeomIntSSTest (aS1, aS2, 2, di);
1640
1641   return 0;
1642 }
1643
1644 #include <ShapeFix_EdgeProjAux.hxx>
1645 static Standard_Integer OCC24370 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
1646 {
1647   if (argc < 5) {
1648     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
1649     return 1;
1650   }
1651
1652   TopoDS_Shape aSh = DBRep::Get(argv[1]);
1653   if (aSh.IsNull()) {
1654     di << argv[0] << " Error: Null input edge\n";
1655     return 1;
1656   }
1657   const TopoDS_Edge& anEdge = TopoDS::Edge (aSh);
1658
1659   Handle(Geom2d_Curve) aC = DrawTrSurf::GetCurve2d(argv[2]);
1660   if (aC.IsNull()) {
1661     di << argv[0] << " Error: Null input curve\n";
1662     return 1;
1663   }
1664
1665   Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(argv[3]);
1666   if (aS.IsNull()) {
1667     di << argv[0] << " Error: Null input surface\n";
1668     return 1;
1669   }
1670
1671   Standard_Real prec = Draw::Atof(argv[4]);
1672
1673   //prepare data
1674   TopoDS_Face aFace;
1675   BRep_Builder aB;
1676   aB.MakeFace (aFace, aS, Precision::Confusion());
1677   aB.UpdateEdge (anEdge, aC, aFace, Precision::Confusion());
1678   aB.Range (anEdge, aFace, aC->FirstParameter(), aC->LastParameter());
1679
1680   //call algorithm
1681   ShapeFix_EdgeProjAux aProj (aFace, anEdge);
1682   aProj.Compute (prec);
1683
1684   Standard_Boolean isfirstdone = aProj.IsFirstDone();
1685   Standard_Boolean islastdone = aProj.IsLastDone();
1686
1687   Standard_Real first = 0.;
1688   Standard_Real last = 0.;
1689   Standard_Integer isfirstdoneInteger = 0;
1690   Standard_Integer islastdoneInteger = 0;
1691
1692
1693   if (isfirstdone) {
1694     first = aProj.FirstParam();
1695     isfirstdoneInteger = 1;
1696   }
1697
1698   if (islastdone) {
1699     last= aProj.LastParam();
1700     islastdoneInteger = 1;
1701   }
1702
1703   di << isfirstdoneInteger << " "<< islastdoneInteger << " "<< first << " "<< last << " \n";
1704
1705   return 0;
1706 }
1707
1708 template<typename T, typename HT>
1709 static void DoIsNull(Draw_Interpretor& di)
1710 {
1711   HT aHandle;
1712   //    QVERIFY (aHandle.IsNull());
1713   QCOMPARE (aHandle.IsNull(), Standard_True);
1714   const T* p = aHandle.Access();
1715 #if OCC_VERSION_HEX > 0x060700
1716   //QVERIFY (!p);
1717   //QVERIFY (p == 0);
1718   QCOMPARE (!p, Standard_True);
1719   QCOMPARE (p == 0, Standard_True);
1720 #endif
1721
1722   aHandle = new T;
1723   //QVERIFY (!aHandle.IsNull());
1724   QCOMPARE (!aHandle.IsNull(), Standard_True);
1725   p = aHandle.Access();
1726   //QVERIFY (p);
1727   //QVERIFY (p != 0);
1728   QCOMPARE (p != NULL, Standard_True);
1729   QCOMPARE (p != 0, Standard_True);
1730 }
1731
1732 //=======================================================================
1733 //function : OCC24533
1734 //purpose  :
1735 //=======================================================================
1736 static Standard_Integer OCC24533 (Draw_Interpretor& di, Standard_Integer n, const char**)
1737 {
1738   if (n != 1) return 1;
1739
1740   DoIsNull<Standard_Transient, Handle(Standard_Transient)>(di);
1741   DoIsNull<Standard_Persistent, Handle(Standard_Persistent)>(di);
1742
1743   return 0;
1744 }
1745
1746 // Dummy class to test interface for compilation issues
1747 class QABugs_HandleClass : public Standard_Transient
1748 {
1749 public:
1750   Standard_Integer HandleProc (Draw_Interpretor& , Standard_Integer  , const char** theArgVec)
1751   {
1752     std::cerr << "QABugs_HandleClass[" << this << "] " << theArgVec[0] << "\n";
1753     return 0;
1754   }
1755   DEFINE_STANDARD_RTTI(QABugs_HandleClass) // Type definition
1756 };
1757 DEFINE_STANDARD_HANDLE    (QABugs_HandleClass, Standard_Transient)
1758 IMPLEMENT_STANDARD_HANDLE (QABugs_HandleClass, Standard_Transient)
1759 IMPLEMENT_STANDARD_RTTIEXT(QABugs_HandleClass, Standard_Transient)
1760
1761 // Dummy class to test interface for compilation issues
1762 struct QABugs_NHandleClass
1763 {
1764   Standard_Integer NHandleProc (Draw_Interpretor& , Standard_Integer  , const char** theArgVec)
1765   {
1766     std::cerr << "QABugs_NHandleClass[" << this << "] " << "" << theArgVec[0] << "\n";
1767     return 0;
1768   }
1769 };
1770
1771 #include <XCAFDoc_ColorTool.hxx>
1772 #include <STEPControl_StepModelType.hxx>
1773 #include <STEPCAFControl_Writer.hxx>
1774 static Standard_Integer OCC23951 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1775 {
1776   if (argc != 1) {
1777     di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
1778     return 1;
1779   }
1780   Handle(TDocStd_Document) aDoc = new TDocStd_Document("dummy");;
1781   TopoDS_Shape s1 = BRepPrimAPI_MakeBox(1,1,1).Shape();
1782   TDF_Label lab1 = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main ())->NewShape();
1783   XCAFDoc_DocumentTool::ShapeTool (aDoc->Main ())->SetShape(lab1, s1);
1784   TDataStd_Name::Set(lab1, "Box1");
1785
1786   Quantity_Color yellow(1,1,0, Quantity_TOC_RGB);
1787   XCAFDoc_DocumentTool::ColorTool (aDoc->Main())->SetColor(lab1, yellow, XCAFDoc_ColorGen);
1788   XCAFDoc_DocumentTool::ColorTool(aDoc->Main())->SetVisibility(lab1, 0);
1789
1790   STEPControl_StepModelType mode = STEPControl_AsIs;
1791   STEPCAFControl_Writer writer;
1792   if ( ! writer.Transfer (aDoc, mode ) )
1793   {
1794     di << "The document cannot be translated or gives no result"  <<  "\n";
1795     return 1;
1796   }
1797
1798   writer.Write("test_box.step");
1799   return 0;
1800 }
1801
1802
1803 //=======================================================================
1804 //function : OCC24622
1805 //purpose  : The command tests sourcing Image_PixMap to AIS_TexturedShape
1806 //=======================================================================
1807 static Standard_Integer OCC24622 (Draw_Interpretor& /*theDi*/, Standard_Integer theArgNb, const char** theArgVec)
1808 {
1809   if (theArgNb != 2)
1810   {
1811     std::cout << "Usage : " << theArgVec[0] << " texture={1D|2D}";
1812     return 1;
1813   }
1814
1815   const Handle(AIS_InteractiveContext)& anAISContext = ViewerTest::GetAISContext();
1816   if (anAISContext.IsNull())
1817   {
1818     std::cout << "Please initialize view with \"vinit\".\n";
1819     return 1;
1820   }
1821
1822   Handle(Image_PixMap) anImage = new Image_PixMap();
1823
1824   static const Image_ColorRGB aBitmap[8] =
1825   {
1826     {{255,   0, 0}}, {{0,  148, 255}}, {{ 0, 148, 255}}, {{255,  94, 0}},
1827     {{255, 121, 0}}, {{76, 255,   0}}, {{76, 255,   0}}, {{255, 202, 0}}
1828   };
1829
1830   TCollection_AsciiString aTextureTypeArg (theArgVec[1]);
1831   aTextureTypeArg.UpperCase();
1832   if (aTextureTypeArg == "1D")
1833   {
1834     anImage->InitWrapper (Image_PixMap::ImgRGB, (Standard_Byte*)aBitmap, 8, 1);
1835   }
1836   else if (aTextureTypeArg == "2D")
1837   {
1838     anImage->InitTrash (Image_PixMap::ImgRGB, 8, 8);
1839     for (Standard_Integer aRow = 0; aRow < 8; ++aRow)
1840     {
1841       for (Standard_Integer aCol = 0; aCol < 8; ++aCol)
1842       {
1843         anImage->ChangeValue<Image_ColorRGB> (aRow, aCol) = aBitmap[aRow];
1844       }
1845     }
1846   }
1847   else
1848   {
1849     std::cout << "Please specify type of texture to test {1D|2D}.\n";
1850     return 1;
1851   }
1852
1853   TopoDS_Shape aBlankShape = BRepPrimAPI_MakeBox (10.0, 10.0, 10.0).Shape();
1854
1855   Handle(AIS_TexturedShape) aTexturedShape = new AIS_TexturedShape (aBlankShape);
1856   aTexturedShape->SetTexturePixMap (anImage);
1857   anAISContext->Display (aTexturedShape, 3, 0);
1858
1859   return 0;
1860 }
1861
1862 //=======================================================================
1863 //function : OCC24667
1864 //purpose  :
1865 //=======================================================================
1866 static Standard_Integer OCC24667 (Draw_Interpretor& di, Standard_Integer n, const char** a)
1867 {
1868   if (n == 1)
1869   {
1870     di << "OCC24667 result Wire_spine Profile [Mode [Approx]]" << "\n";
1871     di << "Mode = 0 - CorrectedFrenet," << "\n";
1872     di << "     = 1 - Frenet," << "\n";
1873     di << "     = 2 - DiscreteTrihedron" << "\n";
1874     di << "Approx - force C1-approximation if result is C0" << "\n";
1875     return 0;
1876   }
1877
1878   if (n > 1 && n < 4) return 1;
1879
1880   TopoDS_Shape Spine = DBRep::Get(a[2],TopAbs_WIRE);
1881   if ( Spine.IsNull()) return 1;
1882
1883   TopoDS_Shape Profile = DBRep::Get(a[3]);
1884   if ( Profile.IsNull()) return 1;
1885
1886   GeomFill_Trihedron Mode = GeomFill_IsCorrectedFrenet;
1887   if (n >= 5)
1888   {
1889     Standard_Integer iMode = atoi(a[4]);
1890     if (iMode == 1)
1891       Mode = GeomFill_IsFrenet;
1892     else if (iMode == 2)
1893       Mode = GeomFill_IsDiscreteTrihedron;
1894   }
1895
1896   Standard_Boolean ForceApproxC1 = Standard_False;
1897   if (n >= 6)
1898     ForceApproxC1 = Standard_True;
1899
1900   BRepOffsetAPI_MakePipe aPipe(TopoDS::Wire(Spine),
1901                                           Profile,
1902                                           Mode,
1903                                           ForceApproxC1);
1904
1905   TopoDS_Shape S = aPipe.Shape();
1906   TopoDS_Shape aSF = aPipe.FirstShape();
1907   TopoDS_Shape aSL = aPipe.LastShape();
1908
1909   DBRep::Set(a[1],S);
1910
1911   TCollection_AsciiString aStrF(a[1], "_f");
1912   TCollection_AsciiString aStrL(a[1], "_l");
1913
1914   DBRep::Set(aStrF.ToCString(), aSF);
1915   DBRep::Set(aStrL.ToCString(), aSL);
1916
1917   return 0;
1918 }
1919
1920 #include <IGESControl_Reader.hxx>
1921 #include <IGESControl_Controller.hxx>
1922 #include <IGESData_IGESEntity.hxx>
1923 #include <BRepCheck_Analyzer.hxx>
1924 #include <PTColStd_TransientPersistentMap.hxx>
1925 #include <PTopoDS_HShape.hxx>
1926 #include <Storage_Data.hxx>
1927 #include <TopExp_Explorer.hxx>
1928 #include <MgtBRep.hxx>
1929 #include <FSD_File.hxx>
1930 #include <ShapeSchema.hxx>
1931 #include <TColStd_HSequenceOfTransient.hxx>
1932 #include <PTColStd_PersistentTransientMap.hxx>
1933 #include <Storage_Root.hxx>
1934
1935 static Standard_Integer OCC24565 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
1936 {
1937   if (argc != 3) {
1938     di << "Usage : " << argv[0] << " FileNameIGS FileNameSTOR";
1939     return 1;
1940   }
1941
1942   Standard_CString sFileNameIGS = argv[1];
1943   Standard_CString sFileNameSTOR = argv[2];
1944
1945   IGESControl_Reader ICReader;
1946
1947   /* * * * * * *
1948    * Read the IGES file and make sure it is valid
1949    *
1950    * * * * * * */
1951   IGESControl_Controller::Init();
1952
1953   if (!ICReader.ReadFile(sFileNameIGS)) {
1954     printf("%s:%d - Error reading '%s'\n",__FUNCTION__,__LINE__,sFileNameIGS);fflush(stdout);
1955     return -1;
1956   }
1957
1958   int nbShapes = ICReader.NbShapes();
1959
1960   printf("%s:%d - nbShapes = '%d'\n",__FUNCTION__,__LINE__,nbShapes);fflush(stdout);
1961
1962   TopoDS_Shape Shape;
1963   if(nbShapes == 0)
1964     {
1965       Handle(TColStd_HSequenceOfTransient) faces=ICReader.GiveList("iges-faces");
1966       Handle(TColStd_HSequenceOfTransient) surfaceList=ICReader.GiveList("xst-transferrable-roots",faces);
1967
1968       if (surfaceList.IsNull())
1969         {
1970           printf("%s:%d - surfaceList.IsNull()\n",__FUNCTION__,__LINE__);fflush(stdout);
1971           return -1;
1972         }
1973       BRep_Builder builder;
1974       TopoDS_Compound* pC = new TopoDS_Compound();
1975       builder.MakeCompound(*pC);
1976
1977       for (int j=1;j<=surfaceList->Length();j++)
1978         {
1979           Handle(IGESData_IGESEntity) igesEntity=Handle(IGESData_IGESEntity)::DownCast(surfaceList->Value(j));
1980           if (igesEntity.IsNull()) continue;
1981           ICReader.ClearShapes();
1982           Standard_Boolean rv;
1983           try {
1984             rv=ICReader.TransferEntity(igesEntity);
1985           }
1986           catch (...) {
1987             rv=Standard_False;
1988           }
1989           if (!rv) {
1990             printf("%s:%d - Error transferring IGES entity\n",__FUNCTION__,__LINE__);fflush(stdout);
1991             printf("%s:%d - FormNumber = %d, TypeNumber = %d\n",__FUNCTION__,__LINE__,igesEntity->FormNumber(),igesEntity->TypeNumber());fflush(stdout);
1992             return -1;
1993           }
1994
1995           TopoDS_Shape S;
1996           try {
1997             S=ICReader.Shape();
1998           }
1999           catch(...) {
2000             printf("%s:%d - Error reading IGES entity\n",__FUNCTION__,__LINE__);fflush(stdout);
2001             printf("%s:%d - FormNumber = %d, TypeNumber = %d\n",__FUNCTION__,__LINE__,igesEntity->FormNumber(),igesEntity->TypeNumber());fflush(stdout);
2002             return -1;
2003           }
2004           if (S.IsNull()) {
2005             printf("%s:%d - NULL Surface encountered\n",__FUNCTION__,__LINE__);
2006             return -1;
2007           }
2008
2009           try
2010             {
2011               builder.Add(*pC,S);
2012             }
2013           catch(...)
2014             {
2015               printf("%s: Exception adding face.\n",__FUNCTION__);
2016             }
2017         }
2018       Shape = TopoDS_Shape(*pC);
2019     }
2020   else
2021     {
2022       Shape = ICReader.OneShape();
2023     }
2024   {
2025     BRepCheck_Analyzer brca(Shape);
2026
2027     if(!brca.IsValid())
2028       {
2029         printf("%s: Invalid shape after reading IGES file.\n",__FUNCTION__);
2030       }
2031   }
2032
2033   /* * * * * * *
2034    * Write the contents of the Shape to a STOR file
2035    *
2036    * * * * * * */
2037   PTColStd_TransientPersistentMap aMapTP;
2038   Handle(PTopoDS_HShape) aPShape_write;
2039   Handle(Storage_Data) d_write=new Storage_Data;
2040   char Name[32];
2041
2042   TopExp_Explorer Ex;
2043   int i;
2044   int max_i = 0;
2045
2046   for (i=0,Ex.Init(Shape,TopAbs_FACE);Ex.More();i++,Ex.Next())
2047     {
2048
2049       max_i = i;
2050       try {
2051         aPShape_write=MgtBRep::Translate(Ex.Current(),aMapTP,MgtBRep_WithoutTriangle);
2052       }
2053       catch (...) {
2054         printf("%s: Error translating surface '%d'\n",__FUNCTION__,i);
2055       }
2056
2057       sprintf(Name,"S%010d",i);
2058
2059       {
2060         BRepCheck_Analyzer brca(Ex.Current());
2061         if(!brca.IsValid())
2062           {
2063             printf("INVALID face '%s' in the shape, which will be written to the STOR file.\n",Name);
2064           }
2065       }
2066       try {
2067         d_write->AddRoot(Name,aPShape_write);
2068       }
2069       catch (...) {
2070         printf("%s: Error adding surface '%d', RootName = '%s'\n",__FUNCTION__,i,Name);
2071       }
2072     }
2073   printf("%s: Going to write %d surfaces.\n",__FUNCTION__,max_i+1);
2074
2075   FSD_File f_write;
2076   if(f_write.Open(sFileNameSTOR, Storage_VSWrite)!=Storage_VSOk)
2077     {
2078       printf("%s: Error opening file: %s\n", __FUNCTION__,sFileNameSTOR);
2079       return -1;
2080     }
2081   Handle(ShapeSchema) s_write=new ShapeSchema;
2082   s_write->Write(f_write,d_write);
2083   f_write.Close();
2084   printf("%s: Wrote to the STOR file.\n",__FUNCTION__);
2085
2086   /* * * * * * *
2087    * Read the contents of the Shape from a STOR file
2088    *
2089    * * * * * * */
2090   FSD_File f_read;
2091   if(f_read.Open(sFileNameSTOR, Storage_VSRead)!=Storage_VSOk)
2092     {
2093       printf("%s: Error opening file: %s\n", __FUNCTION__,sFileNameSTOR);
2094       return -1;
2095     }
2096   Handle(ShapeSchema) s_read=new ShapeSchema;
2097   Handle(Storage_Data) d_read=s_read->Read(f_read);
2098
2099   Handle(Standard_Persistent) p;
2100   Handle(Storage_Root) r;
2101   Handle(PTopoDS_HShape) aPShape_read;
2102   PTColStd_PersistentTransientMap aMapPT;
2103   TopoDS_Shape S_read;
2104
2105   printf("%s: Extracting %d faces from the STOR file.\n",__FUNCTION__,max_i+1);
2106   for(int i = 0; i <= max_i; ++i)
2107     {
2108       sprintf(Name,"S%010d",i);
2109       r=d_read->Find(Name);
2110       if(r.IsNull())
2111         {
2112           printf("%s:%d '%s' IsNull().\n",__FUNCTION__,__LINE__,Name);fflush(stdout);
2113           continue;
2114         }
2115       p=r->Object();
2116       aPShape_read = Handle(PTopoDS_HShape)::DownCast(p);
2117       try {
2118         MgtBRep::Translate(aPShape_read,aMapPT,S_read,MgtBRep_WithoutTriangle);
2119       }
2120       catch (Standard_Failure) {
2121         Handle(Standard_Failure) E=Standard_Failure::Caught();
2122         std::string str;
2123         str="Exception: ";
2124         str+=E->DynamicType()->Name();
2125         str+=" => ";
2126         str+=E->GetMessageString();
2127         printf("%s(1): %s: %s\n",__FUNCTION__,Name,str.c_str());fflush(stdout);
2128       }
2129       catch (...) {
2130         printf("%s(1): Unhandled exception in MgtBRep::Translate\n",__FUNCTION__);
2131       }
2132
2133       BRepCheck_Analyzer brca(S_read);
2134
2135       if(!brca.IsValid())
2136         {
2137           printf("%s: Read INVALID face (%s)!\n",__FUNCTION__,Name);
2138         }
2139     }
2140
2141   printf("Completed.\n");fflush(stdout);
2142
2143   return 0;
2144 }
2145
2146 #include <Handle_BRepTools_NurbsConvertModification.hxx>
2147 #include <BRepPrimAPI_MakeCylinder.hxx>
2148 #include <BRepBuilderAPI_Copy.hxx>
2149 #include <BRepTools_NurbsConvertModification.hxx>
2150 static TopoDS_Shape CreateTestShape (int& theShapeNb)
2151 {
2152   TopoDS_Compound aComp;
2153   BRep_Builder aBuilder;
2154   aBuilder.MakeCompound (aComp);
2155   //NURBS modifier is used to increase footprint of each shape
2156   Handle_BRepTools_NurbsConvertModification aNurbsModif = new BRepTools_NurbsConvertModification;
2157   TopoDS_Shape aRefShape = BRepPrimAPI_MakeCylinder (50., 100.).Solid();
2158   BRepTools_Modifier aModifier (aRefShape, aNurbsModif);
2159   if (aModifier.IsDone()) {
2160     aRefShape = aModifier.ModifiedShape (aRefShape);
2161   }
2162   int aSiblingNb = 0;
2163   for (; theShapeNb > 0; --theShapeNb) {
2164     TopoDS_Shape aShape;
2165     if (++aSiblingNb <= 100) { //number of siblings is limited to avoid long lists
2166                 aShape = BRepBuilderAPI_Copy (aRefShape, Standard_True /*CopyGeom*/).Shape();
2167     } else {
2168       aShape = CreateTestShape (theShapeNb);
2169     }
2170     aBuilder.Add (aComp, aShape);
2171   }
2172   return aComp;
2173 }
2174
2175 #include <AppStd_Application.hxx>
2176 #include <TDataStd_Integer.hxx>
2177 #include <TNaming_Builder.hxx>
2178 static Standard_Integer OCC24931 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
2179 {
2180   if (argc != 1) {
2181     di << "Usage: " << argv[0] << " invalid number of arguments"<<"\n";
2182     return 1;
2183   }
2184   TCollection_ExtendedString aFileName ("testdocument.xml");
2185   PCDM_StoreStatus aSStatus  = PCDM_SS_Failure;
2186
2187   Handle(TDocStd_Application) anApp = new AppStd_Application;
2188   {
2189     Handle(TDocStd_Document) aDoc;
2190     anApp->NewDocument ("XmlOcaf", aDoc);
2191     TDF_Label aLab = aDoc->Main();
2192     TDataStd_Integer::Set (aLab, 0);
2193     int n = 10000; //must be big enough
2194     TopoDS_Shape aShape = CreateTestShape (n);
2195     TNaming_Builder aBuilder (aLab);
2196     aBuilder.Generated (aShape);
2197
2198     aSStatus = anApp->SaveAs (aDoc, aFileName);
2199     anApp->Close (aDoc);
2200   }
2201   QCOMPARE (aSStatus, PCDM_SS_OK);
2202   return 0;
2203 }
2204
2205 #include <AppStdL_Application.hxx>
2206 #include <TDocStd_Application.hxx>
2207 #include <TDataStd_Integer.hxx>
2208 #include <TDF_AttributeIterator.hxx>
2209 //=======================================================================
2210 //function : OCC24755
2211 //purpose  :
2212 //=======================================================================
2213 static Standard_Integer OCC24755 (Draw_Interpretor& di, Standard_Integer n, const char** a)
2214 {
2215   if (n != 1)
2216   {
2217     std::cout << "Usage : " << a[0] << "\n";
2218     return 1;
2219   }
2220
2221   Handle(TDocStd_Application) anApp = new AppStdL_Application;
2222   Handle(TDocStd_Document) aDoc;
2223   anApp->NewDocument ("MDTV-Standard", aDoc);
2224   TDF_Label aLab = aDoc->Main();
2225   TDataStd_Integer::Set (aLab, 0);
2226   TDataStd_Name::Set (aLab, "test");
2227
2228   TDF_AttributeIterator i (aLab);
2229   Handle(TDF_Attribute) anAttr = i.Value();
2230   QCOMPARE (anAttr->IsKind (STANDARD_TYPE (TDataStd_Integer)), Standard_True);
2231   i.Next();
2232   anAttr = i.Value();
2233   QCOMPARE (anAttr->IsKind (STANDARD_TYPE (TDataStd_Name)), Standard_True);
2234
2235   return 0;
2236 }
2237
2238 struct MyStubObject
2239 {
2240   MyStubObject() : ptr(0L) {}
2241   MyStubObject(void* thePtr) : ptr(thePtr) {}
2242   char overhead[40];
2243   void* ptr;
2244 };
2245
2246 //=======================================================================
2247 //function : OCC24834
2248 //purpose  :
2249 //=======================================================================
2250 static Standard_Integer OCC24834 (Draw_Interpretor& di, Standard_Integer n, const char** a)
2251 {
2252   if (n != 1)
2253   {
2254     std::cout << "Usage : " << a[0] << "\n";
2255     return 1;
2256   }
2257
2258   NCollection_List<MyStubObject> aList;
2259   const Standard_Integer aSmallBlockSize = 40;
2260   const Standard_Integer aLargeBlockSize = 1500000;
2261
2262   // quick populate memory with large blocks
2263   try
2264   {
2265     for (;;)
2266     {
2267       aList.Append(MyStubObject(Standard::Allocate(aLargeBlockSize)));
2268     }
2269   }
2270   catch (Standard_Failure)
2271   {
2272     di << "caught out of memory for large blocks: OK\n";
2273   }
2274   catch (...)
2275   {
2276     di << "skept out of memory for large blocks: Error\n";
2277   }
2278
2279   // allocate small blocks
2280   try
2281   {
2282     for (;;)
2283     {
2284       aList.Append(MyStubObject(Standard::Allocate(aSmallBlockSize)));
2285     }
2286   }
2287   catch (Standard_Failure)
2288   {
2289     di << "caught out of memory for small blocks: OK\n";
2290   }
2291   catch (...)
2292   {
2293     di << "skept out of memory for small blocks: Error\n";
2294   }
2295
2296   // release all allocated blocks
2297   for (NCollection_List<MyStubObject>::Iterator it(aList); it.More(); it.Next())
2298   {
2299     Standard::Free(it.Value().ptr);
2300   }
2301   return 0;
2302 }
2303
2304
2305 #include <Geom2dAPI_InterCurveCurve.hxx>
2306 #include <IntRes2d_IntersectionPoint.hxx>
2307 //=======================================================================
2308 //function : OCC24889
2309 //purpose  :
2310 //=======================================================================
2311 static Standard_Integer OCC24889 (Draw_Interpretor& theDI,
2312                                   Standard_Integer /*theNArg*/,
2313                                   const char** /*theArgs*/)
2314 {
2315  // Curves
2316   Handle( Geom2d_Circle ) aCircle1 = new Geom2d_Circle(
2317     gp_Ax22d( gp_Pnt2d( 25, -25 ), gp_Dir2d( 1, 0 ), gp_Dir2d( -0, 1 ) ), 155 );
2318
2319   Handle( Geom2d_Circle ) aCircle2 = new Geom2d_Circle(
2320     gp_Ax22d( gp_Pnt2d( 25, 25 ), gp_Dir2d( 1, 0 ), gp_Dir2d( -0, 1 ) ), 155 );
2321
2322   Handle( Geom2d_TrimmedCurve ) aTrim[2] = {
2323     new Geom2d_TrimmedCurve( aCircle1, 1.57079632679490, 2.97959469729228 ),
2324     new Geom2d_TrimmedCurve( aCircle2, 3.30359060633978, 4.71238898038469 )
2325   };
2326
2327   DrawTrSurf::Set("c_1", aTrim[0]);
2328   DrawTrSurf::Set("c_2", aTrim[1]);
2329
2330   // Intersection
2331   const Standard_Real aTol = Precision::Confusion();
2332   Geom2dAPI_InterCurveCurve aIntTool( aTrim[0], aTrim[1], aTol );
2333
2334   const IntRes2d_IntersectionPoint& aIntPnt =
2335     aIntTool.Intersector().Point( 1 );
2336
2337   gp_Pnt2d aIntRes = aIntTool.Point( 1 );
2338   Standard_Real aPar[2] = {
2339     aIntPnt.ParamOnFirst(),
2340     aIntPnt.ParamOnSecond()
2341   };
2342
2343   //theDI.precision( 5 );
2344   theDI << "Int point: X = " << aIntRes.X() << "; Y = " << aIntRes.Y() << "\n";
2345   for (int i = 0; i < 2; ++i)
2346   {
2347     theDI << "Curve " << i << ": FirstParam = " << aTrim[i]->FirstParameter() <<
2348       "; LastParam = " << aTrim[i]->LastParameter() <<
2349       "; IntParameter = " << aPar[i] << "\n";
2350   }
2351
2352   return 0;
2353 }
2354
2355
2356
2357 void QABugs::Commands_19(Draw_Interpretor& theCommands) {
2358   const char *group = "QABugs";
2359
2360   Handle(QABugs_HandleClass) aClassPtr = new QABugs_HandleClass();
2361   theCommands.Add ("OCC24202_1", "Test Handle-based procedure",
2362                    __FILE__, aClassPtr, &QABugs_HandleClass::HandleProc, group);
2363   NCollection_Handle<QABugs_NHandleClass> aNClassPtr = new QABugs_NHandleClass();
2364   theCommands.Add ("OCC24202_2", "Test NCollection_Handle-based procedure",
2365                    __FILE__, aNClassPtr, &QABugs_NHandleClass::NHandleProc, group);
2366
2367   theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
2368   theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
2369   theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
2370   theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group);
2371   theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
2372   theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
2373   theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
2374   theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
2375   theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
2376   theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
2377   theCommands.Add ("OCC23952sweep", "OCC23952sweep nbupoles shape", __FILE__, OCC23952sweep, group);
2378   theCommands.Add ("OCC23952intersect", "OCC23952intersect nbsol shape1 shape2", __FILE__, OCC23952intersect, group);
2379   theCommands.Add ("test_offset", "test_offset", __FILE__, test_offset, group);
2380   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);
2381   theCommands.Add ("OCC24008", "OCC24008 curve surface", __FILE__, OCC24008, group);
2382   theCommands.Add ("OCC24019", "OCC24019 aShape", __FILE__, OCC24019, group);
2383   theCommands.Add ("OCC11758", "OCC11758", __FILE__, OCC11758, group);
2384   theCommands.Add ("OCC24005", "OCC24005 result", __FILE__, OCC24005, group);
2385   theCommands.Add ("OCC24137", "OCC24137 face vertex U V [N]", __FILE__, OCC24137, group);
2386   theCommands.Add ("OCC24271", "Boolean operations on NCollection_Map", __FILE__, OCC24271, group);
2387   theCommands.Add ("OCC23972", "OCC23972", __FILE__, OCC23972, group);
2388   theCommands.Add ("OCC24370", "OCC24370 edge pcurve surface prec", __FILE__, OCC24370, group);
2389   theCommands.Add ("OCC24533", "OCC24533", __FILE__, OCC24533, group);
2390   theCommands.Add ("OCC24012", "OCC24012 face edge", __FILE__, OCC24012, group);
2391   theCommands.Add ("OCC24051", "OCC24051", __FILE__, OCC24051, group);
2392   theCommands.Add ("OCC24086", "OCC24086 face wire", __FILE__, OCC24086, group);
2393   theCommands.Add ("OCC24622", "OCC24622 texture={1D|2D}\n Tests sourcing of 1D/2D pixmaps for AIS_TexturedShape", __FILE__, OCC24622, group);
2394   theCommands.Add ("OCC24667", "OCC24667 result Wire_spine Profile [Mode [Approx]], no args to get help", __FILE__, OCC24667, group);
2395   theCommands.Add ("OCC24565", "OCC24565 FileNameIGS FileNameSTOR", __FILE__, OCC24565, group);
2396   theCommands.Add ("OCC24755", "OCC24755", __FILE__, OCC24755, group);
2397   theCommands.Add ("OCC24834", "OCC24834", __FILE__, OCC24834, group);
2398   theCommands.Add ("OCC24889", "OCC24889", __FILE__, OCC24889, group);
2399   theCommands.Add ("OCC23951", "OCC23951", __FILE__, OCC23951, group);
2400   theCommands.Add ("OCC24931", "OCC24931", __FILE__, OCC24931, group);
2401   theCommands.Add ("OCC24945", "OCC24945", __FILE__, OCC24945, group);
2402   return;
2403 }