0024931: Stack overflow when writing large shapes to XML
[occt.git] / src / QABugs / QABugs_19.cxx
1 // Created on: 2002-05-21
2 // Created by: QA Admin
3 // Copyright (c) 2002-2014 OPEN CASCADE SAS
4 //
5 // This file is part of Open CASCADE Technology software library.
6 //
7 // This library is free software; you can redistribute it and/or modify it under
8 // the terms of the GNU Lesser General Public License version 2.1 as published
9 // by the Free Software Foundation, with special exception defined in the file
10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
11 // distribution for complete text of the license and disclaimer of any warranty.
12 //
13 // Alternatively, this file may be used under the terms of Open CASCADE
14 // commercial license or contractual agreement.
15
16 #include <QABugs.hxx>
17
18 #include <Draw_Interpretor.hxx>
19 #include <DBRep.hxx>
20 #include <DrawTrSurf.hxx>
21 #include <ViewerTest.hxx>
22 #include <TopoDS_Shape.hxx>
23 #include <AIS_InteractiveContext.hxx>
24 #include <AIS_TexturedShape.hxx>
25 #include <Image_PixMap.hxx>
26 #include <Image_Color.hxx>
27
28 #include <gp_Pnt2d.hxx>
29 #include <gp_Ax1.hxx>
30 #include <GCE2d_MakeSegment.hxx>
31 #include <Geom2d_TrimmedCurve.hxx>
32 #include <DrawTrSurf.hxx>
33
34 #include <Precision.hxx>
35
36 #include <PCollection_HAsciiString.hxx>
37
38 #include <cstdio>
39 #include <cmath>
40 #include <iostream>
41 #include <OSD_PerfMeter.hxx>
42 #include <OSD_Timer.hxx>
43 #include <BRepPrimAPI_MakeBox.hxx>
44 #include <BRepPrimAPI_MakeSphere.hxx>
45 #include <BRepAlgo_Cut.hxx>
46 #include <NCollection_Map.hxx>
47 #include <NCollection_Handle.hxx>
48 #include <TCollection_HAsciiString.hxx>
49 #include <GeomFill_Trihedron.hxx>
50 #include <BRepOffsetAPI_MakePipe.hxx>
51
52 #include <Standard_Version.hxx>
53
54 #define QCOMPARE(val1, val2) \
55   di << "Checking " #val1 " == " #val2 << \
56         ((val1) == (val2) ? ": OK\n" : ": Error\n")
57
58 static Standard_Integer OCC230 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
59 {
60   if ( argc != 4) {
61     di << "ERROR OCC230: Usage : " << argv[0] << " TrimmedCurve Pnt2d Pnt2d" << "\n";
62     return 1;
63   }
64
65   gp_Pnt2d P1, P2;
66   if ( !DrawTrSurf::GetPoint2d(argv[2],P1)) {
67     di << "ERROR OCC230: " << argv[2] << " is not Pnt2d" << "\n";
68     return 1;
69   }
70   if ( !DrawTrSurf::GetPoint2d(argv[3],P2)) {
71     di << "ERROR OCC230: " << argv[3] << " is not Pnt2d" << "\n";
72     return 1;
73   }
74
75   GCE2d_MakeSegment MakeSegment(P1,P2);
76   Handle(Geom2d_TrimmedCurve) TrimmedCurve = MakeSegment.Value();
77   DrawTrSurf::Set(argv[1], TrimmedCurve);
78   return 0;
79 }
80
81 static Standard_Integer OCC142 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
82 {
83   for(Standard_Integer i= 0;i <= 20;i++){
84     Handle(PCollection_HAsciiString) pstr = new PCollection_HAsciiString("TEST");
85     pstr->Clear();
86   }
87   di << "OCC142: OK" << "\n";
88   return 0;
89 }
90
91 static Standard_Integer OCC23361 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
92 {
93   gp_Pnt p(0, 0, 2);
94   
95   gp_Trsf t1, t2;
96   t1.SetRotation(gp_Ax1(p, gp_Dir(0, 1, 0)), -0.49328285294022267);
97   t2.SetRotation(gp_Ax1(p, gp_Dir(0, 0, 1)), 0.87538474718473880);
98
99   gp_Trsf tComp = t2 * t1;
100
101   gp_Pnt p1(10, 3, 4);
102   gp_Pnt p2 = p1.Transformed(tComp);
103   gp_Pnt p3 = p1.Transformed(t1);
104   p3.Transform(t2);
105
106   // points must be equal
107   if ( ! p2.IsEqual(p3, Precision::Confusion()) )
108     di << "ERROR OCC23361: equivalent transformations does not produce equal points" << "\n";
109   else 
110     di << "OCC23361: OK" << "\n";
111
112   return 0;
113 }
114
115 static Standard_Integer OCC23237 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
116 {
117   OSD_PerfMeter aPM("TestMeter",0);
118   OSD_Timer aTM;
119   
120   // run some operation in cycle for about 2 seconds to have good values of times to compare
121   int count = 0;
122   printf("OSD_PerfMeter test.\nRunning Boolean operation on solids in loop.\n");
123   for (; aTM.ElapsedTime() < 2.; count++)
124   {
125     aPM.Start();
126     aTM.Start();
127
128     // do some operation that will take considerable time compared with time of starting / stopping timers
129     BRepPrimAPI_MakeBox aBox (10., 10., 10.);
130     BRepPrimAPI_MakeSphere aSphere (10.);
131     BRepAlgo_Cut aCutter (aBox.Shape(), aSphere.Shape());
132
133     aTM.Stop();
134     aPM.Stop();
135   }
136  
137   int aNbEnters = 0;
138   Standard_Real aPerfMeter_CPUtime = 0., aTimer_ElapsedTime = aTM.ElapsedTime();
139
140   perf_get_meter("TestMeter", &aNbEnters, &aPerfMeter_CPUtime);
141
142   Standard_Real aTimeDiff = (fabs(aTimer_ElapsedTime - aPerfMeter_CPUtime) / aTimer_ElapsedTime);
143
144   printf("\nMeasurement results (%d cycles):\n", count);
145   printf("\nOSD_PerfMeter CPU time: %lf\nOSD_Timer elapsed time: %lf\n", aPerfMeter_CPUtime, aTimer_ElapsedTime);
146   printf("Time delta is: %.3lf %%\n", aTimeDiff * 100);
147
148   if (aTimeDiff > 0.2)
149     di << "OCC23237: Error: too much difference between CPU and elapsed times";
150   else if (aNbEnters != count)
151     di << "OCC23237: Error: counter reported by PerfMeter (" << aNbEnters << ") does not correspond to actual number of cycles";
152   else
153     di << "OCC23237: OK";
154
155   return 0;
156 }
157
158 #ifdef HAVE_TBB
159
160 #include <Standard_Atomic.hxx>
161 #include <tbb/blocked_range.h>
162 #include <tbb/parallel_for.h>
163
164 class IncrementerDecrementer
165 {
166 public:
167     IncrementerDecrementer (Standard_Integer* theVal, Standard_Boolean thePositive) : myVal (theVal), myPositive (thePositive)
168     {}
169     void operator() (const tbb::blocked_range<size_t>& r) const
170     {
171         if (myPositive)
172             for (size_t i = r.begin(); i != r.end(); ++i)
173                 Standard_Atomic_Increment (myVal);
174         else
175             for (size_t i = r.begin(); i != r.end(); ++i)
176                 Standard_Atomic_Decrement (myVal);
177     }
178 private:
179     Standard_Integer*   myVal;
180     Standard_Boolean   myPositive;
181 };
182 #endif
183
184 #ifdef HAVE_TBB
185 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
186 {
187   int aSum = 0;
188
189   //check returned value
190   QCOMPARE (Standard_Atomic_Decrement (&aSum), -1);
191   QCOMPARE (Standard_Atomic_Increment (&aSum), 0);
192   QCOMPARE (Standard_Atomic_Increment (&aSum), 1);
193   QCOMPARE (Standard_Atomic_Increment (&aSum), 2);
194 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 0);
195 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 1);
196
197   //check atomicity 
198   aSum = 0;
199   const int N = 1 << 24; //big enough to ensure concurrency
200
201   //increment
202   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, true));
203   QCOMPARE (aSum, N);
204
205   //decrement
206   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
207   QCOMPARE (aSum, 0);
208
209   return 0;
210 }
211
212 #else /* HAVE_TBB */
213
214 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
215 {
216   di << "Test skipped: command " << argv[0] << " requires TBB library\n";
217   return 0;
218 }
219
220 #endif /* HAVE_TBB */
221
222 #include <TDocStd_Application.hxx>
223 #include <XCAFApp_Application.hxx>
224 #include <TDocStd_Document.hxx>
225 #include <XCAFDoc_ShapeTool.hxx>
226 #include <XCAFDoc_DocumentTool.hxx>
227 #include <TDF_Label.hxx>
228 #include <TDataStd_Name.hxx>
229
230 static Standard_Integer OCC23595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
231 {
232   const Handle(TDocStd_Application)& anApp = XCAFApp_Application::GetApplication();
233   Handle(TDocStd_Document) aDoc;
234   anApp->NewDocument ("XmlXCAF", aDoc);
235   QCOMPARE (!aDoc.IsNull(), Standard_True);
236
237   Handle(XCAFDoc_ShapeTool) aShTool = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main());
238
239   //check default value
240   Standard_Boolean aValue = XCAFDoc_ShapeTool::AutoNaming();
241   QCOMPARE (aValue, Standard_True);
242
243   //true
244   XCAFDoc_ShapeTool::SetAutoNaming (Standard_True);
245   TopoDS_Shape aShape = BRepPrimAPI_MakeBox (100., 200., 300.).Shape();
246   TDF_Label aLabel = aShTool->AddShape (aShape);
247   Handle(TDataStd_Name) anAttr;
248   QCOMPARE (aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
249
250   //false
251   XCAFDoc_ShapeTool::SetAutoNaming (Standard_False);
252   aShape = BRepPrimAPI_MakeBox (300., 200., 100.).Shape();
253   aLabel = aShTool->AddShape (aShape);
254   QCOMPARE (!aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
255
256   //restore
257   XCAFDoc_ShapeTool::SetAutoNaming (aValue);
258
259   return 0;
260 }
261
262 #include <ExprIntrp_GenExp.hxx>
263 Standard_Integer OCC22611 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
264 {
265
266   if (argc != 3) {
267     di << "Usage : " << argv[0] << " string nb\n";
268     return 1;
269   }
270
271   TCollection_AsciiString aToken = argv[1];
272   Standard_Integer aNb = atoi(argv[2]);
273
274   Handle(ExprIntrp_GenExp) aGen = ExprIntrp_GenExp::Create();
275   for (Standard_Integer i=0; i < aNb; i++)
276   {
277     aGen->Process(aToken);
278     Handle(Expr_GeneralExpression) aExpr = aGen->Expression();
279   }
280
281   return 0;
282 }
283
284 Standard_Integer OCC22595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
285 {
286   gp_Mat M0;
287   di << "M0 = "
288   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}"
289   << " {" << M0(2,1) << "} {" << M0(2,2) << "} {" << M0(2,3) <<"}"
290   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}";
291   return 0;
292 }
293
294 #include <TopoDS_Face.hxx>
295 #include <TopoDS_Face.hxx>
296 #include <TopoDS.hxx>
297 #include <BRepBuilderAPI_Transform.hxx>
298 #include <BRepExtrema_DistShapeShape.hxx>
299 #include <BRepTools.hxx>
300
301 static Standard_Boolean OCC23774Test(const TopoDS_Face& grossPlateFace, const TopoDS_Shape& originalWire, Draw_Interpretor& di)
302 {
303   BRepExtrema_DistShapeShape distShapeShape(grossPlateFace,originalWire,Extrema_ExtFlag_MIN);
304   if(!distShapeShape.IsDone()) {
305     di <<"Distance ShapeShape is Not Done\n";
306     return Standard_False;
307   }
308
309   if(distShapeShape.Value() > 0.01) {
310     di << "Wrong Dist = " <<distShapeShape.Value() << "\n";
311     return Standard_False;
312   } else
313     di << "Dist0 = " <<distShapeShape.Value() <<"\n";
314
315   //////////////////////////////////////////////////////////////////////////
316   /// First Flip Y
317   const gp_Pnt2d axis1P1(1474.8199035519228,1249.9995745636970);
318   const gp_Pnt2d axis1P2(1474.8199035519228,1250.9995745636970);
319
320   gp_Vec2d mirrorVector1(axis1P1,axis1P2);
321
322   gp_Trsf2d mirror1;
323   mirror1.SetMirror(gp_Ax2d(axis1P1,mirrorVector1));
324
325   BRepBuilderAPI_Transform transformer1(mirror1);
326   transformer1.Perform(originalWire);
327   if(!transformer1.IsDone()) {
328     di << "Not Done1 " << "\n";
329     return Standard_False;
330   }
331   const TopoDS_Shape& step1ModifiedShape = transformer1.ModifiedShape(originalWire);
332   
333   BRepExtrema_DistShapeShape distShapeShape1(grossPlateFace,step1ModifiedShape,Extrema_ExtFlag_MIN);
334   if(!distShapeShape1.IsDone())
335     return Standard_False;
336   if(distShapeShape1.Value() > 0.01) {
337     di << "Dist = " <<distShapeShape1.Value() <<"\n";
338     return Standard_False;
339   } else
340     di << "Dist1 = " <<distShapeShape1.Value() <<"\n";
341
342   //////////////////////////////////////////////////////////////////////////
343   /// Second flip Y
344   transformer1.Perform(step1ModifiedShape);
345   if(!transformer1.IsDone()) {
346     di << "Not Done1 \n";
347     return Standard_False;
348   }
349   const TopoDS_Shape& step2ModifiedShape = transformer1.ModifiedShape(step1ModifiedShape);
350
351   //This is identity matrix for values but for type is gp_Rotation ?!
352   gp_Trsf2d mirror11 = mirror1;
353   mirror11.PreMultiply(mirror1);
354
355   BRepExtrema_DistShapeShape distShapeShape2(grossPlateFace,step2ModifiedShape);//,Extrema_ExtFlag_MIN);
356   if(!distShapeShape2.IsDone())
357     return Standard_False;
358
359   //This last test case give error (the value is 1008.8822038689706)
360   if(distShapeShape2.Value() > 0.01) {
361     di  << "Wrong Dist2 = " <<distShapeShape2.Value() <<"\n";
362     Standard_Integer N = distShapeShape2.NbSolution();
363     di << "Nb = " <<N <<"\n";
364     for (Standard_Integer i=1;i <= N;i++)
365         di <<"Sol(" <<i<<") = " <<distShapeShape2.PointOnShape1(i).Distance(distShapeShape2.PointOnShape2(i)) <<"\n";
366     return Standard_False;
367   }
368   di << "Distance2 = " <<distShapeShape2.Value() <<"\n";
369  
370   return Standard_True;
371 }
372 static Standard_Integer OCC23774(Draw_Interpretor& di, Standard_Integer n, const char** a)
373
374
375   if (n != 3) {
376         di <<"OCC23774: invalid number of input parameters\n";
377         return 1;
378   }
379
380   const char *ns1 = (a[1]), *ns2 = (a[2]);
381   TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2));
382   if (S1.IsNull() || S2.IsNull()) {
383         di <<"OCC23774: Null input shapes\n";
384         return 1;
385   }
386   const TopoDS_Face& aFace  = TopoDS::Face(S1);
387   if(!OCC23774Test(aFace, S2, di))
388         di << "Something is wrong\n";
389
390  return 0;
391 }
392
393 #include <GeomConvert_ApproxSurface.hxx>
394 #include <Geom_BSplineSurface.hxx>
395 #include <Draw.hxx>
396 #include <OSD_Thread.hxx>
397 static void GeomConvertTest (Draw_Interpretor& di, Standard_Integer theTargetNbUPoles, Standard_CString theFileName)
398 {
399         Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(theFileName);
400         GeomConvert_ApproxSurface aGAS (aSurf, 1e-4, GeomAbs_C1, GeomAbs_C1, 9, 9, 100, 1);
401         if (!aGAS.IsDone()) {
402                 di << "ApproxSurface is not done!" << "\n";
403                 return;
404         }
405         const Handle(Geom_BSplineSurface)& aBSurf = aGAS.Surface();
406         if (aBSurf.IsNull()) {
407                 di << "BSplineSurface is not created!" << "\n";
408                 return;
409         }
410         di << "Number of UPoles:" << aBSurf->NbUPoles() << "\n";
411         QCOMPARE (aBSurf->NbUPoles(), theTargetNbUPoles);
412 }
413
414 struct aData {
415         Draw_Interpretor* di;
416         Standard_Integer nbupoles;
417         Standard_CString filename;
418 };
419
420 Standard_EXPORT Standard_Address convert(Standard_Address data)
421 {
422         aData* info = (aData*) data;
423         GeomConvertTest(*(info->di),info->nbupoles,info->filename);
424         return NULL;
425 }
426
427 static Standard_Integer OCC23952sweep (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
428 {
429         if (argc != 3) {
430                 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
431                 return 1;
432         }
433         struct aData aStorage;
434         aStorage.di = &di;
435         aStorage.nbupoles = Draw::Atoi(argv[1]); 
436         aStorage.filename = argv[2];
437
438         OSD_Thread aThread1(convert);
439         aThread1.Run(&aStorage);
440         GeomConvertTest(di,aStorage.nbupoles,aStorage.filename);
441         cout << "result of thread: " << aThread1.Wait() << endl;
442
443         return 0;
444 }
445
446 #include <GeomInt_IntSS.hxx>
447 static void GeomIntSSTest (Draw_Interpretor& di, Standard_Integer theNbSol, Standard_CString theFileName1, Standard_CString theFileName2)
448 {
449         Handle(Geom_Surface) aSurf1 = DrawTrSurf::GetSurface(theFileName1);
450         Handle(Geom_Surface) aSurf2 = DrawTrSurf::GetSurface(theFileName2);
451         GeomInt_IntSS anInter;
452         anInter.Perform(aSurf1, aSurf2, Precision::Confusion(), Standard_True);
453         if (!anInter.IsDone()) {
454                 di << "An intersection is not done!" << "\n";
455                 return;
456         }
457
458         di << "Number of Lines:" << anInter.NbLines() << "\n";
459         QCOMPARE (anInter.NbLines(), theNbSol);
460 }
461
462 struct aNewData {
463         Draw_Interpretor* di;
464         Standard_Integer nbsol;
465         Standard_CString filename1;
466         Standard_CString filename2;
467 };
468 Standard_EXPORT Standard_Address convert_inter(Standard_Address data)
469 {
470         aNewData* info = (aNewData*) data;
471         GeomIntSSTest(*(info->di),info->nbsol,info->filename1,info->filename2);
472         return NULL;
473 }
474
475 static Standard_Integer OCC23952intersect (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
476 {
477         if (argc != 4) {
478                 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
479                 return 1;
480         }
481         struct aNewData aStorage;
482         aStorage.di = &di;
483         aStorage.nbsol = Draw::Atoi(argv[1]); 
484         aStorage.filename1 = argv[2];
485         aStorage.filename2 = argv[3];
486
487         OSD_Thread aThread1(convert_inter);
488         aThread1.Run(&aStorage);
489         GeomIntSSTest(di,aStorage.nbsol,aStorage.filename1,aStorage.filename2);
490         cout << "result of thread: " << aThread1.Wait() << endl;
491
492         return 0;
493 }
494
495 #include <Geom_SurfaceOfRevolution.hxx> 
496 static Standard_Integer OCC23683 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
497 {
498   if (argc < 2) {
499     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
500     return 1;
501   }
502
503   Standard_Integer ucontinuity = 1;
504   Standard_Integer vcontinuity = 1;
505   Standard_Boolean iscnu = false;
506   Standard_Boolean iscnv = false;
507   
508   Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
509
510   QCOMPARE (aSurf->IsCNu (ucontinuity), iscnu);
511   QCOMPARE (aSurf->IsCNv (vcontinuity), iscnv);
512
513   return 0;
514 }
515
516 #include <gp_Ax1.hxx>
517 #include <gp_Ax22d.hxx>
518 #include <Geom_Plane.hxx>
519 #include <Geom2d_Circle.hxx>
520 #include <Geom2d_TrimmedCurve.hxx>
521 #include <BRepBuilderAPI_MakeEdge.hxx>
522 #include <BRepPrimAPI_MakeRevol.hxx>
523 #include <Geom2d_OffsetCurve.hxx>
524
525 static int test_offset(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
526 {
527   // Check the command arguments
528   if ( argc != 1 )
529   {
530     di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
531     di << "Usage: type help " << argv[0] << "\n";
532     return 1; // TCL_ERROR
533   }
534
535   gp_Ax1 RotoAx( gp::Origin(), gp::DZ() );
536   gp_Ax22d Ax2( gp::Origin2d(), gp::DY2d(), gp::DX2d() );
537   Handle(Geom_Surface) Plane = new Geom_Plane( gp::YOZ() );
538
539   di << "<<<< Preparing sample surface of revolution based on trimmed curve >>>>" << "\n";
540   di << "-----------------------------------------------------------------------" << "\n";
541
542   Handle(Geom2d_Circle) C2d1 = new Geom2d_Circle(Ax2, 1.0);
543   Handle(Geom2d_TrimmedCurve) C2d1Trimmed = new Geom2d_TrimmedCurve(C2d1, 0.0, M_PI/2.0);
544   TopoDS_Edge E1 = BRepBuilderAPI_MakeEdge(C2d1Trimmed, Plane);
545
546   DBRep::Set("e1", E1);
547
548   BRepPrimAPI_MakeRevol aRevolBuilder1(E1, RotoAx);
549   TopoDS_Face F1 = TopoDS::Face( aRevolBuilder1.Shape() );
550
551   DBRep::Set("f1", F1);
552
553   di << "Result: f1" << "\n";
554
555   di << "<<<< Preparing sample surface of revolution based on offset curve  >>>>" << "\n";
556   di << "-----------------------------------------------------------------------" << "\n";
557
558   Handle(Geom2d_OffsetCurve) C2d2Offset = new Geom2d_OffsetCurve(C2d1Trimmed, -0.5);
559   TopoDS_Edge E2 = BRepBuilderAPI_MakeEdge(C2d2Offset, Plane);
560
561   DBRep::Set("e2", E2);
562
563   BRepPrimAPI_MakeRevol aRevolBuilder2(E2, RotoAx);
564   TopoDS_Face F2 = TopoDS::Face( aRevolBuilder2.Shape() );
565
566   DBRep::Set("f2", F2);
567
568   di << "Result: f2" << "\n";
569
570   return 0;
571 }
572
573 #include <Geom_Curve.hxx>
574 #include <Geom_Surface.hxx>
575 #include <Precision.hxx>
576 #include <ShapeConstruct_ProjectCurveOnSurface.hxx>
577 //=======================================================================
578 //function : OCC24008
579 //purpose  : 
580 //=======================================================================
581 static Standard_Integer OCC24008 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
582 {
583   if (argc != 3) {
584     di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
585     return 1;
586   }
587   Handle(Geom_Curve) aCurve = DrawTrSurf::GetCurve(argv[1]);
588   Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[2]);
589   if (aCurve.IsNull()) {
590     di << "Curve was not read" << "\n";
591         return 1;
592   }
593   if (aSurf.IsNull()) {
594         di << "Surface was not read" << "\n";
595         return 1;
596   }
597   ShapeConstruct_ProjectCurveOnSurface aProj;
598   aProj.Init (aSurf, Precision::Confusion());
599   try {
600     Handle(Geom2d_Curve) aPCurve;
601     aProj.Perform (aCurve, aCurve->FirstParameter(), aCurve->LastParameter(), aPCurve);
602     if (aPCurve.IsNull()) {
603           di << "PCurve was not created" << "\n";
604           return 1;
605     }
606   } catch (...) {
607     di << "Exception was caught" << "\n";
608   }
609   return 0;
610 }
611
612 #include <GeomAdaptor_Surface.hxx>
613 #include <Draw.hxx>
614 //=======================================================================
615 //function : OCC23945
616 //purpose  : 
617 //=======================================================================
618
619 static Standard_Integer OCC23945 (Draw_Interpretor& /*di*/,Standard_Integer n, const char** a)
620 {
621   if (n < 5) return 1;
622
623   Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(a[1]);
624   if (aS.IsNull()) return 1;
625
626   GeomAdaptor_Surface GS(aS);
627
628   Standard_Real U = Draw::Atof(a[2]);
629   Standard_Real V = Draw::Atof(a[3]);
630
631   Standard_Boolean DrawPoint = ( n%3 == 2);
632   if ( DrawPoint) n--;
633
634   gp_Pnt P;
635   if (n >= 13) {
636     gp_Vec DU,DV;
637     if (n >= 22) {
638       gp_Vec D2U,D2V,D2UV;
639       GS.D2(U,V,P,DU,DV,D2U,D2V,D2UV);
640       Draw::Set(a[13],D2U.X());
641       Draw::Set(a[14],D2U.Y());
642       Draw::Set(a[15],D2U.Z());
643       Draw::Set(a[16],D2V.X());
644       Draw::Set(a[17],D2V.Y());
645       Draw::Set(a[18],D2V.Z());
646       Draw::Set(a[19],D2UV.X());
647       Draw::Set(a[20],D2UV.Y());
648       Draw::Set(a[21],D2UV.Z());
649     }
650     else
651       GS.D1(U,V,P,DU,DV);
652
653     Draw::Set(a[7],DU.X());
654     Draw::Set(a[8],DU.Y());
655     Draw::Set(a[9],DU.Z());
656     Draw::Set(a[10],DV.X());
657     Draw::Set(a[11],DV.Y());
658     Draw::Set(a[12],DV.Z());
659   }
660   else 
661     GS.D0(U,V,P);
662
663   if ( n > 6) {
664     Draw::Set(a[4],P.X());
665     Draw::Set(a[5],P.Y());
666     Draw::Set(a[6],P.Z());
667   }
668   if ( DrawPoint) {
669     DrawTrSurf::Set(a[n],P);
670   }
671
672   return 0;
673 }
674
675 #include <Voxel_BoolDS.hxx>
676 #include <Voxel_FastConverter.hxx>
677 #include <Voxel_BooleanOperation.hxx>
678 static Standard_Integer OCC24019 (Draw_Interpretor& di, Standard_Integer argc, const char** argv) 
679 {
680   if ( argc != 2 ) {
681     di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
682         return 1;
683   }
684
685   TCollection_AsciiString aFileName = argv[1];
686   TopoDS_Shape aShape;
687   BRep_Builder aBuilder;
688
689   if (!BRepTools::Read(aShape, aFileName.ToCString(), aBuilder)) {
690     di << "Error: Could not read a shape!" << "\n";
691         return 1;
692   }
693   
694   TopoDS_Solid aShape1 = BRepPrimAPI_MakeSphere(gp_Pnt(20,25,35), 7);
695
696   Standard_Real deflection = 0.005;
697   Standard_Integer nbThreads = 1;
698   Standard_Integer nbx = 200, nby = 200, nbz = 200;
699   Voxel_BoolDS theVoxels(0,0,0, 50, 50, 50, nbx, nby, nbz);
700   Voxel_BoolDS theVoxels1(0,0,0, 50, 50, 50, nbx, nby, nbz);
701
702   Standard_Integer progress = 0;
703   Voxel_FastConverter fcp(aShape, theVoxels, deflection, nbx, nby, nbz, nbThreads);
704   fcp.ConvertUsingSAT(progress, 1);
705   fcp.FillInVolume(1);
706
707   Voxel_FastConverter fcp1(aShape1, theVoxels1, deflection, nbx, nby, nbz, nbThreads);
708   fcp1.ConvertUsingSAT(progress, 1);
709   fcp1.FillInVolume(1);
710
711   Voxel_BooleanOperation op;
712   Standard_Boolean result = op.Cut(theVoxels1, theVoxels);
713   if ( result != 1 ) {
714     di << "Error: invalid boolean operation" << "\n";
715   } else {
716         di << "OK: boolean operation is ok" << "\n";
717   }
718
719   return 0;
720 }
721
722 //=======================================================================
723 //function : OCC11758
724 //purpose  : 
725 //=======================================================================
726 static Standard_Integer OCC11758 (Draw_Interpretor& di, Standard_Integer n, const char**)
727 {
728   if (n != 1) return 1;
729
730   const char* theStr = "0123456789";
731   Standard_Integer i, j;
732   for ( i = 0; i < 5; ++i ) {
733     // TCollection_AsciiString(const Standard_CString astring)
734     TCollection_AsciiString a(theStr+i);
735     // IsEqual (const Standard_CString other)const
736     //assert( a == theStr+i );
737     QCOMPARE ( a , theStr+i );
738
739     //TCollection_AsciiString(const Standard_CString astring,const Standard_Integer aLen )
740     TCollection_AsciiString b(theStr+i, 3);
741     //assert( b.Length() == 3 );
742     //assert( strncmp( b.ToCString(), theStr+i, 3 ) == 0 );
743     //assert( strlen( b.ToCString() ) == 3 );
744     QCOMPARE ( b.Length() , 3 );
745     QCOMPARE ( strncmp( b.ToCString() , theStr+i, 3 ) , 0 );
746     QCOMPARE ( b.Length() , 3 );
747
748     //TCollection_AsciiString(const Standard_Integer aValue)
749     TCollection_AsciiString c(i);
750     //assert( c.IsIntegerValue() );
751     //assert( c.IntegerValue() == i );
752     QCOMPARE ( c.IsIntegerValue() , Standard_True );
753     QCOMPARE ( c.IntegerValue() , i );
754
755     //TCollection_AsciiString(const Standard_Real aValue)
756     TCollection_AsciiString d( 0.1*i );
757     //assert( d.IsRealValue() );
758     //assert( TCollection_AsciiString(3.3) == "3.3");
759     QCOMPARE ( d.IsRealValue() , Standard_True );
760     QCOMPARE ( TCollection_AsciiString(3.3) , "3.3" );
761
762     //TCollection_AsciiString(const TCollection_AsciiString& astring)
763     TCollection_AsciiString e(d);
764     //assert( e == d );
765     //assert( e.Length() == d.Length() );
766     //assert( strcmp( e.ToCString(), d.ToCString() ) == 0 );
767     QCOMPARE ( e ,d  );
768     QCOMPARE ( e.Length() , d.Length() );
769     QCOMPARE ( strcmp( e.ToCString(), d.ToCString() ) , 0 );
770
771     // TCollection_AsciiString(const TCollection_AsciiString& astring ,
772     //                         const Standard_Character other )
773     TCollection_AsciiString f(e,'\a');
774     //assert( f.Length() == e.Length() + 1 );
775     //assert( strncmp( f.ToCString(), e.ToCString(), e.Length() ) == 0 );
776     //assert( f.Value( f.Length() ) == '\a');
777     QCOMPARE ( f.Length() , e.Length() + 1 );
778     QCOMPARE ( strncmp( f.ToCString(), e.ToCString(), e.Length() ) , 0 );
779     QCOMPARE ( f.Value( f.Length() ) , '\a' );
780
781     // TCollection_AsciiString(const TCollection_AsciiString& astring ,
782     //                         const Standard_CString other )
783     TCollection_AsciiString g(f, theStr);
784     //assert( g.Length() == f.Length() + strlen( theStr ));
785     //assert( strncmp( g.ToCString(), f.ToCString(), f.Length() ) == 0 );
786     //assert( g.Search( theStr ) == f.Length() + 1 );
787     QCOMPARE ( g.Length() , f.Length() + (Standard_Integer)strlen( theStr ) );
788     QCOMPARE ( strncmp( g.ToCString(), f.ToCString(), f.Length() ) , 0 );
789     QCOMPARE ( g.Search( theStr ) , f.Length() + 1 );
790
791     // TCollection_AsciiString(const TCollection_AsciiString& astring ,
792     //                         const TCollection_AsciiString& other )
793     TCollection_AsciiString h(d,a);
794     //assert( h.Length() == d.Length() + a.Length() );
795     //assert( strncmp( h.ToCString(), d.ToCString(), d.Length() ) == 0 );
796     //assert( strncmp( h.ToCString() + d.Length(), a.ToCString(), a.Length() ) == 0 );
797     QCOMPARE ( h.Length() , d.Length() + a.Length() );
798     QCOMPARE ( strncmp( h.ToCString(), d.ToCString(), d.Length() ) , 0 );
799     QCOMPARE ( strncmp( h.ToCString() + d.Length(), a.ToCString(), a.Length() ) , 0 );
800
801     // AssignCat(const Standard_CString other)
802     c.AssignCat( a.ToCString() );
803     //assert( c.Length() == 1 + a.Length() );
804     //assert( c.Search( a ) == 2 );
805     QCOMPARE ( c.Length() , 1 + a.Length() );
806     QCOMPARE ( c.Search( a ) , 2 );
807
808     // AssignCat(const TCollection_AsciiString& other)
809     Standard_Integer dl = d.Length();
810     d.AssignCat( a );
811     //assert( d.Length() == dl + a.Length() );
812     //assert( d.Search( a ) == dl + 1 );
813     QCOMPARE ( d.Length() , dl + a.Length() );
814     QCOMPARE ( d.Search( a ) , dl + 1 );
815
816     // Capitalize()
817     TCollection_AsciiString capitalize("aBC");
818     capitalize.Capitalize();
819     //assert( capitalize == "Abc" );
820     QCOMPARE ( capitalize , "Abc" );
821
822     // Copy(const Standard_CString fromwhere)
823     d = theStr+i;
824     //assert( d == theStr+i );
825     QCOMPARE ( d , theStr+i );
826
827     // Copy(const TCollection_AsciiString& fromwhere)
828     d = h;
829     // IsEqual (const TCollection_AsciiString& other)const
830     //assert( d == h );
831     QCOMPARE ( d , h );
832
833     // Insert(const Standard_Integer where, const Standard_CString what)
834     dl = d.Length();
835     d.Insert( 2, theStr );
836     //assert( d.Length() == dl + strlen( theStr ));
837     //assert( strncmp( d.ToCString() + 1, theStr, strlen( theStr )) == 0 );
838     QCOMPARE ( d.Length() , dl + (Standard_Integer)strlen( theStr ) );
839     QCOMPARE ( strncmp( d.ToCString() + 1, theStr, strlen( theStr )) , 0 );
840
841     //Insert(const Standard_Integer where,const Standard_Character what)
842     d = theStr;
843     d.Insert( i+1, 'i' );
844     //assert( d.Length() == strlen( theStr ) + 1 );
845     //assert( d.Value( i+1 ) == 'i');
846     //assert( strcmp( d.ToCString() + i + 1, theStr+i ) == 0 );
847     QCOMPARE ( d.Length() , (Standard_Integer)strlen( theStr ) + 1 );
848     QCOMPARE ( d.Value( i+1 ) , 'i' );
849     QCOMPARE ( strcmp( d.ToCString() + i + 1, theStr+i ) , 0 );
850
851     //Insert(const Standard_Integer where,const TCollection_AsciiString& what)
852     d = theStr;
853     d.Insert( i+1, TCollection_AsciiString( "i" ));
854     //assert( d.Length() == strlen( theStr ) + 1 );
855     //assert( d.Value( i+1 ) == 'i');
856     //assert( strcmp( d.ToCString() + i + 1, theStr+i ) == 0 );
857     QCOMPARE ( d.Length() , (Standard_Integer)strlen( theStr ) + 1 );
858     QCOMPARE ( d.Value( i+1 ) , 'i' );
859     QCOMPARE ( strcmp( d.ToCString() + i + 1, theStr+i ) , 0 );
860
861     // IsDifferent (const Standard_CString other)const
862     //assert( d.IsDifferent( theStr ));
863     //assert( d.IsDifferent( "theStr" ));
864     //assert( d.IsDifferent( "" ));
865     //assert( !d.IsDifferent( d.ToCString() ));
866     QCOMPARE ( d.IsDifferent( theStr ) , Standard_True );
867     QCOMPARE ( d.IsDifferent( "theStr" ) , Standard_True );
868     QCOMPARE ( d.IsDifferent( "" ) , Standard_True );
869     QCOMPARE ( !d.IsDifferent( d.ToCString() ) , Standard_True );
870
871     // IsDifferent (const TCollection_AsciiString& other)const
872     //assert( d.IsDifferent( TCollection_AsciiString() ));
873     //assert( d.IsDifferent( a ));
874     //assert( d.IsDifferent( h ));
875     //assert( !d.IsDifferent( d ));
876     QCOMPARE ( d.IsDifferent( TCollection_AsciiString() ) , Standard_True );
877     QCOMPARE ( d.IsDifferent( a ) , Standard_True );
878     QCOMPARE ( d.IsDifferent( h ) , Standard_True );
879     QCOMPARE ( !d.IsDifferent( d ) , Standard_True );
880
881     // IsLess (const Standard_CString other)const
882     //assert( TCollection_AsciiString ("0"). IsLess("1"));
883     //assert( TCollection_AsciiString ("0"). IsLess("00"));
884     //assert( TCollection_AsciiString ("").  IsLess("0"));
885     //assert( !TCollection_AsciiString("1"). IsLess("0"));
886     //assert( !TCollection_AsciiString("00").IsLess("0"));
887     //assert( !TCollection_AsciiString("0"). IsLess(""));
888     //assert( TCollection_AsciiString (theStr+i).IsLess(theStr+i+1));
889     QCOMPARE ( TCollection_AsciiString ("0"). IsLess("1") , Standard_True );
890     QCOMPARE ( TCollection_AsciiString ("0"). IsLess("00") , Standard_True );
891     QCOMPARE ( TCollection_AsciiString ("").  IsLess("0") , Standard_True );
892     QCOMPARE ( !TCollection_AsciiString("1"). IsLess("0"), Standard_True );
893     QCOMPARE ( !TCollection_AsciiString("00").IsLess("0") , Standard_True );
894     QCOMPARE ( !TCollection_AsciiString("0"). IsLess("") , Standard_True );
895     QCOMPARE ( TCollection_AsciiString (theStr+i).IsLess(theStr+i+1) , Standard_True );
896
897     // IsLess (const TCollection_AsciiString& other)const
898     //assert( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("1" )));
899     //assert( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("00")));
900     //assert( TCollection_AsciiString ("").  IsLess(TCollection_AsciiString("0" )));
901     //assert( !TCollection_AsciiString("1"). IsLess(TCollection_AsciiString("0" )));
902     //assert( !TCollection_AsciiString("00").IsLess(TCollection_AsciiString("0" )));
903     //assert( !TCollection_AsciiString("0"). IsLess(TCollection_AsciiString(""  )));
904     //assert( TCollection_AsciiString (theStr+i).IsLess(TCollection_AsciiString(theStr+i+1)));
905     QCOMPARE ( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("1" )) , Standard_True );
906     QCOMPARE ( TCollection_AsciiString ("0"). IsLess(TCollection_AsciiString("00")) , Standard_True );
907     QCOMPARE ( TCollection_AsciiString ("").  IsLess(TCollection_AsciiString("0" )) , Standard_True );
908     QCOMPARE ( !TCollection_AsciiString("1"). IsLess(TCollection_AsciiString("0" )) , Standard_True );
909     QCOMPARE ( !TCollection_AsciiString("00").IsLess(TCollection_AsciiString("0" )) , Standard_True );
910     QCOMPARE ( !TCollection_AsciiString("0"). IsLess(TCollection_AsciiString(""  )) , Standard_True );
911     QCOMPARE ( TCollection_AsciiString (theStr+i).IsLess(TCollection_AsciiString(theStr+i+1)) , Standard_True );
912
913     // IsGreater (const Standard_CString other)const
914     //assert( !TCollection_AsciiString("0"). IsGreater("1"));
915     //assert( !TCollection_AsciiString("0"). IsGreater("00"));
916     //assert( !TCollection_AsciiString("").  IsGreater("0"));
917     //assert( TCollection_AsciiString ("1"). IsGreater("0"));
918     //assert( TCollection_AsciiString ("00").IsGreater("0"));
919     //assert( TCollection_AsciiString ("0"). IsGreater(""));
920     //assert( TCollection_AsciiString (theStr+i+1).IsGreater(theStr+i));
921     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater("1") , Standard_True );
922     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater("00") , Standard_True );
923     QCOMPARE ( !TCollection_AsciiString("").  IsGreater("0") , Standard_True );
924     QCOMPARE ( TCollection_AsciiString ("1"). IsGreater("0") , Standard_True );
925     QCOMPARE ( TCollection_AsciiString ("00").IsGreater("0") , Standard_True );
926     QCOMPARE ( TCollection_AsciiString ("0"). IsGreater("") , Standard_True );
927     QCOMPARE ( TCollection_AsciiString (theStr+i+1).IsGreater(theStr+i) , Standard_True );
928
929     // IsGreater (const TCollection_AsciiString& other)const
930     //assert( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("1" )));
931     //assert( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("00")));
932     //assert( !TCollection_AsciiString("").  IsGreater(TCollection_AsciiString("0" )));
933     //assert( TCollection_AsciiString ("1"). IsGreater(TCollection_AsciiString("0" )));
934     //assert( TCollection_AsciiString ("00").IsGreater(TCollection_AsciiString("0" )));
935     //assert( TCollection_AsciiString ("0"). IsGreater(TCollection_AsciiString(""  )));
936     //assert( TCollection_AsciiString (theStr+i+1).IsGreater(TCollection_AsciiString(theStr+i)));
937     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("1" )) , Standard_True );
938     QCOMPARE ( !TCollection_AsciiString("0"). IsGreater(TCollection_AsciiString("00")) , Standard_True );
939     QCOMPARE ( !TCollection_AsciiString("").  IsGreater(TCollection_AsciiString("0" )) , Standard_True );
940     QCOMPARE ( TCollection_AsciiString ("1"). IsGreater(TCollection_AsciiString("0" )) , Standard_True );
941     QCOMPARE ( TCollection_AsciiString ("00").IsGreater(TCollection_AsciiString("0" )) , Standard_True );
942     QCOMPARE ( TCollection_AsciiString ("0"). IsGreater(TCollection_AsciiString(""  )) , Standard_True );
943     QCOMPARE ( TCollection_AsciiString (theStr+i+1).IsGreater(TCollection_AsciiString(theStr+i)) , Standard_True );
944
945     // void Read(Standard_IStream& astream)
946     std::istringstream is( theStr );
947     e.Read( is );
948     //assert( e == theStr );
949     QCOMPARE ( e , theStr );
950
951     // Standard_Integer SearchFromEnd (const Standard_CString what)const
952     //assert( e.SearchFromEnd( theStr + i ) == i + 1 );
953     QCOMPARE ( e.SearchFromEnd( theStr + i ) , i + 1 );
954
955     // SetValue(const Standard_Integer where, const Standard_CString what)
956     e.SetValue( i+1, "what");
957     //assert( e.Search( "what" ) == i+1 );
958     //assert( e.Length() == strlen( theStr ));
959     QCOMPARE ( e.Search( "what" ) , i+1 );
960     QCOMPARE ( e.Length() , (Standard_Integer)strlen( theStr ) );
961
962     // TCollection_AsciiString Split (const Standard_Integer where)
963     e = theStr;
964     d = e.Split( i+1 );
965     //assert( d.Length() + e.Length() == strlen( theStr ));
966     QCOMPARE ( d.Length() + e.Length() , (Standard_Integer)strlen( theStr ) );
967
968     // TCollection_AsciiString SubString (const Standard_Integer FromIndex,
969     //                                    const Standard_Integer ToIndex) const
970     e = theStr;
971     d = e.SubString( (unsigned int)i+1, (unsigned int)i+3 );
972     //assert( d.Length() == 3 );
973     //assert( d.Value(1) == theStr[ i ]);
974     QCOMPARE ( d.Length() , 3 );
975     QCOMPARE ( d.Value(1) , theStr[ i ] );
976
977     // TCollection_AsciiString Token (const Standard_CString separators,
978     //                                const Standard_Integer whichone) const
979     e = " ";
980     for ( j = 0; j < i; ++j ) {
981       e += TCollection_AsciiString( theStr[j] ) + " ";
982       //assert( e.Token(" ", j+1 ) == TCollection_AsciiString( theStr+j, 1 ));
983       QCOMPARE ( e.Token(" ", j+1 ) , TCollection_AsciiString( theStr+j, 1 ) );
984     }
985   }
986   for ( i = 0; i < 5; ++i )
987   {
988     // TCollection_ExtendedString (const Standard_CString astring, 
989     //                             const Standard_Boolean isMultiByte) 
990     const TCollection_ExtendedString a( theStr+i );
991     //assert( TCollection_AsciiString( a ) == theStr+i );
992     QCOMPARE ( TCollection_AsciiString( a ) , theStr+i );
993
994     //TCollection_ExtendedString (const Standard_ExtString astring)
995     const TCollection_ExtendedString b( a.ToExtString() );
996     //assert( a == b );
997     QCOMPARE ( a , b );
998
999     // TCollection_ExtendedString (const Standard_Integer      length,
1000     //                             const Standard_ExtCharacter filler )
1001     const TCollection_ExtendedString c( i, 1 );
1002     //assert( c.Length() == i );
1003     QCOMPARE ( c.Length() , i );
1004     if ( c.Length() > 0 ) {
1005       //assert( c.Value( i ) == 1 );
1006       QCOMPARE ( c.Value( i ) , 1 );
1007     }
1008
1009     // TCollection_ExtendedString (const Standard_Integer aValue)
1010     TCollection_ExtendedString d( i );
1011     const TCollection_AsciiString da( d );
1012     //assert( da.IsIntegerValue() );
1013     //assert( da.IntegerValue() == i );
1014     QCOMPARE ( da.IsIntegerValue() , Standard_True );
1015     QCOMPARE (  da.IntegerValue(), i );
1016
1017     // TCollection_ExtendedString (const Standard_Real aValue)
1018     const TCollection_ExtendedString e( 0.1 * i );
1019     const TCollection_AsciiString ea( e );
1020     //assert( ea.IsRealValue() );
1021     //assert( Abs( ea.RealValue() - 0.1 * i ) < 1e-10 );
1022     QCOMPARE ( ea.IsRealValue() , Standard_True );
1023     QCOMPARE ( Abs( ea.RealValue() - 0.1 * i ) < 1e-10 , Standard_True );
1024
1025     // TCollection_ExtendedString (const TCollection_ExtendedString& astring)
1026     const TCollection_ExtendedString f(e);
1027     //assert( f.Length() == e.Length());
1028     //assert( f == e );
1029     QCOMPARE ( f.Length() , e.Length() );
1030     QCOMPARE ( f , e );
1031
1032     // TCollection_ExtendedString (const TCollection_AsciiString& astring)
1033     const TCollection_ExtendedString g( ea );
1034     //assert( g.Length() == ea.Length() );
1035     //assert( TCollection_AsciiString( g ) == ea );
1036     QCOMPARE ( g.Length() , ea.Length() );
1037     QCOMPARE ( TCollection_AsciiString( g ) , ea );
1038
1039     // AssignCat (const TCollection_ExtendedString& other)
1040     const TCollection_ExtendedString sep(",");
1041     d.AssignCat( sep );
1042     d.AssignCat( g );
1043     //assert( d.Length() == 2 + g.Length() );
1044     //assert( d.Token( sep.ToExtString(), 1 ) == TCollection_ExtendedString( i ));
1045     //assert( d.Token( sep.ToExtString(), 2 ) == g );
1046     QCOMPARE ( d.Length() , 2 + g.Length() );
1047     QCOMPARE ( d.Token( sep.ToExtString(), 1 ) , TCollection_ExtendedString( i ) );
1048     QCOMPARE ( d.Token( sep.ToExtString(), 2 ) , g );
1049
1050     // TCollection_ExtendedString Cat (const TCollection_ExtendedString& other) const
1051     const TCollection_ExtendedString cat = a.Cat( sep );
1052     //assert( cat.Length() == a.Length() + sep.Length() );
1053     //assert( cat.Search( a ) == 1 );
1054     //assert( cat.Search( sep ) == a.Length() + 1 );
1055     QCOMPARE ( cat.Length() , a.Length() + sep.Length() );
1056     QCOMPARE ( cat.Search( a ) , 1 );
1057     QCOMPARE ( cat.Search( sep ) , a.Length() + 1 );
1058
1059     // Copy (const TCollection_ExtendedString& fromwhere)
1060     d = cat;
1061     //assert( d.Length() == cat.Length() );
1062     //assert( d == cat );
1063     QCOMPARE ( d.Length() , cat.Length() );
1064     QCOMPARE ( d , cat );
1065
1066     // IsEqual (const Standard_ExtString other) const
1067     //assert( d.IsEqual( d.ToExtString() ));
1068     QCOMPARE ( d.IsEqual( d.ToExtString() ) , Standard_True );
1069
1070     // IsDifferent (const Standard_ExtString other ) const
1071     //assert( d.IsDifferent( a.ToExtString() ));
1072     QCOMPARE ( d.IsDifferent( a.ToExtString() ) , Standard_True );
1073
1074     // IsDifferent (const TCollection_ExtendedString& other) const
1075     //assert( d.IsDifferent( a ));
1076     QCOMPARE ( d.IsDifferent( a ) , Standard_True );
1077
1078     // IsLess (const Standard_ExtString other) const
1079     const TCollection_ExtendedString l0("0"), l1("1"), l00("00"), l, ls(theStr+i), ls1(theStr+i+1);
1080     //assert( l0. IsLess( l1.ToExtString() ));
1081     //assert( l0. IsLess( l00.ToExtString() ));
1082     //assert( l.  IsLess( l0.ToExtString() ));
1083     //assert( ! l1. IsLess( l0.ToExtString() ));
1084     //assert( ! l00.IsLess( l0.ToExtString() ));
1085     //assert( ! l0. IsLess( l.ToExtString() ));
1086     //assert( ls.IsLess( ls1.ToExtString() ));
1087     QCOMPARE ( l0. IsLess( l1.ToExtString() ) , Standard_True );
1088     QCOMPARE ( l0. IsLess( l00.ToExtString() ) , Standard_True );
1089     QCOMPARE ( l.  IsLess( l0.ToExtString() ) , Standard_True );
1090     QCOMPARE ( ! l1. IsLess( l0.ToExtString() ) , Standard_True );
1091     QCOMPARE ( ! l00.IsLess( l0.ToExtString() ) , Standard_True );
1092     QCOMPARE ( ! l0. IsLess( l.ToExtString() ) , Standard_True );
1093     QCOMPARE ( ls.IsLess( ls1.ToExtString() ) , Standard_True );
1094
1095     // IsLess (const TCollection_ExtendedString& other) const
1096     //assert( l0. IsLess( l1 ));
1097     //assert( l0. IsLess( l00 ));
1098     //assert( l.  IsLess( l0 ));
1099     //assert( ! l1. IsLess( l0 ));
1100     //assert( ! l00.IsLess( l0 ));
1101     //assert( ! l0. IsLess( l ));
1102     //assert( ls.IsLess( ls1 ));
1103     QCOMPARE ( l0. IsLess( l1 ) , Standard_True );
1104     QCOMPARE ( l0. IsLess( l00 ) , Standard_True );
1105     QCOMPARE ( l.  IsLess( l0 ) , Standard_True );
1106     QCOMPARE ( ! l1. IsLess( l0 ) , Standard_True );
1107     QCOMPARE ( ! l00.IsLess( l0 ) , Standard_True );
1108     QCOMPARE ( ! l0. IsLess( l ) , Standard_True );
1109     QCOMPARE ( ls.IsLess( ls1 ) , Standard_True );
1110
1111     // IsGreater (const Standard_ExtString other) const
1112     //assert( ! l0.IsGreater( l1.ToExtString() ));
1113     //assert( ! l0.IsGreater( l00.ToExtString() ));
1114     //assert( ! l. IsGreater( l0.ToExtString() ));
1115     //assert(  l1. IsGreater( l0.ToExtString() ));
1116     //assert(  l00.IsGreater( l0.ToExtString() ));
1117     //assert(  l0. IsGreater( l.ToExtString() ));
1118     //assert(  ls1.IsGreater( ls.ToExtString() ));
1119     QCOMPARE ( ! l0.IsGreater( l1.ToExtString() ) , Standard_True );
1120     QCOMPARE ( ! l0.IsGreater( l00.ToExtString() ) , Standard_True );
1121     QCOMPARE ( ! l. IsGreater( l0.ToExtString() ) , Standard_True );
1122     QCOMPARE ( l1. IsGreater( l0.ToExtString() ) , Standard_True );
1123     QCOMPARE ( l00.IsGreater( l0.ToExtString() ) , Standard_True );
1124     QCOMPARE ( l0. IsGreater( l.ToExtString() ) , Standard_True );
1125     QCOMPARE ( ls1.IsGreater( ls.ToExtString() ) ,Standard_True  );
1126
1127     // IsGreater (const TCollection_ExtendedString& other) const
1128     //assert( ! l0.IsGreater( l1));
1129     //assert( ! l0.IsGreater( l00));
1130     //assert( ! l. IsGreater( l0));
1131     //assert(  l1. IsGreater( l0));
1132     //assert(  l00.IsGreater( l0));
1133     //assert(  l0. IsGreater( l));
1134     //assert(  ls1.IsGreater( ls));
1135     QCOMPARE ( ! l0.IsGreater( l1) , Standard_True );
1136     QCOMPARE ( ! l0.IsGreater( l00) , Standard_True );
1137     QCOMPARE ( ! l. IsGreater( l0) , Standard_True );
1138     QCOMPARE ( l1. IsGreater( l0) , Standard_True );
1139     QCOMPARE ( l00.IsGreater( l0) , Standard_True );
1140     QCOMPARE ( l0. IsGreater( l) , Standard_True );
1141     QCOMPARE ( ls1.IsGreater( ls) , Standard_True );
1142
1143     // ==========================
1144     //TCollection_HAsciiString::
1145     // ==========================
1146
1147     // IsDifferent(const Handle(TCollection_HAsciiString)& S)
1148     Handle(TCollection_HAsciiString) ha1 = new TCollection_HAsciiString( theStr+i );
1149     Handle(TCollection_HAsciiString) ha2 = new TCollection_HAsciiString( theStr+i+1 );
1150     //assert( ha1->IsDifferent( ha2 ));
1151     //assert( !ha1->IsDifferent( ha1 ));
1152     QCOMPARE ( ha1->IsDifferent( ha2 ) , Standard_True );
1153     QCOMPARE ( !ha1->IsDifferent( ha1 ) , Standard_True );
1154
1155     // IsSameString (const Handle(TCollection_HAsciiString)& S)
1156     //assert( !ha1->IsSameString( ha2 ));
1157     //assert( ha1->IsSameString( ha1 ));
1158     QCOMPARE ( !ha1->IsSameString( ha2 ) , Standard_True );
1159     QCOMPARE ( ha1->IsSameString( ha1 ) , Standard_True );
1160
1161     // IsSameState (const Handle(TCollection_HAsciiString)& other) const
1162     //assert( !ha1->IsSameState( ha2 ));
1163     //assert( ha1->IsSameState( ha1 ));
1164     QCOMPARE ( !ha1->IsSameState( ha2 ) , Standard_True );
1165     QCOMPARE ( ha1->IsSameState( ha1 ) , Standard_True );
1166
1167     // IsSameString (const Handle(TCollection_HAsciiString)& S ,
1168     //               const Standard_Boolean CaseSensitive) const
1169     //assert( !ha1->IsSameString( ha2, true ));
1170     //assert( ha1->IsSameString( ha1, true ));
1171     //assert( !ha1->IsSameString( ha2, false ));
1172     //assert( ha1->IsSameString( ha1, false ));
1173     QCOMPARE ( !ha1->IsSameString( ha2, Standard_True ) , Standard_True );
1174     QCOMPARE ( ha1->IsSameString( ha1, Standard_True ) , Standard_True );
1175     QCOMPARE ( !ha1->IsSameString( ha2, Standard_False ) , Standard_True );
1176     QCOMPARE ( ha1->IsSameString( ha1, Standard_False ) , Standard_True );
1177
1178     ha1->SetValue( 1, "AbC0000000");
1179     ha2->SetValue( 1, "aBc0000000");
1180     //assert( !ha1->IsSameString( ha2, true ));
1181     //assert( ha1->IsSameString( ha2, false ));
1182     QCOMPARE ( !ha1->IsSameString( ha2, Standard_True ) , Standard_True );
1183     QCOMPARE (  ha1->IsSameString( ha2, Standard_False ), Standard_True );
1184   }
1185   return 0;
1186 }
1187
1188 #include <Geom_CylindricalSurface.hxx>
1189 #include <IntTools_FaceFace.hxx>
1190 #include <IntTools_Curve.hxx>
1191 #include <IntTools_PntOn2Faces.hxx>
1192
1193 static Standard_Integer OCC24005 (Draw_Interpretor& theDI, Standard_Integer theNArg, const char** theArgv) 
1194 {
1195   if(theNArg < 2)
1196   {
1197     theDI << "Wrong a number of arguments!\n";
1198     return 1;
1199   }
1200
1201   Handle(Geom_Plane) plane(new Geom_Plane(
1202                                   gp_Ax3( gp_Pnt(-72.948737453424499, 754.30437716359393, 259.52151854671678),
1203                                   gp_Dir(6.2471473085930200e-007, -0.99999999999980493, 0.00000000000000000),
1204                                   gp_Dir(0.99999999999980493, 6.2471473085930200e-007, 0.00000000000000000))));
1205   Handle(Geom_CylindricalSurface) cylinder(
1206                   new Geom_CylindricalSurface(
1207                                   gp_Ax3(gp_Pnt(-6.4812490053250649, 753.39408794522092, 279.16400974257465),
1208                                   gp_Dir(1.0000000000000000, 0.0, 0.00000000000000000),
1209                                   gp_Dir(0.0, 1.0000000000000000, 0.00000000000000000)),
1210                                                                                           19.712534607908712));
1211
1212   DrawTrSurf::Set("pln", plane);
1213   theDI << "pln\n";
1214   DrawTrSurf::Set("cyl", cylinder);
1215   theDI << "cyl\n";
1216
1217   BRep_Builder builder;
1218   TopoDS_Face face1, face2;
1219   builder.MakeFace(face1, plane, Precision::Confusion());
1220   builder.MakeFace(face2, cylinder, Precision::Confusion());
1221   IntTools_FaceFace anInters;
1222   anInters.SetParameters(false, true, true, Precision::Confusion());
1223   anInters.Perform(face1, face2);
1224
1225   if (!anInters.IsDone())
1226   {
1227     theDI<<"No intersections found!"<<"\n";
1228
1229     return 1;
1230   }
1231
1232   //Handle(Geom_Curve) aResult;
1233   //gp_Pnt             aPoint;
1234
1235   const IntTools_SequenceOfCurves& aCvsX=anInters.Lines();
1236   const IntTools_SequenceOfPntOn2Faces& aPntsX=anInters.Points();
1237
1238   char buf[1024];  
1239   Standard_Integer aNbCurves, aNbPoints;
1240
1241   aNbCurves=aCvsX.Length();
1242   aNbPoints=aPntsX.Length();
1243
1244   if (aNbCurves >= 2)
1245   {
1246     for (Standard_Integer i=1; i<=aNbCurves; ++i)
1247     {
1248       Sprintf(buf, "%s_%d",theArgv[1],i);
1249       theDI << buf << " ";
1250       
1251       const IntTools_Curve& aIC = aCvsX(i);
1252       const Handle(Geom_Curve)& aC3D= aIC.Curve();
1253       DrawTrSurf::Set(buf,aC3D);
1254     }
1255   }
1256   else if (aNbCurves == 1)
1257   {
1258     const IntTools_Curve& aIC = aCvsX(1);
1259     const Handle(Geom_Curve)& aC3D= aIC.Curve();
1260     Sprintf(buf, "%s",theArgv[1]);
1261     theDI << buf << " ";
1262     DrawTrSurf::Set(buf,aC3D);
1263   }
1264
1265   for (Standard_Integer i = 1; i<=aNbPoints; ++i)
1266   {
1267     const IntTools_PntOn2Faces& aPi=aPntsX(i);
1268     const gp_Pnt& aP=aPi.P1().Pnt();
1269     
1270     Sprintf(buf,"%s_p_%d",theArgv[1],i);
1271     theDI << buf << " ";
1272     DrawTrSurf::Set(buf, aP);
1273   }
1274
1275   return 0;
1276 }
1277
1278 #include <BRepAlgo_NormalProjection.hxx>
1279 static Standard_Integer OCC24012 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv) 
1280 {
1281         if (argc != 3) {
1282                 di << "Usage : " << argv[0] << " should be 2 arguments (face and edge)";
1283                 return 1;
1284         }
1285         
1286         Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
1287         if(myAISContext.IsNull()) {
1288                 di << "use 'vinit' command before " << argv[0] << "\n";
1289                 return 1;
1290         }
1291
1292         TopoDS_Face m_Face1 = TopoDS::Face(DBRep::Get(argv[1]));
1293         TopoDS_Edge m_Edge = TopoDS::Edge(DBRep::Get(argv[2]));
1294         
1295         BRepAlgo_NormalProjection anormpro(m_Face1);
1296     anormpro.Add(m_Edge);
1297     anormpro.SetDefaultParams();
1298
1299     //anormpro.Compute3d();
1300     //anormpro.SetLimit();
1301
1302     anormpro.Build();
1303
1304     if (anormpro.IsDone())
1305     {
1306         TopoDS_Shape rshape = anormpro.Projection();
1307                 Handle(AIS_InteractiveObject) myShape = new AIS_Shape (rshape);
1308                 myAISContext->SetColor(myShape, Quantity_Color(Quantity_NOC_YELLOW));
1309                 myAISContext->Display(myShape, Standard_True);
1310     }
1311
1312         return 0;
1313 }
1314
1315 #include <Voxel_FastConverter.hxx>
1316 static Standard_Integer OCC24051 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv) 
1317 {
1318         if (argc != 1) {
1319                 di << "Usage : " << argv[0] << " should be one argument (command name only)";
1320                 return 1;
1321         }
1322
1323         TopoDS_Shape shape = BRepPrimAPI_MakeBox(gp_Pnt(5, 10, 10), 10, 20, 30).Shape();
1324         Standard_Integer progress = 0;
1325         Standard_Real deflection = 0.005;
1326         Standard_Integer nbx = 200, nby = 200, nbz = 200;
1327         Voxel_BoolDS theVoxels(-50,-50,-30, 100, 100, 100, nbx, nby, nbz);
1328         Voxel_BoolDS theVoxels1(-50,-50,-30, 100, 100, 100, nbx, nby, nbz);
1329         Standard_Integer nbThreads = 5;
1330         Voxel_FastConverter fcp(shape, theVoxels, deflection, nbx, nby, nbz, nbThreads, Standard_True);
1331         
1332         #ifdef WNT
1333         #pragma omp parallel for
1334         for(int i = 0; i < nbThreads; i++)
1335                         fcp.ConvertUsingSAT(progress, i+1);
1336         #endif
1337         
1338         fcp.ConvertUsingSAT(progress);
1339
1340         return 0;
1341 }
1342
1343 #include <BRepFeat_SplitShape.hxx>
1344 #include <ShapeAnalysis_ShapeContents.hxx>
1345 #include <BRepAlgo.hxx>
1346 static Standard_Integer OCC24086 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv) 
1347 {
1348         if (argc != 3) {
1349                 di << "Usage : " << argv[0] << " should be 2 arguments (face and wire)";
1350                 return 1;
1351         }
1352         
1353         Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
1354         if(myAISContext.IsNull()) {
1355                 di << "use 'vinit' command before " << argv[0] << "\n";
1356                 return 1;
1357         }
1358         
1359         TopoDS_Shape result;
1360         TopoDS_Face face = TopoDS::Face(DBRep::Get(argv[1]));
1361         TopoDS_Wire wire = TopoDS::Wire(DBRep::Get(argv[2]));
1362     
1363         BRepFeat_SplitShape asplit(face);
1364         asplit.Add(wire, face);
1365         asplit.Build();
1366     result = asplit.Shape();
1367     ShapeAnalysis_ShapeContents ana;
1368     ana.Perform(result);
1369     ana.NbFaces();
1370
1371         if (!(BRepAlgo::IsValid(result))) {
1372                 di << "Result was checked and it is INVALID" << "\n";
1373         } else {
1374                 di << "Result was checked and it is VALID" << "\n";
1375         }
1376         
1377         Handle(AIS_InteractiveObject) myShape = new AIS_Shape (result);
1378         myAISContext->Display(myShape, Standard_True);
1379
1380         return 0;
1381 }
1382
1383 #include <Extrema_FuncExtPS.hxx>
1384 #include <math_FunctionSetRoot.hxx>
1385 #include <math_Vector.hxx>
1386 #include <BRepBuilderAPI_MakeVertex.hxx>
1387 static Standard_Integer OCC24137 (Draw_Interpretor& theDI, Standard_Integer theNArg, const char** theArgv) 
1388 {
1389   Standard_Integer anArgIter = 1;
1390   if (theNArg < 5)
1391     {
1392       theDI <<"Usage: " << theArgv[0] << " face vertex U V [N]"<<"\n";
1393       return 1;
1394     }
1395
1396   // get target shape
1397   Standard_CString aFaceName = theArgv[anArgIter++];
1398   Standard_CString aVertName = theArgv[anArgIter++];
1399   const TopoDS_Shape aShapeF = DBRep::Get (aFaceName);
1400   const TopoDS_Shape aShapeV = DBRep::Get (aVertName);
1401   const Standard_Real aUFrom = Atof (theArgv[anArgIter++]);
1402   const Standard_Real aVFrom = Atof (theArgv[anArgIter++]);
1403   const Standard_Integer aNbIts = (anArgIter < theNArg) ? atol (theArgv[anArgIter++]) : 100;
1404   if (aShapeF.IsNull() || aShapeF.ShapeType() != TopAbs_FACE)
1405     {
1406       std::cout << "Error: " << aFaceName << " shape is null / not a face" << std::endl;
1407       return 1;
1408     }
1409   if (aShapeV.IsNull() || aShapeV.ShapeType() != TopAbs_VERTEX)
1410     {
1411       std::cout << "Error: " << aVertName << " shape is null / not a vertex" << std::endl;
1412       return 1;
1413     }
1414   const TopoDS_Face   aFace = TopoDS::Face   (aShapeF);
1415   const TopoDS_Vertex aVert = TopoDS::Vertex (aShapeV);
1416   GeomAdaptor_Surface aSurf (BRep_Tool::Surface (aFace));
1417
1418   gp_Pnt aPnt = BRep_Tool::Pnt (aVert), aRes;
1419
1420   Extrema_FuncExtPS    anExtFunc;
1421   math_FunctionSetRoot aRoot (anExtFunc, aNbIts);
1422
1423   math_Vector aTolUV (1, 2), aUVinf  (1, 2), aUVsup  (1, 2), aFromUV (1, 2);
1424   aTolUV (1) =  Precision::Confusion(); aTolUV (2) =  Precision::Confusion();
1425   aUVinf (1) = -Precision::Infinite();  aUVinf (2) = -Precision::Infinite();
1426   aUVsup (1) =  Precision::Infinite();  aUVsup (2) =  Precision::Infinite();
1427   aFromUV(1) =  aUFrom; aFromUV(2) = aVFrom;
1428
1429   anExtFunc.Initialize (aSurf);
1430   anExtFunc.SetPoint (aPnt);
1431   aRoot.SetTolerance (aTolUV);
1432   aRoot.Perform (anExtFunc, aFromUV, aUVinf, aUVsup);
1433   if (!aRoot.IsDone())
1434     {
1435       std::cerr << "No results!\n";
1436       return 1;
1437     }
1438
1439   theDI << aRoot.Root()(1) << " " << aRoot.Root()(2) << "\n";
1440   
1441   aSurf.D0 (aRoot.Root()(1), aRoot.Root()(2), aRes);
1442   DBRep::Set ("result", BRepBuilderAPI_MakeVertex (aRes));
1443   return 0;
1444 }
1445
1446 //! Check boolean operations on NCollection_Map
1447 static Standard_Integer OCC24271 (Draw_Interpretor& di,
1448                                   Standard_Integer  /*theArgNb*/,
1449                                   const char**      /*theArgVec*/)
1450 {
1451   // input data
1452   const Standard_Integer aLeftLower  = 1;
1453   const Standard_Integer aLeftUpper  = 10;
1454   const Standard_Integer aRightLower = 5;
1455   const Standard_Integer aRightUpper = 15;
1456
1457   // define arguments
1458   NCollection_Map<Standard_Integer> aMapLeft;
1459   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1460   {
1461     aMapLeft.Add (aKeyIter);
1462   }
1463
1464   NCollection_Map<Standard_Integer> aMapRight;
1465   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aRightUpper; ++aKeyIter)
1466   {
1467     aMapRight.Add (aKeyIter);
1468   }
1469
1470   QCOMPARE (aMapLeft .Contains (aMapRight), Standard_False);
1471   QCOMPARE (aMapRight.Contains (aMapLeft),  Standard_False);
1472
1473   // validate Union operation
1474   NCollection_Map<Standard_Integer> aMapUnion;
1475   aMapUnion.Union (aMapLeft, aMapRight);
1476   QCOMPARE (aMapUnion.Extent(), aRightUpper - aLeftLower + 1);
1477   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter <= aRightUpper; ++aKeyIter)
1478   {
1479     QCOMPARE (aMapUnion.Contains (aKeyIter), Standard_True);
1480   }
1481
1482   // validate Intersection operation
1483   NCollection_Map<Standard_Integer> aMapSect;
1484   aMapSect.Intersection (aMapLeft, aMapRight);
1485   QCOMPARE (aMapSect.Extent(), aLeftUpper - aRightLower + 1);
1486   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1487   {
1488     QCOMPARE (aMapSect.Contains (aKeyIter), Standard_True);
1489   }
1490   QCOMPARE (aMapLeft .Contains (aMapSect), Standard_True);
1491   QCOMPARE (aMapRight.Contains (aMapSect), Standard_True);
1492
1493   // validate Substruction operation
1494   NCollection_Map<Standard_Integer> aMapSubsLR;
1495   aMapSubsLR.Subtraction (aMapLeft, aMapRight);
1496   QCOMPARE (aMapSubsLR.Extent(), aRightLower - aLeftLower);
1497   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter < aRightLower; ++aKeyIter)
1498   {
1499     QCOMPARE (aMapSubsLR.Contains (aKeyIter), Standard_True);
1500   }
1501
1502   NCollection_Map<Standard_Integer> aMapSubsRL;
1503   aMapSubsRL.Subtraction (aMapRight, aMapLeft);
1504   QCOMPARE (aMapSubsRL.Extent(), aRightUpper - aLeftUpper);
1505   for (Standard_Integer aKeyIter = aLeftUpper + 1; aKeyIter < aRightUpper; ++aKeyIter)
1506   {
1507     QCOMPARE (aMapSubsRL.Contains (aKeyIter), Standard_True);
1508   }
1509
1510   // validate Difference operation
1511   NCollection_Map<Standard_Integer> aMapDiff;
1512   aMapDiff.Difference (aMapLeft, aMapRight);
1513   QCOMPARE (aMapDiff.Extent(), aRightLower - aLeftLower + aRightUpper - aLeftUpper);
1514   for (Standard_Integer aKeyIter = aLeftLower; aKeyIter < aRightLower; ++aKeyIter)
1515   {
1516     QCOMPARE (aMapDiff.Contains (aKeyIter), Standard_True);
1517   }
1518   for (Standard_Integer aKeyIter = aLeftUpper + 1; aKeyIter < aRightUpper; ++aKeyIter)
1519   {
1520     QCOMPARE (aMapDiff.Contains (aKeyIter), Standard_True);
1521   }
1522
1523   // validate Exchange operation
1524   NCollection_Map<Standard_Integer> aMapSwap;
1525   aMapSwap.Exchange (aMapSect);
1526   for (Standard_Integer aKeyIter = aRightLower; aKeyIter <= aLeftUpper; ++aKeyIter)
1527   {
1528     QCOMPARE (aMapSwap.Contains (aKeyIter), Standard_True);
1529   }
1530   QCOMPARE (aMapSect.IsEmpty(), Standard_True);
1531   aMapSwap.Add (34);
1532   aMapSect.Add (43);
1533
1534   NCollection_Map<Standard_Integer> aMapCopy (aMapSwap);
1535   QCOMPARE (aMapCopy.IsEqual (aMapSwap), Standard_True);
1536   aMapCopy.Remove (34);
1537   aMapCopy.Add    (43);
1538   QCOMPARE (aMapCopy.IsEqual (aMapSwap), Standard_False);
1539
1540   return 0;
1541 }
1542
1543 #define QVERIFY(val1) \
1544   di << "Checking " #val1 " == Standard_True" << \
1545         ((val1) == Standard_True ? ": OK\n" : ": Error\n")
1546
1547 #include <GeomInt_IntSS.hxx>
1548 #include <Geom_ConicalSurface.hxx>
1549 #include <Standard_ErrorHandler.hxx>
1550 //=======================================================================
1551 //function : OCC23972
1552 //purpose  : 
1553 //=======================================================================
1554 static void DoGeomIntSSTest (const Handle(Geom_Surface)& theSurf1,
1555                              const Handle(Geom_Surface)& theSurf2,
1556                              const Standard_Integer theNbSol,
1557                              Draw_Interpretor& di)
1558 {
1559   try {
1560     OCC_CATCH_SIGNALS
1561          GeomInt_IntSS anInter;
1562          anInter.Perform (theSurf1, theSurf2, Precision::Confusion(), Standard_True);
1563          QVERIFY (anInter.IsDone());
1564          QCOMPARE (anInter.NbLines(), theNbSol);
1565   } catch (...) {
1566     QVERIFY (Standard_False);
1567   }
1568 }
1569
1570 namespace {
1571   static Handle(Geom_ConicalSurface) CreateCone (const gp_Pnt& theLoc,
1572                                                  const gp_Dir& theDir,
1573                                                  const gp_Dir& theXDir,
1574                                                  const Standard_Real theRad,
1575                                                  const Standard_Real theSin,
1576                                                  const Standard_Real theCos)
1577   {
1578     const Standard_Real anA = atan (theSin / theCos);
1579     gp_Ax3 anAxis (theLoc, theDir, theXDir);
1580     Handle(Geom_ConicalSurface) aSurf = new Geom_ConicalSurface (anAxis, anA, theRad);
1581     return aSurf;
1582   }
1583 }
1584
1585 static Standard_Integer OCC23972 (Draw_Interpretor& di,Standard_Integer n, const char**)
1586 {
1587   if (n != 1) return 1;
1588
1589   //process specific cones, cannot read them from files because due to rounding the original error
1590   //in math_FunctionRoots gets hidden
1591   Handle(Geom_Surface) aS1 = CreateCone (
1592                                          gp_Pnt (123.694345356663, 789.9, 68.15),
1593                                          gp_Dir (-1, 3.48029791472957e-016, -8.41302743359754e-017),
1594                                          gp_Dir (-3.48029791472957e-016, -1, -3.17572289932207e-016),
1595                                          3.28206830417112,
1596                                          0.780868809443031,
1597                                          0.624695047554424);
1598   Handle(Geom_Surface) aS2 = CreateCone (
1599                                          gp_Pnt (123.694345356663, 784.9, 68.15),
1600                                          gp_Dir (-1, -2.5209507537117e-016, -1.49772808948866e-016),
1601                                          gp_Dir (1.49772808948866e-016, 3.17572289932207e-016, -1),
1602                                          3.28206830417112,
1603                                          0.780868809443031,
1604                                          0.624695047554424);
1605   
1606   DoGeomIntSSTest (aS1, aS2, 2, di);
1607
1608   return 0;
1609 }
1610
1611 #include <ShapeFix_EdgeProjAux.hxx>
1612 static Standard_Integer OCC24370 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
1613 {
1614   if (argc < 5) {
1615     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
1616     return 1;
1617   }
1618
1619   TopoDS_Shape aSh = DBRep::Get(argv[1]);
1620   if (aSh.IsNull()) {
1621     di << argv[0] << " Error: Null input edge\n";
1622     return 1;
1623   }
1624   const TopoDS_Edge& anEdge = TopoDS::Edge (aSh);
1625
1626   Handle(Geom2d_Curve) aC = DrawTrSurf::GetCurve2d(argv[2]);
1627   if (aC.IsNull()) {
1628     di << argv[0] << " Error: Null input curve\n";
1629     return 1;
1630   }
1631
1632   Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(argv[3]);
1633   if (aS.IsNull()) {
1634     di << argv[0] << " Error: Null input surface\n";
1635     return 1;
1636   }
1637
1638   Standard_Real prec = Draw::Atof(argv[4]);
1639   
1640   //prepare data
1641   TopoDS_Face aFace;
1642   BRep_Builder aB;
1643   aB.MakeFace (aFace, aS, Precision::Confusion());
1644   aB.UpdateEdge (anEdge, aC, aFace, Precision::Confusion());
1645   aB.Range (anEdge, aFace, aC->FirstParameter(), aC->LastParameter());
1646
1647   //call algorithm
1648   ShapeFix_EdgeProjAux aProj (aFace, anEdge);
1649   aProj.Compute (prec);
1650   
1651   Standard_Boolean isfirstdone = aProj.IsFirstDone();
1652   Standard_Boolean islastdone = aProj.IsLastDone();
1653
1654   Standard_Real first = 0.;
1655   Standard_Real last = 0.;
1656   Standard_Integer isfirstdoneInteger = 0;
1657   Standard_Integer islastdoneInteger = 0;
1658
1659
1660   if (isfirstdone) {
1661     first = aProj.FirstParam();
1662     isfirstdoneInteger = 1;
1663   }
1664  
1665   if (islastdone) {
1666     last= aProj.LastParam();
1667     islastdoneInteger = 1;
1668   }
1669
1670   di << isfirstdoneInteger << " "<< islastdoneInteger << " "<< first << " "<< last << " \n";
1671
1672   return 0;
1673 }
1674
1675 template<typename T, typename HT>
1676 static void DoIsNull(Draw_Interpretor& di)
1677 {
1678   HT aHandle;
1679   //    QVERIFY (aHandle.IsNull());
1680   QCOMPARE (aHandle.IsNull(), Standard_True);
1681   const T* p = aHandle.Access();
1682 #if OCC_VERSION_HEX > 0x060700
1683   //QVERIFY (!p);
1684   //QVERIFY (p == 0);
1685   QCOMPARE (!p, Standard_True);
1686   QCOMPARE (p == 0, Standard_True);
1687 #endif
1688
1689   aHandle = new T;
1690   //QVERIFY (!aHandle.IsNull());
1691   QCOMPARE (!aHandle.IsNull(), Standard_True);
1692   p = aHandle.Access();
1693   //QVERIFY (p);
1694   //QVERIFY (p != 0);
1695   QCOMPARE (p != NULL, Standard_True);
1696   QCOMPARE (p != 0, Standard_True);
1697 }
1698
1699 //=======================================================================
1700 //function : OCC24533
1701 //purpose  : 
1702 //=======================================================================
1703 static Standard_Integer OCC24533 (Draw_Interpretor& di, Standard_Integer n, const char**)
1704 {
1705   if (n != 1) return 1;
1706
1707   DoIsNull<Standard_Transient, Handle(Standard_Transient)>(di);
1708   DoIsNull<Standard_Persistent, Handle(Standard_Persistent)>(di);
1709
1710   return 0;
1711 }
1712
1713 // Dummy class to test interface for compilation issues
1714 class QABugs_HandleClass : public Standard_Transient
1715 {
1716 public:
1717   Standard_Integer HandleProc (Draw_Interpretor& , Standard_Integer  , const char** theArgVec)
1718   {
1719     std::cerr << "QABugs_HandleClass[" << this << "] " << theArgVec[0] << "\n";
1720     return 0;
1721   }
1722   DEFINE_STANDARD_RTTI(QABugs_HandleClass) // Type definition
1723 };
1724 DEFINE_STANDARD_HANDLE    (QABugs_HandleClass, Standard_Transient)
1725 IMPLEMENT_STANDARD_HANDLE (QABugs_HandleClass, Standard_Transient)
1726 IMPLEMENT_STANDARD_RTTIEXT(QABugs_HandleClass, Standard_Transient)
1727
1728 // Dummy class to test interface for compilation issues
1729 struct QABugs_NHandleClass
1730 {
1731   Standard_Integer NHandleProc (Draw_Interpretor& , Standard_Integer  , const char** theArgVec)
1732   {
1733     std::cerr << "QABugs_NHandleClass[" << this << "] " << "" << theArgVec[0] << "\n";
1734     return 0;
1735   }
1736 };
1737
1738 //=======================================================================
1739 //function : OCC24622
1740 //purpose  : The command tests sourcing Image_PixMap to AIS_TexturedShape
1741 //=======================================================================
1742 static Standard_Integer OCC24622 (Draw_Interpretor& /*theDi*/, Standard_Integer theArgNb, const char** theArgVec)
1743 {
1744   if (theArgNb != 2)
1745   {
1746     std::cout << "Usage : " << theArgVec[0] << " texture={1D|2D}";
1747     return 1;
1748   }
1749
1750   const Handle(AIS_InteractiveContext)& anAISContext = ViewerTest::GetAISContext();
1751   if (anAISContext.IsNull())
1752   {
1753     std::cout << "Please initialize view with \"vinit\".\n";
1754     return 1;
1755   }
1756
1757   Handle(Image_PixMap) anImage = new Image_PixMap();
1758
1759   static const Image_ColorRGB aBitmap[8] =
1760   {
1761     {{255,   0, 0}}, {{0,  148, 255}}, {{ 0, 148, 255}}, {{255,  94, 0}},
1762     {{255, 121, 0}}, {{76, 255,   0}}, {{76, 255,   0}}, {{255, 202, 0}}
1763   };
1764
1765   TCollection_AsciiString aTextureTypeArg (theArgVec[1]);
1766   aTextureTypeArg.UpperCase();
1767   if (aTextureTypeArg == "1D")
1768   {
1769     anImage->InitWrapper (Image_PixMap::ImgRGB, (Standard_Byte*)aBitmap, 8, 1);
1770   }
1771   else if (aTextureTypeArg == "2D")
1772   {
1773     anImage->InitTrash (Image_PixMap::ImgRGB, 8, 8);
1774     for (Standard_Integer aRow = 0; aRow < 8; ++aRow)
1775     {
1776       for (Standard_Integer aCol = 0; aCol < 8; ++aCol)
1777       {
1778         anImage->ChangeValue<Image_ColorRGB> (aRow, aCol) = aBitmap[aRow];
1779       }
1780     }
1781   }
1782   else
1783   {
1784     std::cout << "Please specify type of texture to test {1D|2D}.\n";
1785     return 1;
1786   }
1787
1788   TopoDS_Shape aBlankShape = BRepPrimAPI_MakeBox (10.0, 10.0, 10.0).Shape();
1789
1790   Handle(AIS_TexturedShape) aTexturedShape = new AIS_TexturedShape (aBlankShape);
1791   aTexturedShape->SetTexturePixMap (anImage);
1792   anAISContext->Display (aTexturedShape, 3, 0);
1793
1794   return 0;
1795 }
1796
1797 //=======================================================================
1798 //function : OCC24667
1799 //purpose  : 
1800 //=======================================================================
1801 static Standard_Integer OCC24667 (Draw_Interpretor& di, Standard_Integer n, const char** a)
1802 {
1803   if (n == 1)
1804   {
1805     di << "OCC24667 result Wire_spine Profile [Mode [Approx]]" << "\n";
1806     di << "Mode = 0 - CorrectedFrenet," << "\n";
1807     di << "     = 1 - Frenet," << "\n";
1808     di << "     = 2 - DiscreteTrihedron" << "\n";
1809     di << "Approx - force C1-approximation if result is C0" << "\n";
1810     return 0;
1811   }
1812
1813   if (n > 1 && n < 4) return 1;
1814
1815   TopoDS_Shape Spine = DBRep::Get(a[2],TopAbs_WIRE);
1816   if ( Spine.IsNull()) return 1;
1817
1818   TopoDS_Shape Profile = DBRep::Get(a[3]);
1819   if ( Profile.IsNull()) return 1;
1820
1821   GeomFill_Trihedron Mode = GeomFill_IsCorrectedFrenet;
1822   if (n >= 5)
1823   {
1824     Standard_Integer iMode = atoi(a[4]);
1825     if (iMode == 1)
1826       Mode = GeomFill_IsFrenet;
1827     else if (iMode == 2)
1828       Mode = GeomFill_IsDiscreteTrihedron;
1829   }
1830
1831   Standard_Boolean ForceApproxC1 = Standard_False;
1832   if (n >= 6)
1833     ForceApproxC1 = Standard_True;
1834
1835   BRepOffsetAPI_MakePipe aPipe(TopoDS::Wire(Spine),
1836                                           Profile,
1837                                           Mode,
1838                                           ForceApproxC1);
1839
1840   TopoDS_Shape S = aPipe.Shape();
1841   TopoDS_Shape aSF = aPipe.FirstShape();
1842   TopoDS_Shape aSL = aPipe.LastShape();
1843
1844   DBRep::Set(a[1],S);
1845
1846   TCollection_AsciiString aStrF(a[1], "_f");
1847   TCollection_AsciiString aStrL(a[1], "_l");
1848
1849   DBRep::Set(aStrF.ToCString(), aSF);
1850   DBRep::Set(aStrL.ToCString(), aSL);
1851
1852   return 0;
1853 }
1854
1855 #include <IGESControl_Reader.hxx>
1856 #include <IGESControl_Controller.hxx>
1857 #include <IGESData_IGESEntity.hxx>
1858 #include <BRepCheck_Analyzer.hxx>
1859 #include <PTColStd_TransientPersistentMap.hxx>
1860 #include <PTopoDS_HShape.hxx>
1861 #include <Storage_Data.hxx>
1862 #include <TopExp_Explorer.hxx>
1863 #include <MgtBRep.hxx>
1864 #include <FSD_File.hxx>
1865 #include <ShapeSchema.hxx>
1866 #include <TColStd_HSequenceOfTransient.hxx>
1867 #include <PTColStd_PersistentTransientMap.hxx>
1868 #include <Storage_Root.hxx>
1869
1870 static Standard_Integer OCC24565 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
1871 {
1872   if (argc != 3) {
1873     di << "Usage : " << argv[0] << " FileNameIGS FileNameSTOR";
1874     return 1;
1875   }
1876
1877   Standard_CString sFileNameIGS = argv[1];
1878   Standard_CString sFileNameSTOR = argv[2];
1879
1880   IGESControl_Reader ICReader;
1881
1882   /* * * * * * *
1883    * Read the IGES file and make sure it is valid
1884    *
1885    * * * * * * */
1886   IGESControl_Controller::Init();
1887
1888   if (!ICReader.ReadFile(sFileNameIGS)) {
1889     printf("%s:%d - Error reading '%s'\n",__FUNCTION__,__LINE__,sFileNameIGS);fflush(stdout);
1890     return -1;
1891   }
1892
1893   int nbShapes = ICReader.NbShapes();
1894
1895   printf("%s:%d - nbShapes = '%d'\n",__FUNCTION__,__LINE__,nbShapes);fflush(stdout);
1896
1897   TopoDS_Shape Shape;
1898   if(nbShapes == 0)
1899     {
1900       Handle(TColStd_HSequenceOfTransient) faces=ICReader.GiveList("iges-faces");
1901       Handle(TColStd_HSequenceOfTransient) surfaceList=ICReader.GiveList("xst-transferrable-roots",faces);
1902
1903       if (surfaceList.IsNull())
1904         {
1905           printf("%s:%d - surfaceList.IsNull()\n",__FUNCTION__,__LINE__);fflush(stdout);
1906           return -1;
1907         }
1908       BRep_Builder builder;
1909       TopoDS_Compound* pC = new TopoDS_Compound();
1910       builder.MakeCompound(*pC);
1911
1912       for (int j=1;j<=surfaceList->Length();j++)
1913         {
1914           Handle(IGESData_IGESEntity) igesEntity=Handle(IGESData_IGESEntity)::DownCast(surfaceList->Value(j));
1915           if (igesEntity.IsNull()) continue;
1916           ICReader.ClearShapes();
1917           Standard_Boolean rv;
1918           try {
1919             rv=ICReader.TransferEntity(igesEntity);
1920           }
1921           catch (...) {
1922             rv=Standard_False;
1923           }
1924           if (!rv) {
1925             printf("%s:%d - Error transferring IGES entity\n",__FUNCTION__,__LINE__);fflush(stdout);
1926             printf("%s:%d - FormNumber = %d, TypeNumber = %d\n",__FUNCTION__,__LINE__,igesEntity->FormNumber(),igesEntity->TypeNumber());fflush(stdout);
1927             return -1;
1928           }
1929
1930           TopoDS_Shape S;
1931           try {
1932             S=ICReader.Shape();
1933           }
1934           catch(...) {
1935             printf("%s:%d - Error reading IGES entity\n",__FUNCTION__,__LINE__);fflush(stdout);
1936             printf("%s:%d - FormNumber = %d, TypeNumber = %d\n",__FUNCTION__,__LINE__,igesEntity->FormNumber(),igesEntity->TypeNumber());fflush(stdout);
1937             return -1;
1938           }
1939           if (S.IsNull()) {
1940             printf("%s:%d - NULL Surface encountered\n",__FUNCTION__,__LINE__);
1941             return -1;
1942           }
1943
1944           try
1945             {
1946               builder.Add(*pC,S);
1947             }
1948           catch(...)
1949             {
1950               printf("%s: Exception adding face.\n",__FUNCTION__);
1951             }
1952         }
1953       Shape = TopoDS_Shape(*pC);
1954     }
1955   else
1956     {
1957       Shape = ICReader.OneShape();
1958     }
1959   {
1960     BRepCheck_Analyzer brca(Shape);
1961
1962     if(!brca.IsValid())
1963       {
1964         printf("%s: Invalid shape after reading IGES file.\n",__FUNCTION__);
1965       }
1966   }
1967
1968   /* * * * * * *
1969    * Write the contents of the Shape to a STOR file
1970    *
1971    * * * * * * */
1972   PTColStd_TransientPersistentMap aMapTP;
1973   Handle(PTopoDS_HShape) aPShape_write;
1974   Handle(Storage_Data) d_write=new Storage_Data;
1975   char Name[32];
1976
1977   TopExp_Explorer Ex;
1978   int i;
1979   int max_i = 0;
1980
1981   for (i=0,Ex.Init(Shape,TopAbs_FACE);Ex.More();i++,Ex.Next())
1982     {
1983
1984       max_i = i;
1985       try {
1986         aPShape_write=MgtBRep::Translate(Ex.Current(),aMapTP,MgtBRep_WithoutTriangle);
1987       }
1988       catch (...) {
1989         printf("%s: Error translating surface '%d'\n",__FUNCTION__,i);
1990       }
1991                 
1992       sprintf(Name,"S%010d",i);
1993                 
1994       {
1995         BRepCheck_Analyzer brca(Ex.Current());
1996         if(!brca.IsValid())
1997           {
1998             printf("INVALID face '%s' in the shape, which will be written to the STOR file.\n",Name);
1999           }
2000       }
2001       try {
2002         d_write->AddRoot(Name,aPShape_write);
2003       }
2004       catch (...) {
2005         printf("%s: Error adding surface '%d', RootName = '%s'\n",__FUNCTION__,i,Name);
2006       }
2007     }
2008   printf("%s: Going to write %d surfaces.\n",__FUNCTION__,max_i+1);
2009
2010   FSD_File f_write;
2011   if(f_write.Open(sFileNameSTOR, Storage_VSWrite)!=Storage_VSOk)
2012     {
2013       printf("%s: Error opening file: %s\n", __FUNCTION__,sFileNameSTOR);
2014       return -1;
2015     }
2016   Handle(ShapeSchema) s_write=new ShapeSchema;
2017   s_write->Write(f_write,d_write);
2018   f_write.Close();
2019   printf("%s: Wrote to the STOR file.\n",__FUNCTION__);
2020
2021   /* * * * * * *
2022    * Read the contents of the Shape from a STOR file
2023    *
2024    * * * * * * */
2025   FSD_File f_read;
2026   if(f_read.Open(sFileNameSTOR, Storage_VSRead)!=Storage_VSOk)
2027     {
2028       printf("%s: Error opening file: %s\n", __FUNCTION__,sFileNameSTOR);
2029       return -1;
2030     }
2031   Handle(ShapeSchema) s_read=new ShapeSchema;
2032   Handle(Storage_Data) d_read=s_read->Read(f_read);
2033
2034   Handle(Standard_Persistent) p;
2035   Handle(Storage_Root) r;
2036   Handle(PTopoDS_HShape) aPShape_read;
2037   PTColStd_PersistentTransientMap aMapPT;
2038   TopoDS_Shape S_read;
2039
2040   printf("%s: Extracting %d faces from the STOR file.\n",__FUNCTION__,max_i+1);
2041   for(int i = 0; i <= max_i; ++i)
2042     {
2043       sprintf(Name,"S%010d",i);
2044       r=d_read->Find(Name);
2045       if(r.IsNull())
2046         {
2047           printf("%s:%d '%s' IsNull().\n",__FUNCTION__,__LINE__,Name);fflush(stdout);
2048           continue;
2049         }
2050       p=r->Object();
2051       aPShape_read = Handle(PTopoDS_HShape)::DownCast(p);
2052       try {
2053         MgtBRep::Translate(aPShape_read,aMapPT,S_read,MgtBRep_WithoutTriangle);
2054       }
2055       catch (Standard_Failure) {
2056         Handle(Standard_Failure) E=Standard_Failure::Caught();
2057         std::string str;
2058         str="Exception: ";
2059         str+=E->DynamicType()->Name();
2060         str+=" => ";
2061         str+=E->GetMessageString();
2062         printf("%s(1): %s: %s\n",__FUNCTION__,Name,str.c_str());fflush(stdout);
2063       }
2064       catch (...) {
2065         printf("%s(1): Unhandled exception in MgtBRep::Translate\n",__FUNCTION__);
2066       }
2067
2068       BRepCheck_Analyzer brca(S_read);
2069
2070       if(!brca.IsValid())
2071         {
2072           printf("%s: Read INVALID face (%s)!\n",__FUNCTION__,Name);
2073         }
2074     }
2075
2076   printf("Completed.\n");fflush(stdout);
2077
2078   return 0;
2079 }
2080
2081 #include <Handle_BRepTools_NurbsConvertModification.hxx>
2082 #include <BRepPrimAPI_MakeCylinder.hxx>
2083 #include <BRepBuilderAPI_Copy.hxx>
2084 #include <BRepTools_NurbsConvertModification.hxx>
2085 static TopoDS_Shape CreateTestShape (int& theShapeNb)
2086 {
2087   TopoDS_Compound aComp;
2088   BRep_Builder aBuilder;
2089   aBuilder.MakeCompound (aComp);
2090   //NURBS modifier is used to increase footprint of each shape
2091   Handle_BRepTools_NurbsConvertModification aNurbsModif = new BRepTools_NurbsConvertModification;
2092   TopoDS_Shape aRefShape = BRepPrimAPI_MakeCylinder (50., 100.).Solid();
2093   BRepTools_Modifier aModifier (aRefShape, aNurbsModif);
2094   if (aModifier.IsDone()) {
2095     aRefShape = aModifier.ModifiedShape (aRefShape);
2096   }
2097   int aSiblingNb = 0;
2098   for (; theShapeNb > 0; --theShapeNb) {
2099     TopoDS_Shape aShape;
2100     if (++aSiblingNb <= 100) { //number of siblings is limited to avoid long lists
2101                 aShape = BRepBuilderAPI_Copy (aRefShape, Standard_True /*CopyGeom*/).Shape();
2102     } else {
2103       aShape = CreateTestShape (theShapeNb);
2104     }
2105     aBuilder.Add (aComp, aShape);
2106   }
2107   return aComp;
2108 }
2109
2110 #include <AppStd_Application.hxx>
2111 #include <TDataStd_Integer.hxx>
2112 #include <TNaming_Builder.hxx>
2113 static Standard_Integer OCC24931 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
2114 {
2115   if (argc != 1) {
2116     di << "Usage: " << argv[0] << " invalid number of arguments"<<"\n";
2117     return 1;
2118   }
2119   TCollection_ExtendedString aFileName ("testdocument.xml");
2120   PCDM_StoreStatus aSStatus  = PCDM_SS_Failure;
2121
2122   Handle(TDocStd_Application) anApp = new AppStd_Application;
2123   {
2124     Handle(TDocStd_Document) aDoc;
2125     anApp->NewDocument ("XmlOcaf", aDoc);
2126     TDF_Label aLab = aDoc->Main();
2127     TDataStd_Integer::Set (aLab, 0);
2128     int n = 10000; //must be big enough
2129     TopoDS_Shape aShape = CreateTestShape (n);
2130     TNaming_Builder aBuilder (aLab);
2131     aBuilder.Generated (aShape);
2132
2133     aSStatus = anApp->SaveAs (aDoc, aFileName);
2134     anApp->Close (aDoc);
2135   }
2136   QCOMPARE (aSStatus, PCDM_SS_OK);
2137   return 0;
2138 }
2139
2140 #include <AppStdL_Application.hxx>
2141 #include <TDocStd_Application.hxx>
2142 #include <TDataStd_Integer.hxx>
2143 #include <TDF_AttributeIterator.hxx>
2144 //=======================================================================
2145 //function : OCC24755
2146 //purpose  : 
2147 //=======================================================================
2148 static Standard_Integer OCC24755 (Draw_Interpretor& di, Standard_Integer n, const char** a)
2149 {
2150   if (n != 1)
2151   {
2152     std::cout << "Usage : " << a[0] << "\n";
2153     return 1;
2154   }
2155
2156   Handle(TDocStd_Application) anApp = new AppStdL_Application;
2157   Handle(TDocStd_Document) aDoc;
2158   anApp->NewDocument ("MDTV-Standard", aDoc);
2159   TDF_Label aLab = aDoc->Main();
2160   TDataStd_Integer::Set (aLab, 0);
2161   TDataStd_Name::Set (aLab, "test");
2162
2163   TDF_AttributeIterator i (aLab);
2164   Handle(TDF_Attribute) anAttr = i.Value();
2165   QCOMPARE (anAttr->IsKind (STANDARD_TYPE (TDataStd_Integer)), Standard_True);
2166   i.Next();
2167   anAttr = i.Value();
2168   QCOMPARE (anAttr->IsKind (STANDARD_TYPE (TDataStd_Name)), Standard_True);
2169
2170   return 0;
2171 }
2172
2173 struct MyStubObject
2174 {
2175   MyStubObject() : ptr(0L) {}
2176   MyStubObject(void* thePtr) : ptr(thePtr) {}
2177   char overhead[40];
2178   void* ptr;
2179 };
2180
2181 //=======================================================================
2182 //function : OCC24834
2183 //purpose  : 
2184 //=======================================================================
2185 static Standard_Integer OCC24834 (Draw_Interpretor& di, Standard_Integer n, const char** a)
2186 {
2187   if (n != 1)
2188   {
2189     std::cout << "Usage : " << a[0] << "\n";
2190     return 1;
2191   }
2192
2193   NCollection_List<MyStubObject> aList;
2194   const Standard_Integer aSmallBlockSize = 40;
2195   const Standard_Integer aLargeBlockSize = 1500000;
2196
2197   // quick populate memory with large blocks
2198   try
2199   {
2200     for (;;)
2201     {
2202       aList.Append(MyStubObject(Standard::Allocate(aLargeBlockSize)));
2203     }
2204   }
2205   catch (Standard_Failure)
2206   {
2207     di << "caught out of memory for large blocks: OK\n";
2208   }
2209   catch (...)
2210   {
2211     di << "skept out of memory for large blocks: Error\n";
2212   }
2213
2214   // allocate small blocks
2215   try
2216   {
2217     for (;;)
2218     {
2219       aList.Append(MyStubObject(Standard::Allocate(aSmallBlockSize)));
2220     }
2221   }
2222   catch (Standard_Failure)
2223   {
2224     di << "caught out of memory for small blocks: OK\n";
2225   }
2226   catch (...)
2227   {
2228     di << "skept out of memory for small blocks: Error\n";
2229   }
2230
2231   // release all allocated blocks
2232   for (NCollection_List<MyStubObject>::Iterator it(aList); it.More(); it.Next())
2233   {
2234     Standard::Free(it.Value().ptr);
2235   }
2236   return 0;
2237 }
2238
2239 void QABugs::Commands_19(Draw_Interpretor& theCommands) {
2240   const char *group = "QABugs";
2241
2242   Handle(QABugs_HandleClass) aClassPtr = new QABugs_HandleClass();
2243   theCommands.Add ("OCC24202_1", "Test Handle-based procedure",
2244                    __FILE__, aClassPtr, &QABugs_HandleClass::HandleProc, group);
2245   NCollection_Handle<QABugs_NHandleClass> aNClassPtr = new QABugs_NHandleClass();
2246   theCommands.Add ("OCC24202_2", "Test NCollection_Handle-based procedure",
2247                    __FILE__, aNClassPtr, &QABugs_NHandleClass::NHandleProc, group);
2248
2249   theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
2250   theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
2251   theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
2252   theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group); 
2253   theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
2254   theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
2255   theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
2256   theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
2257   theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
2258   theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
2259   theCommands.Add ("OCC23952sweep", "OCC23952sweep nbupoles shape", __FILE__, OCC23952sweep, group);
2260   theCommands.Add ("OCC23952intersect", "OCC23952intersect nbsol shape1 shape2", __FILE__, OCC23952intersect, group);
2261   theCommands.Add ("test_offset", "test_offset", __FILE__, test_offset, group);
2262   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);
2263   theCommands.Add ("OCC24008", "OCC24008 curve surface", __FILE__, OCC24008, group);
2264   theCommands.Add ("OCC24019", "OCC24019 aShape", __FILE__, OCC24019, group);
2265   theCommands.Add ("OCC11758", "OCC11758", __FILE__, OCC11758, group);
2266   theCommands.Add ("OCC24005", "OCC24005 result", __FILE__, OCC24005, group);
2267   theCommands.Add ("OCC24137", "OCC24137 face vertex U V [N]", __FILE__, OCC24137, group);
2268   theCommands.Add ("OCC24271", "Boolean operations on NCollection_Map", __FILE__, OCC24271, group);
2269   theCommands.Add ("OCC23972", "OCC23972", __FILE__, OCC23972, group);
2270   theCommands.Add ("OCC24370", "OCC24370 edge pcurve surface prec", __FILE__, OCC24370, group);
2271   theCommands.Add ("OCC24533", "OCC24533", __FILE__, OCC24533, group);
2272   theCommands.Add ("OCC24012", "OCC24012 face edge", __FILE__, OCC24012, group);
2273   theCommands.Add ("OCC24051", "OCC24051", __FILE__, OCC24051, group);
2274   theCommands.Add ("OCC24086", "OCC24086 face wire", __FILE__, OCC24086, group);
2275   theCommands.Add ("OCC24622", "OCC24622 texture={1D|2D}\n Tests sourcing of 1D/2D pixmaps for AIS_TexturedShape", __FILE__, OCC24622, group);
2276   theCommands.Add ("OCC24667", "OCC24667 result Wire_spine Profile [Mode [Approx]], no args to get help", __FILE__, OCC24667, group);
2277   theCommands.Add ("OCC24565", "OCC24565 FileNameIGS FileNameSTOR", __FILE__, OCC24565, group);
2278   theCommands.Add ("OCC24755", "OCC24755", __FILE__, OCC24755, group);
2279   theCommands.Add ("OCC24834", "OCC24834", __FILE__, OCC24834, group);
2280   theCommands.Add ("OCC24931", "OCC24931", __FILE__, OCC24931, group);
2281   return;
2282 }