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