0023030: Collect all DRAW commands created for testing bugs into one package
[occt.git] / src / QABugs / QABugs_17.cxx
1 // Created on: 2002-03-18
2 // Created by: QA Admin
3 // Copyright (c) 2002-2012 OPEN CASCADE SAS
4 //
5 // The content of this file is subject to the Open CASCADE Technology Public
6 // License Version 6.5 (the "License"). You may not use the content of this file
7 // except in compliance with the License. Please obtain a copy of the License
8 // at http://www.opencascade.org and read it completely before using this file.
9 //
10 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
11 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
12 //
13 // The Original Code and all software distributed under the License is
14 // distributed on an "AS IS" basis, without warranty of any kind, and the
15 // Initial Developer hereby disclaims all such warranties, including without
16 // limitation, any warranties of merchantability, fitness for a particular
17 // purpose or non-infringement. Please see the License for the specific terms
18 // and conditions governing the rights and limitations under the License.
19
20
21 #if defined(WOKC40)
22 #define V2D
23 #define Viewer2dTest QAViewer2dTest
24 #endif
25
26 #include <stdio.h>
27
28 #include <QABugs.hxx>
29
30 #include <Draw_Interpretor.hxx>
31 #include <DBRep.hxx>
32 #include <DrawTrSurf.hxx>
33 #include <AIS_InteractiveContext.hxx>
34 #include <ViewerTest.hxx>
35 #include <AIS_Shape.hxx>
36 #include <TopoDS_Shape.hxx>
37
38 #include <Geom_Circle.hxx>
39 #include <Geom_Ellipse.hxx>
40 #include <Geom_Plane.hxx>
41 #include <gp_Pln.hxx>
42 #include <Geom2d_Curve.hxx>
43 #include <GeomAPI.hxx>
44 #include <Geom2dAdaptor_Curve.hxx>
45 #include <Geom2dGcc_QualifiedCurve.hxx>
46 #include <Geom2dGcc_Lin2d2Tan.hxx>
47 #include <Geom2d_Line.hxx>
48 #include <BRepBuilderAPI_MakeEdge.hxx>
49 #include <TopoDS_Edge.hxx>
50 #include <Precision.hxx>
51 #include <Geom2d_Circle.hxx>
52 #include <Geom2dGcc_MyQCurve.hxx>
53 #include <GccEnt_QualifiedCirc.hxx>
54 #include <Geom2dGcc_MyL2d2Tan.hxx>
55 #include <BRepBuilderAPI_MakeWire.hxx>
56 #include <TopExp_Explorer.hxx>
57 #include <TopoDS.hxx>
58 #include <TopoDS_Wire.hxx>
59 #include <BRep_Tool.hxx>
60 #include <gp_Circ.hxx>
61 #include <BRepOffsetAPI_MakePipeShell.hxx>
62 #include <AIS_Trihedron.hxx>
63 #include <Geom_Axis2Placement.hxx>
64 #include <V3d_View.hxx>
65 #include <ViewerTest_Tool.hxx>
66 #include <BRepBuilderAPI_MakeVertex.hxx>
67 #include <AIS_LengthDimension.hxx>
68 #include <AIS_RadiusDimension.hxx>
69 #include <BRepPrimAPI_MakeBox.hxx>
70 #include <GeomAdaptor_Surface.hxx>
71 #include <Extrema_ExtPS.hxx>
72 #include <Geom_CartesianPoint.hxx>
73 #include <Geom2d_CartesianPoint.hxx>
74 #include <Geom2dGcc_Circ2d3Tan.hxx>
75 #include <Aspect_Window.hxx>
76
77 #include <TopoDS_Compound.hxx>
78 #include <TopoDS_CompSolid.hxx>
79 #include <StdSelect_ShapeTypeFilter.hxx>
80
81 #include <QABugs_MyText.hxx>
82
83 #if defined(WOKC40)
84 #include <QAViewer2dTest.hxx>
85 #include <QAViewer2dTest_DoubleMapOfInteractiveAndName.hxx>
86 #else
87 #include <Viewer2dTest.hxx>
88 #include <Viewer2dTest_DoubleMapOfInteractiveAndName.hxx>
89 #endif
90
91 #include <Prs3d_Projector.hxx>
92 #include <HLRAlgo_Projector.hxx>
93 #include <AIS2D_ProjShape.hxx>
94 #include <AIS2D_InteractiveContext.hxx>
95 #include <Standard_ErrorHandler.hxx>
96
97 #if ! defined(WNT)
98 #if defined(WOKC40)
99 extern QAViewer2dTest_DoubleMapOfInteractiveAndName& GetMapOfAIS2D();
100 #else
101 extern Viewer2dTest_DoubleMapOfInteractiveAndName& GetMapOfAIS2D();
102 #endif
103 #else
104 Standard_EXPORT Viewer2dTest_DoubleMapOfInteractiveAndName& GetMapOfAIS2D();
105 #endif
106
107 #include <GGraphic2d_SetOfCurves.hxx>
108 #include <Graphic2d_SetOfSegments.hxx>
109
110 #include <Graphic3d_NameOfFont.hxx>
111
112 static Standard_Integer BUC60842 (Draw_Interpretor& di, Standard_Integer /*argc*/,const char ** /*argv*/)
113 {
114   Standard_Character abuf[16];
115   //Standard_Character * abuf = new Standard_Character[16];
116
117   Standard_CString st = abuf;
118
119   Handle(Geom_Circle) cir = new Geom_Circle (gp_Ax2 (gp_Pnt (823.687192, 502.366825, 478.960440), gp_Dir (0.173648, 0.984808, 0.000000), gp_Dir (-0.932169, 0.164367, -0.322560)), 50);
120   Handle(Geom_Ellipse) ell = new Geom_Ellipse (gp_Ax2 (gp_Pnt (1262.224429, 425.040878, 363.609716), gp_Dir (0.173648, 0.984808, 0.000000), gp_Dir (-0.932169, 0.164367, -0.322560)), 150, 100);
121   Handle(Geom_Plane) plne = new Geom_Plane (gp_Ax3 (gp_Ax2 (gp_Pnt (1262.224429, 425.040878, 363.609716), gp_Dir (0.173648, 0.984808, 0.000000), gp_Dir (-0.932169, 0.164367, -0.322560))));
122
123   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
124
125   gp_Pln pln = plne->Pln();
126   Handle(Geom2d_Curve) curve2d = GeomAPI::To2d (ell, pln);
127   sprintf(abuf,"ell");
128   DrawTrSurf::Set(st,curve2d);
129   if(!aContext.IsNull()) {
130     Handle(AIS_Shape) aisp = 
131       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(curve2d, pln)).Edge());
132     aContext->Display (aisp);
133   }
134
135   Handle(Geom2d_Curve) fromcurve2d = GeomAPI::To2d (cir, pln);
136   sprintf(abuf,"cil");
137   DrawTrSurf::Set(st,fromcurve2d);
138   if(!aContext.IsNull()) {
139     Handle(AIS_Shape) aisp = 
140       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(fromcurve2d, pln)).Edge());
141     aContext->Display (aisp);
142   }
143
144   Geom2dAdaptor_Curve acur (curve2d), afromcur (fromcurve2d);
145   Geom2dGcc_QualifiedCurve qcur (acur, GccEnt_outside);
146   Geom2dGcc_QualifiedCurve qfromcur (afromcur, GccEnt_outside);
147   Geom2dGcc_Lin2d2Tan lintan (qcur, qfromcur, 0.1, 0.0, 0.0);
148   Standard_Integer i=0;
149   for(i=0;i<lintan.NbSolutions();i++) {
150     sprintf(abuf,"lintan_%d",i);
151     Handle(Geom2d_Line) glin = new Geom2d_Line(lintan.ThisSolution(i));
152     DrawTrSurf::Set(st,glin);
153     if(!aContext.IsNull()) {
154       Handle(AIS_Shape) aisp = 
155         new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(glin, pln)).Edge());      aContext->Display (aisp);
156     }
157   }
158   di << " Is Done = \n" << (Standard_Integer) lintan.IsDone();
159   return 0;
160 }
161
162 static Standard_Integer BUC60843 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
163 {
164   if (argc != 4 && argc != 5 && argc != 6 && argc != 7) {
165     di << "Usage : " << argv[0] << "result_shape name_of_circle name_of_curve [ par1 [ tol ] ]" << "\n";
166     return 1;
167   }
168
169   Standard_Real par1 = 0.0, par2 = 0.0;
170   Standard_Real tol  = Precision::Angular();
171   if (argc >= 5)
172     par1 = atof(argv[4]);
173   if (argc == 6)
174     par2 = atof(argv[5]);
175   if (argc == 7)
176     tol = atof(argv[6]);
177   Handle(Geom2d_Curve)  aCur2d1 = DrawTrSurf::GetCurve2d(argv[2]);
178   Handle(Geom2d_Curve)  aCur2d2 = DrawTrSurf::GetCurve2d(argv[3]);
179   if (aCur2d1.IsNull() || aCur2d2.IsNull()) {
180     di << "Faulty : not curves" << "\n";
181     return 1;
182   }
183   Standard_Boolean c1IsCircle = Standard_False;
184   Handle(Geom2d_Circle) aCir2d = Handle(Geom2d_Circle)::DownCast(aCur2d1);
185   if (!aCir2d.IsNull()) {
186     c1IsCircle = Standard_True;
187     if (argc == 6) tol = atof(argv[5]);
188   }
189   if(c1IsCircle) {
190     Geom2dAdaptor_Curve acur(aCur2d2);
191     Geom2dGcc_MyQCurve qcur(acur, GccEnt_unqualified);
192     GccEnt_QualifiedCirc qfromcur(aCir2d->Circ2d(), GccEnt_unqualified);
193     Geom2dGcc_MyL2d2Tan lintan(qfromcur, qcur , par1, tol); 
194     if (lintan.IsDone()) {
195       gp_Lin2d lin = lintan.ThisSolution(); 
196       Handle(Geom2d_Line) glin = new Geom2d_Line(lin); 
197       DrawTrSurf::Set(argv[1], glin);
198     } else {
199       di << "Faulty : tangent line is not done." << "\n";
200     }
201   }
202   else {
203     Geom2dAdaptor_Curve acur1(aCur2d1);
204     Geom2dAdaptor_Curve acur2(aCur2d2);
205     Geom2dGcc_MyQCurve qcur1(acur1, GccEnt_unqualified);
206     Geom2dGcc_MyQCurve qcur2(acur2, GccEnt_unqualified);
207     Geom2dGcc_MyL2d2Tan lintan(qcur1, qcur2 , par1, par2, tol); 
208     if (lintan.IsDone()) {
209       gp_Lin2d lin = lintan.ThisSolution(); 
210       Handle(Geom2d_Line) glin = new Geom2d_Line(lin); 
211       DrawTrSurf::Set(argv[1], glin);
212     } else {
213       di << "Faulty : tangent line is not done." << "\n";
214     }
215   }   
216   return 0;
217 }
218
219 static Standard_Integer BUC60970 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
220 {
221   if(argc!=3)
222   {
223     di << "Usage : " << argv[0] << " shape result" <<"\n";
224     return 1;
225   }
226
227   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
228   if(aContext.IsNull()) {
229     di << "use 'vinit' command before " << argv[0] << "\n";
230     return 1;
231   }
232
233   TopoDS_Shape aResult = DBRep::Get(argv[1]);
234
235   BRepBuilderAPI_MakeWire bRepSpineWire; 
236   TopExp_Explorer exp_(aResult, TopAbs_WIRE); 
237
238   Standard_Integer i = 0; 
239   for(;exp_.More();exp_.Next()){     
240     i++; 
241     bRepSpineWire.Add(TopoDS::Wire(exp_.Current())); 
242   } 
243
244   //printf("total no of wires are ............. %d\n", i); 
245   di << "total no of wires are ............. " << i << "\n"; 
246
247   TopoDS_Wire spineWire = bRepSpineWire.Wire(); 
248   aContext->Display(new AIS_Shape(spineWire)); 
249
250   DBRep::Set("slineW",spineWire);
251
252   TopExp_Explorer spineWireExp(spineWire, TopAbs_EDGE); 
253   Standard_Real first, last; 
254   Handle_Geom_Curve curl_ = BRep_Tool::Curve(TopoDS::Edge(spineWireExp.Current()), first, last); 
255   gp_Pnt firstPnt; 
256   gp_Vec tanVec; 
257   curl_->D1(first, firstPnt, tanVec); 
258   tanVec.Normalize(); 
259   gp_Dir tanDir(tanVec.X(), tanVec.Y(), tanVec.Z()); 
260   gp_Ax2 gpAx2(firstPnt, tanDir); 
261   gp_Circ gpCirc(gpAx2, 2.5); 
262   BRepBuilderAPI_MakeWire aMWire(BRepBuilderAPI_MakeEdge(new Geom_Circle(gpCirc)).Edge());
263   TopoDS_Wire topoWire(aMWire); 
264   aContext->Display(new AIS_Shape(topoWire));
265
266   DBRep::Set("topoW",topoWire);
267
268   BRepOffsetAPI_MakePipeShell bRepPipe(spineWire);
269   bRepPipe.Add(topoWire); 
270   bRepPipe.Build(); 
271
272   aContext->Display(new AIS_Shape(bRepPipe.Shape())); 
273
274   DBRep::Set(argv[2],bRepPipe.Shape());
275
276   return 0;
277 }
278
279 static Standard_Integer  BUC60818(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
280 {
281   if(argc!=1)
282   {
283     di << "Usage : " << argv[0] << "\n";
284     return -1;
285   }
286
287   Handle(AIS_InteractiveContext) myAISContext = ViewerTest::GetAISContext();
288   if(myAISContext.IsNull()) 
289   {
290     di << "use 'vinit' command before " << argv[0] << "\n";
291     return -1;
292   }
293
294   // TRIHEDRON
295   Handle(AIS_Trihedron) aTrihedron;
296   Handle(Geom_Axis2Placement) aTrihedronAxis=new Geom_Axis2Placement(gp::XOY());
297   aTrihedron=new AIS_Trihedron(aTrihedronAxis);
298   myAISContext->Display(aTrihedron);
299
300   myAISContext->OpenLocalContext(); 
301   myAISContext->Load(aTrihedron,0); 
302
303   myAISContext->SetAutomaticHilight(  Standard_False );
304
305   Handle(V3d_View) myV3dView = ViewerTest::CurrentView();
306   double Xv,Yv;
307   myV3dView->Project(0,0,0,Xv,Yv);
308
309   Standard_Integer Xp,Yp;
310   myV3dView->Convert(Xv,Yv,Xp,Yp);
311
312   myAISContext->MoveTo( Xp,Yp, myV3dView  );
313   myAISContext->MoveTo( Xp,Yp, myV3dView  );
314
315   if (myAISContext->HasDetected(  )) 
316     di << "has detected shape : OK"   << "\n";
317   else       di << "has detected shape : bugged - Faulty "   << "\n";
318
319   myAISContext->SetAutomaticHilight(  Standard_True );
320
321   return 0;
322 }
323
324 static Standard_Integer BUC60915_1(Draw_Interpretor& di, Standard_Integer argc, const char ** /*argv*/)
325 {
326   if (argc > 1) {
327     di<<"Function don't has parameters"<<"\n";
328     return 1;
329   }
330
331   Handle(AIS_InteractiveContext)   context= ViewerTest_Tool::MakeContext ("buc60915");
332   ViewerTest_Tool::InitViewerTest (context);
333
334   //The following dimesion code has problems regarding arrow_size. The desired effect is not produced.
335   /***************************************/
336   //First view
337   /***************************************/
338   gp_Pnt p1 = gp_Pnt(602.51,50.,0.);
339   gp_Pnt p2 = gp_Pnt(602.51,200.,0.);
340   gp_Pnt p3 = gp_Pnt(102.51,200.,0.);
341   gp_Pnt p4 = gp_Pnt(102.51,170.,0.);
342   gp_Pnt p5 = gp_Pnt(502.51,170.,0.);
343   gp_Pnt p6 = gp_Pnt(502.51,80.,0.);
344   gp_Pnt p7 = gp_Pnt(102.51,80.,0.);
345   gp_Pnt p8 = gp_Pnt(102.51,50.,0.);
346   TopoDS_Vertex V1 = BRepBuilderAPI_MakeVertex(p1);
347   TopoDS_Vertex V2 = BRepBuilderAPI_MakeVertex(p2);
348   TopoDS_Vertex V3 = BRepBuilderAPI_MakeVertex(p3);
349   TopoDS_Vertex V4 = BRepBuilderAPI_MakeVertex(p4);
350   TopoDS_Vertex V5 = BRepBuilderAPI_MakeVertex(p5);
351   TopoDS_Vertex V6 = BRepBuilderAPI_MakeVertex(p6);
352   TopoDS_Vertex V7 = BRepBuilderAPI_MakeVertex(p7);
353   TopoDS_Vertex V8 = BRepBuilderAPI_MakeVertex(p8);
354   gp_Pnt plnpt(0, 0, 0);
355   gp_Dir plndir(0, 0, 1);
356   Handle(Geom_Plane) pln = new Geom_Plane(plnpt,plndir);
357   /***************************************/
358   //dimension "L 502.51"
359   /***************************************/
360   Handle(AIS_LengthDimension) len = new AIS_LengthDimension(V2, V3, pln, 502.51, "502.51");
361   len->SetPosition(gp_Pnt(350, 250, 0));
362   len->SetTypeOfDist(AIS_TOD_Horizontal);
363   len->SetArrowSize(30.0);
364   context->Display(len);
365   /***************************************/
366   //dimension "L 90"
367   /***************************************/
368   Handle(AIS_LengthDimension) len1 = new AIS_LengthDimension(V4, V7, pln, 90, "90");
369   len1->SetPosition(gp_Pnt(70, 120, 0));
370   len1->SetTypeOfDist(AIS_TOD_Vertical);
371   len1->SetArrowSize(100.0);
372   context->Display(len1);
373   /***************************************/
374   //dimension "L 150"
375   /***************************************/
376   Handle(AIS_LengthDimension) len2 = new AIS_LengthDimension(V1, V2, pln, 150, "150",gp_Pnt(650, 120, 0),DsgPrs_AS_BOTHAR,
377     AIS_TOD_Vertical, 100.0 );
378   //len2->SetPosition(gp_Pnt(650, 120, 0));
379   //len2->SetTypeOfDist(AIS_TOD_Vertical);
380   //len1->SetArrowSize(100.0);
381   context->Display(len2);
382   /***************************************/
383   //dimension "R 88.58"
384   /***************************************/
385   gp_Circ cir = gp_Circ(gp_Ax2(gp_Pnt(191.09, -88.58, 0), gp_Dir(0, 0, 1)), 88.58);
386   TopoDS_Edge E1 = BRepBuilderAPI_MakeEdge(cir,gp_Pnt(191.09,0,0.),gp_Pnt(191.09,-177.16,0.) );
387   Handle(AIS_RadiusDimension) dim1 = new AIS_RadiusDimension(E1,88.58, "R 88.58",gp_Pnt(-30.0, -80.0, 0.0),DsgPrs_AS_BOTHAR,
388     100.0 );
389   context->Display(dim1);
390   /***************************************/
391   //dimension "R 43.80"
392   /***************************************/
393   gp_Circ cir1 = gp_Circ(gp_Ax2(gp_Pnt(191.09, -88.58, 0), gp_Dir(0, 0, 1)), 43.80);
394   TopoDS_Edge E_cir1 = BRepBuilderAPI_MakeEdge(cir1);
395   dim1 = new AIS_RadiusDimension(E_cir1,43.80, "R 43.80",gp_Pnt(0.0, -50.0, 0.0),DsgPrs_AS_LASTAR, 60.0 );
396   context->Display(dim1);
397   /***************************************/
398   //dimension "R 17.86"
399   /***************************************/
400   gp_Circ cir2 = gp_Circ(gp_Ax2(gp_Pnt(566.11, -88.58, 0), gp_Dir(0, 0, -1)), 17.86);
401   TopoDS_Edge E_cir2 = BRepBuilderAPI_MakeEdge(cir2);
402   dim1 = new AIS_RadiusDimension(E_cir2,17.86, "R 17.86",gp_Pnt(600.0, -50.0, 0.0),DsgPrs_AS_LASTAR, 40.0 );
403   context->Display(dim1);
404
405   return 0;
406 }
407
408 static Standard_Integer  OCC138 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
409 {
410   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
411   if(aContext.IsNull()) { 
412     di << "use 'vinit' command before " << argv[0] << "\n";
413     return 1;
414   }
415
416   BRepPrimAPI_MakeBox box1(gp_Pnt(0, 0, 0),  gp_Pnt(100, 100, 100));
417   BRepPrimAPI_MakeBox box2(gp_Pnt(120, 120, 120),  gp_Pnt(300, 300,300));
418   BRepPrimAPI_MakeBox box3(gp_Pnt(320, 320, 320),  gp_Pnt(500, 500,500));
419
420   Handle_AIS_Shape ais1 = new AIS_Shape(box1.Shape());
421   Handle_AIS_Shape ais2 = new AIS_Shape(box2.Shape());
422   Handle_AIS_Shape ais3 = new AIS_Shape(box3.Shape());
423
424   aContext->Display(ais1);
425   aContext->Display(ais2);
426   aContext->Display(ais3);
427
428   aContext->AddOrRemoveCurrentObject(ais1);
429   aContext->AddOrRemoveCurrentObject(ais2);
430   aContext->AddOrRemoveCurrentObject(ais3);
431
432   di << "\n No of currents = " << aContext->NbCurrents();
433
434   aContext->InitCurrent();
435
436   int count = 1;
437   while(aContext->MoreCurrent())
438   {
439     di << "\n count is = " << count++;
440     Handle_AIS_InteractiveObject ais = aContext->Current();
441     aContext->AddOrRemoveCurrentObject(ais);
442     aContext->InitCurrent();
443   }
444
445   return 0; 
446 }
447
448 static Standard_Integer BUC60821(Draw_Interpretor& di, Standard_Integer argc,const char ** argv )   
449 {
450
451   if(argc < 3){
452     di << "Usage: " << argv[0] << " TextHight1 TextHight2 TextHight2";
453     return(-1);
454   }
455
456   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
457
458   if(aContext.IsNull()) 
459   { 
460     di << "use 'vinit' command before " << argv[0] << "\n";
461     return -1;
462   }
463
464   Handle(QABugs_MyText) txt1 = new QABugs_MyText("Gosha1",gp_Pnt(0,0,0),Graphic3d_NOF_ASCII_SIMPLEX,Quantity_NOC_RED,atoi(argv[1]));
465   aContext->Display(txt1);
466
467   Handle(QABugs_MyText) txt2 = new QABugs_MyText("Gosha2",gp_Pnt(0,0,100),Graphic3d_NOF_ASCII_SIMPLEX,Quantity_NOC_YELLOW,atoi(argv[2]));
468   aContext->Display(txt2);
469
470   Handle(QABugs_MyText) txt3 = new QABugs_MyText("Gosha3",gp_Pnt(0,100,100),Graphic3d_NOF_ASCII_SIMPLEX,Quantity_NOC_SKYBLUE,atoi(argv[3]));
471   aContext->Display(txt3);
472
473   return 0;
474 }
475
476 static int geom_get_2Dpt_from_3Dpt(const gp_Pnt& pnt3d, const gp_Pln& pln, gp_Pnt2d& pnt2d)
477
478   int ret = 0; 
479   Handle_Geom_Plane gpln = new Geom_Plane(pln); 
480   GeomAdaptor_Surface adsur(gpln); 
481   Extrema_ExtPS extps(pnt3d, adsur, 0.001, 0.001); 
482   if( extps.IsDone() ) 
483   { 
484     int index;
485     for(index=1; index<= extps.NbExt(); index++ ) 
486     { 
487       Extrema_POnSurf pons = extps.Point(index); 
488       Standard_Real U, V; 
489       pons.Parameter(U, V); 
490       pnt2d.SetCoord(U, V); 
491       ret = 1; 
492     } 
493   } 
494   return ret;
495
496
497 static Standard_Integer OCC353 (Draw_Interpretor& di, Standard_Integer , const char ** )
498 {
499   gp_Ax2 ax21( gp_Pnt(100,0,0), gp_Dir(0,0,1) );
500   Handle_Geom_Circle h_cir1 = new Geom_Circle( ax21, 25 );
501
502   gp_Ax2 ax22( gp_Pnt(-100,0,0), gp_Dir(0,0,1) );
503   Handle_Geom_Circle h_cir2 = new Geom_Circle( ax22, 25 );
504
505   gp_Pln refpln( gp_Pnt(0,0,0), gp_Dir(0,0,1) );
506   Handle_Geom2d_Curve cir2d1 = GeomAPI::To2d(h_cir1, refpln);
507   Handle_Geom2d_Curve cir2d2 = GeomAPI::To2d(h_cir2, refpln);
508
509   Geom2dAdaptor_Curve adop1(cir2d1);
510   Geom2dAdaptor_Curve adop2(cir2d2);
511
512   Geom2dGcc_QualifiedCurve qcur1(adop1, GccEnt_enclosing);
513   Geom2dGcc_QualifiedCurve qcur2(adop2, GccEnt_enclosing);
514
515   Handle_Geom_CartesianPoint h_carpt = new Geom_CartesianPoint(0,175,0);
516
517   gp_Pnt pt3d = h_carpt->Pnt();
518   gp_Pnt2d pt2d;
519   geom_get_2Dpt_from_3Dpt(pt3d, refpln, pt2d);
520
521   Standard_CString st;
522   st = "cir2d1";
523   DrawTrSurf::Set(st,cir2d1);
524   st = "cir2d2";
525   DrawTrSurf::Set(st,cir2d2);
526   st = "pt2d";
527   DrawTrSurf::Set(st,pt2d);
528
529   Handle_Geom2d_CartesianPoint pt = new Geom2d_CartesianPoint(pt2d);
530   Geom2dGcc_Circ2d3Tan sol( qcur1, qcur2, pt, 0.001, 0.0, 0.0);
531
532   int res = 0;
533   Standard_Character buf[10];
534   buf[0] = '\0';
535
536   if( sol.IsDone() ) {
537     res = sol.NbSolutions();
538     for(Standard_Integer i=1;i<=res;i++) {
539       Handle(Geom2d_Circle) aC = new Geom2d_Circle(sol.ThisSolution(i));
540       sprintf(buf,"Result_%d",i);
541       st = buf;
542       DrawTrSurf::Set(st,aC);
543     }
544   }
545   else
546     di << "\n Faulty: no solutions found ";
547   //printf("\n Faulty: no solutions found ");
548
549   //printf("\n Num of solutions are %d ", res );
550   di << "\n Num of solutions are " << res << " ";
551
552   return 0;
553 }
554
555 static Standard_Integer OCC280 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
556 {
557   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
558   if(aContext.IsNull()) 
559   { 
560     di << "use 'vinit' command before " << argv[0] << "\n";
561     return -1;
562   }
563   if ( argc != 3) {
564     di << "ERROR : Usage : " << argv[0] << " hlr=0/1 setsurfecedetail=0/1; set perspecrive view" << "\n";
565     return 1;
566   }
567
568   Standard_Integer HLR = atoi(argv[1]);
569   if (HLR != 0) {
570     HLR = 1;
571   }
572
573   Handle(V3d_View) aView = ViewerTest::CurrentView();
574
575   Handle(V3d_Viewer) aViewer = ViewerTest::GetViewerFromContext();
576   if(atoi(argv[2])) {
577     aViewer->SetDefaultSurfaceDetail(V3d_TEX_ALL);
578   }
579   aViewer->SetDefaultTypeOfView(V3d_PERSPECTIVE);
580
581   Handle(Aspect_Window) asp = aView->Window();
582   aViewer->SetViewOff(aView);
583   aView->Remove();
584
585   Handle(V3d_View) aNewView = aViewer->CreateView();
586   ViewerTest::CurrentView(aNewView);
587
588   aNewView->SetWindow(asp);
589   if (!asp->IsMapped()) asp->Map();
590
591   aNewView->Update();
592
593   if (HLR == 1) {
594     di << "HLR" << "\n";
595     aNewView->SetDegenerateModeOff();
596   }
597
598   return 0;
599 }
600
601 static Standard_Integer  OCC232 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
602 {
603   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
604   if(aContext.IsNull()) { 
605     di << "use 'vinit' command before " << argv[0] << "\n";
606     return 1;
607   }
608
609   BRep_Builder builder;
610   TopoDS_Compound comp;
611   TopoDS_CompSolid cs1, cs2;
612
613   builder.MakeCompound(comp);
614   builder.MakeCompSolid(cs1);
615   builder.MakeCompSolid(cs2);
616
617   TopoDS_Shape sh1 = BRepPrimAPI_MakeBox(gp_Pnt(0, 0, 0), 100, 100, 100).Shape();
618   TopoDS_Shape sh2 = BRepPrimAPI_MakeBox(gp_Pnt(100, 0, 0), 100, 100, 100).Shape();
619   TopoDS_Shape sh3 = BRepPrimAPI_MakeBox(gp_Pnt(200, 0, 0), 100, 100, 100).Shape();
620   builder.Add(cs1, sh1);
621   builder.Add(cs1, sh2);
622   builder.Add(cs1, sh3);
623
624   TopoDS_Shape sh4 = BRepPrimAPI_MakeBox(gp_Pnt(0, 500, 0), 100, 100, 100).Shape();
625   TopoDS_Shape sh5 = BRepPrimAPI_MakeBox(gp_Pnt(100, 500, 0), 100, 100, 100).Shape();
626   TopoDS_Shape sh6 = BRepPrimAPI_MakeBox(gp_Pnt(200, 500, 0), 100, 100, 100).Shape();
627
628   builder.Add(cs2, sh4);
629   builder.Add(cs2, sh5);
630   builder.Add(cs2, sh6);
631
632   builder.Add(comp, cs1);
633   builder.Add(comp, cs2);
634
635   Handle_AIS_Shape ais = new AIS_Shape(comp);
636   aContext->Display(ais);
637
638   TopExp_Explorer exp(comp,  TopAbs_COMPSOLID);
639   while(exp.More())
640   {
641     //printf("\n TopAbs_COMPSOLID is  there \n");
642     di << "\n TopAbs_COMPSOLID is  there \n";
643     exp.Next();
644   }
645
646   Handle (StdSelect_ShapeTypeFilter) filt = new StdSelect_ShapeTypeFilter(TopAbs_COMPSOLID);
647   aContext->AddFilter(filt);
648   aContext->CloseAllContexts();
649   aContext->OpenLocalContext(); 
650   aContext->ActivateStandardMode(TopAbs_SOLID);
651
652   return 0; 
653 }
654
655 static Standard_Integer  OCC138LC (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
656 {
657   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
658   if(aContext.IsNull()) { 
659     di << "use 'vinit' command before " << argv[0] << "\n";
660     return 1;
661   }
662
663   BRepPrimAPI_MakeBox box1(gp_Pnt(0, 0, 0),  gp_Pnt(100, 100, 100));
664   BRepPrimAPI_MakeBox box2(gp_Pnt(120, 120, 120),  gp_Pnt(300, 300,300));
665   BRepPrimAPI_MakeBox box3(gp_Pnt(320, 320, 320),  gp_Pnt(500, 500,500));
666
667   Handle_AIS_Shape ais1 = new AIS_Shape(box1.Shape());
668   Handle_AIS_Shape ais2 = new AIS_Shape(box2.Shape());
669   Handle_AIS_Shape ais3 = new AIS_Shape(box3.Shape());
670
671   aContext->Display(ais1);
672   aContext->Display(ais2);
673   aContext->Display(ais3);
674
675   aContext->AddOrRemoveSelected(ais1);
676   aContext->AddOrRemoveSelected(ais2);
677   aContext->AddOrRemoveSelected(ais3);
678
679   di << "\n No of selected = " << aContext->NbSelected();
680
681   aContext->InitSelected();
682
683   int count = 1;
684   while(aContext->MoreSelected())
685   {
686     di << "\n count is = %d" << count++;
687     Handle_AIS_InteractiveObject ais = aContext->SelectedInteractive();
688     aContext->AddOrRemoveSelected(ais);
689     aContext->InitSelected();
690   }
691
692   return 0; 
693 }
694
695 static Standard_Integer  OCC189 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
696 {
697   Handle(AIS_InteractiveContext) aContext1 = ViewerTest::GetAISContext();
698   if(aContext1.IsNull()) { 
699     di << "use 'vinit' command before " << argv[0] << "\n";
700     return 1;
701   }
702   Handle(AIS_InteractiveContext) aContext2 = ViewerTest::GetAISContext();
703   if(aContext2.IsNull()) { 
704     di << "use 'vinit' command before " << argv[0] << "\n";
705     return 1;
706   }
707
708   BRepPrimAPI_MakeBox box1(gp_Pnt(0, 0, 0),  gp_Pnt(100, 100, 100));
709   BRepPrimAPI_MakeBox box2(gp_Pnt(120, 120, 120),  gp_Pnt(300, 300, 300));
710   BRepPrimAPI_MakeBox box3(gp_Pnt(320, 320, 320),  gp_Pnt(500, 500, 500));
711
712   Handle_AIS_Shape ais1 = new AIS_Shape(box1.Shape());
713   Handle_AIS_Shape ais2 = new AIS_Shape(box2.Shape());
714   Handle_AIS_Shape ais3 = new AIS_Shape(box3.Shape());
715
716   aContext1->Display(ais1);
717   aContext1->Display(ais2);
718   aContext1->Display(ais3);
719
720   aContext2->Display(ais1);
721   aContext2->Display(ais2);
722   aContext2->Display(ais3);
723
724   aContext1->AddOrRemoveCurrentObject(ais1);
725   aContext1->AddOrRemoveCurrentObject(ais2);
726   aContext1->AddOrRemoveCurrentObject(ais3);
727
728   di << "\n Stage : 1";
729   di << "\n \t No of currents on aContext1 = " << aContext1->NbCurrents();
730   di << "\n \t No of currents on aContext2 = " << aContext2->NbCurrents() << "\n\n";
731
732   di << "\n aContext1->IsCurrent = " << (Standard_Integer) aContext1->IsCurrent(ais1) << ", aContext2->IsCurrent = " << (Standard_Integer) aContext2->IsCurrent(ais1) << " ";
733
734   aContext2->AddOrRemoveCurrentObject(ais1);
735   aContext2->AddOrRemoveCurrentObject(ais2);
736   aContext2->AddOrRemoveCurrentObject(ais3);
737
738   di << "\n Stage : 2";
739   di << "\n \t No of currents on aContext1 = " << aContext1->NbCurrents();
740   di << "\n \t No of currents on aContext2 = " << aContext2->NbCurrents() << "\n\n";
741
742   aContext1->InitCurrent();
743   int count1 = 1;
744   while(aContext1->MoreCurrent())
745   {
746     di << "\n count1 is = " << count1++;
747     Handle_AIS_InteractiveObject ais = aContext1->Current();
748     aContext1->AddOrRemoveCurrentObject(ais);
749     aContext1->InitCurrent();
750   }
751
752   di << "\n Stage : 3";
753   di << "\n \t No of currents on aContext1 = " << aContext1->NbCurrents();
754   di << "\n \t No of currents on aContext2 = " << aContext2->NbCurrents() << "\n\n";
755
756   aContext2->InitCurrent();
757   int count2 = 1;
758   while(aContext2->MoreCurrent())
759   {
760     di << "\n count2 is = " << count2++;
761     Handle_AIS_InteractiveObject ais = aContext2->Current();
762     aContext2->AddOrRemoveCurrentObject(ais);
763     aContext2->InitCurrent();
764   }
765
766   di << "\n\n Stage : 4";
767   di << "\n \t No of currents on aContext1 = " << aContext1->NbCurrents();
768   di << "\n \t No of currents on aContext2 = " << aContext2->NbCurrents();
769
770   return 0; 
771 }
772
773 static Standard_Integer  OCC389 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
774 {
775   if (argc < 3)
776   {
777     di<<"Usage: " << argv[0] << " name shape1 [shape2] ..."<<"\n";
778     return 1;
779   }
780
781   //if (Viewer2dTest::CurrentView().IsNull())
782   //{ 
783   //  cout << "2D AIS Viewer Initialization ..." << endl;
784   //  Viewer2dTest::ViewerInit(); 
785   //  cout << "Done" << endl;
786   //}
787
788   di<<"Begin!"<<"\n";
789   Handle(AIS2D_InteractiveContext) aContext = Viewer2dTest::GetAIS2DContext();
790   di<<"Check InteractiveContext"<<"\n";
791
792   if(aContext.IsNull()) {
793     di << "use 'v2dinit' command before " << argv[0] << "\n";
794     return 1;
795   }
796
797   di<<"CloseLocalContext"<<"\n";
798   aContext->CloseLocalContext();
799
800   di<<"OpenLocalContext"<<"\n";
801   aContext->OpenLocalContext();
802
803   TCollection_AsciiString name(argv[1]);
804
805   di<<"Found name"<<"\n";
806   Standard_Boolean IsBound = GetMapOfAIS2D().IsBound2(name);
807   if (IsBound)
808   {
809     di<<"Already displayed"<<"\n";
810   }
811   else // Create the AIS2D_ProjShape from a name
812   {
813     di<<"Create the AIS2D_ProjShape from a name; (1)"<<"\n";
814     Prs3d_Projector aPrs3dProjector(Standard_False, 1, -1.0,0.2,0.3, 0.0, 0.0, 0.0, -0.0,0.0,1.0);
815     di<<"Create the AIS2D_ProjShape from a name; (2)"<<"\n";
816     HLRAlgo_Projector aProjector = aPrs3dProjector.Projector();
817     di<<"Create the AIS2D_ProjShape from a name; (3)"<<"\n";
818     Handle(AIS2D_ProjShape) shp_2d =
819       new AIS2D_ProjShape(aProjector, 0, Standard_False, Standard_True);
820     di<<"Create the AIS2D_ProjShape from a name; (4)"<<"\n";
821
822     Standard_Integer i;
823     Standard_Boolean isValid = Standard_False;
824     for (i = 2; i < argc; i++)
825     {
826       di<<"i="<< i-1 <<"\n";
827       TopoDS_Shape aShape = DBRep::Get(argv[i]);
828       if (!aShape.IsNull())
829       {
830         shp_2d->Add( aShape );
831         di<<"Added shape: "<<argv[i]<<"\n";
832         isValid = Standard_True;
833       }
834     }
835
836     if (isValid)
837     {
838       di<<"Now start displaying..."<<"\n";
839       aContext->Display(shp_2d, 1, 1, Standard_True);
840
841       GetMapOfAIS2D().Bind(shp_2d, name);
842     }
843   }
844
845   aContext->UpdateCurrentViewer();
846   return 0; 
847 }
848
849 #include <BRepBndLib.hxx>
850 #include <Draw.hxx>
851 //=======================================================================
852 //function : OCC566
853 //purpose  : 
854 //=======================================================================
855 static Standard_Integer OCC566(Draw_Interpretor& di,Standard_Integer n,const char ** a)
856 {
857   if (n < 2) {
858     di<<"Usage: " << a[0] << " shape [ xmin ymin zmin xmax ymax zmax]"<<"\n";
859     return 1;
860   }
861   TopoDS_Shape S = DBRep::Get(a[1]);
862   if (S.IsNull()) return 1;
863   Bnd_Box B;
864   BRepBndLib::AddClose(S,B);
865   Standard_Real axmin,aymin,azmin,axmax,aymax,azmax;
866   B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
867   di << axmin<<" "<< aymin<<" "<< azmin<<" "<< axmax<<" "<< aymax<<" "<< azmax;
868   if (n >= 8) {
869     Draw::Set(a[2],axmin) ;
870     Draw::Set(a[3],aymin) ;
871     Draw::Set(a[4],azmin) ;
872     Draw::Set(a[5],axmax) ;
873     Draw::Set(a[6],aymax) ;
874     Draw::Set(a[7],azmax) ;
875   }
876   return 0;
877 }
878
879 #include <BRepFilletAPI_MakeFillet.hxx>
880 //=======================================================================
881 //function : OCC570
882 //purpose  : 
883 //=======================================================================
884 static Standard_Integer OCC570 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
885 {
886   if (argc < 2) {
887     di<<"Usage: " << argv[0] << " result"<<"\n";
888     return 1;
889   }
890
891   BRepPrimAPI_MakeBox mkBox(100.,100.,100.);
892   TopoDS_Shape aBox = mkBox.Shape();
893
894   TopExp_Explorer aExp;
895   aExp.Init(aBox,TopAbs_WIRE);
896   if (aExp.More())
897   {
898     TopoDS_Shape aWire = aExp.Current();
899
900     aExp.Init(aWire,TopAbs_EDGE);
901     TopoDS_Edge e1 = TopoDS::Edge(aExp.Current()); aExp.Next();
902     TopoDS_Edge e2 = TopoDS::Edge(aExp.Current()); aExp.Next();
903     TopoDS_Edge e3 = TopoDS::Edge(aExp.Current()); aExp.Next();
904     TopoDS_Edge e4 = TopoDS::Edge(aExp.Current());
905
906     try
907     {
908       OCC_CATCH_SIGNALS
909         BRepFilletAPI_MakeFillet mkFillet(aBox);
910       mkFillet.SetContinuity(GeomAbs_C1,.001);
911
912       // Setup variable fillet data
913       TColgp_Array1OfPnt2d t_pnt(1,4);
914       t_pnt.SetValue(1,gp_Pnt2d(0.0,5.0));
915       t_pnt.SetValue(2,gp_Pnt2d(0.3,15.0));
916       t_pnt.SetValue(3,gp_Pnt2d(0.7,15.0));
917       t_pnt.SetValue(4,gp_Pnt2d(1.0,5.0));
918
919       // HERE:
920       // It is impossible to build fillet if at least one edge
921       // with variable radius is added!!! If all are constant, everything is ok.
922       mkFillet.Add(t_pnt,e1);
923       mkFillet.Add(5.0,e2);
924       mkFillet.Add(t_pnt,e3);
925       mkFillet.Add(5.0,e4);
926
927       mkFillet.Build();
928       TopoDS_Shape aFinalShape = mkFillet.Shape();
929
930       DBRep::Set(argv[1],aFinalShape);
931     }
932     catch(Standard_Failure)
933     {
934       di << argv[0] << ": Exception in fillet"<<"\n";
935       return 2;
936     }
937   }
938
939   return 0;
940 }
941
942 #include <Law_Interpol.hxx>
943
944 static Standard_Real t3d = 1.e-4;
945 static Standard_Real t2d = 1.e-5;
946 static Standard_Real ta  = 1.e-2;
947 static Standard_Real fl  = 1.e-3;
948 static Standard_Real tapp_angle = 1.e-2;
949 static GeomAbs_Shape blend_cont = GeomAbs_C1;
950
951 static BRepFilletAPI_MakeFillet* Rakk = 0;
952 static BRepFilletAPI_MakeFillet* Rake = 0;
953 static char name[100];
954
955 static void printtolblend(Draw_Interpretor& di)
956 {
957   di<<"tolerance ang : "<<ta<<"\n";
958   di<<"tolerance 3d  : "<<t3d<<"\n";
959   di<<"tolerance 2d  : "<<t2d<<"\n";
960   di<<"fleche        : "<<fl<<"\n";
961
962   di<<"tolblend "<<ta<<" "<<t3d<<" "<<t2d<<" "<<fl<<"\n";
963 }
964
965 static Standard_Integer MKEVOL(Draw_Interpretor& di, 
966                                Standard_Integer narg, 
967                                const char ** a)
968 {
969   if(Rake != 0) {delete Rake; Rake = 0;}
970   printtolblend(di);
971   if (narg < 3) return 1;
972   TopoDS_Shape V = DBRep::Get(a[2]);
973   Rake = new BRepFilletAPI_MakeFillet(V);
974   Rake->SetParams(ta,t3d,t2d,t3d,t2d,fl);
975   Rake->SetContinuity(blend_cont, tapp_angle);
976   if (narg == 4) {
977     ChFi3d_FilletShape FSh = ChFi3d_Rational;
978     if (!strcasecmp(a[3], "Q")) {
979       FSh = ChFi3d_QuasiAngular;
980     }
981     else if (!strcasecmp(a[3], "P")) {
982       FSh = ChFi3d_Polynomial;
983     }
984     Rake->SetFilletShape(FSh);
985   }
986   strcpy(name, a[1]);
987   return 0;
988 }
989
990 static Standard_Integer UPDATEVOL(Draw_Interpretor& di, 
991                                   Standard_Integer narg, 
992                                   const char ** a)
993 {
994   if(Rake == 0){
995     di << "MakeFillet non initialise"<<"\n";
996     return 1 ;
997   }
998   if(narg%2 != 0 || narg < 4) return 1;
999   TColgp_Array1OfPnt2d uandr(1,(narg/2)-1);
1000   Standard_Real Rad, Par;
1001   TopoDS_Shape aLocalEdge(DBRep::Get(a[1],TopAbs_EDGE));
1002   TopoDS_Edge E = TopoDS::Edge(aLocalEdge);
1003   for (Standard_Integer ii = 1; ii <= (narg/2)-1; ii++){
1004     Par = atof(a[2*ii]);
1005     Rad = atof(a[2*ii + 1]);
1006     uandr.ChangeValue(ii).SetCoord(Par,Rad);
1007   }
1008   //HELPDESK: Add law creation
1009   const Standard_Boolean aLawMode = !strcmp(a[0],"OCC570evollaw");
1010   if (aLawMode)
1011   {
1012     di<<"INFO: Using law to define variable fillet"<<"\n";
1013     Handle(Law_Interpol) law = new Law_Interpol;
1014     law->Set(uandr);
1015     Rake->Add(law,E);
1016   }
1017   else
1018     Rake->Add(uandr,E);
1019
1020   return 0;
1021 }
1022
1023 static Standard_Integer BUILDEVOL(Draw_Interpretor& di,
1024                                   Standard_Integer, 
1025                                   const char **)
1026 {
1027   if(Rake == 0){
1028     di << "MakeFillet non initialise"<<"\n";
1029     return 1 ;
1030   }
1031   Rake->Build();
1032   if(Rake->IsDone()){
1033     TopoDS_Shape result = Rake->Shape();
1034     DBRep::Set(name,result);
1035     if(Rake != 0) {delete Rake; Rake = 0;}
1036     return 0;
1037   }
1038   if(Rake != 0) {delete Rake; Rake = 0;}
1039   return 1;
1040 }
1041
1042 #include <TColGeom_SequenceOfCurve.hxx>
1043 #include <Geom_TrimmedCurve.hxx>
1044 #include <GeomFill_NSections.hxx>
1045 #include <BRepBuilderAPI_MakeFace.hxx>
1046 //=======================================================================
1047 //function : OCC606
1048 //purpose  : 
1049 //=======================================================================
1050 static Standard_Integer OCC606 ( Draw_Interpretor& di, Standard_Integer n, const char ** a )
1051 {
1052   if (n != 3 && n != 4)
1053   {
1054     di << "Usage : " << a[0] << " result shape [-t]" << "\n";
1055     return 0;
1056   }
1057
1058   Standard_Boolean TrimMode = (n == 4);
1059   if (TrimMode) di<<"INFO: Using trimmed curves..."<<"\n";
1060   else di<<"INFO: Using non trimmed curves..."<<"\n";
1061
1062   TopoDS_Shape S = DBRep::Get(a[2]);
1063
1064   TopExp_Explorer t_exp(S,TopAbs_EDGE);
1065   TColGeom_SequenceOfCurve n_curves1;
1066   TColStd_SequenceOfReal np;
1067   Standard_Real param = 5.0;
1068
1069   for (;t_exp.More();t_exp.Next())
1070   {
1071     Standard_Real f,l;
1072     Handle(Geom_Curve) h_cur = BRep_Tool::Curve(TopoDS::Edge(t_exp.Current()),f,l);
1073     if (!h_cur.IsNull())
1074     {
1075       if (TrimMode) h_cur = new Geom_TrimmedCurve(h_cur,f,l);
1076
1077       n_curves1.Append(h_cur);
1078       np.Append(param);
1079       param += 1.0;
1080     }
1081   }
1082
1083   if (n_curves1.Length()>1 && !np.IsEmpty())
1084   {
1085     try
1086     {
1087       OCC_CATCH_SIGNALS
1088         GeomFill_NSections b_surface1(n_curves1, np);
1089       b_surface1.ComputeSurface();
1090       Handle(Geom_BSplineSurface) result_surf1 = b_surface1.BSplineSurface();
1091       if (!result_surf1.IsNull())
1092       {
1093         BRepBuilderAPI_MakeFace b_face1(result_surf1, Precision::Confusion());
1094         TopoDS_Face bsp_face1 = b_face1.Face();
1095         DBRep::Set(a[1],bsp_face1);
1096       }
1097     }
1098     catch(Standard_Failure)
1099     {
1100       di<<"ERROR: Exception in GeomFill_NSections"<<"\n";
1101     }
1102   }
1103
1104
1105   return 0;
1106 }
1107
1108 //=======================================================================
1109 //function : OCC813
1110 //purpose  : 
1111 //=======================================================================
1112 static Standard_Integer OCC813 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
1113 {
1114   if (argc < 3)
1115   {
1116     di << "Usage : " << argv[0] << " U V" << "\n";
1117     return 1;
1118   }
1119
1120   Standard_CString str;
1121   Standard_Real U = atof(argv[1]);
1122   Standard_Real V = atof(argv[2]);
1123
1124   //Between ellipse and point:
1125
1126   Handle(Geom_Ellipse) ell = new Geom_Ellipse(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
1127     gp_Dir(0.173648,0.984808,0.000000),
1128     gp_Dir(-0.932169,0.164367,-0.322560)), 150, 100);
1129   Handle(Geom_Plane) plne = new Geom_Plane(gp_Ax3(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
1130     gp_Dir(0.173648,0.984808,0.000000),
1131     gp_Dir(-0.932169,0.164367,-0.322560))) );
1132
1133   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1134
1135   gp_Pnt2d pt2d(U,V);
1136   gp_Pln pln = plne->Pln();
1137
1138   str = "OCC813_pnt"; DrawTrSurf::Set(str,pt2d);
1139
1140   Handle(Geom2d_Curve) curve2d = GeomAPI::To2d(ell,pln);
1141   Geom2dAdaptor_Curve acur(curve2d);
1142   Geom2dGcc_QualifiedCurve qcur(acur, GccEnt_outside);
1143
1144   str = "OCC813_ell"; DrawTrSurf::Set(str,curve2d);
1145   if(!aContext.IsNull()) {
1146     Handle(AIS_Shape) aisp = 
1147       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(curve2d, pln)).Edge());
1148     aContext->Display (aisp);
1149   }
1150
1151   //This does not give any solutions.
1152   Geom2dGcc_Lin2d2Tan lintan(qcur, pt2d, 0.1);
1153   di<<"OCC813 nb of solutions = "<<lintan.NbSolutions()<<"\n";
1154
1155   Standard_Character abuf[16];
1156   Standard_CString st = abuf;
1157
1158   Standard_Integer i;
1159   for(i=1;i<=lintan.NbSolutions();i++) {
1160     sprintf(abuf,"lintan_%d",i);
1161     Handle(Geom2d_Line) glin = new Geom2d_Line(lintan.ThisSolution(i));
1162     DrawTrSurf::Set(st,glin);
1163     if(!aContext.IsNull()) {
1164       Handle(AIS_Shape) aisp = 
1165         new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(glin, pln)).Edge());      aContext->Display (aisp);
1166     }
1167   }
1168
1169   return 0;
1170 }
1171
1172 //=======================================================================
1173 //function : OCC814
1174 //purpose  : 
1175 //=======================================================================
1176 static Standard_Integer OCC814 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
1177 {
1178   if (argc > 1)
1179   {
1180     di << "Usage : " << argv[0] << "\n";
1181     return 1;
1182   }
1183
1184   Standard_CString str;
1185
1186   //Between Ellipse and Circle:
1187
1188   Handle(Geom_Circle) cir = new Geom_Circle(gp_Ax2(gp_Pnt(823.687192,502.366825,478.960440),
1189     gp_Dir(0.173648,0.984808,0.000000),
1190     gp_Dir(-0.932169,0.164367,-0.322560)), 50);
1191   Handle(Geom_Ellipse) ell = new Geom_Ellipse(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
1192     gp_Dir(0.173648,0.984808,0.000000),
1193     gp_Dir(-0.932169,0.164367,-0.322560)), 150, 100);
1194   Handle(Geom_Plane) plne = new Geom_Plane(gp_Ax3(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
1195     gp_Dir(0.173648,0.984808,0.000000),
1196     gp_Dir(-0.932169,0.164367,-0.322560))) );
1197
1198   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1199
1200   gp_Pln pln = plne->Pln();
1201   Handle(Geom2d_Curve) curve2d = GeomAPI::To2d(ell, pln);
1202   Handle(Geom2d_Curve) fromcurve2d = GeomAPI::To2d(cir, pln);
1203
1204   str = "OCC814_cir"; DrawTrSurf::Set(str,curve2d);
1205   str = "OCC814_ell"; DrawTrSurf::Set(str,fromcurve2d);
1206   if(!aContext.IsNull()) {
1207     Handle(AIS_Shape) aisp = 
1208       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(curve2d, pln)).Edge());
1209     aContext->Display (aisp);
1210   }
1211   if(!aContext.IsNull()) {
1212     Handle(AIS_Shape) aisp = 
1213       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(fromcurve2d, pln)).Edge());
1214     aContext->Display (aisp);
1215   }
1216
1217   Geom2dAdaptor_Curve acur(curve2d), afromcur(fromcurve2d);
1218
1219   Geom2dGcc_QualifiedCurve qcur(acur, GccEnt_outside) ;
1220   Geom2dGcc_QualifiedCurve qfromcur(afromcur, GccEnt_outside) ;
1221
1222   //This does not give any solutions.
1223   Geom2dGcc_Lin2d2Tan lintan(qcur, qfromcur, 0.1);
1224   di<<"OCC814 nb of solutions = "<<lintan.NbSolutions()<<"\n";
1225
1226   Standard_Character abuf[16];
1227   Standard_CString st = abuf;
1228
1229   Standard_Integer i;
1230   for(i=1;i<=lintan.NbSolutions();i++) {
1231     sprintf(abuf,"lintan_%d",i);
1232     Handle(Geom2d_Line) glin = new Geom2d_Line(lintan.ThisSolution(i));
1233     DrawTrSurf::Set(st,glin);
1234     if(!aContext.IsNull()) {
1235       Handle(AIS_Shape) aisp = 
1236         new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(glin, pln)).Edge());      aContext->Display (aisp);
1237     }
1238   }
1239
1240   return 0;
1241 }
1242
1243 #include <ShapeAnalysis_Wire.hxx>
1244 #include <IntRes2d_SequenceOfIntersectionPoint.hxx>
1245 #include <TColgp_SequenceOfPnt.hxx>
1246 #include <ShapeFix_Wire.hxx>
1247 //=======================================================================
1248 //function : OCC884
1249 //purpose  : 
1250 //=======================================================================
1251 static Standard_Integer OCC884 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1252 {
1253   if (argc < 3 || argc > 5)
1254   {
1255     di << "Usage : " << argv[0] << " result shape [toler [maxtoler]]" << "\n";
1256     return 1;
1257   }
1258
1259   TopoDS_Shape shape = DBRep::Get(argv[2]); //read the shape
1260
1261   TopExp_Explorer exp(shape, TopAbs_EDGE);
1262   if (!exp.More())
1263   {
1264     di<<"Error: no edge found"<<"\n";
1265     return 1;
1266   }
1267
1268   BRep_Builder builder;
1269   TopoDS_Wire wire;
1270   builder.MakeWire(wire);
1271   builder.Add(wire, TopoDS::Edge(exp.Current()));
1272
1273   // HelpDesk: Create planar face if possible
1274   TopoDS_Face face = BRepBuilderAPI_MakeFace(wire,Standard_True);
1275
1276   Handle(ShapeAnalysis_Wire) advWA = new ShapeAnalysis_Wire;
1277   advWA->Load(wire);
1278   advWA->SetFace(face);
1279   advWA->SetPrecision(0.001);
1280
1281   IntRes2d_SequenceOfIntersectionPoint points2d;
1282   TColgp_SequenceOfPnt points3d;
1283
1284   di << "Info: CheckSelfIntersectingEdge = " << (Standard_Integer) advWA->CheckSelfIntersectingEdge(1,points2d,points3d) << "\n";
1285
1286   ShapeExtend_Status status = ShapeExtend_FAIL1;
1287   if (advWA->StatusSelfIntersection(status))
1288     di << "Info: No P Curve found in the edge\n";
1289
1290   status = ShapeExtend_FAIL2; 
1291   if (advWA->StatusSelfIntersection(status))
1292     di << "Info: No Vertices found in the edge\n";
1293
1294   status = ShapeExtend_DONE1; 
1295   if (advWA->StatusSelfIntersection(status))
1296     di << "Info: Self-intersection found in the edge\n";
1297
1298   Standard_Integer i, num = points2d.Length();
1299   di << "Info: No. of self-intersection points : " << num << "\n";
1300
1301   char str[80];
1302   Standard_CString name = str;
1303   for (i = 1; i <= num; ++i)
1304   { 
1305     gp_Pnt pt = points3d(i); 
1306     di << "Info: Intersecting pt : (" << pt.X() << ", " << pt.Y() << ", " << pt.Z() << ")\n";
1307     sprintf(str,"p_%d",i);
1308     DrawTrSurf::Set(name,pt);
1309   }
1310
1311   Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire;
1312   sfw->Load(wire);
1313   sfw->SetFace(face);
1314
1315   if (argc > 3) sfw->SetPrecision(atof(argv[3])/*0.1*/);
1316   if (argc > 4) sfw->SetMaxTolerance(atof(argv[4]));
1317   di << "Info: Precision is set to " << sfw->Precision() << "\n";
1318   di << "Info: MaxTolerance is set to " << sfw->MaxTolerance() << "\n";
1319
1320   sfw->ModifyTopologyMode() = 1;
1321   sfw->ModifyGeometryMode() = 1; 
1322   sfw->FixSelfIntersectingEdgeMode() = 1;
1323   //printf("Info: FixSelfIntersection = %d\n",sfw->FixSelfIntersection());
1324   di << "Info: FixSelfIntersection = " << (Standard_Integer) sfw->FixSelfIntersection() << "\n";
1325
1326   DBRep::Set(argv[1],sfw->Wire());
1327
1328   status = ShapeExtend_OK;
1329   if (sfw->StatusSelfIntersection(status))
1330     di << "Info: ShapeExtend_OK : No intersection found\n";
1331   //printf("Info: ShapeExtend_OK : No intersection found\n");
1332
1333   status = ShapeExtend_FAIL1;
1334   if (sfw->StatusSelfIntersection(status))
1335     di << "Info: ShapeExtend_FAIL1 : analysis failed (edge has no pcurve, or no vertices etc.)\n";
1336
1337   status = ShapeExtend_FAIL2;
1338   if (sfw->StatusSelfIntersection(status))
1339     di << "Info: ShapeExtend_FAIL2 : self-intersection was found, but not fixed because of limit of increasing tolerance (MaxTolerance)\n";
1340
1341   status = ShapeExtend_FAIL3;
1342   if (sfw->StatusSelfIntersection(status))
1343     di << "Info: ShapeExtend_FAIL3 : intercestion of non adjacent edges found, but not fixed because of limit of increasing tolerance (MaxTolerance)\n";
1344
1345   status = ShapeExtend_DONE1;
1346   if (sfw->StatusSelfIntersection(status))
1347     di << "Info: ShapeExtend_DONE1 : tolerance of vertex was increased to fix self-intersection\n";
1348
1349   status = ShapeExtend_DONE2;
1350   if (sfw->StatusSelfIntersection(status))
1351     di << "Info: ShapeExtend_DONE2 : vertex was moved to fix self-intersection\n";
1352
1353   status = ShapeExtend_DONE3;
1354   if (sfw->StatusSelfIntersection(status))
1355     di << "Info: ShapeExtend_DONE3 : some edges were removed because of intersection\n";
1356
1357   status = ShapeExtend_DONE4;
1358   if (sfw->StatusSelfIntersection(status))
1359     di << "Info: ShapeExtend_DONE4 : pcurve(s) was(were) modified\n";
1360
1361   status = ShapeExtend_DONE5;
1362   if (sfw->StatusSelfIntersection(status))
1363     di << "Info: ShapeExtend_DONE5 : non adjacent intersection fixed by increasing tolerance of vertex(vertices)\n";
1364
1365   status = ShapeExtend_DONE6;
1366   if (sfw->StatusSelfIntersection(status))
1367     di << "Info: ShapeExtend_DONE6 : tolerance of edge was increased to hide intersection\n";
1368
1369   status = ShapeExtend_DONE7;
1370   if (sfw->StatusSelfIntersection(status))
1371     di << "Info: ShapeExtend_DONE7 : range of some edges was decreased to avoid intersection\n";
1372
1373   return 0;
1374 }
1375
1376 #include <Aspect_FillMethod.hxx>
1377 //=======================================================================
1378 //function : OCC1188
1379 //purpose  : 
1380 //=======================================================================
1381 static Standard_Integer OCC1188 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1382 {
1383   if (argc < 2 || argc > 3)
1384   {
1385     di << "Usage : " << argv[0] << " imagefile [filltype] : Load image as background" << "\n";
1386     return 1;
1387   }
1388
1389   Handle(AIS_InteractiveContext) AISContext = ViewerTest::GetAISContext();
1390   if(AISContext.IsNull()) 
1391   {
1392     di << "use 'vinit' command before " << argv[0] << "\n";
1393     return 1;
1394   }
1395
1396   Aspect_FillMethod aFillType = Aspect_FM_CENTERED;
1397   if (argc == 3)
1398   {
1399     const char* szType = argv[2];
1400     if      (strcmp(szType, "NONE"    ) == 0) aFillType = Aspect_FM_NONE;
1401     else if (strcmp(szType, "CENTERED") == 0) aFillType = Aspect_FM_CENTERED;
1402     else if (strcmp(szType, "TILED"   ) == 0) aFillType = Aspect_FM_TILED;
1403     else if (strcmp(szType, "STRETCH" ) == 0) aFillType = Aspect_FM_STRETCH;
1404     else
1405     {
1406       di << "Wrong fill type : " << szType << "\n";
1407       di << "Must be one of CENTERED, TILED, STRETCH, NONE" << "\n";
1408       return 1;
1409     }
1410   }
1411
1412   Handle(V3d_View) V3dView = ViewerTest::CurrentView();
1413   V3dView->SetBackgroundImage(argv[1], aFillType, Standard_True);
1414
1415   return 0;
1416 }
1417
1418 #include <AIS_Drawer.hxx>
1419 #include <Graphic3d_MaterialAspect.hxx>
1420 #include <Prs3d_ShadingAspect.hxx>
1421 //=======================================================================
1422 //function : OCC1174_1
1423 //purpose  : 
1424 //=======================================================================
1425 static Standard_Integer OCC1174_1 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1426 {
1427   if (argc != 2)
1428   {
1429     di << "Usage : " << argv[0] << " shape" << "\n";
1430     return 1;
1431   }
1432
1433   Handle(AIS_InteractiveContext) AISContext = ViewerTest::GetAISContext();
1434   if(AISContext.IsNull()) 
1435   {
1436     di << "use 'vinit' command before " << argv[0] << "\n";
1437     return 1;
1438   }
1439
1440   TopoDS_Shape sh = DBRep::Get(argv[1]);
1441
1442   Handle(AIS_Shape) ais = new AIS_Shape(sh);
1443
1444   Quantity_Color colf(0.0, 0.4, 0.0, Quantity_TOC_RGB);
1445   Quantity_Color colb(0.0, 0.0, 0.6, Quantity_TOC_RGB);
1446
1447   Handle(Prs3d_ShadingAspect) sa = ais->Attributes()->ShadingAspect();
1448
1449   Graphic3d_MaterialAspect front = sa->Material(Aspect_TOFM_FRONT_SIDE);
1450   front.SetAmbientColor(colf);
1451   front.SetDiffuseColor(colf);
1452   front.SetSpecularColor(colf);
1453   front.SetEmissiveColor(colf);
1454   front.SetTransparency(0.0);
1455   sa->SetMaterial(front,Aspect_TOFM_FRONT_SIDE);
1456
1457   Graphic3d_MaterialAspect back = sa->Material(Aspect_TOFM_BACK_SIDE);
1458   back.SetAmbientColor(colb);
1459   back.SetDiffuseColor(colb);
1460   back.SetSpecularColor(colb);
1461   back.SetEmissiveColor(colb);
1462   back.SetTransparency(0.0);
1463   sa->SetMaterial(back,Aspect_TOFM_BACK_SIDE);
1464
1465   AISContext->Display(ais,1,0);
1466
1467   Standard_Real r, g, b; 
1468   sa->Color(Aspect_TOFM_FRONT_SIDE).Values(r,g,b, Quantity_TOC_RGB);
1469   di << "Info: color on front side (" << r << "," << g << "," << b << ")\n";
1470   sa->Color(Aspect_TOFM_BACK_SIDE).Values(r,g,b, Quantity_TOC_RGB);
1471   di << "Info: color on back side (" << r << "," << g << "," << b << ")\n";
1472
1473   return 0;
1474 }
1475
1476 //=======================================================================
1477 //function : OCC1174_2
1478 //purpose  : 
1479 //=======================================================================
1480 static Standard_Integer OCC1174_2 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1481 {
1482   if (argc != 2)
1483   {
1484     di << "Usage : " << argv[0] << " shape" << "\n";
1485     return 1;
1486   }
1487
1488   Handle(AIS_InteractiveContext) AISContext = ViewerTest::GetAISContext();
1489   if(AISContext.IsNull()) 
1490   {
1491     di << "use 'vinit' command before " << argv[0] << "\n";
1492     return 1;
1493   }
1494
1495   TopoDS_Shape sh = DBRep::Get(argv[1]);
1496
1497   Handle(AIS_Shape) ais = new AIS_Shape(sh); 
1498   AISContext->Display(ais,1,0); 
1499   AISContext->SetMaterial(ais,Graphic3d_NOM_SHINY_PLASTIC); 
1500
1501   Quantity_Color colf(0.0, 0.4, 0.0, Quantity_TOC_RGB); 
1502   Quantity_Color colb(0.0, 0.0, 0.6, Quantity_TOC_RGB); 
1503   Handle(Prs3d_ShadingAspect) sa = ais->Attributes()->ShadingAspect(); 
1504
1505   Graphic3d_MaterialAspect front = sa->Material(Aspect_TOFM_FRONT_SIDE); 
1506   front.SetAmbientColor(colf); 
1507   front.SetDiffuseColor(colf); 
1508   front.SetSpecularColor(colf); 
1509   front.SetEmissiveColor(colf); 
1510   front.SetTransparency(0.4); 
1511   sa->SetMaterial(front,Aspect_TOFM_FRONT_SIDE); 
1512
1513   Graphic3d_MaterialAspect back = sa->Material(Aspect_TOFM_BACK_SIDE); 
1514   back.SetAmbientColor(colb); 
1515   back.SetDiffuseColor(colb); 
1516   back.SetSpecularColor(colb); 
1517   back.SetEmissiveColor(colb); 
1518   back.SetTransparency(0.2); 
1519   sa->SetMaterial(back,Aspect_TOFM_BACK_SIDE); 
1520
1521   AISContext->Redisplay(ais,1,0);
1522
1523   return 0;
1524 }
1525
1526 #include <Graphic2d_Segment.hxx>
1527 #include <Prs2d_Length.hxx>
1528 #include <Prs2d_AspectLine.hxx>
1529 #include <AIS2D_InteractiveObject.hxx>
1530 //=======================================================================
1531 //  OCC672:
1532 //         OCC672_Length
1533 //         OCC672_Angle
1534 //         OCC672_Diameter
1535 //         OCC672_ShortLength 
1536 //=======================================================================
1537
1538 static Standard_Integer OCC672_Length (Draw_Interpretor& di, Standard_Integer n,const char ** argv)
1539 {
1540   if (n != 17) {
1541     di << " wrong parameters !\n";
1542     di << "must be : OCC672_Length x1 y1 x2 y2 str scale_text length_dist angle_arrow length_arrow x3 y3 x4 y4 TxtAngle TxtPosH TxtPosV\n";
1543     return 1;
1544   }
1545
1546   Standard_Real x1=atof(argv[1]);
1547   Standard_Real y1=atof(argv[2]);
1548   Standard_Real x2=atof(argv[3]);
1549   Standard_Real y2=atof(argv[4]);
1550   Standard_Real x3=atof(argv[10]);
1551   Standard_Real y3=atof(argv[11]);
1552   Standard_Real x4=atof(argv[12]);
1553   Standard_Real y4=atof(argv[13]);
1554   Standard_Real txtAngle=atof(argv[14]);
1555   Standard_Real txtPosH=atof(argv[15]);
1556   Standard_Real txtPosV=atof(argv[16]);
1557
1558   di<<"Begin!\n";
1559   Handle(AIS2D_InteractiveContext) aContext = Viewer2dTest::GetAIS2DContext();
1560   di<<"Check InteractiveContext\n";
1561
1562   if(aContext.IsNull()) {
1563     di << "use 'v2dinit' command before " << argv[0] << "\n";
1564     return 1;
1565   }
1566
1567   di<<"CloseLocalContext\n";
1568   aContext->CloseLocalContext();
1569
1570   di<<"OpenLocalContext\n";
1571   aContext->OpenLocalContext();
1572   //step0 end  
1573   Handle(AIS2D_InteractiveObject) aIO10 = new AIS2D_InteractiveObject();
1574   Handle(Graphic2d_Segment) theSegmA, theSegmB, theSegmC, theSegmD;
1575
1576   theSegmA = new Graphic2d_Segment( aIO10, x1, y1, x2, y2 );
1577   theSegmB = new Graphic2d_Segment( aIO10, x2, y2, x3, y3 );
1578   theSegmC = new Graphic2d_Segment( aIO10, x3, y3, x4, y4 );
1579   theSegmD = new Graphic2d_Segment( aIO10, x4, y4, x1, y1 );
1580
1581   aContext->Display( aIO10, Standard_True ); 
1582
1583   Handle(Prs2d_Length) length = new Prs2d_Length(aIO10, gp_Pnt2d(x1, y1), gp_Pnt2d(x2, y2), argv[5],
1584     atof(argv[6]), atof(argv[7]), Prs2d_TOD_AUTOMATIC, 
1585     atof(argv[8]), atof(argv[9]), Prs2d_TOA_FILLED,
1586     Prs2d_AS_BOTHAR);
1587
1588   length->SetTextRelPos(txtPosH,txtPosV);
1589   length->SetTextRelAngle(txtAngle);
1590
1591   Handle(Prs2d_AspectLine) theAspect = new Prs2d_AspectLine(Quantity_NOC_WHITE,
1592     Aspect_TOL_SOLID, Aspect_WOL_THIN);
1593   aIO10->SetAspect( theAspect,  length);
1594   aContext->Display( aIO10, Standard_True );
1595
1596   Handle(Prs2d_Length) lengthB = new Prs2d_Length(aIO10, gp_Pnt2d(x2, y2), gp_Pnt2d(x3, y3), argv[5],
1597     atof(argv[6]), atof(argv[7]), Prs2d_TOD_AUTOMATIC, 
1598     atof(argv[8]), atof(argv[9]), Prs2d_TOA_FILLED,
1599     Prs2d_AS_BOTHAR);
1600
1601   lengthB->SetTextRelPos(txtPosH,txtPosV);
1602   lengthB->SetTextRelAngle(txtAngle);
1603
1604   aIO10->SetAspect( theAspect,  lengthB);
1605   aContext->Display( aIO10, Standard_True );
1606
1607   Handle(Prs2d_Length) lengthC = new Prs2d_Length(aIO10, gp_Pnt2d(x3, y3), gp_Pnt2d(x4, y4), argv[5],
1608     atof(argv[6]), atof(argv[7]), Prs2d_TOD_AUTOMATIC, 
1609     atof(argv[8]), atof(argv[9]), Prs2d_TOA_FILLED,
1610     Prs2d_AS_BOTHAR);
1611
1612   lengthC->SetTextRelPos(txtPosH,txtPosV);
1613   lengthC->SetTextRelAngle(txtAngle);
1614
1615   aIO10->SetAspect( theAspect,  lengthC);
1616   aContext->Display( aIO10, Standard_True );
1617
1618   Handle(Prs2d_Length) lengthD = new Prs2d_Length(aIO10, gp_Pnt2d(x4, y4), gp_Pnt2d(x1, y1), argv[5],
1619     atof(argv[6]), atof(argv[7]), Prs2d_TOD_AUTOMATIC, 
1620     atof(argv[8]), atof(argv[9]), Prs2d_TOA_FILLED,
1621     Prs2d_AS_BOTHAR);
1622
1623   lengthD->SetTextRelPos(txtPosH,txtPosV);
1624   lengthD->SetTextRelAngle(txtAngle);
1625
1626   aIO10->SetAspect( theAspect,  lengthD);
1627   aContext->Display( aIO10, Standard_True );
1628   return 0;
1629 }
1630
1631 #include <Prs2d_Angle.hxx>
1632 static Standard_Integer OCC672_Angle (Draw_Interpretor& di, Standard_Integer n,const char ** argv)
1633 {
1634   if (n != 15) {
1635     di << "Wrong parameters !\n";
1636     di << "must be : OCC672_Angle x1 y1 x2 y2 x3 y3 aRadius aText aTxtScale anArrAngle anArrLength txtAngle txtPosH txtPosV  \n";
1637     return 1;
1638   }
1639
1640   Standard_Real x1=atof(argv[1]);
1641   Standard_Real y1=atof(argv[2]);
1642   Standard_Real x2=atof(argv[3]);
1643   Standard_Real y2=atof(argv[4]);
1644   Standard_Real x3=atof(argv[5]);
1645   Standard_Real y3=atof(argv[6]);
1646   Standard_Real aRadius=atof(argv[7]);
1647   Standard_Real aTxtScale=atof(argv[9]);
1648   Standard_Real anArrAngle=atof(argv[10]);
1649   Standard_Real anArrLength=atof(argv[11]);
1650   Standard_Real txtAngle=atof(argv[12]);
1651   Standard_Real txtPosH=atof(argv[13]);
1652   Standard_Real txtPosV=atof(argv[14]);
1653
1654   di<<"Begin!\n";
1655   Handle(AIS2D_InteractiveContext) aContext = Viewer2dTest::GetAIS2DContext();
1656   di<<"Check InteractiveContext\n";
1657
1658   if(aContext.IsNull()) {
1659     di << "use 'v2dinit' command before " << argv[0] << "\n";
1660     return 1;
1661   }
1662
1663   di<<"CloseLocalContext\n";
1664   aContext->CloseLocalContext();
1665
1666   di<<"OpenLocalContext\n";
1667   aContext->OpenLocalContext();
1668
1669   Handle(AIS2D_InteractiveObject) aIO10 = new AIS2D_InteractiveObject();
1670   Handle(Graphic2d_Segment) theSegmA, theSegmB;
1671   theSegmA = new Graphic2d_Segment( aIO10, x1, y1, x2, y2 );
1672   theSegmB = new Graphic2d_Segment( aIO10, x1, y1, x3, y3 );
1673   aContext->Display( aIO10, Standard_True ); 
1674
1675   Handle(Prs2d_Angle) angle = new Prs2d_Angle(aIO10, gp_Pnt2d(x1, y1), gp_Pnt2d(x2, y2),gp_Pnt2d(x3,y3),
1676     aRadius,argv[8],aTxtScale,anArrAngle,anArrLength,Prs2d_TOA_FILLED,
1677     Prs2d_AS_BOTHAR);
1678
1679   angle->SetTextRelPos(txtPosH,txtPosV);
1680   angle->SetTextAbsAngle(txtAngle);
1681
1682   aContext->Display( aIO10, Standard_True );
1683   return 0; 
1684 }
1685
1686 #include <Prs2d_Diameter.hxx>
1687 #include <gp_Circ2d.hxx>
1688 #include <Graphic2d_Circle.hxx>
1689 #include <Prs2d_AspectText.hxx>
1690 static Standard_Integer OCC672_Diameter (Draw_Interpretor& di, Standard_Integer n,const char ** argv)
1691 {
1692   if (n != 13) {
1693     di << "Wrong parameters !\n";
1694     di << "must be : OCC672_Diameter x1 y1 cx cy radius aText aTxtScale anArrAngle anArrLength txtAngle txtPosH txtPosV\n";
1695     return 1;
1696   }
1697
1698   Standard_Real x1=atof(argv[1]);
1699   Standard_Real y1=atof(argv[2]);
1700   Standard_Real cx=atof(argv[3]);    //coordinats for circle 
1701   Standard_Real cy=atof(argv[4]);    //
1702   Standard_Real radius=atof(argv[5]);//radius for circle 
1703   Standard_Real aTxtScale=atof(argv[7]);
1704   Standard_Real anArrAngle=atof(argv[8]);
1705   Standard_Real anArrLength=atof(argv[9]);
1706
1707   Standard_Real txtAngle=atof(argv[10]);
1708   Standard_Real txtPosH=atof(argv[11]);
1709   Standard_Real txtPosV=atof(argv[12]);
1710
1711
1712   di<<"Begin!\n";
1713   Handle(AIS2D_InteractiveContext) aContext = Viewer2dTest::GetAIS2DContext();
1714   di<<"Check InteractiveContext\n";
1715
1716   if(aContext.IsNull()) {
1717     di << "use 'v2dinit' command before " << argv[0] << "\n";
1718     return 1;
1719   }
1720
1721   di<<"CloseLocalContext\n";
1722   aContext->CloseLocalContext();
1723
1724   di<<"OpenLocalContext\n";
1725   aContext->OpenLocalContext();
1726   Handle(AIS2D_InteractiveObject) aIO10 = new AIS2D_InteractiveObject(); 
1727
1728   Handle (Graphic2d_Circle) HCirc;
1729   HCirc= new Graphic2d_Circle(aIO10,cx,cy,radius);
1730
1731   gp_Dir2d myDir2d(1,0);
1732   gp_Pnt2d myPnt2d(cx,cy);
1733   gp_Ax2d myAx2d(myPnt2d,myDir2d);
1734   gp_Circ2d myCirc(myAx2d,radius);
1735
1736   Handle(Prs2d_Diameter) diameter = new Prs2d_Diameter(aIO10, gp_Pnt2d(x1, y1),myCirc,
1737     argv[6],aTxtScale,anArrAngle,anArrLength,
1738     Prs2d_TOA_FILLED,
1739     Prs2d_AS_BOTHAR,Standard_True);
1740
1741   diameter->SetTextRelPos(txtPosH,txtPosV);
1742
1743   diameter->SetTextAbsAngle(txtAngle);
1744
1745   aContext->Display( aIO10, Standard_True );
1746
1747   return 0;
1748 }
1749
1750 #include <gp_Lin2d.hxx>
1751 static Standard_Integer OCC672_ShortLength (Draw_Interpretor& di, Standard_Integer n,const char ** argv)
1752 {
1753   if (n != 17) {
1754     di << " wrong parameters !\n";
1755     di << "must be : OCC672_ShortLength x1 y1 x2 y2 str scale_text length_dist angle_arrow length_arrow x3 y3 x4 y4 TxtAngle TxtPosH TxtPosV\n";
1756     return 1;
1757   }
1758
1759   Standard_Real x1=atof(argv[1]);
1760   Standard_Real y1=atof(argv[2]);
1761   Standard_Real x2=atof(argv[3]);
1762   Standard_Real y2=atof(argv[4]);
1763   Standard_Real x3=atof(argv[10]);
1764   Standard_Real y3=atof(argv[11]);
1765   Standard_Real x4=atof(argv[12]);
1766   Standard_Real y4=atof(argv[13]);
1767   Standard_Real txtAngle=atof(argv[14]);
1768   Standard_Real txtPosH=atof(argv[15]);
1769   Standard_Real txtPosV=atof(argv[16]);
1770
1771
1772   di<<"Begin!\n";
1773   Handle(AIS2D_InteractiveContext) aContext = Viewer2dTest::GetAIS2DContext();
1774   di<<"Check InteractiveContext\n";
1775
1776   if(aContext.IsNull()) {
1777     di << "use 'v2dinit' command before " << argv[0] << "\n";
1778     return 1;
1779   }
1780
1781   di<<"CloseLocalContext\n";
1782   aContext->CloseLocalContext();
1783
1784   di<<"OpenLocalContext\n";
1785   aContext->OpenLocalContext();
1786
1787   Handle(AIS2D_InteractiveObject) aIO10 = new AIS2D_InteractiveObject();
1788   Handle(Graphic2d_Segment) theSegmA, theSegmB, theSegmC, theSegmD;
1789   theSegmB = new Graphic2d_Segment( aIO10, x2, y2, x3, y3 );
1790   theSegmD = new Graphic2d_Segment( aIO10, x4, y4, x1, y1 );
1791
1792   aContext->Display( aIO10, Standard_True ); 
1793
1794   gp_Pnt2d p3(x3,y3),p2(x2,y2),p1(x1,y1);
1795   gp_Vec2d theVec(p3,p2);
1796   gp_Dir2d theDir(theVec);
1797   gp_Lin2d theLin1(p1,theDir);
1798   gp_Lin2d theLin2(p3,theDir);
1799
1800   Handle(Prs2d_Length) length = new Prs2d_Length(aIO10, p1, theLin2, argv[5],
1801     atof(argv[6]), 5.0, 20.0, 5.,
1802     Prs2d_TOA_CLOSED,Prs2d_AS_BOTHAR,
1803     Standard_True);
1804
1805
1806   Handle(Prs2d_AspectLine) theAspect = new Prs2d_AspectLine(Quantity_NOC_WHITE,
1807     Aspect_TOL_SOLID, Aspect_WOL_THIN);
1808   length->SetTextRelPos(txtPosH,txtPosV);
1809   length->SetTextRelAngle(txtAngle);
1810
1811   aIO10->SetAspect( theAspect,  length);
1812   aContext->Display( aIO10, Standard_True );
1813
1814   return 0;
1815 }
1816
1817 #include <TopoDS_Solid.hxx>
1818 #include <BRepFeat_MakeDPrism.hxx>
1819 //=======================================================================
1820 //function : OCCN1
1821 //purpose  : FEATURE OPERATION (BOSS and SLOT)
1822 //=======================================================================
1823 static Standard_Integer OCCN1 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1824 {
1825   if (argc != 4)
1826   {
1827     di << "Usage : " << argv[0] << " angle fuse(1 for boss / 0 for slot) length" << "\n";
1828     return 1;
1829   }
1830
1831   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1832   if(aContext.IsNull()) 
1833   {
1834     di << "use 'vinit' command before " << argv[0] << "\n";
1835     return 1;
1836   }
1837   Standard_Real    angle  = atof(argv[1]);
1838   Standard_Integer fuse   = atoi(argv[2]);
1839   Standard_Real    length = atof(argv[3]);
1840
1841   BRepBuilderAPI_MakeEdge edge1(gp_Pnt(0, 0, 0), gp_Pnt(50, 0, 0));
1842   BRepBuilderAPI_MakeEdge edge2(gp_Pnt(50, 0, 0), gp_Pnt(50, 50, 0));
1843   BRepBuilderAPI_MakeEdge edge3(gp_Pnt(50, 50, 0), gp_Pnt(0, 50, 0));
1844   BRepBuilderAPI_MakeEdge edge4(gp_Pnt(0, 50, 0), gp_Pnt(0, 0, 0));
1845   TopoDS_Edge ted1 = edge1.Edge();
1846   TopoDS_Edge ted2 = edge2.Edge();
1847   TopoDS_Edge ted3 = edge3.Edge();
1848   TopoDS_Edge ted4 = edge4.Edge();
1849
1850   BRepBuilderAPI_MakeWire wire(ted1, ted2, ted3, ted4);
1851   TopoDS_Wire twire = wire.Wire();
1852
1853   BRepBuilderAPI_MakeFace face(twire);
1854   TopoDS_Face tface = face.Face();
1855   ////////Handle_AIS_Shape face_ais = new AIS_Shape( tface );
1856   ////////aContext->Display(face_ais);
1857
1858   TopoDS_Solid box;
1859   BRepPrimAPI_MakeBox block(gp_Pnt(-50, -50, -150), 150, 150, 100);
1860   box = block.Solid();
1861
1862   TopoDS_Face top_face;
1863   top_face = block.TopFace();
1864
1865   BRepFeat_MakeDPrism slot(box, tface, top_face, angle, fuse,
1866     Standard_True);
1867
1868   slot.Perform(length);
1869   TopoDS_Shape S2 = slot.Shape();
1870
1871   DBRep::Set("OCCN1_tface",tface);
1872   DBRep::Set("OCCN1_slot",S2);
1873
1874   return 0;
1875 }
1876
1877 #include <BRepPrimAPI_MakeCylinder.hxx>
1878 #include <BRepPrimAPI_MakeSphere.hxx>
1879 #include <BRepAlgoAPI_Section.hxx>
1880 #include <BRepAlgo_Section.hxx>
1881 //=======================================================================
1882 //function : OCCN2
1883 //purpose  : BOOLEAN OPERATION
1884 //=======================================================================
1885 static Standard_Integer OCCN2 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1886 {
1887   if (argc > 2) {
1888     di << "Usage : " << argv[0] << " [BRepAlgoAPI/BRepAlgo = 1/0]" << "\n";
1889     return 1;
1890   }
1891   Standard_Boolean IsBRepAlgoAPI = Standard_True;
1892   if (argc == 2) {
1893     Standard_Integer IsB = atoi(argv[1]);
1894     if (IsB != 1) {
1895       IsBRepAlgoAPI = Standard_False;
1896 #if ! defined(BRepAlgo_def04)
1897 #endif
1898     }
1899   }
1900
1901   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1902   if(aContext.IsNull()) 
1903   {
1904     di << "use 'vinit' command before " << argv[0] << "\n";
1905     return 1;
1906   }
1907
1908   BRepPrimAPI_MakeCylinder cylinder(50,200);
1909   TopoDS_Shape cylinder_sh = cylinder.Shape();
1910
1911   BRepPrimAPI_MakeSphere sphere(gp_Pnt(60,0,100),50);
1912   TopoDS_Shape sphere_sh = sphere.Shape();
1913
1914   //BRepAlgoAPI_Section section(cylinder_sh, sphere_sh);
1915   //TopoDS_Shape shape = section.Shape();
1916
1917   TopoDS_Shape shape;
1918   if (IsBRepAlgoAPI) {
1919     di << "BRepAlgoAPI_Section section(cylinder_sh, sphere_sh)" <<"\n";
1920     BRepAlgoAPI_Section section(cylinder_sh, sphere_sh);
1921     section.Build();
1922     if(!section.IsDone()){
1923       di << "Error performing intersection: not done." << "\n";
1924     }
1925     shape = section.Shape();
1926   } else {
1927     di << "BRepAlgo_Section section(cylinder_sh, sphere_sh)" <<"\n";
1928     BRepAlgo_Section section(cylinder_sh, sphere_sh);
1929     section.Build();
1930     if(!section.IsDone()){
1931       di << "Error performing intersection: not done." << "\n";
1932     }
1933     shape = section.Shape();
1934   }
1935
1936   DBRep::Set("OCCN2_cylinder",cylinder_sh);
1937   DBRep::Set("OCCN2_sphere",sphere_sh);
1938   DBRep::Set("OCCN2_section",shape);
1939
1940   return 0;
1941 }
1942
1943 #include <TColgp_Array1OfPnt.hxx>
1944 #include <Geom_BezierCurve.hxx>
1945
1946 static Standard_Integer OCC2569 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1947 {
1948   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1949   if(aContext.IsNull()) 
1950   {
1951     di << "use 'vinit' command before " << argv[0] << "\n";
1952     return 1;
1953   }
1954   if(argc != 3) {
1955     di << "Usage : " << argv[0] << " nbpoles result" << "\n";
1956     return 1;
1957   }
1958
1959   int poles=atoi(argv[1]); 
1960
1961   TColgp_Array1OfPnt arr(1, poles); 
1962   for(int i=1; i<=poles; i++) 
1963     arr.SetValue(i, gp_Pnt(i+10, i*2+20, i*3+45)); 
1964
1965   Handle_Geom_BezierCurve bez = new Geom_BezierCurve(arr); 
1966   if(bez.IsNull()) {
1967     di << "\n The curve is not created.\n";
1968   } else {
1969     di << "\n Degree = " << bez->Degree() << "\n";   
1970   }
1971   TopoDS_Edge sh = BRepBuilderAPI_MakeEdge(bez).Edge(); 
1972   Handle_AIS_Shape ais = new AIS_Shape(sh); 
1973   aContext->Display(ais); 
1974   DrawTrSurf::Set(argv[2],bez);
1975   return 0;
1976 }
1977
1978 #include <TopTools_IndexedMapOfShape.hxx>
1979 #include <TopExp.hxx>
1980
1981 static Standard_Integer OCC1642 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1982 {
1983   if (argc != 10) {
1984     di<<"Usage: " << argv[0] << " FinalWare FinalFace InitWare InitFace shape FixReorder FixDegenerated FixConnected FixSelfIntersection"<<"\n";
1985     return 1;
1986   }
1987
1988   TopoDS_Shape shape = DBRep::Get(argv[5]);
1989
1990   TopExp_Explorer exp(shape, TopAbs_WIRE);
1991   TopoDS_Shape wire = exp.Current();
1992
1993   Handle(ShapeAnalysis_Wire) advWA = new ShapeAnalysis_Wire;
1994   advWA->Load(TopoDS::Wire(wire));
1995
1996   DBRep::Set(argv[3],wire);
1997
1998   TopoDS_Face face =
1999     BRepBuilderAPI_MakeFace(TopoDS::Wire(wire),Standard_True);
2000
2001   DBRep::Set(argv[4],face);
2002
2003   advWA->SetFace(face);
2004   float precision_to_ana = 0.0001;
2005   advWA->SetPrecision(precision_to_ana);
2006
2007   TopTools_IndexedMapOfShape M;
2008   TopExp::MapShapes(wire, TopAbs_EDGE, M);
2009
2010   int j= 1;
2011   for(j= 1; j<=M.Extent(); ++j) {
2012     Standard_Integer num=1;
2013     IntRes2d_SequenceOfIntersectionPoint points2d;
2014     TColgp_SequenceOfPnt points3d;
2015
2016     cout << "\n j =" << j << ",  CheckSelfIntersectingEdge = " <<
2017       advWA->CheckSelfIntersectingEdge(j, points2d, points3d);
2018
2019
2020     ShapeExtend_Status status = ShapeExtend_FAIL1;
2021     if(advWA->StatusSelfIntersection(status))
2022
2023       status = ShapeExtend_FAIL2;
2024     if(advWA->StatusSelfIntersection(status))
2025       di << "\n No Vertices found in the edge";
2026
2027     status = ShapeExtend_DONE1;
2028     if(advWA->StatusSelfIntersection(status))
2029       di << "\n Self-intersection found in the edge";
2030
2031     num = points2d.Length();
2032     di << "\n No. of self-intersecting edges : " << num;
2033
2034     for(int i=1; i<=num; ++i) {
2035       gp_Pnt pt = points3d(i);
2036       di << "\n Intersecting pt : (" << pt.X() << ", " << pt.Y()<< ", " << pt.Z() << ")";
2037     }
2038
2039   }
2040
2041   Handle_ShapeFix_Wire sfw = new ShapeFix_Wire;
2042   sfw->Load(TopoDS::Wire(wire));
2043   sfw->SetFace(face);
2044   sfw->SetPrecision(Precision::Confusion());
2045
2046   sfw->FixReorderMode() = 1;
2047   sfw->ClosedWireMode() = 1;
2048   sfw->FixDegeneratedMode() = 1;
2049   sfw->FixConnectedMode() = 1;
2050   sfw->ModifyTopologyMode() = 1;
2051   sfw->ModifyGeometryMode() = 1;
2052   sfw->FixSelfIntersectionMode() = 1;
2053   sfw->FixSelfIntersectingEdgeMode() = 1;
2054   sfw->FixIntersectingEdgesMode() = 1;
2055   sfw->FixNonAdjacentIntersectingEdgesMode() = 1;
2056   sfw->FixEdgeCurvesMode() = 1;
2057
2058   sfw->ModifyRemoveLoopMode() = 1;
2059   sfw->SetMaxTolerance(1.0);
2060
2061   di << "\n FixReorder = " << (Standard_Integer) sfw->FixReorder();
2062   di << "\n FixDegenerated = " << (Standard_Integer) sfw->FixDegenerated();
2063   di << "\n FixConnected = " << (Standard_Integer) sfw->FixConnected();
2064   di << "\n";
2065   di << "\n FixSelfIntersection = " << (Standard_Integer) sfw->FixSelfIntersection();
2066   di << "\n";
2067
2068   Draw::Set(argv[6],sfw->FixReorder()) ;
2069   Draw::Set(argv[7],sfw->FixDegenerated()) ;
2070   Draw::Set(argv[8],sfw->FixConnected()) ;
2071   Draw::Set(argv[9],sfw->FixSelfIntersection()) ;
2072
2073   ShapeExtend_Status status = ShapeExtend_OK;
2074   if(sfw->StatusSelfIntersection(status))
2075     di << "\n ShapeExtend_OK : No intersection found";
2076
2077   status = ShapeExtend_FAIL1;
2078   if(sfw->StatusSelfIntersection(status))
2079     di << "\n ShapeExtend_FAIL1 : analysis failed (edge has no pcurve,or no vertices etc.)";
2080
2081   status = ShapeExtend_FAIL2;
2082   if(sfw->StatusSelfIntersection(status))
2083     di << "\n ShapeExtend_FAIL2 : self-intersection was found, but not fixed because of limit of increasing tolerance (MaxTolerance)";
2084
2085   status = ShapeExtend_FAIL3;
2086   if(sfw->StatusSelfIntersection(status))
2087     di << "\n ShapeExtend_FAIL3 : intercestion of non adjacent edges found, but not fixed because of limit of increasing tolerance (MaxTolerance)";
2088
2089   status = ShapeExtend_DONE1;
2090   if(sfw->StatusSelfIntersection(status))
2091     di << "\n ShapeExtend_DONE1 : tolerance of vertex was increased to fix self-intersection";
2092
2093   status = ShapeExtend_DONE2;
2094   if(sfw->StatusSelfIntersection(status))
2095     di << "\n ShapeExtend_DONE2 : vertex was moved to fix self-intersection";
2096
2097   status = ShapeExtend_DONE3;
2098   if(sfw->StatusSelfIntersection(status))
2099     di << "\n ShapeExtend_DONE3 : some edges were removed because of intersection";
2100
2101   status = ShapeExtend_DONE4;
2102   if(sfw->StatusSelfIntersection(status))
2103     di << "\n ShapeExtend_DONE4 : pcurve(s) was(were) modified";
2104
2105   status = ShapeExtend_DONE5;
2106   if(sfw->StatusSelfIntersection(status))
2107     di << "\n ShapeExtend_DONE5 : non adjacent intersection fixed by increasing tolerance of vertex(vertices)";
2108
2109   status = ShapeExtend_DONE6;
2110   if(sfw->StatusSelfIntersection(status))
2111     di << "\n ShapeExtend_DONE6 : tolerance of edge was increased to hide intersection";
2112
2113   status = ShapeExtend_DONE7;
2114   if(sfw->StatusSelfIntersection(status))
2115     di << "\n ShapeExtend_DONE7 : range of some edges was decreased to avoid intersection";
2116
2117   TopoDS_Wire finalwire = sfw->Wire();
2118
2119   DBRep::Set(argv[1],finalwire);
2120
2121
2122   advWA->Load(TopoDS::Wire(finalwire));
2123
2124   TopoDS_Face fface =
2125     BRepBuilderAPI_MakeFace(TopoDS::Wire(finalwire),Standard_True);
2126
2127   DBRep::Set(argv[2],fface);
2128
2129   advWA->SetFace(fface);
2130
2131   TopTools_IndexedMapOfShape fM;
2132   TopExp::MapShapes(finalwire, TopAbs_EDGE, fM);
2133
2134   for(j= 1; j<=fM.Extent(); ++j) {
2135     Standard_Integer num=1;
2136     IntRes2d_SequenceOfIntersectionPoint points2d;
2137     TColgp_SequenceOfPnt points3d;
2138
2139     di << "\n j =" << j << ",  CheckSelfIntersectingEdge = " <<
2140       (Standard_Integer) advWA->CheckSelfIntersectingEdge(j, points2d, points3d);
2141
2142
2143     ShapeExtend_Status status = ShapeExtend_FAIL1;
2144     if(advWA->StatusSelfIntersection(status))
2145
2146       status = ShapeExtend_FAIL2;
2147     if(advWA->StatusSelfIntersection(status))
2148       di << "\n No Vertices found in the edge";
2149
2150     status = ShapeExtend_DONE1;
2151     if(advWA->StatusSelfIntersection(status))
2152       di << "\n Self-intersection found in the edge";
2153
2154     num = points2d.Length();
2155     di << "\n No. of self-intersecting edges : " << num;
2156
2157     for(int i=1; i<=num; ++i) {
2158       gp_Pnt pt = points3d(i);
2159       di << "\n Intersecting pt : (" << pt.X() << ", " << pt.Y() << ", " << pt.Z() << ")";
2160     }
2161
2162   }
2163
2164   di << "\n";
2165
2166   return 0;
2167 }
2168
2169 Standard_Integer OCC17480(Draw_Interpretor& di,   Standard_Integer n,   const char ** a)
2170 {
2171   if ( n < 2 ) {
2172     cout << "Usage: OCC17480 basic_shape <mode[0:1]>" <<endl;
2173     return 0;    
2174   }
2175
2176   TopoDS_Shape brep_pipe =  DBRep::Get(a[1]);
2177
2178   Standard_Boolean mode = Standard_False;
2179   if(n == 3 && atoi(a[2]) == 1)
2180     mode = Standard_True;
2181
2182   Handle_AIS_Shape ais_pipe = new AIS_Shape( brep_pipe );
2183   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
2184   Handle(AIS2D_InteractiveContext) aContext_2d = Viewer2dTest::GetAIS2DContext();
2185
2186   if(!aContext.IsNull())
2187     aContext->Display(ais_pipe);      
2188
2189   if(!aContext_2d.IsNull())
2190     aContext_2d->EraseAll(Standard_True);
2191
2192   gp_Ax3 ax3(gp_Pnt(0, 0, 0), gp_Dir(0, 0, 1), gp_Dir(0, 1, 0));
2193   gp_Pln plane(ax3);
2194   gp_Ax2 ax2 = ax3.Ax2();
2195
2196   HLRAlgo_Projector aProjector = HLRAlgo_Projector(ax2);
2197   Handle_AIS2D_ProjShape myDisplayableShape = 
2198     new AIS2D_ProjShape(aProjector, 0, mode, Standard_True);
2199
2200   myDisplayableShape->Add( brep_pipe );
2201
2202   di << "\n Total number of primitives : " << myDisplayableShape->Length();
2203
2204   Handle(Graphic2d_Line) mySectProjLines = myDisplayableShape->GetPrimitives();
2205   Handle(Graphic2d_Line) mySectProjLinesHLR = myDisplayableShape->GetHLPrimitives();
2206
2207   Handle_Standard_Type sty = mySectProjLines->DynamicType();
2208   Handle_Standard_Type sty1 = mySectProjLinesHLR->DynamicType();
2209
2210   di << "\n Class Name = " << (char*)sty->Name();
2211   di << "\n HLR Class Name = " << (char*)sty1->Name();
2212   di << "\n mySectProjLines.IsNull = " << (Standard_Integer) mySectProjLines.IsNull();
2213   di << "\n mySectProjLinesHLR.IsNull = " << (Standard_Integer) mySectProjLinesHLR.IsNull();
2214
2215   Handle(GGraphic2d_SetOfCurves) setCurves;
2216   Handle(Graphic2d_SetOfSegments) setSegments;
2217
2218   if(mySectProjLines->IsKind(STANDARD_TYPE(GGraphic2d_SetOfCurves)))
2219   {
2220     setCurves = Handle(GGraphic2d_SetOfCurves)::DownCast(mySectProjLines);
2221     di << "\n Number of Curves in set = " << setCurves->Length();
2222   }
2223   if(mySectProjLines->IsKind(STANDARD_TYPE(Graphic2d_SetOfSegments)))    
2224   {
2225     setSegments = Handle(Graphic2d_SetOfSegments)::DownCast(mySectProjLines);
2226     di << "\n Number of Curves in set = " << setSegments->Length();
2227   }
2228
2229   Handle(GGraphic2d_SetOfCurves) setCurvesHLR;
2230   Handle(Graphic2d_SetOfSegments) setSegmentsHLR;
2231
2232   if(mySectProjLinesHLR->IsKind(STANDARD_TYPE(GGraphic2d_SetOfCurves)))
2233   {
2234     setCurvesHLR = Handle(GGraphic2d_SetOfCurves)::DownCast(mySectProjLinesHLR);
2235     di << "\n HLR Number of Curves in set = " << setCurvesHLR->Length();
2236   }
2237   if(mySectProjLinesHLR->IsKind(STANDARD_TYPE(Graphic2d_SetOfSegments)))
2238   {
2239     setSegmentsHLR = Handle(Graphic2d_SetOfSegments)::DownCast(mySectProjLinesHLR);
2240     di << "\n HLR Number of Curves in set = " << setSegmentsHLR->Length();
2241   }
2242
2243   aContext_2d->Display( myDisplayableShape,Standard_True );
2244   return 0;
2245 }
2246
2247 void QABugs::Commands_17(Draw_Interpretor& theCommands) {
2248   char *group = "QABugs";
2249
2250   theCommands.Add ("BUC60842", "BUC60842", __FILE__, BUC60842, group);
2251   theCommands.Add ("BUC60843", "BUC60843 result_shape name_of_circle name_of_curve [ par1 [ tol ] ]", __FILE__, BUC60843, group);
2252   theCommands.Add ("BUC60970", "BUC60970 shape result", __FILE__, BUC60970, group);
2253   theCommands.Add ("BUC60818", "BUC60818", __FILE__, BUC60818, group);
2254   theCommands.Add ("BUC60915", "BUC60915", __FILE__, BUC60915_1, group);
2255   theCommands.Add ("OCC138", "OCC138", __FILE__, OCC138, group);
2256   theCommands.Add ("BUC60821","BUC60821",__FILE__,BUC60821,group);
2257   theCommands.Add ("OCC353","OCC353",__FILE__,OCC353,group);
2258   theCommands.Add ("OCC280","OCC280 hlr=0/1 setsurfecedetail=0/1; set perspecrive view",__FILE__,OCC280,group);
2259   theCommands.Add ("OCC232", "OCC232", __FILE__, OCC232 , group);
2260   theCommands.Add ("OCC138LC", "OCC138LC", __FILE__, OCC138LC, group);
2261   theCommands.Add ("OCC189", "OCC189", __FILE__, OCC189, group);
2262   theCommands.Add ("OCC389", "OCC389 name shape1 [shape2] ...", __FILE__, OCC389, group);
2263
2264   theCommands.Add ("OCC566", "OCC566 shape [ xmin ymin zmin xmax ymax zmax] ; print bounding box", __FILE__, OCC566, group);
2265   theCommands.Add ("OCC570", "OCC570 result", __FILE__, OCC570, group);
2266
2267   theCommands.Add("OCC570mkevol",
2268     "OCC570mkevol result object (then use updatevol) [R/Q/P]; mkevol",
2269     __FILE__, MKEVOL,group);
2270   theCommands.Add("OCC570updatevol",
2271     "OCC570updatevol edge u1 rad1 u2 rad2 ...; updatevol",
2272     __FILE__, UPDATEVOL,group);
2273   theCommands.Add("OCC570updatevollaw",
2274     "OCC570updatevollaw edge u1 rad1 u2 rad2 ...; updatevollaw",
2275     __FILE__, UPDATEVOL,group);
2276   theCommands.Add("OCC570buildevol",
2277     "OCC570buildevol; end of the evol fillet computation",
2278     __FILE__, BUILDEVOL,group);
2279
2280   theCommands.Add ("OCC606", "OCC606 result shape [-t]", __FILE__, OCC606, group);
2281
2282   theCommands.Add ("OCC813", "OCC813 U V", __FILE__, OCC813, group);
2283   theCommands.Add ("OCC814", "OCC814", __FILE__, OCC814, group);
2284   theCommands.Add ("OCC884", "OCC884 result shape [toler [maxtoler]]", __FILE__, OCC884, group);
2285   theCommands.Add ("OCC1188", "OCC1188 imagefile [filltype] : Load image as background", __FILE__, OCC1188, group);
2286
2287   theCommands.Add ("OCC1174_1", "OCC1174_1 shape", __FILE__, OCC1174_1, group);
2288   theCommands.Add ("OCC1174_2", "OCC1174_2 shape", __FILE__, OCC1174_2, group);
2289
2290   theCommands.Add("OCC672_Length"," OCC672_Length x1 y1 x2 y2 str scale_text length_dist angle_arrow length_arrow x3 y3 x4 y4 textAngle textPosH textPosV",__FILE__,OCC672_Length,group);
2291   theCommands.Add("OCC672_Angle","OCC672_Angle x1 y1 x2 y2 x3 y3 aRadius aText aTxtScale anArrAngle anArrLength txtAngle txtPosH txtPosV",__FILE__,OCC672_Angle,group);
2292   theCommands.Add("OCC672_Diameter","OCC672_Diameter x1 y1 cx cy radius aText aTxtScale anArrAngle anArrLength txtAngle txtPosH txtPosV",__FILE__,OCC672_Diameter,group);
2293   theCommands.Add("OCC672_ShortLength","",__FILE__,OCC672_ShortLength,group);
2294
2295   theCommands.Add ("OCCN1", "OCCN1 angle fuse(1 for boss / 0 for slot) length", __FILE__, OCCN1, group);
2296   theCommands.Add ("OCCN2", "OCCN2", __FILE__, OCCN2, group);
2297
2298   theCommands.Add ("OCC2569", "OCC2569 nbpoles result", __FILE__, OCC2569, group);
2299
2300   theCommands.Add ("OCC1642", "OCC1642 FinalWare FinalFace InitWare InitFace shape FixReorder FixDegenerated FixConnected FixSelfIntersection", __FILE__, OCC1642, group);
2301
2302   theCommands.Add ("OCC17480", "OCC17480 Shape", __FILE__, OCC17480, group);
2303
2304   return;
2305 }