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