0023663: Removing 2D viewer library
[occt.git] / src / QABugs / QABugs_11.cxx
1 // Created on: 2002-03-20
2 // Created by: QA Admin
3 // Copyright (c) 2002-2012 OPEN CASCADE SAS
4 //
5 // The content of this file is subject to the Open CASCADE Technology Public
6 // License Version 6.5 (the "License"). You may not use the content of this file
7 // except in compliance with the License. Please obtain a copy of the License
8 // at http://www.opencascade.org and read it completely before using this file.
9 //
10 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
11 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
12 //
13 // The Original Code and all software distributed under the License is
14 // distributed on an "AS IS" basis, without warranty of any kind, and the
15 // Initial Developer hereby disclaims all such warranties, including without
16 // limitation, any warranties of merchantability, fitness for a particular
17 // purpose or non-infringement. Please see the License for the specific terms
18 // and conditions governing the rights and limitations under the License.
19
20
21
22 #include <stdio.h>
23
24 #include <QABugs.hxx>
25
26 #include <Draw_Interpretor.hxx>
27 #include <DBRep.hxx>
28 #include <DrawTrSurf.hxx>
29 #include <AIS_InteractiveContext.hxx>
30 #include <ViewerTest.hxx>
31 #include <AIS_Shape.hxx>
32 #include <TopoDS_Shape.hxx>
33
34 #include <Geom_Axis2Placement.hxx>
35 #include <gp.hxx>
36 #include <gp_Trsf.hxx>
37 #include <AIS_Trihedron.hxx>
38 #include <BRepPrimAPI_MakeBox.hxx>
39 #include <Graphic3d_MaterialAspect.hxx>
40 #include <ViewerTest_DoubleMapOfInteractiveAndName.hxx>
41 #include <TopoDS_Solid.hxx>
42 #include <BRepPrimAPI_MakeSphere.hxx>
43 #include <BRepPrimAPI_MakeCone.hxx>
44 #include <BRepPrimAPI_MakeCylinder.hxx>
45 #include <IGESToBRep_Reader.hxx>
46 #include <TopExp_Explorer.hxx>
47 #include <TopoDS.hxx>
48 #include <GCPnts_UniformDeflection.hxx>
49 #include <BRepAdaptor_Curve.hxx>
50 #include <IGESToBRep.hxx>
51 #include <V3d_Viewer.hxx>
52 #include <BRepAdaptor_CompCurve.hxx>
53 #include <GCPnts_AbscissaPoint.hxx>
54 #include <Standard_ErrorHandler.hxx>
55 #include <Standard_Overflow.hxx>
56 #include <Standard_Underflow.hxx>
57 #include <Standard_DivideByZero.hxx>
58 #include <OSD_SIGSEGV.hxx>
59 #include <OSD_Exception_ACCESS_VIOLATION.hxx>
60 #include <OSD_Exception_STACK_OVERFLOW.hxx>
61 #include <OSD.hxx>
62 #include <STEPCAFControl_Writer.hxx>
63 #include <STEPControl_StepModelType.hxx>
64 #include <Interface_Static.hxx>
65 #include <IFSelect_ReturnStatus.hxx>
66 #include <Standard_Failure.hxx>
67 #include <TColgp_HArray1OfPnt2d.hxx>
68 #include <Geom2dAPI_Interpolate.hxx>
69 #include <Geom2d_BSplineCurve.hxx>
70 #include <Geom2dConvert_BSplineCurveToBezierCurve.hxx>
71 #include <Geom2d_BezierCurve.hxx>
72 #include <BRep_Tool.hxx>
73 #include <GeomProjLib.hxx>
74 #include <Geom2dAPI_InterCurveCurve.hxx>
75 #include <IntRes2d_IntersectionSegment.hxx>
76 #include <TDataStd_RealArray.hxx>
77 #include <TDF_CopyLabel.hxx>
78 #include <NCollection_Vector.hxx>
79 #include <TColStd_Array1OfInteger.hxx>
80 #include <Geom_BSplineCurve.hxx>
81 #include <TColgp_Array1OfPnt.hxx>
82 #include <AIS_ListOfInteractive.hxx>
83 #include <AIS_ListIteratorOfListOfInteractive.hxx>
84 #include <ViewerTest_DoubleMapOfInteractiveAndName.hxx>
85 #include <ViewerTest_DoubleMapIteratorOfDoubleMapOfInteractiveAndName.hxx>
86 #include <BRepBuilderAPI_MakePolygon.hxx>
87 #include <gp_GTrsf.hxx>
88 #include <Poly_Triangulation.hxx>
89 #include <IGESControl_Reader.hxx>
90 #include <IGESData_IGESModel.hxx>
91 #include <IGESData_IGESEntity.hxx>
92 #include <V3d_View.hxx>
93 #include <BRepFeat_SplitShape.hxx>
94 #include <BRepAlgoAPI_Section.hxx>
95
96 #include <tcl.h>
97
98 #if ! defined(WNT)
99 extern ViewerTest_DoubleMapOfInteractiveAndName& GetMapOfAIS();
100 #else
101 Standard_EXPORT ViewerTest_DoubleMapOfInteractiveAndName& GetMapOfAIS();
102 #endif
103
104 static Standard_Integer  OCC128 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
105 {
106   Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
107   if(myAISContext.IsNull()) {
108     di << "use 'vinit' command before " << argv[0] ;
109     return 1;
110   }
111
112   Handle(Geom_Axis2Placement) aTrihedronAxis = new Geom_Axis2Placement(gp::XOY());
113
114   gp_Trsf trsf1;
115   trsf1.SetTranslation(gp_Vec(100, 100, 0));
116   aTrihedronAxis->Transform(trsf1);
117   Handle(AIS_Trihedron) myTrihedron = new AIS_Trihedron(aTrihedronAxis);
118   myTrihedron->SetColor(Quantity_NOC_LIGHTSTEELBLUE4);
119   myTrihedron->SetSize(100);
120   myAISContext->Display(myTrihedron, Standard_True);
121
122 //  TopoDS_Shape shape1 = (TopoDS_Shape) BRepPrimAPI_MakeBox(50,50,50);
123   TopoDS_Shape shape1 = BRepPrimAPI_MakeBox(50,50,50).Shape();
124   Handle(AIS_Shape) AS = new AIS_Shape(shape1);
125   AS->SetDisplayMode(1);
126   Graphic3d_MaterialAspect mat(Graphic3d_NOM_PLASTIC);
127   AS->SetMaterial(mat);
128   AS->SetColor(Quantity_NOC_RED);
129   myAISContext->Display(AS);
130
131   gp_Trsf TouchTrsf;
132   TouchTrsf.SetTranslation(gp_Vec(20, 20, 0));
133
134   myAISContext->ResetLocation(AS);
135   myAISContext->SetLocation(AS , TouchTrsf) ;
136   myAISContext->Redisplay(AS, Standard_True);
137
138  return 0;
139 }
140
141   // Remove as bad version of QAAddOrRemoveSelected from QADraw
142 //static Standard_Integer OCC129 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
143 //{
144 //  if( argc != 3) {
145 //    di<<"Usage : " << argv[0] << " shape islocal\n";
146 //    return 1;
147 //  }
148 //  //get AIS_Shape:
149 //  Handle(AIS_InteractiveContext) anAISCtx = ViewerTest::GetAISContext();
150 //
151 // //   ViewerTest_DoubleMapOfInteractiveAndName& aMap =
152 // //                          ViewerTest::GetDataMapOfAIS ();
153 //  ViewerTest_DoubleMapOfInteractiveAndName& aMap = GetMapOfAIS();
154 //
155 //  TCollection_AsciiString aName(argv[1]);
156 //  Handle(AIS_InteractiveObject) AISObj;
157 //
158 //  if(aMap.IsBound2(aName)) {
159 //    AISObj = aMap.Find2(aName);
160 //    if(AISObj.IsNull()){
161 //      di<<"No interactive object \n";
162 //      return 1;
163 //    }
164 //
165 //    Standard_Integer aNum = -1;
166 //
167 //    if(atoi(argv[2])) {
168 //      aNum = anAISCtx->OpenLocalContext();
169 //    }
170 //
171 //    if(anAISCtx->HasOpenedContext()){
172 //      anAISCtx->InitSelected();
173 //      anAISCtx->AddOrRemoveSelected(AISObj);
174 //    }
175 //    else {
176 //      anAISCtx->InitCurrent();
177 //      anAISCtx->AddOrRemoveCurrentObject(AISObj);
178 //    }
179 //
180 //    if(aNum >= 0) {
181 // //      anAISCtx->CloseLocalContext(aNum);
182 //    }
183 //
184 //    return 0;
185 //  }
186 //  //select this shape:
187 //  else {
188 //    di<<"Use 'vdisplay' before";
189 //    return 1;
190 //  }
191 //}
192
193 static Standard_Integer OCC136 (Draw_Interpretor& di, Standard_Integer argc, const char ** /*argv*/)
194 {
195   if(argc > 1){
196     di<<"Usage: OCC136\n";
197     return 1;
198   }
199
200   //create some primitives:
201   // Two basic points:
202   Standard_Real Size=100;
203   gp_Pnt P0(0,0,0), P1(Size,Size,Size);
204   //box
205   TopoDS_Solid aBox = BRepPrimAPI_MakeBox(P0,P1);
206   //sphere
207   TopoDS_Solid aSphere = BRepPrimAPI_MakeSphere (Size*0.5);
208   //cone
209   gp_Ax2 anAx2(P1, gp_Dir(1,1,1));
210   TopoDS_Solid aCone = BRepPrimAPI_MakeCone(anAx2, Size*0.7, Size*0.3, Size);
211   //cylinder
212   anAx2.SetLocation(gp_Pnt(Size,0,0));
213   anAx2.SetDirection(gp_Dir(-1,-1,1));
214   TopoDS_Solid aCyl = BRepPrimAPI_MakeCylinder(anAx2, Size*0.5, Size);
215
216   Handle(AIS_InteractiveContext) anAISCtx = ViewerTest::GetAISContext();
217   if(anAISCtx.IsNull()){
218     di<<"Null interactive context. Use 'vinit' at first.\n";
219     return 1;
220   }
221
222   if(anAISCtx->HasOpenedContext()) anAISCtx->CloseAllContexts();
223   anAISCtx->EraseAll();
224
225   //load primitives to context
226   Handle(AIS_Shape) aSh1 = new AIS_Shape(aBox);
227   anAISCtx->Display(aSh1);
228
229   Handle(AIS_Shape) aSh2 = new AIS_Shape(aSphere);
230   anAISCtx->Display(aSh2);
231
232   Handle(AIS_Shape) aSh3 = new AIS_Shape(aCone);
233   anAISCtx->Display(aSh3);
234
235   Handle(AIS_Shape) aSh4 = new AIS_Shape(aCyl);
236   anAISCtx->Display(aSh4);
237
238   //set selected
239   anAISCtx->InitCurrent();
240   anAISCtx->AddOrRemoveCurrentObject(aSh1);
241   anAISCtx->AddOrRemoveCurrentObject(aSh2);
242   anAISCtx->AddOrRemoveCurrentObject(aSh3);
243   anAISCtx->AddOrRemoveCurrentObject(aSh4);
244
245   //remove all this objects from context
246   anAISCtx->Clear(aSh1, Standard_False);
247   anAISCtx->Clear(aSh2, Standard_False);
248   anAISCtx->Clear(aSh3, Standard_False);
249   anAISCtx->Clear(aSh4, Standard_False);
250   return 0;
251 }
252
253 static int BUC60610(Draw_Interpretor& di, Standard_Integer argc, const char ** argv) {
254   if(argc < 2){
255     printf("Usage: %s  iges_input [name]\n",argv[0]);
256     return(1);
257   }
258   Standard_Character *Ch;
259
260   if(argc > 2) {
261     Ch = new Standard_Character[strlen(argv[2])+3];
262   }
263   IGESToBRep_Reader IR;
264   IR.LoadFile (argv[1]);
265   IR.Clear();
266   IR.TransferRoots();
267   TopoDS_Shape aTopShape = IR.OneShape();
268   TopExp_Explorer ex(aTopShape, TopAbs_EDGE);
269   Standard_Integer i=0;
270   for( ; ex.More(); ex.Next()){
271     const TopoDS_Edge &E = TopoDS::Edge(ex.Current());
272     BRepAdaptor_Curve aCurve(E);
273     GCPnts_UniformDeflection plin(aCurve, 0.1);
274     di << "Num points = " << plin.NbPoints() << "\n";
275     if(argc > 2) {
276       i++;
277       sprintf(Ch,"%s_%i",argv[2],1);
278       DBRep::Set(Ch,E);
279     }
280   }
281   return (1);
282 }
283
284 static Standard_Integer BUC60661(Draw_Interpretor& di, Standard_Integer argc, const char ** a)
285 {
286   if(argc!=2)
287     {
288       di << "Usage : " << a[0] << " file.igs" << "\n";
289       return -1;
290     }
291
292   Handle(AIS_InteractiveContext) myContext = ViewerTest::GetAISContext();
293
294   if(myContext.IsNull()) {
295     di << "use 'vinit' command before " << a[0] << "\n";
296     return -1;
297   }
298
299   // MKV 30.03.05
300 #if ((TCL_MAJOR_VERSION > 8) || ((TCL_MAJOR_VERSION == 8) && (TCL_MINOR_VERSION >= 4))) && !defined(USE_NON_CONST)
301   const Standard_Character *DD = Tcl_GetVar(di.Interp(),"Draw_DataDir",TCL_GLOBAL_ONLY);
302 #else
303   Standard_Character *DD = Tcl_GetVar(di.Interp(),"Draw_DataDir",TCL_GLOBAL_ONLY);
304 #endif
305
306   Standard_Character  *file1 = new Standard_Character [strlen(DD)+strlen(a[1])+2];
307   sprintf(file1,"%s/%s",DD,a[1]);
308
309   IGESToBRep_Reader reader;
310   Standard_Integer status = reader.LoadFile(file1);
311   if( !status ) {
312     IGESToBRep::Init();
313       reader.TransferRoots();
314       TopoDS_Shape shape = reader.OneShape();
315
316       Handle(AIS_Shape) importedShape = new AIS_Shape(shape);
317
318       Handle(V3d_Viewer) myPView = myContext->CurrentViewer();
319
320       if( ! myPView.IsNull() && (myPView->DefaultVisualization() == V3d_WIREFRAME) )
321         importedShape->SetDisplayMode(AIS_WireFrame);
322       else importedShape->SetDisplayMode(AIS_Shaded);
323       myContext->Display(importedShape);
324     }
325
326   printf("\n End of my  IGES to 3D-viewer *****************>\n");
327   return 0;
328 }
329
330
331 //====================================================
332 //
333 // Following code is inserted from
334 // /dn03/KAS/dev/QAopt/src/QADraw/QADraw_TOPOLOGY.cxx
335 // ( 75455 Apr 16 18:59)
336 //
337 //====================================================
338
339 //OCC105
340 #include <BRepTools_WireExplorer.hxx>
341 #include <BRep_Tool.hxx>
342 #include <GCPnts_UniformAbscissa.hxx>
343 #include <TopExp.hxx>
344
345 //
346 // usage : OCC105 shape
347 //
348 // comments:
349 //GCPnts_UniformAbscissa returns bad end point foe first edge. Its value is
350
351 //Value Pnt = -338.556216693211 -394.465571897208 0
352 //should be
353 //Value Pnt = -307.47165394 -340.18073533 0
354
355 static int OCC105(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
356 {
357   if (argc != 2){
358     di<<"Usage : OCC105 shape\n";
359     return 1;
360   }
361 //  TopoDS_Wire myTopoDSWire = TopoDS::Wire(DBRep::Get("aa.brep"));
362   TopoDS_Wire myTopoDSWire = TopoDS::Wire(DBRep::Get(argv[1]));
363   Standard_Real l = 0.5; //atof(argv[2]);
364   // Find the first vertex of the wire
365   BRepTools_WireExplorer wire_exp(myTopoDSWire);
366   TopoDS_Vertex vlast;
367   {
368     TopoDS_Vertex vw1, vw2;
369     TopExp::Vertices(myTopoDSWire,vw1,vw2);
370     TopoDS_Vertex ve1, ve2;
371     TopoDS_Edge edge = TopoDS::Edge(wire_exp.Current());
372     TopExp::Vertices(edge,ve1,ve2);
373     if (vw1.IsSame(ve1) || vw1.IsSame(ve2))
374       vlast = vw1;
375     else {
376 //      assert(vw2.IsSame(ve1) || vw2.IsSame(ve2));
377       vlast = vw2;
378     }
379   }
380   for ( ; wire_exp.More(); wire_exp.Next())
381     {
382       di << "\n\n New Edge \n"   << "\n";
383       Standard_Real newufirst, newulast;
384       TopoDS_Edge edge = TopoDS::Edge(wire_exp.Current());
385       Standard_Real ufirst, ulast;
386       Handle(Geom_Curve) acurve;
387       TopoDS_Vertex ve1, ve2;
388       TopExp::Vertices(edge,ve1,ve2);
389       if (ve1.IsSame(vlast))
390         {
391           acurve = BRep_Tool::Curve(edge, ufirst, ulast);
392           newufirst = ufirst;
393           newulast  = ulast;
394           vlast = ve2;
395         }
396       else
397         {
398 //          assert(ve2.IsSame(vlast));
399 //          assert ( wire_exp.Orientation( ) == TopAbs_REVERSED );
400           acurve = BRep_Tool::Curve( edge, ufirst, ulast );
401           newufirst = acurve->ReversedParameter( ufirst );
402           newulast  = acurve->ReversedParameter( ulast );
403           acurve = acurve->Reversed( );
404           vlast = ve1;
405         }
406
407       GeomAdaptor_Curve   curve;
408       GCPnts_UniformAbscissa  algo;
409       curve.Load(acurve);
410       algo.Initialize( curve, l, newufirst, newulast );
411       if (!algo.IsDone())
412         di << "Not Done!!!"   << "\n";
413       for (Standard_Integer Index = 1; Index<=algo.NbPoints();Index++) {
414         Standard_Real t = algo.Parameter(Index);
415         gp_Pnt      pt3 = curve.Value(t);
416         di << "Parameter t = " << t   << "\n";
417         di << "Value Pnt = " << pt3.X()<<" " <<pt3.Y()<<" " << pt3.Z()  << "\n";
418       }
419     }
420   return 0;
421
422 }
423
424 #include <TColStd_SequenceOfTransient.hxx>
425 #include <GeomFill_Pipe.hxx>
426 static int pipe_OCC9 (Draw_Interpretor& di,
427                       Standard_Integer n, const char ** a)
428 {
429   if (n < 6) {
430     di << "Usage: " << a[0] << " result path cur1 cur2 radius [tolerance]" << "\n";
431     return 1;
432   }
433
434   TColStd_SequenceOfTransient aCurveSeq;
435   Standard_Integer i;
436   for (i=2 ; i<=4; i++) {
437     Handle(Geom_Curve) aC = Handle(Geom_Curve)::DownCast( DrawTrSurf::Get(a[i]) );
438     if (aC.IsNull()) {
439       di << a[i] << " is not a curve" << "\n";
440       return 1;
441     }
442     aCurveSeq.Append(aC);
443   }
444
445   GeomFill_Pipe aPipe(Handle(Geom_Curve)::DownCast( aCurveSeq(1) ),
446                       Handle(Geom_Curve)::DownCast( aCurveSeq(2) ),
447                       Handle(Geom_Curve)::DownCast( aCurveSeq(3) ),
448                       atof (a[5]) );
449
450   if (n == 7) {
451     aPipe.Perform(atof (a[6]), Standard_True);
452   } else {
453     aPipe.Perform(Standard_True/*, Standard_True*/);
454   }
455
456   Handle(Geom_Surface) aSurf = aPipe.Surface();
457
458   DrawTrSurf::Set(a[1], aSurf);
459   return 0;
460 }
461
462 //======================================================================
463 // OCC125
464 // usage : OCC125 shell
465 //======================================================================
466 #include <ShapeFix_Shell.hxx>
467
468 Standard_Integer  OCC125(Draw_Interpretor& di ,
469                          Standard_Integer n,
470                          const char ** a)
471 {
472   if (n!=2) {
473     di<<" Use OCC125 shell";
474     return 1;
475   }
476
477   TopoDS_Shape S = DBRep::Get(a[1]);
478
479   if (S.IsNull()) {
480     di<<" Null shape is not allowed";
481     return 1;
482   }
483
484   TopAbs_ShapeEnum aT;
485   aT=S.ShapeType();
486   if (aT!=TopAbs_SHELL) {
487     di<<" Shape Type must be SHELL";
488     return 1;
489   }
490
491   const TopoDS_Shell& aShell = TopoDS::Shell(S);
492   //
493   Standard_Boolean isAccountMultiConex, bNonManifold, bResult;
494
495   isAccountMultiConex=Standard_True;
496   bNonManifold=Standard_False;
497
498   Handle (ShapeFix_Shell) aFix = new ShapeFix_Shell(aShell);
499   bResult=aFix->FixFaceOrientation(aShell, isAccountMultiConex, bNonManifold);
500
501   di<<"bResult="<<(Standard_Integer)bResult;
502
503   TopoDS_Shape aShape;
504   aShape=aFix->Shape();
505
506   TCollection_AsciiString aName(a[1]), aDef("_sh"), aRName;
507   aRName=aName;
508   aRName=aRName+aDef;
509   DBRep::Set (aRName.ToCString(), aShape);
510   di<<aRName.ToCString();
511   //
512   return 0;
513 }
514
515 #include <BRepLib_FindSurface.hxx>
516 #include <BRepBuilderAPI_MakeFace.hxx>
517 Standard_Integer  OCC157(Draw_Interpretor& di,
518                          Standard_Integer n,
519                          const char ** a)
520 //static Standard_Integer findplanarsurface(Draw_Interpretor&, Standard_Integer n, const char ** a)
521 {
522   if (n<3) {
523     di << "bad number of arguments" <<"\n";
524     return 1;
525   }
526
527   // try to read a shape:
528   TopoDS_Shape inputShape=DBRep::Get(a[2]);
529   if (inputShape.IsNull() || inputShape.ShapeType() != TopAbs_WIRE) {
530     di << "Invalid input shape"<< "\n";
531     return 1;
532   }
533   Standard_Real toler = atof(a[3]);
534   TopoDS_Wire aWire = TopoDS::Wire(inputShape);
535   BRepLib_FindSurface FS(aWire, toler, Standard_True);
536   if(FS.Found()) {
537     di<<"OCC157: OK; Planar surface is found"<<"\n";
538     Handle(Geom_Surface) aSurf = FS.Surface();
539     BRepBuilderAPI_MakeFace aMakeFace(aSurf,aWire,Standard_True);
540     if(aMakeFace.IsDone()) {
541       TopoDS_Face aFace = aMakeFace.Face();
542       DBRep::Set(a[1],aFace);
543     }
544   }
545   else di<<"OCC157: ERROR; Planar surface is not found with toler = "<<toler <<"\n";
546   return 0;
547
548 }
549
550 // #include <MyCommandsCMD.h>
551 #include <ShapeFix_Shape.hxx>
552 #include <BRepOffset_MakeOffset.hxx>
553 #include <BRepOffsetAPI_MakeOffset.hxx>
554 #include <BRepOffset_Mode.hxx>
555 #include <GeomAbs_JoinType.hxx>
556 #include <AIS_Shape.hxx>
557
558 #include <BRepTools.hxx>
559
560 Standard_Integer  OCC165(Draw_Interpretor& di ,
561                          Standard_Integer n,
562                          const char ** a)
563
564
565 //=======================================================================
566
567 // static int YOffset (Draw_Interpretor& di, Standard_Integer argc, const char ** argv);
568
569 // void MyOffsets_Commands(Draw_Interpretor& theCommands)
570 // {
571 //      theCommands.Add("yoffset" , "yoffset" , __FILE__, YOffset, " Offset on Z Direction");
572 // }
573
574 //=======================================================================
575
576 // static int YOffset (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
577  {
578    if (n > 2)
579      {
580        di <<"Usage : " << a[0] << " [file]"<<"\n";
581        return 1;
582      }
583         di.Eval ("axo");
584
585 #define _OFFSET_TELCO_
586 #ifdef _OFFSET_TELCO_
587
588         Standard_CString file = a[1];
589
590         BRep_Builder aBuilder;
591         TopoDS_Shape theShape;
592         //BRepTools::Read(theShape, Standard_CString("/dn02/users_SUN/inv/3/OCC165/2d_tr_line.brep"), aBuilder);
593         BRepTools::Read(theShape, file, aBuilder);
594         DBRep::Set("shape", theShape);
595
596         TopoDS_Wire theWire = TopoDS::Wire(theShape);
597
598         Standard_Real anOffset = 1.5;
599
600 #else
601
602         Standard_Real xA = 0.0, xB = 200.0, xC = 200.0, xD = 0.0,
603                 yA = 0.0, yB = 0.0, yC = 200.0, yD = 200.0,
604                 zA = 0.0, zB = 0.0, zC = 0.0, zD = 0.0;
605
606         BRepBuilderAPI_MakePolygon theSquare;
607         TopoDS_Vertex theA = BRepBuilderAPI_MakeVertex(gp_Pnt(xA, yA, zA));
608         theSquare.Add(theA);
609         TopoDS_Vertex theB = BRepBuilderAPI_MakeVertex(gp_Pnt(xB, yB, zB));
610         theSquare.Add(theB);
611         TopoDS_Vertex theC = BRepBuilderAPI_MakeVertex(gp_Pnt(xC, yC, zC));
612         theSquare.Add(theC);
613         TopoDS_Vertex theD = BRepBuilderAPI_MakeVertex(gp_Pnt(xD, yD, zD));
614         theSquare.Add(theD);
615
616         theSquare.Close();
617         TopoDS_Wire theWire = theSquare.Wire();
618
619         Standard_Real anOffset = 10;
620
621
622 #endif /* _OFFSET_TELCO_ */
623
624
625         TopoDS_Face theFace = BRepBuilderAPI_MakeFace(theWire).Face();
626         DBRep::Set("face", theFace);
627
628
629         Standard_Real anAlt = 0.;
630         GeomAbs_JoinType theJoin = GeomAbs_Intersection;
631 //GeomAbs_Intersection; //GeomAbs_Arc;
632         BRepOffsetAPI_MakeOffset aMakeOffset(theFace, theJoin);
633         aMakeOffset.AddWire(theWire);
634
635         aMakeOffset.Perform(anOffset, anAlt);
636
637         TopoDS_Shape theOffsetShapePos = aMakeOffset.Shape();
638         DBRep::Set("offset", theOffsetShapePos);
639         return 0;
640 //      return TCL_OK;
641 }
642
643 #include<BRepAlgoAPI_Cut.hxx>
644 #include<BRepAlgo_Cut.hxx>
645
646 #include<BRepPrimAPI_MakeHalfSpace.hxx>
647 #include<Handle_Geom_CartesianPoint.hxx>
648 #include<Geom_CartesianPoint.hxx>
649 #include<AIS_Point.hxx>
650
651 #include <BRepBuilderAPI_MakeEdge.hxx>
652 #include <BRepBuilderAPI_MakeWire.hxx>
653 #include <BRepBuilderAPI_MakeFace.hxx>
654
655 static Standard_Integer OCC297 (Draw_Interpretor& di,Standard_Integer /*argc*/, const char ** argv )
656
657 {
658
659   Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
660   if(myAISContext.IsNull()) {
661     di << "use 'vinit' command before " << argv[0] << "\n";
662     return -1;
663   }
664
665
666     gp_Pnt pt1_(250.,250.,0.);
667     gp_Pnt pt2_(-250.,250.,0.);
668     gp_Pnt pt3_(-250.,-250.,0.);
669     gp_Pnt pt4_(250.,-250.,0.);
670     BRepBuilderAPI_MakeEdge edg1_(pt1_, pt2_);
671     BRepBuilderAPI_MakeEdge edg2_(pt2_, pt3_);
672     BRepBuilderAPI_MakeEdge edg3_(pt3_, pt4_);
673     BRepBuilderAPI_MakeEdge edg4_(pt4_, pt1_);
674
675     BRepBuilderAPI_MakeWire wire_(edg1_, edg2_, edg3_, edg4_);
676     BRepBuilderAPI_MakeFace face_(wire_);
677     TopoDS_Face sh_ = face_.Face();
678
679     int up = 1;
680
681     gp_Pnt g_pnt;
682     if(up)
683         g_pnt = gp_Pnt(0,0,-100);
684     else
685           g_pnt = gp_Pnt(0,0,100);
686
687     myAISContext->EraseAll(Standard_False);
688         Handle(Geom_CartesianPoint) GEOMPoint = new Geom_CartesianPoint(g_pnt);
689         Handle(AIS_Point) AISPoint = new AIS_Point(GEOMPoint);
690     myAISContext->Display(AISPoint);
691
692     BRepPrimAPI_MakeHalfSpace half_(sh_, g_pnt);
693     TopoDS_Solid sol1_ = half_.Solid();
694 //         Handle(AIS_Shape) AISHalf = new AIS_Shape(sol1_);
695 //         AISHalf->SetColor(Quantity_NOC_GREEN);
696 //     myAISContext->Display(AISHalf);
697
698   DBRep::Set("Face",sol1_);
699
700     gp_Ax1 ax1_(gp_Pnt(0., 0., -100.), gp_Dir(0., 0., 1.));
701
702     Standard_Real x=0., y=0., z=-80.;
703
704     BRepPrimAPI_MakeBox box(gp_Pnt(x,y,z),gp_Pnt(x+150,y+200,z+200));
705 //         Handle(AIS_Shape) AISBox = new AIS_Shape(box);
706 //         AISBox->SetColor(Quantity_NOC_BLUE1);
707 //     myAISContext->Display(AISBox);
708
709  DBRep::Set("Box",box.Shape());
710
711 //     BRepAlgoAPI_Cut cut( sol1_, box.Shape() );
712 //         //BRepAlgoAPI_Cut cut(  box.Shape(), sol1_ );
713 //     TopoDS_Shape sh1_ = cut.Shape();
714 //         Handle(AIS_Shape) AISCut = new AIS_Shape(sh1_);
715 //         AISCut->SetColor(Quantity_NOC_RED);
716 //     myAISContext->Display(AISCut);
717
718 //  DBRep::Set("Cut",sh1_);
719
720   return 0;
721
722 }
723 #include<GProp_GProps.hxx>
724 #include<BRepGProp.hxx>
725
726 static Standard_Integer OCC305 (Draw_Interpretor& di,Standard_Integer argc, const char ** argv )
727
728 {
729   if (argc != 2)
730   {
731     di <<"Usage : " << argv[0] << " file"<<"\n";
732     return 1;
733   }
734   Standard_CString file = argv[1];
735
736   Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
737   if(myAISContext.IsNull()) {
738     di << "use 'vinit' command before " << argv[0] << "\n";
739     return -1;
740   }
741
742 TopoDS_Shape sh;
743 BRep_Builder builder;
744 //BRepTools::Read(sh, "/dn02/users_SUN/inv/3/OCC305/testc.brep", builder);
745 BRepTools::Read(sh, file, builder);
746
747 TopoDS_Wire wire;
748 builder.MakeWire(wire);
749 TopoDS_Edge ed;
750 TopoDS_Vertex vt1, vt2;
751 TopExp_Explorer wex(sh, TopAbs_EDGE);
752 for(;wex.More();wex.Next())
753 {
754     ed = TopoDS::Edge(wex.Current());
755     TopExp::Vertices(ed, vt1, vt2);
756     builder.UpdateVertex(vt1, 0.01);
757     builder.UpdateVertex(vt2, 0.01);
758     builder.UpdateEdge(ed, 0.01);
759     builder.Add(wire, ed);
760
761
762     GProp_GProps lprop;
763     BRepGProp::LinearProperties(ed, lprop);
764     printf("\n length = %f", lprop.Mass());
765 }
766  DBRep::Set("Wire",wire);
767 // Handle(AIS_Shape) res = new AIS_Shape( wire );
768 // aContext->SetColor( res, Quantity_NOC_RED );
769 // aContext->Display( res );
770
771 // BRepOffsetAPI_MakeOffset off(wire, GeomAbs_Arc);
772 // off.Perform(0.5, 0);
773
774 // printf("\n IsDone = %d", off.IsDone());
775 // sh = off.Shape();
776 // res = new AIS_Shape( sh );
777 // aContext->SetColor( res, Quantity_NOC_GREEN );
778 // aContext->Display( res );
779
780   return 0;
781
782 }
783
784 static Standard_Integer OCC166 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
785 {
786
787   Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
788   if(myAISContext.IsNull()) {
789     di << "use 'vinit' command before " << argv[0];
790     return 1;
791   }
792
793   BRepPrimAPI_MakeBox aBox(gp_Pnt(0, 0, 0), 100, 100, 100);
794   Handle(AIS_Shape) anAISBox = new AIS_Shape(aBox.Shape());
795   myAISContext->Display(anAISBox, 1);
796   anAISBox->SetSelectionMode(-1);
797   Standard_Integer myLocContInd = myAISContext->OpenLocalContext();
798   myAISContext->CloseLocalContext(myLocContInd);
799   Standard_Integer aSelMode = ((Handle(AIS_InteractiveObject)) anAISBox)->SelectionMode();
800   if(aSelMode != -1)
801     return 1;
802
803   return 0;
804 }
805
806 #include <TDocStd_Document.hxx>
807 #include <DDocStd.hxx>
808 #include <PCDM_StoreStatus.hxx>
809 #include <TDocStd_Application.hxx>
810
811 static Standard_Integer OCC381_Save (Draw_Interpretor& di, Standard_Integer nb, const char ** a)
812 {
813   if (nb != 2) {
814     di << "Usage: " << a[0] << " Doc" << "\n";
815     return 1;
816   }
817
818   Handle(TDocStd_Document) D;
819   if (!DDocStd::GetDocument(a[1],D)) return 1;
820
821   Handle(TDocStd_Application) A;
822   if (!DDocStd::Find(A)) return 1;
823
824   TCollection_ExtendedString theStatusMessage;
825   if (!D->IsSaved()) {
826     di << "this document has never been saved" << "\n";
827     return 0;
828   }
829   PCDM_StoreStatus theStatus = A->Save(D, theStatusMessage);
830   if (theStatus != PCDM_SS_OK ) {
831     switch ( theStatus ) {
832     case PCDM_SS_DriverFailure: {
833       di << " Could not store , no driver found to make it " <<"\n" ;
834       break ;
835     }
836     case PCDM_SS_WriteFailure: {
837       di << " Write access failure " << "\n" ;
838       break;
839     }
840     case PCDM_SS_Failure: {
841       di << " Write failure " << "\n" ;
842     }
843     }
844     return 1;
845   }
846   return 0;
847 }
848
849 static Standard_Integer OCC381_SaveAs (Draw_Interpretor& di, Standard_Integer nb, const char ** a)
850 {
851   if (nb != 3) {
852     di << "Usage: " << a[0] << " Doc Path" << "\n";
853     return 1;
854   }
855
856   Handle(TDocStd_Document) D;
857   if (!DDocStd::GetDocument(a[1],D)) return 1;
858
859   TCollection_ExtendedString path (a[2]);
860   Handle(TDocStd_Application) A;
861   if (!DDocStd::Find(A)) return 1;
862
863   TCollection_ExtendedString theStatusMessage;
864   PCDM_StoreStatus theStatus = A->SaveAs(D,path, theStatusMessage);
865   if (theStatus != PCDM_SS_OK ) {
866     switch ( theStatus ) {
867     case PCDM_SS_DriverFailure: {
868       di << " Could not store , no driver found to make it " <<"\n" ;
869       break ;
870     }
871     case PCDM_SS_WriteFailure: {
872       di << " Write access failure " << "\n" ;
873       break;
874     }
875     case PCDM_SS_Failure: {
876       di << " Write failure " << "\n" ;
877     }
878     }
879     return 1;
880   }
881
882   return 0;
883 }
884
885 #include <BRepClass3d_SolidClassifier.hxx>
886
887 Standard_Integer OCC299bug (Draw_Interpretor& di,
888                             Standard_Integer n,
889                             const char ** a)
890 {
891   char sbf[512];
892
893   if (n < 3) {
894     di << "Usage : " << a[0] << " Solid Point [Tolerance=1.e-7]" << "\n";
895     return -1;
896   }
897
898   TopoDS_Shape aS = DBRep::Get(a[1]);
899   if (aS.IsNull()) {
900     sprintf(sbf, " Null Shape is not allowed here\n");
901     di<<sbf;
902     return 1;
903   }
904
905   if (aS.ShapeType()!=TopAbs_SOLID) {
906     sprintf(sbf, " Shape type must be SOLID\n");
907     di<<sbf;
908     return 1;
909   }
910   //
911   Standard_Real aTol=1.e-7;
912   TCollection_AsciiString sIN("IN"), sOUT("OUT of"), sON("ON"), sUNKNOWN("UNKNOWN");
913   TopAbs_State aState = TopAbs_UNKNOWN;
914   gp_Pnt aP(8., 9., 10.);
915
916   if (!DrawTrSurf::GetPoint(a[2], aP) ) {
917     sprintf(sbf, " Null Point is not allowed here\n");
918     di<<sbf;
919     return 1;
920   }
921
922   aTol=1.e-7;
923   if (n==4) {
924     aTol=atof(a[3]);
925   }
926   //
927   BRepClass3d_SolidClassifier aSC(aS);
928   aSC.Perform(aP,aTol);
929   //
930   aState = aSC.State();
931   //
932   sprintf(sbf, "The point is "); di<<sbf;
933   //
934   switch (aState) {
935   case TopAbs_IN:
936     sprintf(sbf, sIN.ToCString());
937     break;
938   case TopAbs_OUT:
939     sprintf(sbf, sOUT.ToCString());
940     break;
941   case TopAbs_ON:
942     sprintf(sbf, sON.ToCString());
943     break;
944   case TopAbs_UNKNOWN:
945     sprintf(sbf, sUNKNOWN.ToCString());
946     break;
947   default:
948     sprintf(sbf, sUNKNOWN.ToCString());
949     break;
950   }
951   di<<sbf;
952         //
953   sprintf(sbf, " shape\n");
954   di<<sbf;
955
956   return 0;
957 }
958
959 #include <OSD_Process.hxx>
960 #include <OSD_Path.hxx>
961
962 static Standard_Integer OCC309bug (Draw_Interpretor& di, Standard_Integer nb, const char ** a)
963 {
964   if (nb != 1) {
965     di << "Usage: " << a[0] << "\n";
966     return 1;
967   }
968   OSD_Process p;
969   OSD_Path d = p.CurrentDirectory();
970   TCollection_AsciiString s;
971   d.SystemName(s);
972   di << "*" <<  s.ToCString() << "*" << "\n";
973   d.UpTrek();
974   d.SystemName(s);
975   di <<  "*" <<  s.ToCString() <<  "*" <<"\n";
976   return 0;
977 }
978
979 static Standard_Integer OCC310bug (Draw_Interpretor& di, Standard_Integer nb, const char ** a)
980 {
981   if (nb != 1) {
982     di << "Usage: " << a[0] << "\n";
983     return 1;
984   }
985   OSD_Path p("/where/you/want/tmp/qwerty/tmp/");
986   di << p.Trek().ToCString() << "\n";
987   p.UpTrek();
988   di << p.Trek().ToCString() << "\n";
989   return 0;
990 }
991
992 #include <BRepAlgoAPI_Fuse.hxx>
993 #include <BRepAlgoAPI_Common.hxx>
994 #include <BRepAlgo_Fuse.hxx>
995 #include <BRepAlgo_Common.hxx>
996
997 static Standard_Integer OCC277bug (Draw_Interpretor& di, Standard_Integer nb, const char ** a)
998 {
999   //if (nb != 1) {
1000   //  cout << "Usage: " << a[0] << endl;
1001   //  return 1;
1002   //}
1003   if(nb < 1 || nb > 2) {
1004     di << "Usage : " << a[0] << " [BRepAlgoAPI/BRepAlgo = 1/0]" << "\n";
1005     return 1;
1006   }
1007   Standard_Boolean IsBRepAlgoAPI = Standard_True;
1008   if (nb == 2) {
1009     Standard_Integer IsB = atoi(a[1]);
1010     if (IsB != 1) {
1011       IsBRepAlgoAPI = Standard_False;
1012     }
1013   }
1014
1015   BRepPrimAPI_MakeBox box1( 100, 100, 100 );
1016   BRepPrimAPI_MakeBox box2( gp_Pnt( 50, 50,50 ), 200, 200, 200 );
1017
1018   TopoDS_Shape shape1 = box1.Shape();
1019   TopoDS_Shape shape2 = box2.Shape();
1020
1021 //#if ! defined(BRepAlgoAPI_def01)
1022 //  TopoDS_Shape fuse = BRepAlgoAPI_Fuse( shape1, shape2 );
1023 //  TopoDS_Shape comm = BRepAlgoAPI_Common( shape1, shape2 );
1024 //#else
1025 //  TopoDS_Shape fuse = BRepAlgo_Fuse( shape1, shape2 );
1026 //  TopoDS_Shape comm = BRepAlgo_Common( shape1, shape2 );
1027 //#endif
1028   TopoDS_Shape fuse,comm;
1029   if (IsBRepAlgoAPI) {
1030     di << "fuse = BRepAlgoAPI_Fuse( shape1, shape2 )" <<"\n";
1031     di << "comm = BRepAlgoAPI_Common( shape1, shape2 )" <<"\n";
1032     fuse = BRepAlgoAPI_Fuse( shape1, shape2 );
1033     comm = BRepAlgoAPI_Common( shape1, shape2 );
1034   } else {
1035     di << "fuse = BRepAlgo_Fuse( shape1, shape2 )" <<"\n";
1036     fuse = BRepAlgo_Fuse( shape1, shape2 );
1037     di << "comm = BRepAlgo_Common( shape1, shape2 )" <<"\n";
1038     comm = BRepAlgo_Common( shape1, shape2 );
1039   }
1040
1041   return 0;
1042 }
1043
1044 #include <ShapeAnalysis_Edge.hxx>
1045
1046 static Standard_Integer OCC333bug (Draw_Interpretor& di, Standard_Integer n, const char ** a)
1047 {
1048   if( n < 3) {
1049     di<<"-1"<<"\n";
1050     di << "Usage: " << a[0] << " edge1 edge2 [toler domaindist]" << "\n";
1051     return 1;
1052   }
1053   TopoDS_Shape Sh1 = DBRep::Get(a[1]);
1054   TopoDS_Shape Sh2 = DBRep::Get(a[2]);
1055   if(Sh1.IsNull() || Sh2.IsNull()) {
1056     di<<"-2"<<"\n";
1057     di<<"Invalid arguments"<<"\n";
1058     return 1;
1059   }
1060   TopoDS_Edge e1 = TopoDS::Edge(Sh1);
1061   TopoDS_Edge e2 = TopoDS::Edge(Sh2);
1062   if(e1.IsNull() || e2.IsNull()) {
1063     di<<"-3"<<"\n";
1064     di<<"Invalid type of arguments"<<"\n";
1065     return 1;
1066   }
1067   Standard_Real aTol = Precision::Confusion();
1068   Standard_Real aDistDomain = 0.0;
1069   Standard_Integer k = 3;
1070   if(k < n)
1071     aTol = atof(a[k++]);
1072   if(k < n)
1073     aDistDomain = atof(a[k++]);
1074
1075   ShapeAnalysis_Edge sae;
1076   if(sae.CheckOverlapping(e1,e2,aTol,aDistDomain)) {
1077     if(aDistDomain ==0.0) {
1078       di<<"1"<<"\n";
1079       di<<"Edges is overlaping comletly"<<"\n";
1080     } else {
1081       di<<"2"<<"\n";
1082       di<<"Edges is overlaped"<<"\n";
1083       di<<"with tolerance = "<<aTol<<"\n";
1084       di<<"on segment length = "<<aDistDomain<<"\n";
1085     }
1086   } else {
1087     di<<"3"<<"\n";
1088     di<<"Edges is not overlaped"<<"\n";
1089   }
1090   return 0;
1091 }
1092
1093
1094 #include <DDocStd_DrawDocument.hxx>
1095 #include <TDataStd_Name.hxx>
1096 #include <Draw.hxx>
1097 #include <XCAFDoc_ShapeTool.hxx>
1098 #include <XCAFDoc_DocumentTool.hxx>
1099 #include <TDF_LabelSequence.hxx>
1100 #include <TPrsStd_AISPresentation.hxx>
1101 #include <TDF_Data.hxx>
1102 #include <TDF_Label.hxx>
1103 #include <XCAFPrs_Driver.hxx>
1104
1105 //------------------------------------------------------------------------------------------
1106 // name    : OCC363
1107 // Purpose :
1108 //------------------------------------------------------------------------------------------
1109 static Standard_Integer OCC363 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1110 {
1111   try
1112   {
1113     OCC_CATCH_SIGNALS
1114     // 1. Verufy amount of arguments
1115     if(argc < 3) { di <<"Error OCC363 : Use : OCC363 document filename\n"; return 1; }
1116
1117     // 2. Retrieve DDocStd application
1118     Handle(TDocStd_Application) App;
1119     if (!DDocStd::Find(App)) { di << "Error OCC363 : There is no current DDocStd application\n";return 1;}
1120
1121     // 3. Open document
1122     TCollection_ExtendedString name(argv[2]);
1123     Handle(TDocStd_Document) Doc;
1124     if(App->Open(name, Doc) != PCDM_RS_OK) { di << "Error OCC363 : document was not opened successfully\n"; return 1;}
1125     Handle(DDocStd_DrawDocument) DD = new DDocStd_DrawDocument(Doc);
1126     TDataStd_Name::Set(Doc->GetData()->Root(),argv[1]);
1127     Draw::Set(argv[1],DD);
1128
1129     // 4. Create prsentations
1130     Handle(XCAFDoc_ShapeTool) shapes = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
1131     TDF_LabelSequence seq;
1132     shapes->GetFreeShapes ( seq );
1133     Handle(TPrsStd_AISPresentation) prs;
1134     for ( Standard_Integer i=1; i <= seq.Length(); i++ )
1135       if ( ! seq.Value(i).FindAttribute ( TPrsStd_AISPresentation::GetID(), prs ) )
1136         prs = TPrsStd_AISPresentation::Set(seq.Value(i),XCAFPrs_Driver::GetID());
1137   }
1138   catch(Standard_Failure) { di << "FAULTY OCC363 : Exception during reading document.\n";return 0;}
1139
1140   di << "OCC363 OK\n";
1141   return 0;
1142 }
1143
1144 // Must use OCC299
1145 ////======================================================================================
1146 //// Function : OCC372
1147 //// Purpose  :
1148 ////======================================================================================
1149 //static Standard_Integer OCC372 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1150 //{
1151 //  try
1152 //  {
1153 //    OCC_CATCH_SIGNALS
1154 //    // 1. Verufy amount of arguments
1155 //    if(argc < 2) {di << "OCC372 FAULTY. Use : OCC372 brep-file";return 0;}
1156 //
1157 //    // 2. Read solid
1158 //    BRep_Builder B;
1159 //    TopoDS_Shape Ref;
1160 //    BRepTools::Read(Ref, argv[1], B);
1161 //
1162 //    // 3. Calculate location of aP3d in relation to the solid
1163 //    gp_Pnt aP3d(6311.4862583184, -2841.3092756034, 16.461053497188);
1164 //    BRepClass3d_SolidClassifier SC(Ref);
1165 //    SC.Perform(aP3d, 1e-7);
1166 //
1167 //    // 4. Check returned state. The point must be inside the solid.
1168 //    TopAbs_State aState=SC.State();
1169 //    switch (aState)
1170 //    {
1171 //    case TopAbs_OUT:
1172 //      di<<"OCC372 FAULTY. aState = TopAbs_OUT";
1173 //      return 0;
1174 //    case TopAbs_ON:
1175 //      di<<"OCC372 FAULTY. aState = TopAbs_ON";
1176 //     return 0;
1177 //    case TopAbs_IN:
1178 //      di<<"OCC372 OK. aState = TopAbs_IN" ;
1179 //      return 0;
1180 //    default:
1181 //      di<<"OCC372 FAULTY. aState = UNKNOWN";
1182 //      return 0;
1183 //    }
1184 //  }
1185 //  catch (Standard_Failure) { di<<"OCC372 FAULTY. Exception raised"; }
1186 //
1187 //  return 0;
1188 //}
1189
1190 #include <BRepTopAdaptor_FClass2d.hxx>
1191
1192 //======================================================================================
1193 // Function : OCC377
1194 // Purpose  :
1195 //======================================================================================
1196 static Standard_Integer OCC377 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1197 {
1198   try
1199   {
1200     OCC_CATCH_SIGNALS
1201     // 1. Verify validity of arguments
1202     if ( argc < 1 ) {di << "Error OCC377. Use  OCC377 file x y precuv \n";return 0;}
1203
1204     // 2. Initialize parameters
1205     gp_Pnt2d p2d;
1206     p2d.SetX ( atof(argv[2]) );
1207     p2d.SetY ( atof(argv[3]) );
1208     Standard_Real precuv = atof (argv[4] );
1209
1210     // 3. Read shape
1211     BRep_Builder B;
1212     TopoDS_Shape Shape;
1213     BRepTools::Read ( Shape, argv[1], B );
1214
1215     // 4. Verify whether enrtry point is on wire and reversed ones (indeed results of veridying must be same)
1216     TopExp_Explorer exp;
1217     Standard_Integer i=1;
1218     for (exp.Init(Shape.Oriented(TopAbs_FORWARD),TopAbs_WIRE); exp.More(); exp.Next(), i++)
1219     {
1220       // 4.1. Verify whether enrtry point is on wire
1221       const TopoDS_Wire& wir = TopoDS::Wire(exp.Current());
1222       TopoDS_Face newFace = TopoDS::Face(Shape.EmptyCopied());
1223
1224       TopAbs_Orientation orWire = wir.Orientation();
1225       newFace.Orientation(TopAbs_FORWARD);
1226       B.Add(newFace,wir);
1227
1228       BRepTopAdaptor_FClass2d FClass2d1(newFace,precuv);
1229       TopAbs_State stat1 = FClass2d1.PerformInfinitePoint();
1230       //di << "Wire " << i << ": Infinite point is " <<
1231       //  ( stat1 == TopAbs_IN ? "IN" : stat1 == TopAbs_OUT ? "OUT" : stat1 == TopAbs_ON ? "ON" : "UNKNOWN" ) << "\n";
1232
1233       TCollection_AsciiString TmpString;
1234       stat1 == TopAbs_IN ? TmpString.AssignCat("IN") : stat1 == TopAbs_OUT ? TmpString.AssignCat("OUT") : stat1 == TopAbs_ON ? TmpString.AssignCat("ON") : TmpString.AssignCat("UNKNOWN");
1235       di << "Wire " << i << ": Infinite point is " << TmpString.ToCString() << "\n";
1236
1237       stat1 = FClass2d1.Perform(p2d);
1238       //di << "Wire " << i << ": point ( " << p2d.X() << ", " << p2d.Y() << " ) is " <<
1239       //  ( stat1 == TopAbs_IN ? "IN" : stat1 == TopAbs_OUT ? "OUT" : stat1 == TopAbs_ON ? "ON" : "UNKNOWN" ) << "\n";
1240
1241       TmpString.Clear();
1242       stat1 == TopAbs_IN ? TmpString.AssignCat("IN") : stat1 == TopAbs_OUT ? TmpString.AssignCat("OUT") : stat1 == TopAbs_ON ? TmpString.AssignCat("ON") : TmpString.AssignCat("UNKNOWN");
1243       di << "Wire " << i << ": point ( " << p2d.X() << ", " << p2d.Y() << " ) is " << TmpString.ToCString() << "\n";
1244
1245       // 4.2. Verify whether enrtry point is on reversed wire
1246       newFace = TopoDS::Face(Shape.EmptyCopied());
1247       newFace.Orientation(TopAbs_FORWARD);
1248       orWire = TopAbs::Reverse(orWire);
1249       B.Add(newFace,wir.Oriented(orWire));
1250       BRepTopAdaptor_FClass2d FClass2d2(newFace,precuv);
1251       TopAbs_State stat2 = FClass2d2.PerformInfinitePoint();
1252       //di << "Reversed Wire " << i << ": Infinite point is " <<
1253       //  ( stat2 == TopAbs_IN ? "IN" : stat2 == TopAbs_OUT ? "OUT" : stat2 == TopAbs_ON ? "ON" : "UNKNOWN" ) << "\n";
1254
1255       TmpString.Clear();
1256       stat2 == TopAbs_IN ? TmpString.AssignCat("IN") : stat2 == TopAbs_OUT ? TmpString.AssignCat("OUT") : stat2 == TopAbs_ON ? TmpString.AssignCat("ON") : TmpString.AssignCat("UNKNOWN");
1257       di << "Reversed Wire " << i << ": Infinite point is " << TmpString.ToCString() << "\n";
1258
1259       stat2 = FClass2d2.Perform(p2d);
1260       //di << "Reversed Wire " << i << ": point ( " << p2d.X() << ", " << p2d.Y() << " ) is " <<
1261       //  ( stat2 == TopAbs_IN ? "IN" : stat2 == TopAbs_OUT ? "OUT" : stat2 == TopAbs_ON ? "ON" : "UNKNOWN" ) << "\n";
1262
1263       TmpString.Clear();
1264       stat2 == TopAbs_IN ? TmpString.AssignCat("IN") : stat2 == TopAbs_OUT ? TmpString.AssignCat("OUT") : stat2 == TopAbs_ON ? TmpString.AssignCat("ON") : TmpString.AssignCat("UNKNOWN");
1265       di << "Reversed Wire " << i << ": point ( " << p2d.X() << ", " << p2d.Y() << " ) is " << TmpString.ToCString() << "\n";
1266
1267       // 4.3. Compare results (they must be same)
1268       if(stat1 ==stat2) di << "OCC377 OK" << "\n";
1269       else {di << "OCC377 FAULTY" << "\n"; return 0;}
1270     }
1271   }
1272   catch(Standard_Failure)
1273   {
1274     di << "OCC377 Exception";
1275   }
1276
1277   return 0;
1278 }
1279
1280 #include <ShapeUpgrade_ShapeDivideAngle.hxx>
1281 #include <ShapeBuild_ReShape.hxx>
1282
1283 //=======================================================================
1284 //function : OCC22
1285 //purpose  :
1286 //=======================================================================
1287 static Standard_Integer OCC22 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1288 {
1289   try
1290   {
1291     OCC_CATCH_SIGNALS
1292     // 1. Verify arguments of the command
1293     if (argc < 5) { di << "OCC22 FAULTY. Use : OCC22 Result Shape CompoundOfSubshapesToBeDivided ConsiderLocation"; return 0;}
1294
1295     Standard_Boolean aConsiderLocation;
1296     if(strcmp(argv[4], "0")==0) aConsiderLocation = Standard_False;
1297     else aConsiderLocation = Standard_True;
1298
1299     // 2. Iniitialize aShapeUpgrade
1300     ShapeUpgrade_ShapeDivideAngle aShapeUpgrade(M_PI/2.);
1301     // precision
1302     aShapeUpgrade.SetPrecision (Precision::Confusion());
1303     // tolerance
1304     aShapeUpgrade.SetMaxTolerance(0.1);
1305     // subshapes to be divided
1306     TopoDS_Shape aSubShapesToBeDivided = DBRep::Get(argv[3]);
1307     if(aSubShapesToBeDivided.IsNull()) {di << "OCC22 FAULTY. Compound of subshapes to be divided is not exist. Please, verify input values. \n";return 0;}
1308     aShapeUpgrade.Init(aSubShapesToBeDivided);
1309     // context
1310     Handle(ShapeBuild_ReShape) aReshape = new ShapeBuild_ReShape;
1311     aShapeUpgrade.SetContext(aReshape);
1312     if(aConsiderLocation) aReshape->ModeConsiderLocation() = Standard_True;
1313
1314     // 3. Perform splitting
1315     if (aShapeUpgrade.Perform (Standard_False))         di << "Upgrade_SplitRevolution_Done \n";
1316     else if (aShapeUpgrade.Status (ShapeExtend_OK))     di << "Upgrade_SplitRevolution_OK \n";
1317     else if (aShapeUpgrade.Status (ShapeExtend_FAIL)) { di << "OCC22 FAULTY. Operation failed. Angle was not divided\n";return 0;}
1318
1319     // 4. Perform rebuilding shape
1320     // 4.1. Retrieve Shape
1321     TopoDS_Shape anInitShape = DBRep::Get(argv[2]);
1322     if(anInitShape.IsNull()) { di << "OCC22 FAULTY. Initial shape is not exist. Please verify input values \n"; return 0;}
1323     // 4.2 Rebuid retrieved shape
1324     TopoDS_Shape aResultShape = aReshape->Apply(anInitShape);
1325     // 4.3. Create result Draw shape
1326     DBRep::Set(argv[1], aResultShape);
1327   }
1328   catch (Standard_Failure) {di << "OCC22 Exception \n" ;return 0;}
1329
1330   return 0;
1331 }
1332
1333
1334 #include <ShapeProcess_OperLibrary.hxx>
1335 #include <ShapeProcess_ShapeContext.hxx>
1336 #include <ShapeProcess.hxx>
1337
1338 #include <TopTools_DataMapIteratorOfDataMapOfShapeShape.hxx>
1339 #include <BRepMesh_IncrementalMesh.hxx>
1340
1341 //=======================================================================
1342 //function : OCC24
1343 //purpose  :
1344 //=======================================================================
1345 static Standard_Integer OCC24 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1346 {
1347   try
1348   {
1349     OCC_CATCH_SIGNALS
1350     // 1. Verify amount of arguments of the command
1351     if (argc < 6) { di << "OCC24 FAULTY. Use : OCC22 Result Shape CompoundOfSubshapes ResourceFileName SequenceName"; return 0;}
1352
1353     // 2. Retrieve parameters
1354     // initial shape
1355     TopoDS_Shape anInitShape = DBRep::Get(argv[2]);
1356     if(anInitShape.IsNull()) { di << "OCC24 FAULTY. Initial shape is not exist. Please verify input values \n"; return 0;}
1357     // compound of subshapes
1358     TopoDS_Shape aSubShapes = DBRep::Get(argv[3]);
1359     if(aSubShapes.IsNull()) {di << "OCC24 FAULTY. Compound of subshapes is not exist. Please, verify input values. \n";return 0;}
1360     // name of resource file
1361     const char* aResourceFile = argv[4];
1362     // name of sequence from resource file to be executed
1363     const char* aSequenceName = argv[5];
1364
1365     // 3. Initialize ShapeContext and perform sequence of operation specified with resource file
1366     ShapeProcess_OperLibrary::Init();
1367     Handle(ShapeProcess_ShapeContext) aShapeContext = new ShapeProcess_ShapeContext (aSubShapes, aResourceFile);
1368     aShapeContext->SetDetalisation (TopAbs_EDGE);
1369     ShapeProcess::Perform (aShapeContext, aSequenceName);
1370
1371     // 4. Rebuild initil shape in accordance with performed operation
1372     Handle(ShapeBuild_ReShape) aReshape = new ShapeBuild_ReShape;
1373     TopTools_DataMapIteratorOfDataMapOfShapeShape anIter (aShapeContext->Map());
1374     for (; anIter.More(); anIter.Next())
1375       aReshape->Replace(anIter.Key(), anIter.Value());
1376     TopoDS_Shape aResultShape = aReshape->Apply(anInitShape);
1377
1378     // 5 Create resultant Draw shape
1379     DBRep::Set(argv[1], aResultShape);
1380
1381   }
1382   catch (Standard_Failure) {di << "OCC24 Exception \n" ;return 0;}
1383
1384   return 0;
1385 }
1386
1387 //=======================================================================
1388 //function : OCC369
1389 //purpose  : Verify whether exception occurs during building mesh
1390 //=======================================================================
1391 static Standard_Integer OCC369(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1392 {
1393   try
1394   {
1395     OCC_CATCH_SIGNALS
1396     // 1. Verify amount of arguments of the command
1397     if (argc < 2) { di << "OCC369 FAULTY. Use : OCC369 Shape \n"; return 0;}
1398
1399     // 2. Retrieve shape
1400     TopoDS_Shape aShape = DBRep::Get(argv[1]);
1401     if(aShape.IsNull()) {di << "OCC369 FAULTY. Entry shape is NULL \n"; return 0;}
1402
1403     // 3. Build mesh
1404     BRepMesh_IncrementalMesh aMesh(aShape, 0.2, Standard_True, M_PI/6);
1405
1406   }
1407   catch (Standard_Failure) {di << "OCC369 Exception \n" ;return 0;}
1408
1409   di << "OCC369 OK \n";
1410   return 0;
1411 }
1412
1413 #include <math_Vector.hxx>
1414 #include <math_Matrix.hxx>
1415 static Standard_Integer OCC524 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1416 {
1417   if(argc != 9){
1418     di<<"Usage : " << argv[0] << " LowerVector UpperVector InitialValueVector LowerRowMatrix UpperRowMatrix LowerColMatrix UpperColMatrix InitialValueMatrix\n";
1419     return 1;
1420   }
1421   Standard_Integer LowerVector = atoi(argv[1]);
1422   Standard_Integer UpperVector = atoi(argv[2]);
1423   Standard_Real InitialValueVector = atof(argv[3]);
1424   Standard_Integer LowerRowMatrix = atoi(argv[4]);
1425   Standard_Integer UpperRowMatrix = atoi(argv[5]);
1426   Standard_Integer LowerColMatrix = atoi(argv[6]);
1427   Standard_Integer UpperColMatrix = atoi(argv[7]);
1428   Standard_Real InitialValueMatrix = atof(argv[8]);
1429
1430   math_Vector Vector1(LowerVector, UpperVector);
1431   math_Vector Vector2(LowerVector, UpperVector);
1432
1433   math_Vector Vector(LowerVector, UpperVector, InitialValueVector);
1434   math_Matrix Matrix(LowerRowMatrix, UpperRowMatrix, LowerColMatrix, UpperColMatrix, InitialValueMatrix);
1435
1436   //Vector.Dump(cout);
1437   //cout<<endl;
1438
1439   //Matrix.Dump(cout);
1440   //cout<<endl;
1441
1442   Vector1.Multiply(Vector, Matrix);
1443
1444   //Vector1.Dump(cout);
1445   Standard_SStream aSStream1;
1446   Vector1.Dump(aSStream1);
1447   di << aSStream1;
1448   di<<"\n";
1449
1450   Vector2.TMultiply(Vector, Matrix);
1451
1452   //Vector2.Dump(cout);
1453   Standard_SStream aSStream2;
1454   Vector2.Dump(aSStream2);
1455   di << aSStream2;
1456   di<<"\n";
1457
1458   return 0;
1459 }
1460
1461 #include <GeomPlate_BuildPlateSurface.hxx>
1462 //=======================================================================
1463 //function : OCC525
1464 //purpose  :
1465 //=======================================================================
1466 static Standard_Integer OCC525(Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
1467 {
1468   try
1469   {
1470     OCC_CATCH_SIGNALS
1471     GeomPlate_BuildPlateSurface aBuilder;
1472     aBuilder.Perform();
1473   }
1474   catch (Standard_RangeError) { di << "OCC525 Exception \n" ;return 0; }
1475   //catch (...) { di << "OCC525 Exception \n" ;return 0; }
1476
1477   di << "OCC525 OK \n";
1478   return 0;
1479 }
1480
1481 #include <BRepPrimAPI_MakeWedge.hxx>
1482 #include <gce_MakeRotation.hxx>
1483 #include <gce_MakeTranslation.hxx>
1484 #include <BRepBuilderAPI_Transform.hxx>
1485 #include <BRepPrimAPI_MakeWedge.hxx>
1486 #include <BRepAlgoAPI_Fuse.hxx>
1487 #include <BRepAlgoAPI_Cut.hxx>
1488 #include <BRepAlgo_Fuse.hxx>
1489 #include <BRepAlgo_Cut.hxx>
1490 //=======================================================================
1491 //function :  OCC578
1492 //purpose  :
1493 //=======================================================================
1494 static Standard_Integer OCC578 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1495 {
1496   //if (argc!=4) {
1497   //  di<<"Usage : " << argv[0] << " shape1 shape2 shape3\n";
1498   //  return 1;
1499   //}
1500   if(argc < 4 || argc > 5) {
1501     di << "Usage : " << argv[0] << " shape1 shape2 shape3 [BRepAlgoAPI/BRepAlgo = 1/0]" << "\n";
1502     return 1;
1503   }
1504   Standard_Boolean IsBRepAlgoAPI = Standard_True;
1505   if (argc == 5) {
1506     Standard_Integer IsB = atoi(argv[4]);
1507     if (IsB != 1) {
1508       IsBRepAlgoAPI = Standard_False;
1509 //      di << "Error: There is not BRepAlgo_Fuse class" << "\n";
1510 //      return 1;
1511 //      di << "Error: There is not BRepAlgo_Cut class" << "\n";
1512 //      return 1;
1513     }
1514   }
1515
1516   gp_Pnt P0(0,0,0.0);
1517   double xperiod = 1.0;
1518   double yperiod = 1.0;
1519   double sub_thick = 0.5;
1520
1521   // mask_substrate
1522   //TopoDS_Shape substrate = BRepPrimAPI_MakeBox( P0, xperiod, yperiod, sub_thick );
1523   TopoDS_Shape substrate = BRepPrimAPI_MakeBox( P0, xperiod, yperiod, sub_thick ).Shape();
1524
1525   // --------------------------------------------------------------------
1526
1527   // wedge
1528   //TopoDS_Shape wedge1 = BRepPrimAPI_MakeWedge(0.5, 0.05, 0.5,
1529         //                                    0.1,  0.1  , 0.4, 0.4 );
1530   TopoDS_Shape wedge1 = BRepPrimAPI_MakeWedge(0.5, 0.05, 0.5,
1531                                               0.1,  0.1  , 0.4, 0.4 ).Shape();
1532
1533   gp_Trsf rotate = gce_MakeRotation ( gp_Pnt(0.0,0.0,0.0),
1534                                      gp_Dir(1.0,0.0,0.0),
1535                                      1.570795 );
1536
1537   gp_Trsf translate = gce_MakeTranslation(gp_Pnt( 0.0, -0.5, 0.0),
1538                                           gp_Pnt( 0.25, 0.25, 0.5)
1539                                           );
1540
1541   rotate.PreMultiply( translate );
1542
1543   TopoDS_Shape wedge1a = BRepBuilderAPI_Transform( wedge1, rotate );
1544
1545   if (wedge1a.IsNull()) {
1546     di<<" Null shape1 is not allowed\n";
1547     return 1;
1548   }
1549   DBRep::Set(argv[1], wedge1a);
1550
1551   // --------------------------------------------------------------------
1552
1553   // wedge top
1554   //TopoDS_Shape wedge2 = BRepPrimAPI_MakeWedge(0.5, 0.3, 0.5,
1555         //                                    0.1,  0.1  , 0.4, 0.4 );
1556   TopoDS_Shape wedge2 = BRepPrimAPI_MakeWedge(0.5, 0.3, 0.5,
1557                                               0.1,  0.1  , 0.4, 0.4 ).Shape();
1558
1559   gp_Trsf rotate2 = gce_MakeRotation ( gp_Pnt(0.0,0.0,0.0),
1560                                       gp_Dir(1.0,0.0,0.0),
1561                                       1.570795 * 3.0 );
1562
1563   gp_Trsf translate2 = gce_MakeTranslation(gp_Pnt( 0.0, 0.0, 0.0),
1564                                            gp_Pnt( 0.25, 0.25, 0.5)
1565                                            );
1566
1567   rotate2.PreMultiply( translate2 );
1568
1569   TopoDS_Shape wedge2a = BRepBuilderAPI_Transform( wedge2, rotate2 );
1570
1571   if (wedge2a.IsNull()) {
1572     di<<" Null shape2 is not allowed\n";
1573     return 1;
1574   }
1575   DBRep::Set(argv[2], wedge2a);
1576
1577
1578   // combine wedges
1579 //#if ! defined(BRepAlgoAPI_def01)
1580 //  TopoDS_Shape wedge_common = BRepAlgoAPI_Fuse(wedge1a , wedge2a);
1581 //#else
1582 //  TopoDS_Shape wedge_common = BRepAlgo_Fuse(wedge1a , wedge2a);
1583 //#endif
1584   TopoDS_Shape wedge_common;
1585   if (IsBRepAlgoAPI) {
1586     di << "wedge_common = BRepAlgoAPI_Fuse(wedge1a , wedge2a)" <<"\n";
1587     wedge_common = BRepAlgoAPI_Fuse(wedge1a , wedge2a);
1588   } else {
1589     di << "wedge_common = BRepAlgo_Fuse(wedge1a , wedge2a)" <<"\n";
1590     wedge_common = BRepAlgo_Fuse(wedge1a , wedge2a);
1591   }
1592
1593   // remove wedge area from substrate
1594 //#if ! defined(BRepAlgoAPI_def01)
1595 //  TopoDS_Shape sub_etch1 = BRepAlgoAPI_Cut(substrate, wedge_common);
1596 //#else
1597 //  TopoDS_Shape sub_etch1 = BRepAlgo_Cut(substrate, wedge_common);
1598 //#endif
1599   TopoDS_Shape sub_etch1;
1600   if (IsBRepAlgoAPI) {
1601     di << "sub_etch1 = BRepAlgoAPI_Cut(substrate, wedge_common)" <<"\n";
1602     sub_etch1 = BRepAlgoAPI_Cut(substrate, wedge_common);
1603   } else {
1604     di << "sub_etch1 = BRepAlgo_Cut(substrate, wedge_common)" <<"\n";
1605     sub_etch1 = BRepAlgo_Cut(substrate, wedge_common);
1606   }
1607
1608   if (sub_etch1.IsNull()) {
1609     di<<" Null shape3 is not allowed\n";
1610     return 1;
1611   }
1612   DBRep::Set(argv[3], sub_etch1);
1613
1614   return 0;
1615 }
1616
1617 #include <Standard_GUID.hxx>
1618 //=======================================================================
1619 //function :  OCC669
1620 //purpose  :
1621 //=======================================================================
1622 static Standard_Integer OCC669 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1623 {
1624   if(argc != 2){
1625     di<<"Usage : " << argv[0] << " GUID\n";
1626     return -1;
1627   }
1628   Standard_GUID guid(argv[1]);
1629   //guid.ShallowDump(cout);
1630   Standard_SStream aSStream;
1631   guid.ShallowDump(aSStream);
1632   di << aSStream;
1633   di<<"\n";
1634   return 0;
1635 }
1636
1637 #include <XCAFDoc.hxx>
1638 //=======================================================================
1639 //function :  OCC738_ShapeRef
1640 //purpose  :
1641 //=======================================================================
1642 static Standard_Integer OCC738_ShapeRef (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1643 {
1644   if(argc != 1){
1645     di<<"Usage : " << argv[0] << "\n";
1646     return -1;
1647   }
1648   const Standard_GUID& guid = XCAFDoc::ShapeRefGUID ();
1649   //guid.ShallowDump(cout);
1650   Standard_SStream aSStream;
1651   guid.ShallowDump(aSStream);
1652   di << aSStream;
1653   return 0;
1654 }
1655
1656 //=======================================================================
1657 //function :  OCC738_Assembly
1658 //purpose  : 
1659 //=======================================================================
1660 static Standard_Integer OCC738_Assembly (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1661 {
1662   if(argc != 1){
1663     di<<"Usage : " << argv[0] << "\n";
1664     return -1;
1665   }
1666   const Standard_GUID& guid = XCAFDoc::AssemblyGUID ();
1667   //guid.ShallowDump(cout);
1668   Standard_SStream aSStream;
1669   guid.ShallowDump(aSStream);
1670   di << aSStream;
1671   return 0;
1672 }
1673
1674 #if defined(DDataStd_def01)
1675 #include <DDataStd_DrawPresentation.hxx>
1676 //=======================================================================
1677 //function :  OCC739_DrawPresentation
1678 //purpose  : 
1679 //=======================================================================
1680 static Standard_Integer OCC739_DrawPresentation (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1681 {
1682   if(argc != 1){
1683     di<<"Usage : " << argv[0] << "\n";
1684     return -1;
1685   }
1686   const Standard_GUID& guid = DDataStd_DrawPresentation::GetID() ;
1687   //guid.ShallowDump(cout);
1688   Standard_SStream aSStream;
1689   guid.ShallowDump(aSStream);
1690   di << aSStream;
1691   return 0;
1692 }
1693 #endif
1694
1695 //=======================================================================
1696 //function :  OCC708
1697 //purpose  : 
1698 //=======================================================================
1699 static Standard_Integer OCC708 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1700 {
1701   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1702   if(aContext.IsNull()) { 
1703     di << argv[0] << "ERROR : use 'vinit' command before " << "\n";
1704     return 1;
1705   }
1706
1707   if ( argc != 2) {
1708     di << "ERROR : Usage : " << argv[0] << " shape ; Deactivate the current transformation" << "\n";
1709     return 1;
1710   }
1711   
1712   Standard_Boolean updateviewer = Standard_True, PutInCollector = Standard_True;
1713
1714   ViewerTest_DoubleMapOfInteractiveAndName& aMap = GetMapOfAIS();
1715   
1716   TCollection_AsciiString aName(argv[1]);
1717   Handle(AIS_InteractiveObject) AISObj;
1718
1719   if(!aMap.IsBound2(aName)) {
1720     di << "Use 'vdisplay' before" << "\n";
1721     return 1;
1722   } else {
1723     AISObj = Handle(AIS_InteractiveObject)::DownCast(aMap.Find2(aName));
1724     if(AISObj.IsNull()){
1725       di << argv[1] << " : No interactive object" << "\n";
1726       return 1;
1727     } 
1728     AISObj->UnsetTransformation();
1729     if (!aContext->HasOpenedContext()) {
1730       aContext->OpenLocalContext();
1731     }
1732     aContext->Erase(AISObj, updateviewer, PutInCollector);
1733     aContext->UpdateCurrentViewer();
1734     aContext->Display(AISObj, updateviewer);
1735     aContext->UpdateCurrentViewer();
1736   }
1737   return 0;
1738 }
1739
1740 //=======================================================================
1741 //function :  OCC670
1742 //purpose  :
1743 //=======================================================================
1744 #include <TColStd_Array2OfInteger.hxx>
1745 static Standard_Integer OCC670 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1746 {
1747   if(argc != 1){
1748     di<<"Usage : " << argv[0] << "\n";
1749     return -1;
1750   }
1751   TColStd_Array2OfInteger Array2OfInteger(1,1,1,1);
1752   Array2OfInteger.SetValue(5,5,55);
1753   return 0;
1754 }
1755
1756 #include <GeomAPI_ProjectPointOnSurf.hxx>
1757 //=======================================================================
1758 //function :  OCC867
1759 //purpose  : 
1760 //=======================================================================
1761 static Standard_Integer OCC867(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1762 {
1763   if (argc!=7)  
1764   {
1765     di<<"Usage : " << argv[0] << " Point Surface Umin Usup Vmin Vsup \n";
1766     return 1;   
1767   }
1768   
1769   gp_Pnt aPoint3d;        
1770   DrawTrSurf::GetPoint(argv[1],aPoint3d);
1771   Handle (Geom_Surface) aSurface=DrawTrSurf::GetSurface(argv[2]);
1772   Standard_Real             Umin=atof(argv[3]);
1773   Standard_Real             Usup=atof(argv[4]);
1774   Standard_Real             Vmin=atof(argv[5]);
1775   Standard_Real             Vsup=atof(argv[6]);
1776  
1777   if (aSurface.IsNull()) {
1778     di << argv[2] << " Null surface \n" ;
1779     return 1;
1780   }
1781   
1782   GeomAPI_ProjectPointOnSurf PonSurf;
1783   PonSurf.Init(aSurface, Umin, Usup, Vmin, Vsup);
1784   PonSurf.Perform(aPoint3d);
1785
1786   return 0; 
1787 }
1788
1789 //=======================================================================
1790 //function :  OCC909
1791 //purpose  : 
1792 //=======================================================================
1793 static Standard_Integer OCC909 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1794 {
1795   if (argc!=3)  
1796   {
1797     di<<"Usage : " << argv[0] << " wire face\n";
1798     return 1;   
1799   }
1800   
1801   TopoDS_Wire awire = TopoDS::Wire(DBRep::Get(argv[1])); //read the wire
1802   TopoDS_Face aface = TopoDS::Face(DBRep::Get(argv[2])); //read the face
1803   if (awire.IsNull() || aface.IsNull()) {
1804     di << "Null object" << "\n";
1805     return 1;
1806   }
1807
1808   Standard_Integer count = 0;
1809   TopExp_Explorer TE(awire, TopAbs_VERTEX);
1810   if ( TE.More()) {
1811     BRepTools_WireExplorer WE;
1812     for ( WE.Init(awire,aface); WE.More(); WE.Next()) {
1813       TopoDS_Edge E = WE.Current();
1814       count++;
1815     }
1816   }
1817   di << "Count = " << count << "\n";
1818
1819   return 0; 
1820 }
1821
1822 //=======================================================================
1823 //function :  OCC921
1824 //purpose  : 
1825 //=======================================================================
1826 static Standard_Integer OCC921 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1827 {
1828   if (argc != 2)
1829   {
1830     di <<"Usage : " << argv[0] << " face"<<"\n";
1831     return 1;
1832   }
1833   Standard_Real u1, u2, v1, v2;
1834   TopoDS_Face F = TopoDS::Face( DBRep::Get(argv[1]) ); //read the shape
1835   if (F.IsNull())
1836     return 1;
1837   BRepTools::UVBounds(F, u1, u2, v1, v2);
1838   di << "Bounds: " << u1 << "   " << u2 << "   " << v1 << "   " << v2 << "\n";
1839   return 0;
1840 }
1841
1842 #include <Expr_NamedUnknown.hxx>
1843 #include <Expr_GeneralExpression.hxx>
1844 #include <Expr_Exponential.hxx>
1845 //=======================================================================
1846 //function :  OCC902
1847 //purpose  : 
1848 //=======================================================================
1849 static Standard_Integer OCC902(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1850 {
1851   if (argc != 2)
1852   {
1853     di <<"Usage : " << argv[0] << " expression"<<"\n";
1854     return 1;
1855   }
1856
1857  TCollection_AsciiString  myStr(argv[1]);
1858
1859  Handle (Expr_NamedUnknown)      myNamed = new Expr_NamedUnknown(myStr);
1860  Handle (Expr_Exponential)       oldExpr = new Expr_Exponential(myNamed); 
1861  Handle (Expr_GeneralExpression) newExpr = oldExpr->Derivative(myNamed);
1862
1863  
1864  TCollection_AsciiString  res        = newExpr->String();
1865  Standard_CString         resStr     = res.ToCString();
1866  TCollection_AsciiString  res_old    = oldExpr->String();
1867  Standard_CString         res_oldStr = res_old.ToCString();
1868  
1869
1870  di << "X = " << argv[1] << "\n";
1871  di << "Y = " << res_oldStr << "\n";
1872  di << "Y' = " << resStr  << "\n";
1873
1874  return 0;
1875 }
1876
1877 #include <DDF.hxx>
1878 #include <TPrsStd_AISViewer.hxx>
1879 #include <TDF_Label.hxx>
1880 #include <TPrsStd_AISPresentation.hxx>
1881 //=======================================================================
1882 //function : OCC1029_AISTransparency 
1883 //purpose  : OCC1029_AISTransparency  (DOC,entry,[real])
1884 //=======================================================================
1885
1886 static Standard_Integer OCC1029_AISTransparency (Draw_Interpretor& di,
1887                                              Standard_Integer nb, 
1888                                              const char ** arg) 
1889 {
1890   if (nb >= 3 ) {     
1891     Handle(TDocStd_Document) D;
1892     if (!DDocStd::GetDocument(arg[1],D)) return 1;  
1893     TDF_Label L;
1894     if (!DDF::FindLabel(D->GetData(),arg[2],L)) return 1;  
1895
1896     Handle(TPrsStd_AISViewer) viewer;
1897     if( !TPrsStd_AISViewer::Find(L, viewer) ) return 1;  
1898
1899     Handle(TPrsStd_AISPresentation) prs;
1900     if(L.FindAttribute( TPrsStd_AISPresentation::GetID(), prs) ) {   
1901       if( nb == 4 ) {
1902         prs->SetTransparency(atof(arg[3]));
1903         TPrsStd_AISViewer::Update(L);
1904       }
1905       else {
1906          di << "Transparency = " << prs->Transparency() << "\n";
1907       }
1908       return 0;
1909     }
1910   }
1911   di << arg[0] << " : Error" << "\n";
1912   return 1;
1913 }
1914
1915 //=======================================================================
1916 //function : OCC1030_AISColor 
1917 //purpose  : OCC1030_AISColor (DOC,entry,[color])
1918 //=======================================================================
1919
1920 static Standard_Integer OCC1030_AISColor (Draw_Interpretor& di,
1921                                       Standard_Integer nb, 
1922                                       const char ** arg) 
1923 {
1924   if (nb >= 3) {     
1925     Handle(TDocStd_Document) D;
1926     if (!DDocStd::GetDocument(arg[1],D)) return 1;  
1927     TDF_Label L;
1928     if (!DDF::FindLabel(D->GetData(),arg[2],L)) return 1;  
1929
1930     Handle(TPrsStd_AISViewer) viewer;
1931     if( !TPrsStd_AISViewer::Find(L, viewer) ) return 1;  
1932
1933     Handle(TPrsStd_AISPresentation) prs;
1934     if(L.FindAttribute( TPrsStd_AISPresentation::GetID(), prs) ) {   
1935       if( nb == 4 ) {
1936         prs->SetColor((Quantity_NameOfColor)atoi(arg[3]));
1937         TPrsStd_AISViewer::Update(L);
1938       }
1939       else
1940          di << "Color = " << prs->Color() << "\n";
1941       return 0; 
1942     }
1943   }
1944   di << arg[0] << " : Error" << "\n";
1945   return 1;
1946 }
1947
1948 //=======================================================================
1949 //function : OCC1031_AISMaterial
1950 //purpose  : OCC1031_AISMaterial (DOC,entry,[material])
1951 //=======================================================================
1952
1953 static Standard_Integer OCC1031_AISMaterial (Draw_Interpretor& di,
1954                                          Standard_Integer nb,
1955                                          const char ** arg)
1956 {
1957   if (nb >= 3) {     
1958     Handle(TDocStd_Document) D;
1959     if (!DDocStd::GetDocument(arg[1],D)) return 1;  
1960     TDF_Label L;
1961     if (!DDF::FindLabel(D->GetData(),arg[2],L)) return 1;  
1962
1963     Handle(TPrsStd_AISViewer) viewer;
1964     if( !TPrsStd_AISViewer::Find(L, viewer) ) return 1;  
1965
1966     Handle(TPrsStd_AISPresentation) prs;
1967     if(L.FindAttribute( TPrsStd_AISPresentation::GetID(), prs) ) {   
1968       if( nb == 4 ) {
1969         prs->SetMaterial((Graphic3d_NameOfMaterial)atoi(arg[3]));
1970         TPrsStd_AISViewer::Update(L);
1971       }
1972       else {
1973          di << "Material = " << prs->Material() << "\n";
1974       }
1975       return 0;
1976     }
1977   }
1978   di << arg[0] << " : Error" << "\n";
1979   return 1;
1980 }
1981
1982 //=======================================================================
1983 //function : OCC1032_AISWidth
1984 //purpose  : OCC1032_AISWidth (DOC,entry,[width])
1985 //=======================================================================
1986
1987 static Standard_Integer OCC1032_AISWidth (Draw_Interpretor& di,
1988                                       Standard_Integer nb, 
1989                                       const char ** arg) 
1990 {
1991   if (nb >= 3) {     
1992     Handle(TDocStd_Document) D;
1993     if (!DDocStd::GetDocument(arg[1],D)) return 1;  
1994     TDF_Label L;
1995     if (!DDF::FindLabel(D->GetData(),arg[2],L)) return 1;  
1996
1997     Handle(TPrsStd_AISViewer) viewer;
1998     if( !TPrsStd_AISViewer::Find(L, viewer) ) return 1;  
1999
2000     Handle(TPrsStd_AISPresentation) prs;
2001     if(L.FindAttribute( TPrsStd_AISPresentation::GetID(), prs) ) {   
2002       if( nb == 4 ) {
2003         prs->SetWidth(atof(arg[3]));
2004         TPrsStd_AISViewer::Update(L);
2005       }
2006       else {
2007          di << "Width = " << prs->Width() << "\n";
2008       }
2009       return 0;
2010     }
2011   }
2012   di << arg[0] << " : Error" << "\n";
2013   return 1;
2014 }
2015
2016 //=======================================================================
2017 //function : OCC1033_AISMode
2018 //purpose  : OCC1033_AISMode (DOC,entry,[mode])
2019 //=======================================================================
2020
2021 static Standard_Integer OCC1033_AISMode (Draw_Interpretor& di,
2022                                      Standard_Integer nb, 
2023                                      const char ** arg) 
2024 {
2025   if (nb >= 3) {     
2026     Handle(TDocStd_Document) D;
2027     if (!DDocStd::GetDocument(arg[1],D)) return 1;  
2028     TDF_Label L;
2029     if (!DDF::FindLabel(D->GetData(),arg[2],L)) return 1;  
2030
2031     Handle(TPrsStd_AISViewer) viewer;
2032     if( !TPrsStd_AISViewer::Find(L, viewer) ) return 1;  
2033
2034     Handle(TPrsStd_AISPresentation) prs;
2035     if(L.FindAttribute( TPrsStd_AISPresentation::GetID(), prs) ) {   
2036       if( nb == 4 ) {
2037         prs->SetMode(atoi(arg[3]));
2038         TPrsStd_AISViewer::Update(L);
2039       }
2040       else {
2041          di << "Mode = " << prs->Mode() << "\n";
2042       }
2043       return 0;
2044     }
2045   }
2046   di << arg[0] << " : Error" << "\n";
2047   return 1;
2048 }
2049
2050 //=======================================================================
2051 //function : OCC1034_AISSelectionMode
2052 //purpose  : OCC1034_AISSelectionMode (DOC,entry,[selectionmode])
2053 //=======================================================================
2054
2055 static Standard_Integer OCC1034_AISSelectionMode (Draw_Interpretor& di,
2056                                               Standard_Integer nb, 
2057                                               const char ** arg) 
2058 {
2059   if (nb >= 3) {     
2060     Handle(TDocStd_Document) D;
2061     if (!DDocStd::GetDocument(arg[1],D)) return 1;  
2062     TDF_Label L;
2063     if (!DDF::FindLabel(D->GetData(),arg[2],L)) return 1;  
2064
2065     Handle(TPrsStd_AISViewer) viewer;
2066     if( !TPrsStd_AISViewer::Find(L, viewer) ) return 1;  
2067
2068     Handle(TPrsStd_AISPresentation) prs;
2069     if(L.FindAttribute( TPrsStd_AISPresentation::GetID(), prs) ) {   
2070       if( nb == 4 ) {
2071         prs->SetSelectionMode(atoi(arg[3]));
2072         TPrsStd_AISViewer::Update(L);
2073       }
2074       else {
2075          di << "SelectionMode = " << prs->SelectionMode() << "\n";
2076       }
2077       return 0;
2078     }
2079   }
2080   di << arg[0] << " : Error" << "\n";
2081   return 1;
2082 }
2083
2084 #include<BRepAlgoAPI_Cut.hxx>
2085 #include<BRepAlgo_Cut.hxx>
2086 //=======================================================================
2087 //function :  OCC1487
2088 //purpose  :
2089 //=======================================================================
2090 static Standard_Integer OCC1487 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2091 {
2092   //if(argc != 5) {
2093   //  cerr << "Usage : " << argv[0] << " CylinderVariant(=1/2) cylinder1 cylinder2 cutshape" << endl;
2094   //  return -1;
2095   //}
2096   if(argc < 5 || argc > 6) {
2097     di << "Usage : " << argv[0] << " CylinderVariant(=1/2) cylinder1 cylinder2 cutshape [BRepAlgoAPI/BRepAlgo = 1/0]" << "\n";
2098     return 1;
2099   }
2100   Standard_Boolean IsBRepAlgoAPI = Standard_True;
2101   if (argc == 6) {
2102     Standard_Integer IsB = atoi(argv[5]);
2103     if (IsB != 1) {
2104       IsBRepAlgoAPI = Standard_False;
2105 //      di << "Error: There is not BRepAlgo_Cut class" << "\n";
2106 //      return 1;
2107     }
2108   }
2109
2110   Standard_Integer CaseNumber = atoi(argv[1]);
2111
2112   //BRepPrimAPI_MakeCylinder o_mc1 (gp_Ax2 (gp_Pnt(0,-50,140), gp_Dir(1,0,0)), 50,1000);
2113   gp_Dir myDir(1,0,0);
2114   gp_Pnt myPnt(0,-50,140);
2115   gp_Ax2 myAx2(myPnt, myDir);
2116   BRepPrimAPI_MakeCylinder o_mc1 (myAx2, 50,1000);
2117
2118   TopoDS_Shape cyl1 = o_mc1.Shape();
2119
2120   TopoDS_Shape cyl2;
2121   TopoDS_Shape o_cut_shape;
2122   if (CaseNumber == 1) {
2123     //BRepPrimAPI_MakeCylinder o_mc2 (gp_Ax2 (gp_Pnt(21.65064, -50.0, 127.5),gp_Dir(-sin(M_PI/3), 0.0, 0.5)), 5, 150);
2124     gp_Dir myDir_mc2(-sin(M_PI/3), 0.0, 0.5);
2125     gp_Pnt myPnt_mc2(21.65064, -50.0, 127.5);
2126     gp_Ax2 myAx2_mc2(myPnt_mc2, myDir_mc2);
2127     BRepPrimAPI_MakeCylinder o_mc2 (myAx2_mc2, 5, 150);
2128
2129     cyl2 = o_mc2.Shape();
2130 //#if ! defined(BRepAlgoAPI_def01)
2131 //    o_cut_shape = BRepAlgoAPI_Cut (o_mc1.Solid (), o_mc2.Solid ());
2132 //#else
2133 //    o_cut_shape = BRepAlgo_Cut (o_mc1.Solid (), o_mc2.Solid ());
2134 //#endif
2135     if (IsBRepAlgoAPI) {
2136       di << "o_cut_shape = BRepAlgoAPI_Cut (o_mc1.Solid (), o_mc2.Solid ())" <<"\n";
2137       o_cut_shape = BRepAlgoAPI_Cut (o_mc1.Solid (), o_mc2.Solid ());
2138     } else {
2139       di << "o_cut_shape = BRepAlgo_Cut (o_mc1.Solid (), o_mc2.Solid ())" <<"\n";
2140       o_cut_shape = BRepAlgo_Cut (o_mc1.Solid (), o_mc2.Solid ());
2141     }
2142   } else {
2143     //BRepPrimAPI_MakeCylinder o_mc2 (gp_Ax2 (gp_Pnt(978.34936, -50.0, 127.5),gp_Dir(sin(M_PI/3), 0.0, 0.5)), 5, 150);
2144     gp_Dir myDir_mc2(sin(M_PI/3), 0.0, 0.5);
2145     gp_Pnt myPnt_mc2(978.34936, -50.0, 127.5);
2146     gp_Ax2 myAx2_mc2(myPnt_mc2, myDir_mc2);
2147     BRepPrimAPI_MakeCylinder o_mc2 (myAx2_mc2, 5, 150);
2148
2149     cyl2 = o_mc2.Shape();
2150 //#if ! defined(BRepAlgoAPI_def01)
2151 //    o_cut_shape = BRepAlgoAPI_Cut (o_mc1.Solid (), o_mc2.Solid ());
2152 //#else
2153 //    o_cut_shape = BRepAlgo_Cut (o_mc1.Solid (), o_mc2.Solid ());
2154 //#endif
2155     if (IsBRepAlgoAPI) {
2156       di << "o_cut_shape = BRepAlgoAPI_Cut (o_mc1.Solid (), o_mc2.Solid ())" <<"\n";
2157       o_cut_shape = BRepAlgoAPI_Cut (o_mc1.Solid (), o_mc2.Solid ());
2158     } else {
2159       di << "o_cut_shape = BRepAlgo_Cut (o_mc1.Solid (), o_mc2.Solid ())" <<"\n";
2160       o_cut_shape = BRepAlgo_Cut (o_mc1.Solid (), o_mc2.Solid ());
2161     }
2162   }
2163
2164   DBRep::Set(argv[2],cyl1);
2165   DBRep::Set(argv[3],cyl2);
2166   DBRep::Set(argv[4],o_cut_shape);
2167
2168   return 0;
2169 }
2170
2171 #include<TopTools_ListIteratorOfListOfShape.hxx>
2172 #include<BRepFilletAPI_MakeFillet.hxx>
2173 //=======================================================================
2174 //function :  OCC1077
2175 //purpose  :
2176 //=======================================================================
2177 TopoDS_Shape OCC1077_boolbl(BRepAlgoAPI_BooleanOperation& aBoolenaOperation,const Standard_Real aRadius)
2178 {
2179   Standard_Real t3d = 1.e-4;
2180   Standard_Real t2d = 1.e-5;
2181   Standard_Real ta  = 1.e-2;
2182   Standard_Real fl  = 1.e-3;
2183   Standard_Real tapp_angle = 1.e-2;
2184   GeomAbs_Shape blend_cont = GeomAbs_C1;
2185
2186   TopoDS_Shape ShapeCut = aBoolenaOperation.Shape();
2187
2188 //#ifdef OCC40 
2189 //  Handle_TopOpeBRepBuild_HBuilder build = aBoolenaOperation.Builder();
2190 //#endif 
2191   TopTools_ListIteratorOfListOfShape its;
2192
2193   TopoDS_Compound result;
2194   BRep_Builder B;
2195   B.MakeCompound(result);
2196
2197   TopExp_Explorer ex;
2198   for (ex.Init(ShapeCut, TopAbs_SOLID); ex.More(); ex.Next())
2199     {
2200       const TopoDS_Shape& cutsol = ex.Current();
2201
2202       BRepFilletAPI_MakeFillet fill(cutsol);
2203       fill.SetParams(ta, t3d, t2d, t3d, t2d, fl);
2204       fill.SetContinuity(blend_cont, tapp_angle);
2205 //#ifdef OCC40
2206 //      its = build->Section();
2207 //#else //OCC40DEV
2208       its = aBoolenaOperation.SectionEdges();
2209 //#endif
2210       while (its.More())
2211         {
2212           TopoDS_Edge E = TopoDS::Edge(its.Value());
2213           fill.Add(aRadius, E);
2214           its.Next();
2215         }
2216
2217       fill.Build();
2218       if (fill.IsDone())
2219         {
2220           B.Add(result, fill.Shape());
2221         }
2222       else
2223         {
2224           B.Add(result, cutsol);
2225         }
2226     }
2227   return result;
2228 }
2229
2230 TopoDS_Shape OCC1077_cut_blend(const TopoDS_Shape& aShapeToCut, const TopoDS_Shape& aTool, const Standard_Real aRadius)
2231 {
2232   //return OCC1077_boolbl(BRepAlgoAPI_Cut(aShapeToCut, aTool),aRadius);
2233   BRepAlgoAPI_Cut aCut(aShapeToCut, aTool);
2234   return OCC1077_boolbl(aCut,aRadius);
2235 }
2236
2237 //TopoDS_Shape OCC1077_common_blend(const TopoDS_Shape& aShape1, const TopoDS_Shape& aShape2, const Standard_Real aRadius)
2238 //{
2239 //  return OCC1077_boolbl(BRepAlgoAPI_Common(aShape1, aShape2),aRadius);
2240 //}
2241
2242 TopoDS_Shape OCC1077_Bug()
2243 {
2244   TopoDS_Shape theBox = BRepPrimAPI_MakeBox(gp_Pnt(-5, - 5, - 5), 10, 10, 10).Shape();
2245   TopoDS_Shape theSphere = BRepPrimAPI_MakeSphere(7).Shape();
2246
2247   TopoDS_Shape theCommon = BRepAlgoAPI_Common(theBox,theSphere);
2248   TopoDS_Shape theCylinder1 = BRepPrimAPI_MakeCylinder(gp_Ax2(gp_Pnt(0, 0, - 10),
2249                                                               gp_Dir(0, 0, 1)), 3, 20).Shape();
2250   TopoDS_Shape theCylinder2 = BRepPrimAPI_MakeCylinder(gp_Ax2(gp_Pnt(-10, 0, 0),
2251                                                               gp_Dir(1, 0, 0)), 3, 20).Shape();
2252   TopoDS_Shape theCylinder3 = BRepPrimAPI_MakeCylinder(gp_Ax2(gp_Pnt(0, - 10, 0),
2253                                                               gp_Dir(0, 1, 0)), 3, 20).Shape();
2254   TopoDS_Shape theTmp1 = OCC1077_cut_blend(theCommon,theCylinder1,0.7);
2255   Handle_ShapeFix_Shape fixer = new ShapeFix_Shape(theTmp1);
2256   fixer->Perform();
2257   theTmp1 = fixer->Shape();
2258   TopoDS_Shape theTmp2 = OCC1077_cut_blend(theTmp1,theCylinder2,0.7);
2259   fixer->Init(theTmp2);
2260   fixer->Perform();
2261   theTmp2 = fixer->Shape();
2262   TopoDS_Shape theResult = OCC1077_cut_blend(theTmp2,theCylinder3,0.7);
2263   fixer->Init(theResult);
2264   fixer->Perform();
2265   theResult = fixer->Shape();
2266   return theResult;
2267 }
2268
2269 static Standard_Integer OCC1077 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2270 {
2271   if(argc < 1 || argc > 2) {
2272     di << "Usage : " << argv[0] << " result" << "\n";
2273     return 1;
2274   }
2275
2276   TopoDS_Shape S = OCC1077_Bug();
2277   DBRep::Set(argv[1],S);
2278
2279   return 0;
2280 }
2281
2282 //////////////////////////////////////////////////////////////
2283 /*!
2284  * Compute uniform distribution of points using GCPnts_UniformAbscissa
2285  */
2286 //////////////////////////////////////////////////////////////
2287 static Standard_Integer OCC5739_UniAbs (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2288 {
2289   if (argc < 4)
2290   {
2291     di << "Usage : " << argv[0] << " name shape step" << "\n";
2292     return 1;
2293   }
2294   const char *name = argv[1];
2295   Adaptor3d_Curve *adapCurve=NULL;
2296   Handle(Geom_Curve) curve = DrawTrSurf::GetCurve(argv[2]);
2297   if (!curve.IsNull())
2298     adapCurve = new GeomAdaptor_Curve(curve);
2299   else
2300   {
2301     TopoDS_Shape wire = DBRep::Get(argv[2]);
2302     if (wire.IsNull() || wire.ShapeType() != TopAbs_WIRE)
2303     {
2304       di << argv[0] <<" Faulty : incorrect 1st parameter, curve or wire expected"<<"\n";
2305       return 1;
2306     }
2307     adapCurve = new BRepAdaptor_CompCurve(TopoDS::Wire(wire));
2308   }
2309   double step = atof(argv[3]);
2310   GCPnts_UniformAbscissa aUni(*adapCurve, step);
2311   int res;
2312   if (!aUni.IsDone())
2313   {
2314     di << argv[0] <<" : fail"<<"\n";
2315     res = 1;
2316   }
2317   else
2318   {
2319     int i, np = aUni.NbPoints();
2320     for (i=0; i < np; i++)
2321     {
2322       double par = aUni.Parameter(i+1);
2323       gp_Pnt p = adapCurve->Value(par);
2324       char n[20], *pname=n;
2325       sprintf(n,"%s_%d",name,i+1);
2326       DrawTrSurf::Set(pname,p);
2327       di<<pname<<" ";
2328     }
2329     res = 0;
2330   }
2331   delete adapCurve;
2332   return res;
2333 }
2334
2335 static Standard_Integer OCC6046 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2336 {
2337   if (argc != 3)
2338   {
2339     di << "Usage : " << argv[0] << " nb_of_vectors size" << "\n";
2340     return 1;
2341   }
2342
2343   Standard_Integer nb = atoi(argv[1]);
2344   Standard_Integer sz = atoi(argv[2]);
2345   Standard_Real val = 10;
2346   math_Vector **pv = new math_Vector *[nb];
2347
2348   di<<"creating "<<nb<<" vectors "<<sz<<" elements each..."<<"\n";
2349   Standard_Integer i;
2350   for (i=0; i < nb; i++) {
2351     pv[i] = new math_Vector (1, sz, val);
2352     if ((i % (nb/10)) == 0) {
2353       di<<" "<<i;
2354       //cout.flush();
2355       di<<"\n";
2356     }
2357   }
2358   di<<" done"<<"\n";
2359   di<<"deleting them ..."<<"\n";
2360   for (i=0; i < nb; i++) {
2361     delete pv[i];
2362     if ((i % (nb/10)) == 0) {
2363       di<<" "<<i;
2364       //cout.flush();
2365       di<<"\n";
2366     }
2367   }
2368   di<<" done"<<"\n";
2369
2370   delete [] pv;
2371
2372   return 0;
2373 }
2374
2375 static Standard_Integer OCC5698 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2376 {
2377   if (argc != 2)
2378   {
2379     di << "Usage : " << argv[0] << " wire" << "\n";
2380     return 1;
2381   }
2382   TopoDS_Shape shape = DBRep::Get(argv[1],TopAbs_WIRE);
2383   if (shape.IsNull())
2384     return 1;
2385   TopoDS_Wire wire = TopoDS::Wire(shape);
2386   // create curve parameterised by curvilinear distance
2387   BRepAdaptor_CompCurve curve(wire,Standard_True);
2388   Standard_Real length = curve.LastParameter();
2389   Standard_Real need_length = length/2;
2390   gp_Pnt pnt;
2391   curve.D0(need_length,pnt);
2392   // create check_curve parameterised in a general way
2393   BRepAdaptor_CompCurve check_curve(wire);
2394   Standard_Real check_par =
2395     GCPnts_AbscissaPoint(check_curve, need_length, 0).Parameter();
2396   gp_Pnt check_pnt;
2397   check_curve.D0(check_par,check_pnt);
2398   // check that points are coinsiding
2399   Standard_Real error_dist = pnt.Distance(check_pnt);
2400   if (error_dist > Precision::Confusion()) {
2401     //cout.precision(3);
2402     di<<"error_dist = "<<error_dist<<
2403       "  ( "<<error_dist/need_length*100<<" %)"<<"\n";
2404     return 0;
2405   }
2406   di<<"OK"<<"\n";
2407   return 0;
2408 }
2409
2410 static char sarr[2000];
2411 static int si=1;
2412 static int StackOverflow(int i = -1)
2413 {
2414   char arr[2000];
2415   if (si == 1) {
2416     si = 0;
2417     memcpy(arr,sarr,2000);
2418     arr[1999]=0;
2419     int n = strlen(arr), s=0;
2420     while (n--)
2421       s += StackOverflow(i-1);
2422     return i + s + StackOverflow(i-1);
2423   }
2424   else if (i != 0) {
2425     return i + StackOverflow(i-1);
2426   }
2427   si = 1;
2428   return i;
2429 }
2430
2431 #ifdef WNT
2432 // this code does not work with optimize mode on Windows
2433 #pragma optimize( "", off )
2434 #endif
2435 static Standard_Integer OCC6143 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2436 {
2437   if (argc != 1)
2438     {
2439       di << "Usage : " << argv[0] << "\n";
2440       return 1;
2441     }
2442   Standard_Boolean Succes;
2443   
2444   Succes = Standard_True;
2445   //OSD::SetSignal();
2446
2447   {//==== Test Divide ByZero (Integer) ========================================
2448     try{
2449       OCC_CATCH_SIGNALS
2450       di << "(Integer) Divide By Zero...";
2451       //cout.flush();
2452       di << "\n";
2453       Standard_Integer res, a =4, b = 0 ;
2454       res = a / b;
2455       di << " 4 / 0 = " << res << "  Does not Caught... KO"<< "\n";
2456       Succes = Standard_False;
2457     }
2458 #if defined(SOLARIS) || defined(WNT)
2459     catch(Standard_DivideByZero)
2460 #else
2461     catch(Standard_NumericError)
2462 #endif
2463     {
2464       di << " Ok"<< "\n";
2465     }
2466     catch(Standard_Failure) {
2467       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2468       di << " Caught (";
2469       di << Standard_Failure::Caught()->GetMessageString();
2470       di << ")... KO" << "\n";
2471       Succes = Standard_False;
2472     }
2473 #ifndef NO_CXX_EXCEPTION
2474     // this case tests if (...) supersedes (Standard_*),
2475     // the normal behaviour is not
2476     catch(...) {
2477       di<<" unknown exception... (But) Ok"<<"\n";
2478     }
2479 #endif
2480   }
2481
2482   {//==== Test Divide ByZero (Real) ===========================================
2483     try{
2484       OCC_CATCH_SIGNALS
2485       di << "(Real) Divide By Zero...";
2486       //cout.flush();
2487       di << "\n";
2488       Standard_Real res, a= 4.0, b=0.0;
2489       res = a / b;
2490       di << " 4.0 / 0.0 = " << res << "  Does not Caught... KO"<< "\n";
2491       Succes = Standard_False;
2492     }
2493 #if defined(SOLARIS) || defined(WNT)
2494     catch(Standard_DivideByZero)
2495 #else
2496     catch(Standard_NumericError)
2497 #endif
2498     {
2499       di << " Ok"<< "\n";
2500     }
2501     catch(Standard_Failure) {
2502       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2503       di << " Caught (";
2504       di << Standard_Failure::Caught()->GetMessageString();
2505       di << ")... KO" << "\n";
2506       Succes = Standard_False;
2507     }
2508   }
2509
2510   {//==== Test Overflow (Integer) =============================================
2511     try{
2512       OCC_CATCH_SIGNALS
2513       di << "(Integer) Overflow...";
2514       //cout.flush();
2515       di << "\n";
2516       Standard_Integer res, i=IntegerLast();
2517       res = i + 1;
2518       //++++ cout << " -- "<<res<<"="<<i<<"+1   Does not Caught... KO"<< endl;
2519       //++++ Succes = Standard_False;
2520       di << " "<<res<<"="<<i<<"+1  Does not Caught... (But) Ok"<< "\n";
2521     }
2522     catch(Standard_Overflow) {
2523       di << " Ok"<< "\n";
2524     }
2525     catch(Standard_Failure) {
2526       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2527       di << " Caught (";
2528       di << Standard_Failure::Caught()->GetMessageString();
2529       di << ")... KO" << "\n";
2530       Succes = Standard_False;
2531     }
2532   }
2533
2534   {//==== Test Overflow (Real) ================================================ 
2535     try{
2536       OCC_CATCH_SIGNALS
2537       di << "(Real) Overflow...";
2538       //cout.flush();
2539       di << "\n";
2540       Standard_Real res, r=RealLast();
2541       res = r * r;
2542       
2543       sin(1.); //this function tests FPU flags and raises signal (tested on LINUX).
2544
2545       di << "-- "<<res<<"="<<r<<"*"<<r<<"   Does not Caught... KO"<< "\n";
2546       Succes = Standard_False;
2547     }
2548 #if defined(SOLARIS) || defined(WNT)
2549     catch(Standard_Overflow)
2550 #else
2551     catch(Standard_NumericError)
2552 #endif
2553     {
2554       di << " Ok"<< "\n";
2555     }
2556     catch(Standard_Failure) {
2557       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2558       di << " Caught (";
2559       di << Standard_Failure::Caught()->GetMessageString();
2560       di << ")... KO" << "\n";
2561       Succes = Standard_False;
2562     }
2563   }
2564
2565   {//==== Test Underflow (Real) ===============================================
2566     try{
2567       OCC_CATCH_SIGNALS
2568       di << "(Real) Underflow";
2569       //cout.flush();
2570       di << "\n";
2571       Standard_Real res, r=1.0e-308;
2572       res = r * r;
2573       //res = res + 1.;
2574       //++++ cout<<"-- "<<res<<"="<<r<<"*"<<r<<"   Does not Caught... KO"<<endl;
2575       //++++ Succes = Standard_False;
2576       di<<" -- "<<res<<"="<<r<<"*"<<r<<"   Does not Caught... (But) Ok"<<"\n";
2577     }
2578 #if defined(SOLARIS) || defined(WNT)
2579     catch(Standard_Underflow)
2580 #else
2581     catch(Standard_NumericError)
2582 #endif
2583     {
2584       di << " Ok"<< "\n";
2585     }
2586     catch(Standard_Failure) {
2587       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2588       di << " Caught (";
2589       di << Standard_Failure::Caught()->GetMessageString();
2590       di << ")... KO" << "\n";
2591       Succes = Standard_False;
2592     }
2593   }
2594
2595   {//==== Test Invalid Operation (Real) ===============================================
2596     try{
2597       OCC_CATCH_SIGNALS
2598       di << "(Real) Invalid Operation...";
2599       //cout.flush();
2600       di << "\n";
2601       Standard_Real res, r=-1;
2602       res = sqrt(r);
2603       di<<" "<<res<<"=sqrt("<<r<<")  Does not Caught... KO"<<"\n";
2604       Succes = Standard_False;
2605     }
2606     catch(Standard_NumericError) {
2607       di << " Ok"<< "\n";
2608     }
2609     catch(Standard_Failure) {
2610       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2611       di << " Caught (";
2612       di << Standard_Failure::Caught()->GetMessageString();
2613       di << ")... KO" << "\n";
2614       Succes = Standard_False;
2615     }
2616   }
2617
2618   {//==== Test Access Violation ===============================================
2619     try {
2620       OCC_CATCH_SIGNALS
2621       di << "Segmentation Fault...";
2622       //cout.flush();
2623       di << "\n";
2624       int* pint=NULL;
2625       *pint = 4;
2626       di << "  Does not Caught... KO"<<"\n";
2627       Succes = Standard_False;
2628     }
2629 #ifdef WNT
2630     catch(OSD_Exception_ACCESS_VIOLATION)
2631 #else
2632     catch(OSD_SIGSEGV)
2633 #endif
2634     {
2635       di << " Ok"<< "\n";
2636     } catch(Standard_Failure) {
2637       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2638       di << " Caught (";
2639       di << Standard_Failure::Caught()->GetMessageString();
2640       di << ")... KO" << "\n";
2641       Succes = Standard_False;
2642     }
2643   }
2644
2645 #ifdef WNT
2646   {//==== Test Stack Overflow ===============================================
2647     try {
2648       OCC_CATCH_SIGNALS
2649       di << "Stack Overflow...";
2650       //cout.flush();
2651       di << "\n";
2652       StackOverflow();
2653       di << "  Does not Caught... KO"<<"\n";
2654       Succes = Standard_False;
2655     }
2656     catch(OSD_Exception_STACK_OVERFLOW) {
2657       di << " Ok"<< "\n";
2658     }
2659     catch(Standard_Failure) {
2660       //cout << " Caught (" << Standard_Failure::Caught() << ")... KO" << endl;
2661       di << " Caught (";
2662       di << Standard_Failure::Caught()->GetMessageString();
2663       di << ")... KO" << "\n";
2664       Succes = Standard_False;
2665     }
2666   }
2667 #endif
2668
2669  if(Succes) {
2670    di << "TestExcept: Successfull completion" << "\n";
2671  } else {
2672    di << "TestExcept: failure" << "\n";
2673  }
2674
2675   return 0;
2676 }
2677 #ifdef WNT
2678 #pragma optimize( "", on )
2679 #endif
2680
2681 static TopoDS_Compound AddTestStructure(int nCount_)
2682 {
2683   BRep_Builder B;
2684   int nCount=nCount_;
2685   TopoDS_Compound C;
2686   B.MakeCompound(C);
2687   BRepPrimAPI_MakeBox mkBox(1.0, 2.0, 3.0);
2688   for (int i=0; i<nCount; i++) {
2689     for (int j=0; j<nCount; j++) {
2690       gp_Trsf trsf;
2691       trsf.SetTranslationPart(gp_Vec(5.0*i, 05.0*j, 0.0));
2692       TopLoc_Location topLoc(trsf);
2693       TopoDS_Shape tempShape=mkBox.Shape().Located(topLoc);
2694       B.Add(C, tempShape);
2695     }
2696   }
2697   return C;
2698 }
2699
2700 static Standard_Integer OCC7141 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2701 {
2702   int nCount = 10;
2703
2704   if (argc > 2)
2705     {
2706       di << "Usage : " << argv[0] << " [nCount]" << "\n";
2707       return 1;
2708     }
2709
2710   if (argc > 1)
2711     nCount = atoi(argv[1]);
2712   STEPCAFControl_Writer writer;
2713   Handle_TDocStd_Document document;
2714   document = new TDocStd_Document("Pace Test-StepExporter-");
2715   Handle_XCAFDoc_ShapeTool shapeTool;
2716   shapeTool = XCAFDoc_DocumentTool::ShapeTool(document->Main());
2717   shapeTool->AddShape(AddTestStructure(nCount), Standard_True);
2718   STEPControl_StepModelType mode = STEPControl_AsIs;
2719   if (!Interface_Static::SetIVal("write.step.assembly",1)) { //assembly mode
2720     di << "Failed to set assembly mode for step data\n" << "\n";
2721     return 0;
2722   }
2723   try {
2724     OCC_CATCH_SIGNALS
2725     if( writer.Transfer(document, mode)) {
2726       IFSelect_ReturnStatus stat =
2727         writer.Write("TestExportStructure.step");
2728     }
2729   }
2730   catch(OSD_Exception_STACK_OVERFLOW) {
2731     di << "Failed : STACK OVERFLOW\n" << "\n";
2732   }
2733   catch (Standard_Failure) {
2734     di << "Failed :\n" << "\n";
2735     //cout << Standard_Failure::Caught() << endl;
2736     di << Standard_Failure::Caught()->GetMessageString();
2737   }
2738   di << argv[0] << " : Finish" << "\n";
2739   
2740   return 0;
2741 }
2742
2743 static Standard_Integer OCC7372 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2744 {
2745   if (argc != 1)
2746     {
2747       di << "Usage : " << argv[0] << "\n";
2748       return 1;
2749     }
2750   
2751   // 1. Create an array of points
2752   Handle(TColgp_HArray1OfPnt2d) ap = new TColgp_HArray1OfPnt2d(1,5);
2753   ap->SetValue(1,gp_Pnt2d(100.0,0.0));
2754   ap->SetValue(2,gp_Pnt2d(100.0,100.0));
2755   ap->SetValue(3,gp_Pnt2d(0.0,100.0));
2756   ap->SetValue(4,gp_Pnt2d(0.0,0.0));
2757   ap->SetValue(5,gp_Pnt2d(50.0,-50.0));
2758
2759   // 2. Create a periodic bspline through these 5 points
2760   Geom2dAPI_Interpolate intp(ap,Standard_True,1e-6);
2761   intp.Perform();
2762   Handle(Geom2d_BSplineCurve) bspline1 = intp.Curve();
2763
2764   // 3. Increase degree of curve from 3 to 8
2765   bspline1->IncreaseDegree(8); // Increase degree to demonstrate the error
2766   Standard_CString CString1 = "BSplineCurve";
2767   DrawTrSurf::Set(CString1,bspline1);
2768
2769   // 4. Convers BSpline curve to Bezier segments
2770   Geom2dConvert_BSplineCurveToBezierCurve bc(bspline1);
2771
2772   // 5. Test the result of conversion
2773   TCollection_AsciiString aRName;
2774   for(Standard_Integer i = 1; i <= bc.NbArcs(); i++) {
2775     Handle(Geom2d_BezierCurve) arc = bc.Arc(i);
2776     aRName="segment_";
2777     aRName=aRName+TCollection_AsciiString(i);
2778     Standard_CString aRNameStr = aRName.ToCString();
2779     DrawTrSurf::Set(aRNameStr,arc);
2780     di << aRNameStr << " ";
2781   }
2782
2783   return 0;
2784 }
2785
2786 static Standard_Integer OCC8169 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2787 {
2788   if (argc != 4)
2789   {
2790     di << "Usage : " << argv[0] << " edge1 edge2 plane" << "\n";
2791     return 1;
2792   }
2793   TopoDS_Edge theEdge1 = TopoDS::Edge(DBRep::Get(argv[1],TopAbs_EDGE));
2794   if (theEdge1.IsNull()) {
2795     di << "Invalid input shape " << argv[1] << "\n";
2796     return 1;
2797   }
2798   TopoDS_Edge theEdge2 = TopoDS::Edge(DBRep::Get(argv[2],TopAbs_EDGE));
2799   if (theEdge2.IsNull()) {
2800     di << "Invalid input shape " << argv[2] << "\n";
2801     return 1;
2802   }
2803   TopoDS_Face theFace = TopoDS::Face(DBRep::Get(argv[3],TopAbs_FACE));
2804   if (theFace.IsNull()) {
2805     di << "Invalid input shape " << argv[3] << "\n";
2806     return 1;
2807   }
2808
2809   Handle(Geom_Surface) thePlane = BRep_Tool::Surface(theFace);
2810
2811   Standard_Real aConfusion = Precision::Confusion();
2812   Standard_Real aP1first, aP1last, aP2first, aP2last;
2813
2814   Handle(Geom_Curve) aCurve1 = BRep_Tool::Curve(theEdge1, aP1first, aP1last);
2815   Handle(Geom_Curve) aCurve2 = BRep_Tool::Curve(theEdge2, aP2first, aP2last);
2816   Handle(Geom2d_Curve) aCurve2d1 = GeomProjLib::Curve2d(aCurve1, aP1first, aP1last, thePlane);
2817   Handle(Geom2d_Curve) aCurve2d2 = GeomProjLib::Curve2d(aCurve2, aP2first, aP2last, thePlane);
2818
2819   Geom2dAPI_InterCurveCurve anInter(aCurve2d1, aCurve2d2, aConfusion);
2820
2821   Standard_Integer NbPoints = anInter.NbPoints();
2822
2823   di << "NbPoints = " << NbPoints << "\n" ;
2824
2825   if (NbPoints > 0) {
2826     Standard_Integer i;
2827     for (i=1; i<=NbPoints; i++) {
2828       gp_Pnt2d aPi = anInter.Point(i);
2829       di << "Point.X(" << i << ") = " << aPi.X() << "   " << "Point.Y(" << i << ") = " << aPi.Y() << "\n" ;
2830     }
2831   }
2832
2833   Standard_Integer NbSegments = anInter.NbSegments();
2834
2835   di << "\nNbSegments = " << NbSegments << "\n" ;
2836
2837   if (NbSegments > 0) {
2838     IntRes2d_IntersectionSegment aSegment = anInter.Intersector().Segment(1);
2839
2840     gp_Pnt2d aP1 = aCurve2d1->Value(aSegment.FirstPoint().ParamOnFirst());
2841     gp_Pnt2d aP2 = aCurve2d2->Value(aSegment.FirstPoint().ParamOnSecond());
2842   
2843     Standard_Real aDist = aP1.Distance(aP2);
2844   
2845     di << "aP1.X() = " << aP1.X() << "   " << "aP1.Y() = " << aP1.Y() << "\n" ;
2846     di << "aP2.X() = " << aP2.X() << "   " << "aP2.Y() = " << aP2.Y() << "\n" ;
2847
2848     di << "Distance = " << aDist << "\n" ;
2849
2850     di << "Confusion = " << aConfusion << "\n" ;
2851
2852     if (aDist > aConfusion) {
2853       di << "\n" << argv[0] << " Faulty" << "\n" ;
2854     } else {
2855       di << "\n" << argv[0] << " OK" << "\n" ;
2856     }
2857   } else {
2858     di << "\n" << argv[0] << " OK" << "\n" ;
2859   }
2860
2861   return 0;
2862 }
2863 static Standard_Integer OCC10138 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
2864 {
2865   if (argc != 3)
2866   {
2867     di << "Usage : " << argv[0] << " lower upper" << "\n";
2868     return 1;
2869   }
2870
2871   Standard_Integer LOWER = atoi(argv[1]);
2872   Standard_Integer UPPER = atoi(argv[2]);
2873
2874   //! 0. Create an empty document with several test labels
2875   Handle(TDocStd_Document) doc = new TDocStd_Document("XmlOcaf");
2876   doc->SetUndoLimit(100);
2877   TDF_Label main_label = doc->Main();
2878   TDF_Label label1 = main_label.FindChild(1, Standard_True);
2879   TDF_Label label2 = main_label.FindChild(2, Standard_True);
2880   
2881   //! 1. Set/Get OCAF attribute
2882   doc->OpenCommand();
2883   TDataStd_RealArray::Set(label1, LOWER, UPPER);
2884   Handle(TDataStd_RealArray) array;
2885   if (label1.FindAttribute(TDataStd_RealArray::GetID(), array) && 
2886       array->Lower() == LOWER && array->Upper() == UPPER)
2887     cout<<"1: OK"<<endl;
2888   else
2889   {
2890     cout<<"1: Failed.."<<endl;
2891     return 1;
2892   }
2893   doc->CommitCommand();
2894
2895   //! 2. Set/Get value
2896   doc->OpenCommand();
2897   Standard_Integer i;
2898   for (i = LOWER; i <= UPPER; i++)
2899     array->SetValue(i, i);
2900   for (i = LOWER; i <= UPPER; i++)
2901   {  
2902     if (array->Value(i) != i)
2903     {
2904       cout<<"2: Failed.."<<endl;
2905       return 2;
2906     }
2907   }
2908   cout<<"2: OK"<<endl;
2909   doc->CommitCommand();
2910
2911   //! 3. Re-init the array
2912   doc->OpenCommand();
2913   array->Init(LOWER + 2, UPPER + 4);
2914   if (array->Lower() != LOWER + 2 && array->Upper() != UPPER + 4)
2915   {
2916     cout<<"3: Failed.."<<endl;
2917     return 3;
2918   }
2919   for (i = LOWER + 2; i <= UPPER + 4; i++)
2920     array->SetValue(i, i);
2921   for (i = LOWER + 2; i <= UPPER + 4; i++)
2922   {  
2923     if (array->Value(i) != i)
2924     {
2925       cout<<"3: Failed.."<<endl;
2926       return 3;
2927     }
2928   }
2929   cout<<"3: OK"<<endl;
2930   doc->CommitCommand();
2931
2932   //! 4. Change array
2933   doc->OpenCommand();
2934   Handle(TColStd_HArray1OfReal) arr = new TColStd_HArray1OfReal(LOWER + 5, UPPER + 5);
2935   for (i = LOWER + 5; i <= UPPER + 5; i++)
2936     arr->SetValue(i, i);
2937   array->ChangeArray(arr);
2938   for (i = LOWER + 5; i <= UPPER + 5; i++)
2939   {  
2940     if (array->Value(i) != i)
2941     {
2942       cout<<"4: Failed.."<<endl;
2943       return 4;
2944     }
2945   }
2946   cout<<"4: OK"<<endl;
2947   doc->CommitCommand();
2948
2949   //! 5. Copy the array
2950   doc->OpenCommand();
2951   TDF_CopyLabel copier(label1, label2);
2952   copier.Perform();
2953   if (!copier.IsDone())
2954   {
2955     cout<<"5: Failed.."<<endl;
2956     return 5;
2957   }
2958   Handle(TDataStd_RealArray) array2;
2959   if (!label2.FindAttribute(TDataStd_RealArray::GetID(), array2))
2960   {
2961     cout<<"5: Failed.."<<endl;
2962     return 5;
2963   }
2964   for (i = LOWER + 5; i <= UPPER + 5; i++)
2965   {  
2966     if (array->Value(i) != i)
2967     {
2968       cout<<"5: Failed.."<<endl;
2969       return 5;
2970     }
2971   }
2972   cout<<"5: OK"<<endl;
2973   doc->CommitCommand();
2974
2975   //! 6. Undo/Redo
2976   //! 6.a: undoes the 5th action: the copied array should disappear
2977   doc->Undo();
2978   if (!label1.FindAttribute(TDataStd_RealArray::GetID(), array) ||
2979       label2.FindAttribute(TDataStd_RealArray::GetID(), array2))
2980   {
2981     cout<<"6.a: Failed.."<<endl;
2982     return 6;
2983   }
2984   //! 6.b: undoes the 4th action: the array should be changed to (lower+2,upper+4)
2985   doc->Undo();
2986   if (!label1.FindAttribute(TDataStd_RealArray::GetID(), array) || 
2987       array->Lower() != LOWER + 2 ||
2988       array->Upper() != UPPER + 4)
2989   {
2990     cout<<"6.b: Failed.."<<endl;
2991     return 6;
2992   }
2993   for (i = LOWER + 2; i <= UPPER + 4; i++)
2994   {
2995     if (array->Value(i) != i)
2996     {
2997       cout<<"6.b: Failed.."<<endl;
2998       return 6;
2999     }
3000   }
3001   //! 6.c: undoes the 3d action: the array should be changed to (lower,upper)
3002   doc->Undo();
3003   if (!label1.FindAttribute(TDataStd_RealArray::GetID(), array) || 
3004       array->Lower() != LOWER ||
3005       array->Upper() != UPPER)
3006   {
3007     cout<<"6.c: Failed.."<<endl;
3008     return 6;
3009   }
3010   for (i = LOWER; i <= UPPER; i++)
3011   {
3012     if (array->Value(i) != i)
3013     {
3014       cout<<"6.c: Failed.."<<endl;
3015       return 6;
3016     }
3017   }
3018   //! 6.d: undoes and redoes the 2nd action: no change is expected.
3019   doc->Undo();
3020   doc->Redo();
3021   if (!label1.FindAttribute(TDataStd_RealArray::GetID(), array) || 
3022       array->Lower() != LOWER ||
3023       array->Upper() != UPPER)
3024   {
3025     cout<<"6.d: Failed.."<<endl;
3026     return 6;
3027   }
3028   for (i = LOWER; i <= UPPER; i++)
3029   {
3030     if (array->Value(i) != i)
3031     {
3032       cout<<"6.d: Failed.."<<endl;
3033       return 6;
3034     }
3035   }
3036   cout<<"6: OK"<<endl;
3037
3038   //! 7. Re-set the array
3039   doc->OpenCommand();
3040   array = TDataStd_RealArray::Set(label1, LOWER + 1, UPPER + 1);
3041   if (array->Lower() != LOWER + 1 && array->Upper() != UPPER + 1)
3042   {
3043     cout<<"7: Failed.."<<endl;
3044     return 7;
3045   }
3046   for (i = LOWER + 1; i <= UPPER + 1; i++)
3047     array->SetValue(i, i);
3048   for (i = LOWER + 1; i <= UPPER + 1; i++)
3049   {  
3050     if (array->Value(i) != i)
3051     {
3052       cout<<"7: Failed.."<<endl;
3053       return 7;
3054     }
3055   }
3056   cout<<"7: OK"<<endl;
3057   doc->CommitCommand();
3058
3059   //! 8.Test of speed: set LOWER and UPPER equal to great integer number and 
3060   //! measure the time spent by this test.
3061   //! Good luck!
3062
3063   return 0;
3064 }
3065
3066 static Standard_Integer OCC7639 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3067 {
3068   Standard_Boolean IsEvenArgc =  Standard_True;
3069   if (argc % 2 == 0) {
3070     IsEvenArgc =  Standard_True;
3071   } else {
3072     IsEvenArgc =  Standard_False;
3073   }
3074
3075   if (argc < 3 || IsEvenArgc)
3076     {
3077       di << "Usage : " << argv[0] << " index1 value1 ... [indexN valueN]" << "\n";
3078       return 1;
3079     }
3080
3081   Standard_Integer i, aValue, aPosition;
3082   NCollection_Vector<int> vec;
3083   for (i = 0; i < argc - 1; i++) {
3084     i++;
3085     aValue = atoi(argv[i]);
3086     aPosition = atoi(argv[i+1]);
3087     vec.SetValue(aValue, aPosition);
3088   }
3089   NCollection_Vector<int>::Iterator it(vec);
3090   Standard_Integer j;
3091   for (j = 0; it.More(); it.Next(), j++) {
3092     //di << it.Value() << "\n";                                               
3093     di << j << " " << it.Value() << "\n";                                               
3094   }
3095
3096   return 0;
3097 }
3098
3099 static Standard_Integer OCC8797 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3100 {
3101   if (argc != 1) {
3102     di << "Usage : " << argv[0] << "\n";
3103     return 1;
3104   }
3105
3106   gp_Pnt point(0.0,0.0,0.0);
3107
3108   TColgp_Array1OfPnt poles(0,6);
3109   poles(0)=point;
3110
3111   point.SetCoord(1.0,1.0,0.0);
3112   poles(1)=point;
3113
3114   point.SetCoord(2.0,1.0,0.0);
3115   poles(2)=point;
3116
3117   point.SetCoord(3.0,0.0,0.0);
3118   poles(3)=point;
3119
3120   point.SetCoord(4.0,1.0,0.0);
3121   poles(4)=point;
3122
3123   point.SetCoord(5.0,1.0,0.0);
3124   poles(5)=point;
3125
3126   point.SetCoord(6.0,0.0,0.0);
3127   poles(6)=point;
3128
3129   TColStd_Array1OfReal knots(0,2);
3130   knots(0)=0.0;
3131   knots(1)=0.5;
3132   knots(2)=1.0;
3133
3134   TColStd_Array1OfInteger multi(0,2);
3135   multi(0)=4;
3136   multi(1)=3;
3137   multi(2)=4;
3138
3139   Handle(Geom_BSplineCurve) spline = new Geom_BSplineCurve(poles,knots,multi,3);
3140
3141   //length!! 1.
3142   Standard_Real l_abcissa,l_gprop;
3143   GeomAdaptor_Curve adaptor_spline(spline);
3144   GCPnts_AbscissaPoint temp;
3145   l_abcissa=temp.Length(adaptor_spline);
3146   cout<<"Length Spline(abcissa_Pnt): "<<l_abcissa<<endl;
3147
3148   //length!! 2.
3149   TopoDS_Edge edge = BRepBuilderAPI_MakeEdge (spline);
3150   GProp_GProps prop;
3151   BRepGProp::LinearProperties(edge,prop);
3152   l_gprop=prop.Mass();
3153   cout<<"Length Spline(GProp_GProps): "<<l_gprop<<endl;
3154
3155   cout<<"Difference (abcissa_Pnt<->GProp_GProps): "<<l_gprop-l_abcissa<<endl;
3156
3157   return 0;
3158 }
3159
3160 static Standard_Integer OCC7068 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3161 {
3162   if (argc != 1)
3163   {
3164     di << "Usage : " << argv[0] << "\n";
3165     return 1;
3166   }
3167
3168   Handle(AIS_InteractiveContext) AISContext = ViewerTest::GetAISContext();
3169   if(AISContext.IsNull())
3170   {
3171     di << "use 'vinit' command before " << argv[0] << "\n";
3172     return 1;
3173   }
3174
3175   // ObjectsInside
3176   AIS_ListOfInteractive ListOfIO_1;
3177   AISContext->ObjectsInside(ListOfIO_1);
3178   di<< "ObjectsInside = " << ListOfIO_1.Extent() <<"\n";
3179   if (!ListOfIO_1.IsEmpty() ) {
3180     AIS_ListIteratorOfListOfInteractive iter;
3181     for (iter.Initialize(ListOfIO_1); iter.More() ; iter.Next() ) {
3182       Handle(AIS_InteractiveObject) aIO=iter.Value();
3183       di<< GetMapOfAIS().Find1(aIO).ToCString() <<"\n";
3184     }
3185   }
3186
3187   // ObjectsInCollector
3188   AIS_ListOfInteractive ListOfIO_2;
3189   AISContext->ObjectsInCollector(ListOfIO_2);
3190   di<< "ObjectsInCollector = " << ListOfIO_2.Extent() <<"\n";
3191   if (!ListOfIO_2.IsEmpty() ) {
3192     AIS_ListIteratorOfListOfInteractive iter;
3193     for (iter.Initialize(ListOfIO_2); iter.More() ; iter.Next() ) {
3194       Handle(AIS_InteractiveObject) aIO=iter.Value();
3195       di<< GetMapOfAIS().Find1(aIO).ToCString() <<"\n";
3196     }
3197   }
3198
3199   return 0;
3200 }
3201
3202 static Standard_Integer OCC11457 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3203 {
3204   if ((argc < 9) || (((argc-3) % 3) != 0))
3205   {
3206     di << "Usage : " << argv[0] << "polygon lastedge x1 y1 z1 x2 y2 z2 ...\n";
3207     return 1;
3208   }
3209   Standard_Integer i, j, np = (argc-3) / 3;
3210   BRepBuilderAPI_MakePolygon W;
3211   j = 3;
3212   for (i = 1; i <= np; i ++) {
3213     W.Add(gp_Pnt(atof(argv[j]),atof(argv[j+1]),atof(argv[j+2])));
3214     j += 3;
3215   }
3216   W.Close();
3217   DBRep::Set(argv[1],W.Wire());
3218   DBRep::Set(argv[2],W.Edge());
3219   return 0;
3220 }
3221
3222 static Standard_Integer OCC13963 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3223 {
3224   if (argc < 5) {
3225     di << "Usage : " << argv[0] << " ratio origin_x origin_y origin_z\n";
3226     return 1;
3227   }
3228   gp_Ax2 aPln (gp_Pnt(0.,0.,0.),
3229                gp_Dir(1., -1., 0.));
3230   gp_GTrsf aTrf;
3231   aTrf.SetAffinity (aPln, atof(argv[4]));
3232   gp_XYZ aOrigin (atof(argv[1]),atof(argv[2]),atof(argv[3]));
3233   gp_XYZ aResult (aOrigin);
3234   aTrf.Transforms(aResult);
3235   char sbf[512];
3236   sprintf(sbf, "( %8.3f %8.3f %8.3f ) => ( %8.3f %8.3f %8.3f )\n",
3237           aOrigin.X(), aOrigin.Y(), aOrigin.Z(),
3238           aResult.X(), aResult.Y(), aResult.Z());
3239   di<<sbf;
3240   return 0;
3241 }
3242
3243 Standard_Integer OCC14376(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3244 {
3245   if (argc < 2) {
3246     di << "Usage : " << argv[0] << " shape [deflection]\n";
3247     return 1;
3248   }
3249
3250   TopoDS_Shape aShape = DBRep::Get(argv[1]);
3251
3252   if (aShape.IsNull()) {
3253     di<<" Null shape is not allowed";
3254     return 1;
3255   }
3256
3257   Standard_Real aDeflection = 0.45110277533;
3258   if (argc > 2) {
3259     aDeflection = atof(argv[2]);
3260   }
3261   di<<"deflection="<< aDeflection << "\n";
3262
3263   BRepMesh_IncrementalMesh aIMesh(aShape, aDeflection, Standard_False, M_PI/9.);
3264   TopLoc_Location aLocation;
3265   Handle(Poly_Triangulation) aTriang = BRep_Tool::Triangulation(TopoDS::Face(aShape), aLocation);
3266
3267   if(aTriang.IsNull()) {
3268     di << argv[0] << " : Faulty\n" ;
3269   } else {
3270     di << argv[0] << " : OK\n" ;
3271     di<<"NbNodes="<< aTriang->NbNodes()<< "\n";
3272     di<<"NbTriangles="<< aTriang->NbTriangles()<< "\n";
3273   }
3274   return 0;
3275 }
3276
3277 static Standard_Integer OCC15489 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3278 {
3279   if (argc != 4) {
3280     di << "Usage : " << argv[0] << " A B C\n";
3281     return 1;
3282   }
3283   try
3284     {
3285       gp_Lin2d aLin2d (atof(argv[1]),atof(argv[2]),atof(argv[3]));
3286       gp_Pnt2d anOrigin = aLin2d.Location();
3287       di << "X_0 = " << anOrigin.X() << "   Y_0 = " << anOrigin.Y() << "\n" ;
3288     }
3289   catch(Standard_ConstructionError)
3290     {
3291       di << argv[0] << " Exception: Sqrt(A*A + B*B) <= Resolution from gp\n";
3292     }
3293   return 0;
3294 }
3295
3296 static Standard_Integer OCC15755 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
3297 {
3298   if (argc != 3) {
3299     di << "Usage : " << argv[0] << " file shape\n";
3300     return 1;
3301   }
3302
3303   IGESControl_Reader aReader;
3304   aReader.ReadFile(argv[1]);
3305   aReader.SetReadVisible(Standard_True);
3306   aReader.TransferRoots();
3307
3308   Handle(IGESData_IGESModel) model = aReader.IGESModel();
3309   if (model.IsNull()) {
3310     di << "model.IsNull()\n";
3311     return 1;
3312   }
3313   Standard_Integer nb = model->NbEntities();
3314   for (Standard_Integer i = 1; i <= nb; i ++) {
3315     Handle(IGESData_IGESEntity) ent = model->Entity(i);
3316     Handle(TCollection_HAsciiString) name;
3317     name = ent->NameValue();
3318     Standard_CString aStr = name->ToCString();
3319     di << "NameValue = " << aStr << "\n";
3320   }
3321
3322   TopoDS_Shape shape = aReader.OneShape();
3323   DBRep::Set(argv[2],shape);
3324   return 0;
3325 }
3326
3327 // For OCC16782 testing
3328 #include <AppStd_Application.hxx>
3329 #include <TDF_Tool.hxx>
3330 #include <TColStd_HArray1OfInteger.hxx>
3331 // Iterators
3332 #include <TColStd_ListIteratorOfListOfInteger.hxx>
3333 #include <TColStd_ListIteratorOfListOfReal.hxx>
3334 #include <TDataStd_ListIteratorOfListOfExtendedString.hxx>
3335 #include <TDataStd_ListIteratorOfListOfByte.hxx>
3336 #include <TDF_ListIteratorOfLabelList.hxx>
3337 // Attributes
3338 #include <TDataStd_Tick.hxx>
3339 #include <TDataStd_IntegerList.hxx>
3340 #include <TDataStd_RealList.hxx>
3341 #include <TDataStd_ExtStringList.hxx>
3342 #include <TDataStd_BooleanList.hxx>
3343 #include <TDataStd_ReferenceList.hxx>
3344 #include <TDataStd_BooleanArray.hxx>
3345 #include <TDataStd_ReferenceArray.hxx>
3346 #include <TDataStd_ByteArray.hxx>
3347 #include <TDataStd_NamedData.hxx>
3348 #include <TDF_Reference.hxx>
3349 //
3350 Handle(AppStd_Application) app;
3351 int TestSetGet(const Handle(TDocStd_Document)& doc)
3352 {
3353   // TDataStd_Tick:
3354   // Set
3355   TDataStd_Tick::Set(doc->Main());
3356   // Get
3357   Handle(TDataStd_Tick) tick;
3358   if (!doc->Main().FindAttribute(TDataStd_Tick::GetID(), tick))
3359     return 1;
3360   // Forget
3361   doc->Main().ForgetAttribute(TDataStd_Tick::GetID());
3362   if (doc->Main().IsAttribute(TDataStd_Tick::GetID()))
3363     return 2;
3364   doc->Main().ResumeAttribute(tick);
3365   if (!doc->Main().IsAttribute(TDataStd_Tick::GetID()))
3366     return 3;
3367   // Forget
3368   doc->Main().ForgetAttribute(TDataStd_Tick::GetID());
3369   if (doc->Main().IsAttribute(TDataStd_Tick::GetID()))
3370     return 2;
3371
3372   // TDataStd_IntegerList:
3373   // Set
3374   Handle(TDataStd_IntegerList) setintlist = TDataStd_IntegerList::Set(doc->Main());
3375   setintlist->Append(2);
3376   setintlist->Prepend(1);
3377   setintlist->InsertAfter(3, 2);
3378   setintlist->InsertBefore(0, 1);
3379   setintlist->Append(200);
3380   setintlist->Remove(0);
3381   setintlist->Remove(200);
3382   // Get
3383   Handle(TDataStd_IntegerList) getintlist;
3384   if (!doc->Main().FindAttribute(TDataStd_IntegerList::GetID(), getintlist))
3385     return 1;
3386   if (getintlist->First() != 1)
3387     return 2;
3388   if (getintlist->Last() != 3)
3389     return 3;
3390   const TColStd_ListOfInteger& intlist = getintlist->List();
3391   TColStd_ListIteratorOfListOfInteger itr_intlist(intlist);
3392   for (; itr_intlist.More(); itr_intlist.Next())
3393   {
3394     if (itr_intlist.Value() != 1 &&
3395         itr_intlist.Value() != 2 &&
3396         itr_intlist.Value() != 3)
3397     {
3398       return 4;
3399     }
3400   }
3401   getintlist->Clear();
3402
3403   // TDataStd_RealList:
3404   // Set
3405   Handle(TDataStd_RealList) setdbllist = TDataStd_RealList::Set(doc->Main());
3406   setdbllist->Append(2.5);
3407   setdbllist->Prepend(1.5);
3408   setdbllist->InsertAfter(3.5, 2.5);
3409   setdbllist->InsertBefore(0.5, 1.5);
3410   setdbllist->Append(200.5);
3411   setdbllist->Remove(0.5);
3412   setdbllist->Remove(200.5);
3413   // Get
3414   Handle(TDataStd_RealList) getdbllist;
3415   if (!doc->Main().FindAttribute(TDataStd_RealList::GetID(), getdbllist))
3416     return 1;
3417   if (getdbllist->First() != 1.5)
3418     return 2;
3419   if (getdbllist->Last() != 3.5)
3420     return 3;
3421   const TColStd_ListOfReal& dbllist = getdbllist->List();
3422   TColStd_ListIteratorOfListOfReal itr_dbllist(dbllist);
3423   for (; itr_dbllist.More(); itr_dbllist.Next())
3424   {
3425     if (itr_dbllist.Value() != 1.5 &&
3426         itr_dbllist.Value() != 2.5 &&
3427         itr_dbllist.Value() != 3.5)
3428     {
3429       return 4;
3430     }
3431   }
3432   getdbllist->Clear();
3433
3434   // TDataStd_ExtStringList:
3435   // Set
3436   Handle(TDataStd_ExtStringList) setstrlist = TDataStd_ExtStringList::Set(doc->Main());
3437   setstrlist->Append("Hello");
3438   setstrlist->Prepend("Guten Tag");
3439   setstrlist->InsertAfter("Bonjour", "Guten Tag");
3440   setstrlist->InsertBefore("Bonsoir", "Hello");
3441   setstrlist->Append("Good bye");
3442   setstrlist->Remove("Bonsoir");
3443   setstrlist->Remove("Good bye");
3444   // Get
3445   Handle(TDataStd_ExtStringList) getstrlist;
3446   if (!doc->Main().FindAttribute(TDataStd_ExtStringList::GetID(), getstrlist))
3447     return 1;
3448   if (getstrlist->First() != "Guten Tag")
3449     return 2;
3450   if (getstrlist->Last() != "Hello")
3451     return 3;
3452   const TDataStd_ListOfExtendedString& strlist = getstrlist->List();
3453   TDataStd_ListIteratorOfListOfExtendedString itr_strlist(strlist);
3454   for (; itr_strlist.More(); itr_strlist.Next())
3455   {
3456     if (itr_strlist.Value() != "Guten Tag" &&
3457         itr_strlist.Value() != "Bonjour" &&
3458         itr_strlist.Value() != "Hello")
3459     {
3460       return 4;
3461     }
3462   }
3463   getstrlist->Clear();
3464
3465   // TDataStd_BooleanList:
3466   // Set
3467   Handle(TDataStd_BooleanList) setboollist = TDataStd_BooleanList::Set(doc->Main());
3468   setboollist->Append(Standard_True);
3469   setboollist->Prepend(Standard_False);
3470   // Get
3471   Handle(TDataStd_BooleanList) getboollist;
3472   if (!doc->Main().FindAttribute(TDataStd_BooleanList::GetID(), getboollist))
3473     return 1;
3474   if (getboollist->First() != Standard_False)
3475     return 2;
3476   if (getboollist->Last() != Standard_True)
3477     return 3;
3478   const TDataStd_ListOfByte& boollist = getboollist->List();
3479   TDataStd_ListIteratorOfListOfByte itr_boollist(boollist);
3480   for (; itr_boollist.More(); itr_boollist.Next())
3481   {
3482     if (itr_boollist.Value() != Standard_True &&
3483         itr_boollist.Value() != Standard_False)
3484     {
3485       return 4;
3486     }
3487   }
3488   getboollist->Clear();
3489
3490   // TDataStd_ReferenceList:
3491   TDF_Label L1 = doc->Main().FindChild(100);
3492   TDF_Label L2 = doc->Main().FindChild(101);
3493   TDF_Label L3 = doc->Main().FindChild(102);
3494   TDF_Label L4 = doc->Main().FindChild(103);
3495   TDF_Label L5 = doc->Main().FindChild(104);
3496   // Set
3497   Handle(TDataStd_ReferenceList) setreflist = TDataStd_ReferenceList::Set(doc->Main());
3498   setreflist->Append(L1);
3499   setreflist->Prepend(L2);
3500   setreflist->InsertAfter(L3, L2);
3501   setreflist->InsertBefore(L4, L1);
3502   setreflist->Append(L5);
3503   setreflist->Remove(L4);
3504   setreflist->Remove(L5);
3505   // Get
3506   Handle(TDataStd_ReferenceList) getreflist;
3507   if (!doc->Main().FindAttribute(TDataStd_ReferenceList::GetID(), getreflist))
3508     return 1;
3509   if (getreflist->First() != L2)
3510     return 2;
3511   if (getreflist->Last() != L1)
3512     return 3;
3513   const TDF_LabelList& reflist = getreflist->List();
3514   TDF_ListIteratorOfLabelList itr_reflist(reflist);
3515   for (; itr_reflist.More(); itr_reflist.Next())
3516   {
3517     if (itr_reflist.Value() != L1 &&
3518         itr_reflist.Value() != L2 &&
3519         itr_reflist.Value() != L3)
3520     {
3521       return 4;
3522     }
3523   }
3524   getreflist->Clear();
3525
3526   // TDataStd_BooleanArray:
3527   // Set
3528   Handle(TDataStd_BooleanArray) setboolarr = TDataStd_BooleanArray::Set(doc->Main(), 12, 16);
3529   setboolarr->SetValue(12, Standard_True);
3530   setboolarr->SetValue(13, Standard_False);
3531   setboolarr->SetValue(14, Standard_False);
3532   setboolarr->SetValue(15, Standard_False);
3533   setboolarr->SetValue(16, Standard_True);
3534   setboolarr->SetValue(14, Standard_True);
3535   // Get
3536   Handle(TDataStd_BooleanArray) getboolarr;
3537   if (!doc->Main().FindAttribute(TDataStd_BooleanArray::GetID(), getboolarr))
3538     return 1;
3539   if (getboolarr->Value(12) != Standard_True)
3540     return 2;
3541   if (getboolarr->Value(13) != Standard_False)
3542     return 2;
3543   if (getboolarr->Value(14) != Standard_True)
3544     return 2;
3545   if (getboolarr->Value(15) != Standard_False)
3546     return 2;
3547   if (getboolarr->Value(16) != Standard_True)
3548     return 2;
3549
3550   // TDataStd_ReferenceArray:
3551   // Set
3552   Handle(TDataStd_ReferenceArray) setrefarr = TDataStd_ReferenceArray::Set(doc->Main(), 0, 4);
3553   setrefarr->SetValue(0, L1);
3554   setrefarr->SetValue(1, L2);
3555   setrefarr->SetValue(2, L3);
3556   setrefarr->SetValue(3, L4);
3557   setrefarr->SetValue(4, L5);
3558   // Get
3559   Handle(TDataStd_ReferenceArray) getrefarr;
3560   if (!doc->Main().FindAttribute(TDataStd_ReferenceArray::GetID(), getrefarr))
3561     return 1;
3562   if (getrefarr->Value(0) != L1)
3563     return 2;
3564   if (getrefarr->Value(1) != L2)
3565     return 2;
3566   if (getrefarr->Value(2) != L3)
3567     return 2;
3568   if (getrefarr->Value(3) != L4)
3569     return 2;
3570   if (getrefarr->Value(4) != L5)
3571     return 2;
3572
3573   // TDataStd_ByteArray:
3574   // Set
3575   Handle(TDataStd_ByteArray) setbytearr = TDataStd_ByteArray::Set(doc->Main(), 12, 16);
3576   setbytearr->SetValue(12, 0);
3577   setbytearr->SetValue(13, 1);
3578   setbytearr->SetValue(14, 2);
3579   setbytearr->SetValue(15, 3);
3580   setbytearr->SetValue(16, 255);
3581   // Get
3582   Handle(TDataStd_ByteArray) getbytearr;
3583   if (!doc->Main().FindAttribute(TDataStd_ByteArray::GetID(), getbytearr))
3584     return 1;
3585   if (getbytearr->Value(12) != 0)
3586     return 2;
3587   if (getbytearr->Value(13) != 1)
3588     return 2;
3589   if (getbytearr->Value(14) != 2)
3590     return 2;
3591   if (getbytearr->Value(15) != 3)
3592     return 2;
3593   if (getbytearr->Value(16) != 255)
3594     return 2;
3595   
3596   // TDataStd_NamedData:
3597   // Set:
3598   Handle(TDataStd_NamedData) setnd = TDataStd_NamedData::Set(doc->Main());
3599   setnd->SetInteger("Integer1", 1);
3600   setnd->SetInteger("Integer2", 2);
3601   setnd->SetInteger("Integer3", 8);
3602   setnd->SetInteger("Integer3", 3);
3603   // Get:
3604   Handle(TDataStd_NamedData) getnd;
3605   if (!doc->Main().FindAttribute(TDataStd_NamedData::GetID(), getnd))
3606     return 1;
3607   if (!getnd->HasIntegers())
3608     return 2;
3609   if (!getnd->HasInteger("Integer1"))
3610     return 3;
3611   if (getnd->GetInteger("Integer2") != 2)
3612     return 4;
3613   if (getnd->GetInteger("Integer3") != 3)
3614     return 4;
3615
3616   return 0;
3617 }
3618
3619 int TestUndoRedo(const Handle(TDocStd_Document)& doc)
3620 {
3621   // TDataStd_Tick:
3622   doc->OpenCommand();
3623   Handle(TDataStd_Tick) tick = TDataStd_Tick::Set(doc->Main());
3624   doc->CommitCommand();
3625   if (!doc->Main().IsAttribute(TDataStd_Tick::GetID()))
3626     return 1;
3627   doc->Undo();
3628   if (doc->Main().IsAttribute(TDataStd_Tick::GetID()))
3629     return 2;
3630   doc->Redo();
3631   if (!doc->Main().IsAttribute(TDataStd_Tick::GetID()))
3632     return 3;
3633
3634   // TDataStd_IntegerList:
3635   doc->OpenCommand();
3636   Handle(TDataStd_IntegerList) intlist = TDataStd_IntegerList::Set(doc->Main());
3637   intlist->Append(2);
3638   intlist->Prepend(1);
3639   intlist->InsertBefore(0, 1);
3640   intlist->InsertAfter(3, 2);
3641   doc->CommitCommand();
3642   if (!doc->Main().IsAttribute(TDataStd_IntegerList::GetID()))
3643     return 1;
3644   doc->Undo();
3645   if (!intlist->IsEmpty())
3646     return 2;
3647   doc->Redo();
3648   if (!intlist->Extent())
3649     return 3;
3650   if (intlist->First() != 0)
3651     return 4;
3652   if (intlist->Last() != 3)
3653     return 5;
3654   intlist->Clear();
3655
3656   // TDataStd_RealList:
3657   doc->OpenCommand();
3658   Handle(TDataStd_RealList) dbllist = TDataStd_RealList::Set(doc->Main());
3659   dbllist->Append(2.5);
3660   dbllist->Prepend(1.5);
3661   dbllist->InsertBefore(0.5, 1.5);
3662   dbllist->InsertAfter(3.5, 2.5);
3663   doc->CommitCommand();
3664   if (!doc->Main().IsAttribute(TDataStd_RealList::GetID()))
3665     return 1;
3666   doc->Undo();
3667   if (!dbllist->IsEmpty())
3668     return 2;
3669   doc->Redo();
3670   if (!dbllist->Extent())
3671     return 3;
3672   if (dbllist->First() != 0.5)
3673     return 4;
3674   if (dbllist->Last() != 3.5)
3675     return 5;
3676   dbllist->Clear();
3677
3678   // TDataStd_ExtStringList:
3679   doc->OpenCommand();
3680   Handle(TDataStd_ExtStringList) strlist = TDataStd_ExtStringList::Set(doc->Main());
3681   strlist->Append("Hello");
3682   strlist->Prepend("Guten Tag");
3683   strlist->InsertAfter("Bonjour", "Guten Tag");
3684   strlist->InsertBefore("Bonsoir", "Hello");
3685   doc->CommitCommand();
3686   if (!doc->Main().IsAttribute(TDataStd_ExtStringList::GetID()))
3687     return 1;
3688   doc->Undo();
3689   if (!strlist->IsEmpty())
3690     return 2;
3691   doc->Redo();
3692   if (!strlist->Extent())
3693     return 3;
3694   if (strlist->First() != "Guten Tag")
3695     return 4;
3696   if (strlist->Last() != "Hello")
3697     return 5;
3698   strlist->Clear();
3699
3700   // TDataStd_BooleanList:
3701   doc->OpenCommand();
3702   Handle(TDataStd_BooleanList) boollist = TDataStd_BooleanList::Set(doc->Main());
3703   boollist->Append(Standard_True);
3704   boollist->Prepend(Standard_False);
3705   doc->CommitCommand();
3706   if (!doc->Main().IsAttribute(TDataStd_BooleanList::GetID()))
3707     return 1;
3708   doc->Undo();
3709   if (!boollist->IsEmpty())
3710     return 2;
3711   doc->Redo();
3712   if (!boollist->Extent())
3713     return 3;
3714   if (boollist->First() != Standard_False)
3715     return 4;
3716   if (boollist->Last() != Standard_True)
3717     return 5;
3718   boollist->Clear();
3719
3720   // TDataStd_ReferenceList:
3721   TDF_Label L1 = doc->Main().FindChild(100);
3722   TDF_Label L2 = doc->Main().FindChild(101);
3723   TDF_Label L3 = doc->Main().FindChild(102);
3724   TDF_Label L4 = doc->Main().FindChild(103);
3725   doc->OpenCommand();
3726   Handle(TDataStd_ReferenceList) reflist = TDataStd_ReferenceList::Set(doc->Main());
3727   reflist->Append(L1);
3728   reflist->Prepend(L2);
3729   reflist->InsertBefore(L3, L1);
3730   reflist->InsertAfter(L4, L2);
3731   doc->CommitCommand();
3732   if (!doc->Main().IsAttribute(TDataStd_ReferenceList::GetID()))
3733     return 1;
3734   doc->Undo();
3735   if (!reflist->IsEmpty())
3736     return 2;
3737   doc->Redo();
3738   if (!reflist->Extent())
3739     return 3;
3740   if (reflist->First() != L2)
3741     return 4;
3742   if (reflist->Last() != L1)
3743     return 5;
3744   reflist->Clear();
3745
3746   // TDataStd_BooleanArray:
3747   doc->OpenCommand();
3748   Handle(TDataStd_BooleanArray) boolarr = TDataStd_BooleanArray::Set(doc->Main(), 23, 25);
3749   boolarr->SetValue(23, Standard_True);
3750   boolarr->SetValue(25, Standard_True);
3751   doc->CommitCommand();
3752   doc->OpenCommand();
3753   boolarr = TDataStd_BooleanArray::Set(doc->Main(), 230, 250);
3754   boolarr->SetValue(230, Standard_True);
3755   boolarr->SetValue(250, Standard_True);  
3756   doc->CommitCommand();
3757   doc->Undo();
3758   if (boolarr->Value(23) != Standard_True)
3759     return 2;
3760   if (boolarr->Value(24) != Standard_False)
3761     return 2;
3762   if (boolarr->Value(25) != Standard_True)
3763     return 2;
3764   doc->Redo();
3765   if (boolarr->Value(230) != Standard_True)
3766     return 3;
3767   if (boolarr->Value(240) != Standard_False)
3768     return 3;
3769   if (boolarr->Value(250) != Standard_True)
3770     return 3;
3771
3772   // TDataStd_ReferenceArray:
3773   doc->OpenCommand();
3774   Handle(TDataStd_ReferenceArray) refarr = TDataStd_ReferenceArray::Set(doc->Main(), 5, 8);
3775   refarr->SetValue(5, L1);
3776   refarr->SetValue(6, L2);
3777   refarr->SetValue(7, L3);
3778   refarr->SetValue(8, L4);
3779   doc->CommitCommand();
3780   if (!doc->Main().IsAttribute(TDataStd_ReferenceArray::GetID()))
3781     return 1;
3782   doc->Undo();
3783   doc->Redo();
3784   if (refarr->Value(5) != L1)
3785     return 4;
3786   if (refarr->Value(6) != L2)
3787     return 4;
3788   if (refarr->Value(7) != L3)
3789     return 4;
3790   if (refarr->Value(8) != L4)
3791     return 4;
3792
3793   // TDataStd_ByteArray:
3794   doc->OpenCommand();
3795   Handle(TDataStd_ByteArray) bytearr = TDataStd_ByteArray::Set(doc->Main(), 23, 25);
3796   bytearr->SetValue(23, 23);
3797   bytearr->SetValue(25, 25);
3798   doc->CommitCommand();
3799   doc->OpenCommand();
3800   bytearr = TDataStd_ByteArray::Set(doc->Main(), 230, 250);
3801   bytearr->SetValue(230, 230);
3802   bytearr->SetValue(250, 250);  
3803   doc->CommitCommand();
3804   doc->Undo();
3805   if (bytearr->Value(23) != 23)
3806     return 2;
3807   if (bytearr->Value(25) != 25)
3808     return 2;
3809   doc->Redo();
3810   if (bytearr->Value(230) != 230)
3811     return 3;
3812   if (bytearr->Value(250) != 250)
3813     return 3;
3814
3815   // TDataStd_NamedData:
3816   doc->OpenCommand();
3817   Hand