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