0032155: Tests - remove broken command OCC280
[occt.git] / src / QABugs / QABugs_17.cxx
1 // Created on: 2002-03-18
2 // Created by: QA Admin
3 // Copyright (c) 2002-2014 OPEN CASCADE SAS
4 //
5 // This file is part of Open CASCADE Technology software library.
6 //
7 // This library is free software; you can redistribute it and/or modify it under
8 // the terms of the GNU Lesser General Public License version 2.1 as published
9 // by the Free Software Foundation, with special exception defined in the file
10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
11 // distribution for complete text of the license and disclaimer of any warranty.
12 //
13 // Alternatively, this file may be used under the terms of Open CASCADE
14 // commercial license or contractual agreement.
15
16 #include <stdio.h>
17
18 #include <QABugs.hxx>
19
20 #include <Draw.hxx>
21 #include <Draw_Interpretor.hxx>
22 #include <DBRep.hxx>
23 #include <DrawTrSurf.hxx>
24 #include <AIS_InteractiveContext.hxx>
25 #include <ViewerTest.hxx>
26 #include <ViewerTest_EventManager.hxx>
27 #include <AIS_Shape.hxx>
28 #include <TopoDS_Shape.hxx>
29
30 #include <Geom_Circle.hxx>
31 #include <Geom_Ellipse.hxx>
32 #include <Geom_Plane.hxx>
33 #include <Geom_BSplineSurface.hxx>
34 #include <gp_Pln.hxx>
35 #include <Geom2d_Curve.hxx>
36 #include <GeomAPI.hxx>
37 #include <Geom2dAdaptor_Curve.hxx>
38 #include <Geom2dGcc_QualifiedCurve.hxx>
39 #include <Geom2dGcc_Lin2d2Tan.hxx>
40 #include <Geom2d_Line.hxx>
41 #include <BRepBuilderAPI_MakeEdge.hxx>
42 #include <TopoDS_Edge.hxx>
43 #include <Precision.hxx>
44 #include <Geom2d_Circle.hxx>
45 #include <Geom2dGcc_QCurve.hxx>
46 #include <GccEnt_QualifiedCirc.hxx>
47 #include <Geom2dGcc_Lin2d2TanIter.hxx>
48 #include <BRepBuilderAPI_MakeWire.hxx>
49 #include <TopExp_Explorer.hxx>
50 #include <TopoDS.hxx>
51 #include <TopoDS_Wire.hxx>
52 #include <BRep_Tool.hxx>
53 #include <gp_Circ.hxx>
54 #include <BRepOffsetAPI_MakePipeShell.hxx>
55 #include <AIS_Trihedron.hxx>
56 #include <Geom_Axis2Placement.hxx>
57 #include <V3d_View.hxx>
58 #include <BRepBuilderAPI_MakeVertex.hxx>
59 #include <PrsDim_LengthDimension.hxx>
60 #include <PrsDim_RadiusDimension.hxx>
61 #include <BRepPrimAPI_MakeBox.hxx>
62 #include <GeomAdaptor_Surface.hxx>
63 #include <Extrema_ExtPS.hxx>
64 #include <Geom_CartesianPoint.hxx>
65 #include <Geom2d_CartesianPoint.hxx>
66 #include <Geom2dGcc_Circ2d3Tan.hxx>
67 #include <Aspect_Window.hxx>
68
69 #include <TopoDS_Compound.hxx>
70 #include <TopoDS_CompSolid.hxx>
71 #include <HLRAlgo_Projector.hxx>
72 #include <Standard_ErrorHandler.hxx>
73 #include <Font_NameOfFont.hxx>
74
75 static Standard_Integer BUC60842 (Draw_Interpretor& di, Standard_Integer /*argc*/,const char ** /*argv*/)
76 {
77   Standard_Character abuf[16];
78   //Standard_Character * abuf = new Standard_Character[16];
79
80   Standard_CString st = abuf;
81
82   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);
83   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);
84   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))));
85
86   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
87
88   gp_Pln pln = plne->Pln();
89   Handle(Geom2d_Curve) curve2d = GeomAPI::To2d (ell, pln);
90   Sprintf(abuf,"ell");
91   DrawTrSurf::Set(st,curve2d);
92   if(!aContext.IsNull()) {
93     Handle(AIS_Shape) aisp = 
94       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(curve2d, pln)).Edge());
95     aContext->Display (aisp, Standard_False);
96   }
97
98   Handle(Geom2d_Curve) fromcurve2d = GeomAPI::To2d (cir, pln);
99   Sprintf(abuf,"cil");
100   DrawTrSurf::Set(st,fromcurve2d);
101   if(!aContext.IsNull()) {
102     Handle(AIS_Shape) aisp = 
103       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(fromcurve2d, pln)).Edge());
104     aContext->Display (aisp, Standard_False);
105   }
106
107   Geom2dAdaptor_Curve acur (curve2d), afromcur (fromcurve2d);
108   Geom2dGcc_QualifiedCurve qcur (acur, GccEnt_outside);
109   Geom2dGcc_QualifiedCurve qfromcur (afromcur, GccEnt_outside);
110   Geom2dGcc_Lin2d2Tan lintan (qcur, qfromcur, 0.1, 0.0, 0.0);
111   Standard_Integer i=0;
112   for(i=0;i<lintan.NbSolutions();i++) {
113     Sprintf(abuf,"lintan_%d",i);
114     Handle(Geom2d_Line) glin = new Geom2d_Line(lintan.ThisSolution(i));
115     DrawTrSurf::Set(st,glin);
116     if(!aContext.IsNull()) {
117       Handle(AIS_Shape) aisp = 
118         new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(glin, pln)).Edge());
119       aContext->Display (aisp, Standard_False);
120     }
121   }
122   aContext->UpdateCurrentViewer();
123   di << " Is Done = \n" << (Standard_Integer) lintan.IsDone();
124   return 0;
125 }
126
127 static Standard_Integer BUC60843 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
128 {
129   if (argc != 4 && argc != 5 && argc != 6 && argc != 7) {
130     di << "Usage : " << argv[0] << "result_shape name_of_circle name_of_curve [ par1 [ tol ] ]\n";
131     return 1;
132   }
133
134   Standard_Real par1 = 0.0, par2 = 0.0;
135   Standard_Real tol  = Precision::Angular();
136   if (argc >= 5)
137     par1 = Draw::Atof(argv[4]);
138   if (argc == 6)
139     par2 = Draw::Atof(argv[5]);
140   if (argc == 7)
141     tol = Draw::Atof(argv[6]);
142   Handle(Geom2d_Curve)  aCur2d1 = DrawTrSurf::GetCurve2d(argv[2]);
143   Handle(Geom2d_Curve)  aCur2d2 = DrawTrSurf::GetCurve2d(argv[3]);
144   if (aCur2d1.IsNull() || aCur2d2.IsNull()) {
145     di << "Faulty : not curves\n";
146     return 1;
147   }
148   Standard_Boolean c1IsCircle = Standard_False;
149   Handle(Geom2d_Circle) aCir2d = Handle(Geom2d_Circle)::DownCast(aCur2d1);
150   if (!aCir2d.IsNull()) {
151     c1IsCircle = Standard_True;
152     if (argc == 6) tol = Draw::Atof(argv[5]);
153   }
154   if(c1IsCircle) {
155     Geom2dAdaptor_Curve acur(aCur2d2);
156     Geom2dGcc_QCurve qcur(acur, GccEnt_unqualified);
157     GccEnt_QualifiedCirc qfromcur(aCir2d->Circ2d(), GccEnt_unqualified);
158     Geom2dGcc_Lin2d2TanIter lintan(qfromcur, qcur , par1, tol); 
159     if (lintan.IsDone()) {
160       gp_Lin2d lin = lintan.ThisSolution(); 
161       Handle(Geom2d_Line) glin = new Geom2d_Line(lin); 
162       DrawTrSurf::Set(argv[1], glin);
163     } else {
164       di << "Faulty : tangent line is not done.\n";
165     }
166   }
167   else {
168     Geom2dAdaptor_Curve acur1(aCur2d1);
169     Geom2dAdaptor_Curve acur2(aCur2d2);
170     Geom2dGcc_QCurve qcur1(acur1, GccEnt_unqualified);
171     Geom2dGcc_QCurve qcur2(acur2, GccEnt_unqualified);
172     Geom2dGcc_Lin2d2TanIter lintan(qcur1, qcur2 , par1, par2, tol); 
173     if (lintan.IsDone()) {
174       gp_Lin2d lin = lintan.ThisSolution(); 
175       Handle(Geom2d_Line) glin = new Geom2d_Line(lin); 
176       DrawTrSurf::Set(argv[1], glin);
177     } else {
178       di << "Faulty : tangent line is not done.\n";
179     }
180   }   
181   return 0;
182 }
183
184 static Standard_Integer BUC60970 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
185 {
186   if(argc!=3)
187   {
188     di << "Usage : " << argv[0] << " shape result\n";
189     return 1;
190   }
191
192   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
193   if(aContext.IsNull()) {
194     di << "use 'vinit' command before " << argv[0] << "\n";
195     return 1;
196   }
197
198   TopoDS_Shape aResult = DBRep::Get(argv[1]);
199
200   BRepBuilderAPI_MakeWire bRepSpineWire; 
201   TopExp_Explorer exp_(aResult, TopAbs_WIRE); 
202
203   Standard_Integer i = 0; 
204   for(;exp_.More();exp_.Next()){     
205     i++; 
206     bRepSpineWire.Add(TopoDS::Wire(exp_.Current())); 
207   } 
208
209   //printf("total no of wires are ............. %d\n", i); 
210   di << "total no of wires are ............. " << i << "\n"; 
211
212   TopoDS_Wire spineWire = bRepSpineWire.Wire(); 
213   aContext->Display (new AIS_Shape(spineWire), Standard_False);
214
215   DBRep::Set("slineW",spineWire);
216
217   TopExp_Explorer spineWireExp(spineWire, TopAbs_EDGE); 
218   Standard_Real first, last; 
219   Handle(Geom_Curve) curl_ = BRep_Tool::Curve(TopoDS::Edge(spineWireExp.Current()), first, last); 
220   gp_Pnt firstPnt; 
221   gp_Vec tanVec; 
222   curl_->D1(first, firstPnt, tanVec); 
223   tanVec.Normalize(); 
224   gp_Dir tanDir(tanVec.X(), tanVec.Y(), tanVec.Z()); 
225   gp_Ax2 gpAx2(firstPnt, tanDir); 
226   gp_Circ gpCirc(gpAx2, 2.5); 
227   BRepBuilderAPI_MakeWire aMWire(BRepBuilderAPI_MakeEdge(new Geom_Circle(gpCirc)).Edge());
228   TopoDS_Wire topoWire(aMWire); 
229   aContext->Display (new AIS_Shape(topoWire), Standard_False);
230
231   DBRep::Set("topoW",topoWire);
232
233   BRepOffsetAPI_MakePipeShell bRepPipe(spineWire);
234   bRepPipe.Add(topoWire); 
235   bRepPipe.Build(); 
236
237   aContext->Display (new AIS_Shape(bRepPipe.Shape()), Standard_True);
238
239   DBRep::Set(argv[2],bRepPipe.Shape());
240
241   return 0;
242 }
243
244 static Standard_Integer BUC60915_1(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
245 {
246   if (argc > 1) {
247     di<<"Function don't has parameters\n";
248     return 1;
249   }
250
251   if(ViewerTest::GetAISContext().IsNull())
252   {
253     di << "View was not created. Use 'vinit' command before "
254        << argv[0] << "\n";
255     return 1;
256   }
257   Handle(AIS_InteractiveContext) context = ViewerTest::GetAISContext();
258
259   //The following dimesion code has problems regarding arrow_size. The desired effect is not produced.
260   /***************************************/
261   //First view
262   /***************************************/
263   gp_Pnt p1 = gp_Pnt(602.51,50.,0.);
264   gp_Pnt p2 = gp_Pnt(602.51,200.,0.);
265   gp_Pnt p3 = gp_Pnt(102.51,200.,0.);
266   gp_Pnt p4 = gp_Pnt(102.51,170.,0.);
267   gp_Pnt p5 = gp_Pnt(502.51,170.,0.);
268   gp_Pnt p6 = gp_Pnt(502.51,80.,0.);
269   gp_Pnt p7 = gp_Pnt(102.51,80.,0.);
270   gp_Pnt p8 = gp_Pnt(102.51,50.,0.);
271   TopoDS_Vertex V1 = BRepBuilderAPI_MakeVertex(p1);
272   TopoDS_Vertex V2 = BRepBuilderAPI_MakeVertex(p2);
273   TopoDS_Vertex V3 = BRepBuilderAPI_MakeVertex(p3);
274   TopoDS_Vertex V4 = BRepBuilderAPI_MakeVertex(p4);
275   TopoDS_Vertex V5 = BRepBuilderAPI_MakeVertex(p5);
276   TopoDS_Vertex V6 = BRepBuilderAPI_MakeVertex(p6);
277   TopoDS_Vertex V7 = BRepBuilderAPI_MakeVertex(p7);
278   TopoDS_Vertex V8 = BRepBuilderAPI_MakeVertex(p8);
279   gp_Pnt plnpt(0, 0, 0);
280   gp_Dir plndir(0, 0, 1);
281   Handle(Geom_Plane) pln = new Geom_Plane(plnpt,plndir);
282   Handle(Prs3d_DimensionAspect) anAspect = new Prs3d_DimensionAspect();
283   anAspect->MakeArrows3d (Standard_True);
284   anAspect->ArrowAspect()->SetAngle (M_PI/180.0 * 10.0);
285   /***************************************/
286   //dimension "L 502.51"
287   /***************************************/
288   Handle(PrsDim_LengthDimension) len = new PrsDim_LengthDimension(V2, V3, pln->Pln());
289   anAspect->ArrowAspect()->SetLength (30.0);
290   len->SetDimensionAspect (anAspect);
291   context->Display (len, Standard_False);
292   /***************************************/
293   //dimension "L 90"
294   /***************************************/
295   Handle(PrsDim_LengthDimension) len1 = new PrsDim_LengthDimension(V7, V4, pln->Pln());
296   len1->SetDimensionAspect (anAspect);
297   len1->SetFlyout (30.0);
298   anAspect->ArrowAspect()->SetLength (100.0);
299   context->Display (len1, Standard_False);
300   /***************************************/
301   //dimension "L 150"
302   /***************************************/
303   Handle(PrsDim_LengthDimension) len2 = new PrsDim_LengthDimension(V1, V2, pln->Pln());
304   len2->SetDimensionAspect (anAspect);
305   context->Display (len2, Standard_False);
306   /***************************************/
307   //dimension "R 88.58"
308   /***************************************/
309   gp_Circ cir = gp_Circ(gp_Ax2(gp_Pnt(191.09, -88.58, 0), gp_Dir(0, 0, 1)), 88.58);
310   TopoDS_Edge E1 = BRepBuilderAPI_MakeEdge(cir,gp_Pnt(191.09,0,0.),gp_Pnt(191.09,-177.16,0.) );
311   Handle(PrsDim_RadiusDimension) dim1 = new PrsDim_RadiusDimension(E1);
312   dim1->SetDimensionAspect (anAspect);
313   context->Display (dim1, Standard_False);
314   /***************************************/
315   //dimension "R 43.80"
316   /***************************************/
317   gp_Circ cir1 = gp_Circ(gp_Ax2(gp_Pnt(191.09, -88.58, 0), gp_Dir(0, 0, 1)), 43.80);
318   TopoDS_Edge E_cir1 = BRepBuilderAPI_MakeEdge(cir1);
319   dim1 = new PrsDim_RadiusDimension(E_cir1);
320   anAspect->ArrowAspect()->SetLength (60.0);
321   dim1->SetDimensionAspect (anAspect);
322   context->Display (dim1, Standard_False);
323   /***************************************/
324   //dimension "R 17.86"
325   /***************************************/
326   gp_Circ cir2 = gp_Circ(gp_Ax2(gp_Pnt(566.11, -88.58, 0), gp_Dir(0, 0, -1)), 17.86);
327   TopoDS_Edge E_cir2 = BRepBuilderAPI_MakeEdge(cir2);
328   dim1 = new PrsDim_RadiusDimension(E_cir2);
329   anAspect->ArrowAspect()->SetLength (40.0);
330   dim1->SetDimensionAspect (anAspect);
331   context->Display (dim1, Standard_True);
332
333   return 0;
334 }
335
336 static Standard_Integer  OCC138 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
337 {
338   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
339   if(aContext.IsNull()) { 
340     di << "use 'vinit' command before " << argv[0] << "\n";
341     return 1;
342   }
343
344   BRepPrimAPI_MakeBox box1(gp_Pnt(0, 0, 0),  gp_Pnt(100, 100, 100));
345   BRepPrimAPI_MakeBox box2(gp_Pnt(120, 120, 120),  gp_Pnt(300, 300,300));
346   BRepPrimAPI_MakeBox box3(gp_Pnt(320, 320, 320),  gp_Pnt(500, 500,500));
347
348   Handle(AIS_InteractiveObject) ais1 = new AIS_Shape(box1.Shape());
349   Handle(AIS_InteractiveObject) ais2 = new AIS_Shape(box2.Shape());
350   Handle(AIS_InteractiveObject) ais3 = new AIS_Shape(box3.Shape());
351
352   aContext->Display (ais1, Standard_False);
353   aContext->Display (ais2, Standard_False);
354   aContext->Display (ais3, Standard_False);
355
356   aContext->AddOrRemoveSelected (ais1, Standard_False);
357   aContext->AddOrRemoveSelected (ais2, Standard_False);
358   aContext->AddOrRemoveSelected (ais3, Standard_False);
359
360   di << "\n No of currents = " << aContext->NbSelected();
361
362   aContext->InitSelected();
363
364   int count = 1;
365   while(aContext->MoreSelected())
366   {
367     di << "\n count is = " << count++;
368     Handle(AIS_InteractiveObject) ais = aContext->SelectedInteractive();
369     aContext->AddOrRemoveSelected (ais, Standard_False);
370     aContext->InitSelected();
371   }
372
373   aContext->UpdateCurrentViewer();
374
375   return 0; 
376 }
377
378 static int geom_get_2Dpt_from_3Dpt(const gp_Pnt& pnt3d, const gp_Pln& pln, gp_Pnt2d& pnt2d)
379
380   int ret = 0; 
381   Handle(Geom_Plane) gpln = new Geom_Plane(pln); 
382   GeomAdaptor_Surface adsur(gpln); 
383   Extrema_ExtPS extps(pnt3d, adsur, 0.001, 0.001); 
384   if( extps.IsDone() ) 
385   { 
386     int index;
387     for(index=1; index<= extps.NbExt(); index++ ) 
388     { 
389       Extrema_POnSurf pons = extps.Point(index); 
390       Standard_Real U, V; 
391       pons.Parameter(U, V); 
392       pnt2d.SetCoord(U, V); 
393       ret = 1; 
394     } 
395   } 
396   return ret;
397
398
399 static Standard_Integer OCC353 (Draw_Interpretor& di, Standard_Integer , const char ** )
400 {
401   gp_Ax2 ax21( gp_Pnt(100,0,0), gp_Dir(0,0,1) );
402   Handle(Geom_Circle) h_cir1 = new Geom_Circle( ax21, 25 );
403
404   gp_Ax2 ax22( gp_Pnt(-100,0,0), gp_Dir(0,0,1) );
405   Handle(Geom_Circle) h_cir2 = new Geom_Circle( ax22, 25 );
406
407   gp_Pln refpln( gp_Pnt(0,0,0), gp_Dir(0,0,1) );
408   Handle(Geom2d_Curve) cir2d1 = GeomAPI::To2d(h_cir1, refpln);
409   Handle(Geom2d_Curve) cir2d2 = GeomAPI::To2d(h_cir2, refpln);
410
411   Geom2dAdaptor_Curve adop1(cir2d1);
412   Geom2dAdaptor_Curve adop2(cir2d2);
413
414   Geom2dGcc_QualifiedCurve qcur1(adop1, GccEnt_enclosing);
415   Geom2dGcc_QualifiedCurve qcur2(adop2, GccEnt_enclosing);
416
417   Handle(Geom_CartesianPoint) h_carpt = new Geom_CartesianPoint(0,175,0);
418
419   gp_Pnt pt3d = h_carpt->Pnt();
420   gp_Pnt2d pt2d;
421   geom_get_2Dpt_from_3Dpt(pt3d, refpln, pt2d);
422
423   Standard_CString st;
424   st = "cir2d1";
425   DrawTrSurf::Set(st,cir2d1);
426   st = "cir2d2";
427   DrawTrSurf::Set(st,cir2d2);
428   st = "pt2d";
429   DrawTrSurf::Set(st,pt2d);
430
431   Handle(Geom2d_CartesianPoint) pt = new Geom2d_CartesianPoint(pt2d);
432   Geom2dGcc_Circ2d3Tan sol( qcur1, qcur2, pt, 0.001, 0.0, 0.0);
433
434   int res = 0;
435   Standard_Character buf[10];
436   buf[0] = '\0';
437
438   if( sol.IsDone() ) {
439     res = sol.NbSolutions();
440     for(Standard_Integer i=1;i<=res;i++) {
441       Handle(Geom2d_Circle) aC = new Geom2d_Circle(sol.ThisSolution(i));
442       Sprintf(buf,"Result_%d",i);
443       st = buf;
444       DrawTrSurf::Set(st,aC);
445     }
446   }
447   else
448     di << "\n Faulty: no solutions found ";
449   //printf("\n Faulty: no solutions found ");
450
451   //printf("\n Num of solutions are %d ", res );
452   di << "\n Num of solutions are " << res << " ";
453
454   return 0;
455 }
456
457 static Standard_Integer  OCC138LC (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** argv)
458 {
459   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
460   if(aContext.IsNull()) { 
461     di << "use 'vinit' command before " << argv[0] << "\n";
462     return 1;
463   }
464
465   BRepPrimAPI_MakeBox box1(gp_Pnt(0, 0, 0),  gp_Pnt(100, 100, 100));
466   BRepPrimAPI_MakeBox box2(gp_Pnt(120, 120, 120),  gp_Pnt(300, 300,300));
467   BRepPrimAPI_MakeBox box3(gp_Pnt(320, 320, 320),  gp_Pnt(500, 500,500));
468
469   Handle(AIS_InteractiveObject) ais1 = new AIS_Shape(box1.Shape());
470   Handle(AIS_InteractiveObject) ais2 = new AIS_Shape(box2.Shape());
471   Handle(AIS_InteractiveObject) ais3 = new AIS_Shape(box3.Shape());
472
473   aContext->Display (ais1, Standard_False);
474   aContext->Display (ais2, Standard_False);
475   aContext->Display (ais3, Standard_False);
476
477   aContext->AddOrRemoveSelected (ais1, Standard_False);
478   aContext->AddOrRemoveSelected (ais2, Standard_False);
479   aContext->AddOrRemoveSelected (ais3, Standard_False);
480
481   di << "\n No of selected = " << aContext->NbSelected();
482
483   aContext->InitSelected();
484
485   int count = 1;
486   while(aContext->MoreSelected())
487   {
488     di << "\n count is = %d" << count++;
489     Handle(AIS_InteractiveObject) ais = aContext->SelectedInteractive();
490     aContext->AddOrRemoveSelected (ais, Standard_False);
491     aContext->InitSelected();
492   }
493
494   aContext->UpdateCurrentViewer();
495
496   return 0; 
497 }
498
499 #include <BRepBndLib.hxx>
500 #include <Draw.hxx>
501 //=======================================================================
502 //function : OCC566
503 //purpose  : 
504 //=======================================================================
505 static Standard_Integer OCC566(Draw_Interpretor& di,Standard_Integer n,const char ** a)
506 {
507   if (n < 2) {
508     di<<"Usage: " << a[0] << " shape [ xmin ymin zmin xmax ymax zmax]\n";
509     return 1;
510   }
511   TopoDS_Shape S = DBRep::Get(a[1]);
512   if (S.IsNull()) return 1;
513   Bnd_Box B;
514   BRepBndLib::AddClose(S,B);
515   Standard_Real axmin,aymin,azmin,axmax,aymax,azmax;
516   B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
517   di << axmin<<" "<< aymin<<" "<< azmin<<" "<< axmax<<" "<< aymax<<" "<< azmax;
518   if (n >= 8) {
519     Draw::Set(a[2],axmin) ;
520     Draw::Set(a[3],aymin) ;
521     Draw::Set(a[4],azmin) ;
522     Draw::Set(a[5],axmax) ;
523     Draw::Set(a[6],aymax) ;
524     Draw::Set(a[7],azmax) ;
525   }
526   return 0;
527 }
528
529 #include <BRepFilletAPI_MakeFillet.hxx>
530 //=======================================================================
531 //function : OCC570
532 //purpose  : 
533 //=======================================================================
534 static Standard_Integer OCC570 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
535 {
536   if (argc < 2) {
537     di<<"Usage: " << argv[0] << " result\n";
538     return 1;
539   }
540
541   BRepPrimAPI_MakeBox mkBox(100.,100.,100.);
542   TopoDS_Shape aBox = mkBox.Shape();
543
544   TopExp_Explorer aExp;
545   aExp.Init(aBox,TopAbs_WIRE);
546   if (aExp.More())
547   {
548     TopoDS_Shape aWire = aExp.Current();
549
550     aExp.Init(aWire,TopAbs_EDGE);
551     TopoDS_Edge e1 = TopoDS::Edge(aExp.Current()); aExp.Next();
552     TopoDS_Edge e2 = TopoDS::Edge(aExp.Current()); aExp.Next();
553     TopoDS_Edge e3 = TopoDS::Edge(aExp.Current()); aExp.Next();
554     TopoDS_Edge e4 = TopoDS::Edge(aExp.Current());
555
556     try
557     {
558       OCC_CATCH_SIGNALS
559         BRepFilletAPI_MakeFillet mkFillet(aBox);
560       mkFillet.SetContinuity(GeomAbs_C1,.001);
561
562       // Setup variable fillet data
563       TColgp_Array1OfPnt2d t_pnt(1,4);
564       t_pnt.SetValue(1,gp_Pnt2d(0.0,5.0));
565       t_pnt.SetValue(2,gp_Pnt2d(0.3,15.0));
566       t_pnt.SetValue(3,gp_Pnt2d(0.7,15.0));
567       t_pnt.SetValue(4,gp_Pnt2d(1.0,5.0));
568
569       // HERE:
570       // It is impossible to build fillet if at least one edge
571       // with variable radius is added!!! If all are constant, everything is ok.
572       mkFillet.Add(t_pnt,e1);
573       mkFillet.Add(5.0,e2);
574       mkFillet.Add(t_pnt,e3);
575       mkFillet.Add(5.0,e4);
576
577       mkFillet.Build();
578       TopoDS_Shape aFinalShape = mkFillet.Shape();
579
580       DBRep::Set(argv[1],aFinalShape);
581     }
582     catch(Standard_Failure const&)
583     {
584       di << argv[0] << ": Exception in fillet\n";
585       return 2;
586     }
587   }
588
589   return 0;
590 }
591
592 #include <Law_Interpol.hxx>
593
594 static Standard_Real t3d = 1.e-4;
595 static Standard_Real t2d = 1.e-5;
596 static Standard_Real ta  = 1.e-2;
597 static Standard_Real fl  = 1.e-3;
598 static Standard_Real tapp_angle = 1.e-2;
599 static GeomAbs_Shape blend_cont = GeomAbs_C1;
600
601 static BRepFilletAPI_MakeFillet* Rake = 0;
602 static char name[100];
603
604 static void printtolblend(Draw_Interpretor& di)
605 {
606   di<<"tolerance ang : "<<ta<<"\n";
607   di<<"tolerance 3d  : "<<t3d<<"\n";
608   di<<"tolerance 2d  : "<<t2d<<"\n";
609   di<<"fleche        : "<<fl<<"\n";
610
611   di<<"tolblend "<<ta<<" "<<t3d<<" "<<t2d<<" "<<fl<<"\n";
612 }
613
614 static Standard_Integer MKEVOL(Draw_Interpretor& di, 
615                                Standard_Integer narg, 
616                                const char ** a)
617 {
618   if(Rake != 0) {delete Rake; Rake = 0;}
619   printtolblend(di);
620   if (narg < 3) return 1;
621   TopoDS_Shape V = DBRep::Get(a[2]);
622   Rake = new BRepFilletAPI_MakeFillet(V);
623   Rake->SetParams(ta,t3d,t2d,t3d,t2d,fl);
624   Rake->SetContinuity(blend_cont, tapp_angle);
625   if (narg == 4) {
626     ChFi3d_FilletShape FSh = ChFi3d_Rational;
627     if (!strcasecmp(a[3], "Q")) {
628       FSh = ChFi3d_QuasiAngular;
629     }
630     else if (!strcasecmp(a[3], "P")) {
631       FSh = ChFi3d_Polynomial;
632     }
633     Rake->SetFilletShape(FSh);
634   }
635   strcpy(name, a[1]);
636   return 0;
637 }
638
639 static Standard_Integer UPDATEVOL(Draw_Interpretor& di, 
640                                   Standard_Integer narg, 
641                                   const char ** a)
642 {
643   if(Rake == 0){
644     di << "MakeFillet non initialise\n";
645     return 1 ;
646   }
647   if(narg%2 != 0 || narg < 4) return 1;
648   TColgp_Array1OfPnt2d uandr(1,(narg/2)-1);
649   Standard_Real Rad, Par;
650   TopoDS_Shape aLocalEdge(DBRep::Get(a[1],TopAbs_EDGE));
651   TopoDS_Edge E = TopoDS::Edge(aLocalEdge);
652   for (Standard_Integer ii = 1; ii <= (narg/2)-1; ii++){
653     Par = Draw::Atof(a[2*ii]);
654     Rad = Draw::Atof(a[2*ii + 1]);
655     uandr.ChangeValue(ii).SetCoord(Par,Rad);
656   }
657   //HELPDESK: Add law creation
658   const Standard_Boolean aLawMode = !strcmp(a[0],"OCC570evollaw");
659   if (aLawMode)
660   {
661     di<<"INFO: Using law to define variable fillet\n";
662     Handle(Law_Interpol) law = new Law_Interpol;
663     law->Set(uandr);
664     Rake->Add(law,E);
665   }
666   else
667     Rake->Add(uandr,E);
668
669   return 0;
670 }
671
672 static Standard_Integer BUILDEVOL(Draw_Interpretor& di,
673                                   Standard_Integer, 
674                                   const char **)
675 {
676   if(Rake == 0){
677     di << "MakeFillet non initialise\n";
678     return 1 ;
679   }
680   Rake->Build();
681   if(Rake->IsDone()){
682     TopoDS_Shape result = Rake->Shape();
683     DBRep::Set(name,result);
684     if(Rake != 0) {delete Rake; Rake = 0;}
685     return 0;
686   }
687   if(Rake != 0) {delete Rake; Rake = 0;}
688   return 1;
689 }
690
691 #include <TColGeom_SequenceOfCurve.hxx>
692 #include <Geom_TrimmedCurve.hxx>
693 #include <GeomFill_NSections.hxx>
694 #include <BRepBuilderAPI_MakeFace.hxx>
695 //=======================================================================
696 //function : OCC606
697 //purpose  : 
698 //=======================================================================
699 static Standard_Integer OCC606 ( Draw_Interpretor& di, Standard_Integer n, const char ** a )
700 {
701   if (n != 3 && n != 4)
702   {
703     di << "Usage : " << a[0] << " result shape [-t]\n";
704     return 0;
705   }
706
707   Standard_Boolean TrimMode = (n == 4);
708   if (TrimMode) di<<"INFO: Using trimmed curves...\n";
709   else di<<"INFO: Using non trimmed curves...\n";
710
711   TopoDS_Shape S = DBRep::Get(a[2]);
712
713   TopExp_Explorer t_exp(S,TopAbs_EDGE);
714   TColGeom_SequenceOfCurve n_curves1;
715   TColStd_SequenceOfReal np;
716   Standard_Real param = 5.0;
717
718   for (;t_exp.More();t_exp.Next())
719   {
720     Standard_Real f,l;
721     Handle(Geom_Curve) h_cur = BRep_Tool::Curve(TopoDS::Edge(t_exp.Current()),f,l);
722     if (!h_cur.IsNull())
723     {
724       if (TrimMode) h_cur = new Geom_TrimmedCurve(h_cur,f,l);
725
726       n_curves1.Append(h_cur);
727       np.Append(param);
728       param += 1.0;
729     }
730   }
731
732   if (n_curves1.Length()>1 && !np.IsEmpty())
733   {
734     try
735     {
736       OCC_CATCH_SIGNALS
737         GeomFill_NSections b_surface1(n_curves1, np);
738       Handle(Geom_BSplineSurface) result_surf1 = b_surface1.BSplineSurface();
739       if (!result_surf1.IsNull())
740       {
741         BRepBuilderAPI_MakeFace b_face1(result_surf1, Precision::Confusion());
742         TopoDS_Face bsp_face1 = b_face1.Face();
743         DBRep::Set(a[1],bsp_face1);
744       }
745     }
746     catch(Standard_Failure const&)
747     {
748       di<<"ERROR: Exception in GeomFill_NSections\n";
749     }
750   }
751
752
753   return 0;
754 }
755
756 //=======================================================================
757 //function : OCC813
758 //purpose  : 
759 //=======================================================================
760 static Standard_Integer OCC813 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
761 {
762   if (argc < 3)
763   {
764     di << "Usage : " << argv[0] << " U V\n";
765     return 1;
766   }
767
768   Standard_CString str;
769   Standard_Real U = Draw::Atof(argv[1]);
770   Standard_Real V = Draw::Atof(argv[2]);
771
772   //Between ellipse and point:
773
774   Handle(Geom_Ellipse) ell = new Geom_Ellipse(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
775     gp_Dir(0.173648,0.984808,0.000000),
776     gp_Dir(-0.932169,0.164367,-0.322560)), 150, 100);
777   Handle(Geom_Plane) plne = new Geom_Plane(gp_Ax3(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
778     gp_Dir(0.173648,0.984808,0.000000),
779     gp_Dir(-0.932169,0.164367,-0.322560))) );
780
781   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
782
783   gp_Pnt2d pt2d(U,V);
784   gp_Pln pln = plne->Pln();
785
786   str = "OCC813_pnt"; DrawTrSurf::Set(str,pt2d);
787
788   Handle(Geom2d_Curve) curve2d = GeomAPI::To2d(ell,pln);
789   Geom2dAdaptor_Curve acur(curve2d);
790   Geom2dGcc_QualifiedCurve qcur(acur, GccEnt_outside);
791
792   str = "OCC813_ell"; DrawTrSurf::Set(str,curve2d);
793   if(!aContext.IsNull()) {
794     Handle(AIS_Shape) aisp = 
795       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(curve2d, pln)).Edge());
796     aContext->Display (aisp, Standard_False);
797   }
798
799   //This does not give any solutions.
800   Geom2dGcc_Lin2d2Tan lintan(qcur, pt2d, 0.1);
801   di<<"OCC813 nb of solutions = "<<lintan.NbSolutions()<<"\n";
802
803   Standard_Character abuf[16];
804   Standard_CString st = abuf;
805
806   Standard_Integer i;
807   for(i=1;i<=lintan.NbSolutions();i++) {
808     Sprintf(abuf,"lintan_%d",i);
809     Handle(Geom2d_Line) glin = new Geom2d_Line(lintan.ThisSolution(i));
810     DrawTrSurf::Set(st,glin);
811     if(!aContext.IsNull()) {
812       Handle(AIS_Shape) aisp = 
813         new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(glin, pln)).Edge());
814       aContext->Display (aisp, Standard_False);
815     }
816   }
817
818   if (!aContext.IsNull())
819   {
820     aContext->UpdateCurrentViewer();
821   }
822
823   return 0;
824 }
825
826 //=======================================================================
827 //function : OCC814
828 //purpose  : 
829 //=======================================================================
830 static Standard_Integer OCC814 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
831 {
832   if (argc > 1)
833   {
834     di << "Usage : " << argv[0] << "\n";
835     return 1;
836   }
837
838   Standard_CString str;
839
840   //Between Ellipse and Circle:
841
842   Handle(Geom_Circle) cir = new Geom_Circle(gp_Ax2(gp_Pnt(823.687192,502.366825,478.960440),
843     gp_Dir(0.173648,0.984808,0.000000),
844     gp_Dir(-0.932169,0.164367,-0.322560)), 50);
845   Handle(Geom_Ellipse) ell = new Geom_Ellipse(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
846     gp_Dir(0.173648,0.984808,0.000000),
847     gp_Dir(-0.932169,0.164367,-0.322560)), 150, 100);
848   Handle(Geom_Plane) plne = new Geom_Plane(gp_Ax3(gp_Ax2(gp_Pnt(1262.224429,425.040878,363.609716),
849     gp_Dir(0.173648,0.984808,0.000000),
850     gp_Dir(-0.932169,0.164367,-0.322560))) );
851
852   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
853
854   gp_Pln pln = plne->Pln();
855   Handle(Geom2d_Curve) curve2d = GeomAPI::To2d(ell, pln);
856   Handle(Geom2d_Curve) fromcurve2d = GeomAPI::To2d(cir, pln);
857
858   str = "OCC814_cir"; DrawTrSurf::Set(str,curve2d);
859   str = "OCC814_ell"; DrawTrSurf::Set(str,fromcurve2d);
860   if(!aContext.IsNull()) {
861     Handle(AIS_Shape) aisp = 
862       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(curve2d, pln)).Edge());
863     aContext->Display (aisp, Standard_False);
864   }
865   if(!aContext.IsNull()) {
866     Handle(AIS_Shape) aisp = 
867       new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(fromcurve2d, pln)).Edge());
868     aContext->Display (aisp, Standard_False);
869   }
870
871   Geom2dAdaptor_Curve acur(curve2d), afromcur(fromcurve2d);
872
873   Geom2dGcc_QualifiedCurve qcur(acur, GccEnt_outside) ;
874   Geom2dGcc_QualifiedCurve qfromcur(afromcur, GccEnt_outside) ;
875
876   //This does not give any solutions.
877   Geom2dGcc_Lin2d2Tan lintan(qcur, qfromcur, 0.1);
878   di<<"OCC814 nb of solutions = "<<lintan.NbSolutions()<<"\n";
879
880   Standard_Character abuf[16];
881   Standard_CString st = abuf;
882
883   Standard_Integer i;
884   for(i=1;i<=lintan.NbSolutions();i++) {
885     Sprintf(abuf,"lintan_%d",i);
886     Handle(Geom2d_Line) glin = new Geom2d_Line(lintan.ThisSolution(i));
887     DrawTrSurf::Set(st,glin);
888     if(!aContext.IsNull()) {
889       Handle(AIS_Shape) aisp = 
890         new AIS_Shape (BRepBuilderAPI_MakeEdge(GeomAPI::To3d(glin, pln)).Edge());
891       aContext->Display (aisp, Standard_False);
892     }
893   }
894
895   if (!aContext.IsNull())
896   {
897     aContext->UpdateCurrentViewer();
898   }
899
900   return 0;
901 }
902
903 #include <ShapeAnalysis_Wire.hxx>
904 #include <IntRes2d_SequenceOfIntersectionPoint.hxx>
905 #include <TColgp_SequenceOfPnt.hxx>
906 #include <ShapeFix_Wire.hxx>
907 //=======================================================================
908 //function : OCC884
909 //purpose  : 
910 //=======================================================================
911 static Standard_Integer OCC884 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
912 {
913   if (argc < 3 || argc > 5)
914   {
915     di << "Usage : " << argv[0] << " result shape [toler [maxtoler]]\n";
916     return 1;
917   }
918
919   TopoDS_Shape shape = DBRep::Get(argv[2]); //read the shape
920
921   TopExp_Explorer exp(shape, TopAbs_EDGE);
922   if (!exp.More())
923   {
924     di<<"Error: no edge found\n";
925     return 1;
926   }
927
928   BRep_Builder builder;
929   TopoDS_Wire wire;
930   builder.MakeWire(wire);
931   builder.Add(wire, TopoDS::Edge(exp.Current()));
932
933   // HelpDesk: Create planar face if possible
934   TopoDS_Face face = BRepBuilderAPI_MakeFace(wire,Standard_True);
935
936   Handle(ShapeAnalysis_Wire) advWA = new ShapeAnalysis_Wire;
937   advWA->Load(wire);
938   advWA->SetFace(face);
939   advWA->SetPrecision(0.001);
940
941   IntRes2d_SequenceOfIntersectionPoint points2d;
942   TColgp_SequenceOfPnt points3d;
943
944   di << "Info: CheckSelfIntersectingEdge = " << (Standard_Integer) advWA->CheckSelfIntersectingEdge(1,points2d,points3d) << "\n";
945
946   ShapeExtend_Status status = ShapeExtend_FAIL1;
947   if (advWA->StatusSelfIntersection(status))
948     di << "Info: No P Curve found in the edge\n";
949
950   status = ShapeExtend_FAIL2; 
951   if (advWA->StatusSelfIntersection(status))
952     di << "Info: No Vertices found in the edge\n";
953
954   status = ShapeExtend_DONE1; 
955   if (advWA->StatusSelfIntersection(status))
956     di << "Info: Self-intersection found in the edge\n";
957
958   Standard_Integer i, num = points2d.Length();
959   di << "Info: No. of self-intersection points : " << num << "\n";
960
961   char str[80];
962   Standard_CString aName = str;
963   for (i = 1; i <= num; ++i)
964   { 
965     gp_Pnt pt = points3d(i); 
966     di << "Info: Intersecting pt : (" << pt.X() << ", " << pt.Y() << ", " << pt.Z() << ")\n";
967     Sprintf(str,"p_%d",i);
968     DrawTrSurf::Set(aName,pt);
969   }
970
971   Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire;
972   sfw->Load(wire);
973   sfw->SetFace(face);
974
975   if (argc > 3) sfw->SetPrecision(Draw::Atof(argv[3])/*0.1*/);
976   if (argc > 4) sfw->SetMaxTolerance(Draw::Atof(argv[4]));
977   di << "Info: Precision is set to " << sfw->Precision() << "\n";
978   di << "Info: MaxTolerance is set to " << sfw->MaxTolerance() << "\n";
979
980   sfw->ModifyTopologyMode() = 1;
981   sfw->ModifyGeometryMode() = 1; 
982   sfw->FixSelfIntersectingEdgeMode() = 1;
983   //printf("Info: FixSelfIntersection = %d\n",sfw->FixSelfIntersection());
984   di << "Info: FixSelfIntersection = " << (Standard_Integer) sfw->FixSelfIntersection() << "\n";
985
986   DBRep::Set(argv[1],sfw->Wire());
987
988   status = ShapeExtend_OK;
989   if (sfw->StatusSelfIntersection(status))
990     di << "Info: ShapeExtend_OK : No intersection found\n";
991   //printf("Info: ShapeExtend_OK : No intersection found\n");
992
993   status = ShapeExtend_FAIL1;
994   if (sfw->StatusSelfIntersection(status))
995     di << "Info: ShapeExtend_FAIL1 : analysis failed (edge has no pcurve, or no vertices etc.)\n";
996
997   status = ShapeExtend_FAIL2;
998   if (sfw->StatusSelfIntersection(status))
999     di << "Info: ShapeExtend_FAIL2 : self-intersection was found, but not fixed because of limit of increasing tolerance (MaxTolerance)\n";
1000
1001   status = ShapeExtend_FAIL3;
1002   if (sfw->StatusSelfIntersection(status))
1003     di << "Info: ShapeExtend_FAIL3 : intercestion of non adjacent edges found, but not fixed because of limit of increasing tolerance (MaxTolerance)\n";
1004
1005   status = ShapeExtend_DONE1;
1006   if (sfw->StatusSelfIntersection(status))
1007     di << "Info: ShapeExtend_DONE1 : tolerance of vertex was increased to fix self-intersection\n";
1008
1009   status = ShapeExtend_DONE2;
1010   if (sfw->StatusSelfIntersection(status))
1011     di << "Info: ShapeExtend_DONE2 : vertex was moved to fix self-intersection\n";
1012
1013   status = ShapeExtend_DONE3;
1014   if (sfw->StatusSelfIntersection(status))
1015     di << "Info: ShapeExtend_DONE3 : some edges were removed because of intersection\n";
1016
1017   status = ShapeExtend_DONE4;
1018   if (sfw->StatusSelfIntersection(status))
1019     di << "Info: ShapeExtend_DONE4 : pcurve(s) was(were) modified\n";
1020
1021   status = ShapeExtend_DONE5;
1022   if (sfw->StatusSelfIntersection(status))
1023     di << "Info: ShapeExtend_DONE5 : non adjacent intersection fixed by increasing tolerance of vertex(vertices)\n";
1024
1025   status = ShapeExtend_DONE6;
1026   if (sfw->StatusSelfIntersection(status))
1027     di << "Info: ShapeExtend_DONE6 : tolerance of edge was increased to hide intersection\n";
1028
1029   status = ShapeExtend_DONE7;
1030   if (sfw->StatusSelfIntersection(status))
1031     di << "Info: ShapeExtend_DONE7 : range of some edges was decreased to avoid intersection\n";
1032
1033   return 0;
1034 }
1035
1036 #include <Graphic3d_MaterialAspect.hxx>
1037 #include <Prs3d_Drawer.hxx>
1038 #include <Prs3d_ShadingAspect.hxx>
1039 //=======================================================================
1040 //function : OCC1174_1
1041 //purpose  : 
1042 //=======================================================================
1043 static Standard_Integer OCC1174_1 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1044 {
1045   if (argc != 2)
1046   {
1047     di << "Usage : " << argv[0] << " shape\n";
1048     return 1;
1049   }
1050
1051   Handle(AIS_InteractiveContext) anAISContext = ViewerTest::GetAISContext();
1052   if(anAISContext.IsNull())
1053   {
1054     di << "use 'vinit' command before " << argv[0] << "\n";
1055     return 1;
1056   }
1057
1058   TopoDS_Shape aShape = DBRep::Get(argv[1]);
1059
1060   Handle(AIS_Shape) anAisIO = new AIS_Shape(aShape);
1061
1062   Quantity_Color aColF(0.0, 0.4, 0.0, Quantity_TOC_sRGB);
1063   Quantity_Color aColB(0.0, 0.0, 0.6, Quantity_TOC_sRGB);
1064
1065   Handle(Prs3d_Drawer) aDrawer = anAisIO->Attributes();
1066   Handle(Prs3d_ShadingAspect) aShadingAspect = aDrawer->ShadingAspect();
1067
1068   Graphic3d_MaterialAspect aFront = aShadingAspect->Material(Aspect_TOFM_FRONT_SIDE);
1069   aFront.SetAmbientColor(aColF);
1070   aFront.SetDiffuseColor(aColF);
1071   aFront.SetSpecularColor(aColF);
1072   aFront.SetEmissiveColor(Quantity_NOC_BLACK);
1073   aFront.SetTransparency(0.0);
1074   aShadingAspect->SetMaterial(aFront,Aspect_TOFM_FRONT_SIDE);
1075
1076   Graphic3d_MaterialAspect aBack = aShadingAspect->Material(Aspect_TOFM_BACK_SIDE);
1077   aBack.SetAmbientColor(aColB);
1078   aBack.SetDiffuseColor(aColB);
1079   aBack.SetSpecularColor(aColB);
1080   aBack.SetEmissiveColor(Quantity_NOC_BLACK);
1081   aBack.SetTransparency(0.0);
1082   aShadingAspect->SetMaterial(aBack,Aspect_TOFM_BACK_SIDE);
1083
1084   aDrawer->SetShadingAspect (aShadingAspect);
1085
1086   anAISContext->Display (anAisIO, 1, 0, Standard_True);
1087
1088   Standard_Real r, g, b; 
1089   aShadingAspect->Color(Aspect_TOFM_FRONT_SIDE).Values(r,g,b, Quantity_TOC_sRGB);
1090   di << "Info: color on front side (" << r << "," << g << "," << b << ")\n";
1091   aShadingAspect->Color(Aspect_TOFM_BACK_SIDE).Values(r,g,b, Quantity_TOC_sRGB);
1092   di << "Info: color on back side (" << r << "," << g << "," << b << ")\n";
1093
1094   return 0;
1095 }
1096
1097 //=======================================================================
1098 //function : OCC1174_2
1099 //purpose  : 
1100 //=======================================================================
1101 static Standard_Integer OCC1174_2 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1102 {
1103   if (argc != 2)
1104   {
1105     di << "Usage : " << argv[0] << " shape\n";
1106     return 1;
1107   }
1108
1109   Handle(AIS_InteractiveContext) AISContext = ViewerTest::GetAISContext();
1110   if(AISContext.IsNull()) 
1111   {
1112     di << "use 'vinit' command before " << argv[0] << "\n";
1113     return 1;
1114   }
1115
1116   TopoDS_Shape sh = DBRep::Get(argv[1]);
1117
1118   Handle(AIS_Shape) ais = new AIS_Shape(sh); 
1119   AISContext->Display (ais, 1, 0, Standard_False);
1120   AISContext->SetMaterial (ais, Graphic3d_NameOfMaterial_ShinyPlastified, Standard_False);
1121
1122   Quantity_Color colf(0.0, 0.4, 0.0, Quantity_TOC_sRGB);
1123   Quantity_Color colb(0.0, 0.0, 0.6, Quantity_TOC_sRGB);
1124   Handle(Prs3d_ShadingAspect) sa = ais->Attributes()->ShadingAspect(); 
1125
1126   Graphic3d_MaterialAspect front = sa->Material(Aspect_TOFM_FRONT_SIDE); 
1127   front.SetAmbientColor(colf); 
1128   front.SetDiffuseColor(colf); 
1129   front.SetSpecularColor(colf); 
1130   front.SetEmissiveColor(Quantity_NOC_BLACK);
1131   front.SetTransparency (0.4f);
1132   sa->SetMaterial(front,Aspect_TOFM_FRONT_SIDE); 
1133
1134   Graphic3d_MaterialAspect back = sa->Material(Aspect_TOFM_BACK_SIDE); 
1135   back.SetAmbientColor(colb); 
1136   back.SetDiffuseColor(colb); 
1137   back.SetSpecularColor(colb); 
1138   back.SetEmissiveColor(Quantity_NOC_BLACK);
1139   back.SetTransparency (0.2f);
1140   sa->SetMaterial(back,Aspect_TOFM_BACK_SIDE); 
1141
1142   AISContext->Redisplay (ais, 1, 0);
1143
1144   return 0;
1145 }
1146
1147 #include <TopoDS_Solid.hxx>
1148 #include <BRepFeat_MakeDPrism.hxx>
1149 //=======================================================================
1150 //function : OCCN1
1151 //purpose  : FEATURE OPERATION (BOSS and SLOT)
1152 //=======================================================================
1153 static Standard_Integer OCCN1 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1154 {
1155   if (argc != 4)
1156   {
1157     di << "Usage : " << argv[0] << " angle fuse(1 for boss / 0 for slot) length\n";
1158     return 1;
1159   }
1160
1161   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1162   if(aContext.IsNull()) 
1163   {
1164     di << "use 'vinit' command before " << argv[0] << "\n";
1165     return 1;
1166   }
1167   Standard_Real    angle  = Draw::Atof(argv[1]);
1168   Standard_Integer fuse   = Draw::Atoi(argv[2]);
1169   Standard_Real    length = Draw::Atof(argv[3]);
1170
1171   BRepBuilderAPI_MakeEdge edge1(gp_Pnt(0, 0, 0), gp_Pnt(50, 0, 0));
1172   BRepBuilderAPI_MakeEdge edge2(gp_Pnt(50, 0, 0), gp_Pnt(50, 50, 0));
1173   BRepBuilderAPI_MakeEdge edge3(gp_Pnt(50, 50, 0), gp_Pnt(0, 50, 0));
1174   BRepBuilderAPI_MakeEdge edge4(gp_Pnt(0, 50, 0), gp_Pnt(0, 0, 0));
1175   TopoDS_Edge ted1 = edge1.Edge();
1176   TopoDS_Edge ted2 = edge2.Edge();
1177   TopoDS_Edge ted3 = edge3.Edge();
1178   TopoDS_Edge ted4 = edge4.Edge();
1179
1180   BRepBuilderAPI_MakeWire wire(ted1, ted2, ted3, ted4);
1181   TopoDS_Wire twire = wire.Wire();
1182
1183   BRepBuilderAPI_MakeFace face(twire);
1184   TopoDS_Face tface = face.Face();
1185   ////////Handle(AIS_Shape) face_ais = new AIS_Shape( tface );
1186   ////////aContext->Display(face_ais);
1187
1188   TopoDS_Solid box;
1189   BRepPrimAPI_MakeBox block(gp_Pnt(-50, -50, -150), 150, 150, 100);
1190   box = block.Solid();
1191
1192   TopoDS_Face top_face;
1193   top_face = block.TopFace();
1194
1195   BRepFeat_MakeDPrism slot(box, tface, top_face, angle, fuse,
1196     Standard_True);
1197
1198   slot.Perform(length);
1199   TopoDS_Shape S2 = slot.Shape();
1200
1201   DBRep::Set("OCCN1_tface",tface);
1202   DBRep::Set("OCCN1_slot",S2);
1203
1204   return 0;
1205 }
1206
1207 #include <BRepPrimAPI_MakeCylinder.hxx>
1208 #include <BRepPrimAPI_MakeSphere.hxx>
1209 #include <BRepAlgoAPI_Section.hxx>
1210 //=======================================================================
1211 //function : OCCN2
1212 //purpose  : BOOLEAN OPERATION
1213 //=======================================================================
1214 static Standard_Integer OCCN2 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1215 {
1216   if (argc > 2) {
1217     di << "Usage : " << argv[0] << "\n";
1218     return 1;
1219   }
1220
1221   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1222   if(aContext.IsNull()) 
1223   {
1224     di << "use 'vinit' command before " << argv[0] << "\n";
1225     return 1;
1226   }
1227
1228   BRepPrimAPI_MakeCylinder cylinder(50,200);
1229   TopoDS_Shape cylinder_sh = cylinder.Shape();
1230
1231   BRepPrimAPI_MakeSphere sphere(gp_Pnt(60,0,100),50);
1232   TopoDS_Shape sphere_sh = sphere.Shape();
1233
1234   di << "BRepAlgoAPI_Section section(cylinder_sh, sphere_sh)\n";
1235   BRepAlgoAPI_Section section(cylinder_sh, sphere_sh);
1236   if (!section.IsDone()){
1237     di << "Error performing intersection: not done.\n";
1238   }
1239   const TopoDS_Shape& shape = section.Shape();
1240
1241   DBRep::Set("OCCN2_cylinder",cylinder_sh);
1242   DBRep::Set("OCCN2_sphere",sphere_sh);
1243   DBRep::Set("OCCN2_section",shape);
1244
1245   return 0;
1246 }
1247
1248 #include <TColgp_Array1OfPnt.hxx>
1249 #include <Geom_BezierCurve.hxx>
1250
1251 static Standard_Integer OCC2569 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1252 {
1253   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
1254   if(aContext.IsNull()) 
1255   {
1256     di << "use 'vinit' command before " << argv[0] << "\n";
1257     return 1;
1258   }
1259   if(argc != 3) {
1260     di << "Usage : " << argv[0] << " nbpoles result\n";
1261     return 1;
1262   }
1263
1264   int poles=Draw::Atoi(argv[1]); 
1265
1266   TColgp_Array1OfPnt arr(1, poles); 
1267   for(int i=1; i<=poles; i++) 
1268     arr.SetValue(i, gp_Pnt(i+10, i*2+20, i*3+45)); 
1269
1270   Handle(Geom_BezierCurve) bez = new Geom_BezierCurve(arr); 
1271   if(bez.IsNull()) {
1272     di << "\n The curve is not created.\n";
1273   } else {
1274     di << "\n Degree = " << bez->Degree() << "\n";   
1275   }
1276   TopoDS_Edge sh = BRepBuilderAPI_MakeEdge(bez).Edge(); 
1277   Handle(AIS_Shape) ais = new AIS_Shape(sh); 
1278   aContext->Display (ais, Standard_True);
1279   DrawTrSurf::Set(argv[2],bez);
1280   return 0;
1281 }
1282
1283 #include <TopTools_IndexedMapOfShape.hxx>
1284 #include <TopExp.hxx>
1285
1286 static Standard_Integer OCC1642 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
1287 {
1288   if (argc != 10) {
1289     di<<"Usage: " << argv[0] << " FinalWare FinalFace InitWare InitFace shape FixReorder FixDegenerated FixConnected FixSelfIntersection\n";
1290     return 1;
1291   }
1292
1293   TopoDS_Shape shape = DBRep::Get(argv[5]);
1294
1295   TopExp_Explorer exp(shape, TopAbs_WIRE);
1296   TopoDS_Shape wire = exp.Current();
1297
1298   Handle(ShapeAnalysis_Wire) advWA = new ShapeAnalysis_Wire;
1299   advWA->Load(TopoDS::Wire(wire));
1300
1301   DBRep::Set(argv[3],wire);
1302
1303   TopoDS_Face face =
1304     BRepBuilderAPI_MakeFace(TopoDS::Wire(wire),Standard_True);
1305
1306   DBRep::Set(argv[4],face);
1307
1308   advWA->SetFace(face);
1309   Standard_Real precision_to_ana = 0.0001;
1310   advWA->SetPrecision(precision_to_ana);
1311
1312   TopTools_IndexedMapOfShape M;
1313   TopExp::MapShapes(wire, TopAbs_EDGE, M);
1314
1315   int j= 1;
1316   for(j= 1; j<=M.Extent(); ++j) {
1317     Standard_Integer num=1;
1318     IntRes2d_SequenceOfIntersectionPoint points2d;
1319     TColgp_SequenceOfPnt points3d;
1320
1321     std::cout << "\n j =" << j << ",  CheckSelfIntersectingEdge = " <<
1322       advWA->CheckSelfIntersectingEdge(j, points2d, points3d);
1323
1324
1325     ShapeExtend_Status status = ShapeExtend_FAIL1;
1326     if(advWA->StatusSelfIntersection(status))
1327
1328       status = ShapeExtend_FAIL2;
1329     if(advWA->StatusSelfIntersection(status))
1330       di << "\n No Vertices found in the edge";
1331
1332     status = ShapeExtend_DONE1;
1333     if(advWA->StatusSelfIntersection(status))
1334       di << "\n Self-intersection found in the edge";
1335
1336     num = points2d.Length();
1337     di << "\n No. of self-intersecting edges : " << num;
1338
1339     for(int i=1; i<=num; ++i) {
1340       gp_Pnt pt = points3d(i);
1341       di << "\n Intersecting pt : (" << pt.X() << ", " << pt.Y()<< ", " << pt.Z() << ")";
1342     }
1343
1344   }
1345
1346   Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire;
1347   sfw->Load(TopoDS::Wire(wire));
1348   sfw->SetFace(face);
1349   sfw->SetPrecision(Precision::Confusion());
1350
1351   sfw->FixReorderMode() = 1;
1352   sfw->ClosedWireMode() = 1;
1353   sfw->FixDegeneratedMode() = 1;
1354   sfw->FixConnectedMode() = 1;
1355   sfw->ModifyTopologyMode() = 1;
1356   sfw->ModifyGeometryMode() = 1;
1357   sfw->FixSelfIntersectionMode() = 1;
1358   sfw->FixSelfIntersectingEdgeMode() = 1;
1359   sfw->FixIntersectingEdgesMode() = 1;
1360   sfw->FixNonAdjacentIntersectingEdgesMode() = 1;
1361   sfw->FixEdgeCurvesMode() = 1;
1362
1363   sfw->ModifyRemoveLoopMode() = 1;
1364   sfw->SetMaxTolerance(1.0);
1365
1366   di << "\n FixReorder = " << (Standard_Integer) sfw->FixReorder();
1367   di << "\n FixDegenerated = " << (Standard_Integer) sfw->FixDegenerated();
1368   di << "\n FixConnected = " << (Standard_Integer) sfw->FixConnected();
1369   di << "\n";
1370   di << "\n FixSelfIntersection = " << (Standard_Integer) sfw->FixSelfIntersection();
1371   di << "\n";
1372
1373   Draw::Set(argv[6],sfw->FixReorder()) ;
1374   Draw::Set(argv[7],sfw->FixDegenerated()) ;
1375   Draw::Set(argv[8],sfw->FixConnected()) ;
1376   Draw::Set(argv[9],sfw->FixSelfIntersection()) ;
1377
1378   ShapeExtend_Status status = ShapeExtend_OK;
1379   if(sfw->StatusSelfIntersection(status))
1380     di << "\n ShapeExtend_OK : No intersection found";
1381
1382   status = ShapeExtend_FAIL1;
1383   if(sfw->StatusSelfIntersection(status))
1384     di << "\n ShapeExtend_FAIL1 : analysis failed (edge has no pcurve,or no vertices etc.)";
1385
1386   status = ShapeExtend_FAIL2;
1387   if(sfw->StatusSelfIntersection(status))
1388     di << "\n ShapeExtend_FAIL2 : self-intersection was found, but not fixed because of limit of increasing tolerance (MaxTolerance)";
1389
1390   status = ShapeExtend_FAIL3;
1391   if(sfw->StatusSelfIntersection(status))
1392     di << "\n ShapeExtend_FAIL3 : intercestion of non adjacent edges found, but not fixed because of limit of increasing tolerance (MaxTolerance)";
1393
1394   status = ShapeExtend_DONE1;
1395   if(sfw->StatusSelfIntersection(status))
1396     di << "\n ShapeExtend_DONE1 : tolerance of vertex was increased to fix self-intersection";
1397
1398   status = ShapeExtend_DONE2;
1399   if(sfw->StatusSelfIntersection(status))
1400     di << "\n ShapeExtend_DONE2 : vertex was moved to fix self-intersection";
1401
1402   status = ShapeExtend_DONE3;
1403   if(sfw->StatusSelfIntersection(status))
1404     di << "\n ShapeExtend_DONE3 : some edges were removed because of intersection";
1405
1406   status = ShapeExtend_DONE4;
1407   if(sfw->StatusSelfIntersection(status))
1408     di << "\n ShapeExtend_DONE4 : pcurve(s) was(were) modified";
1409
1410   status = ShapeExtend_DONE5;
1411   if(sfw->StatusSelfIntersection(status))
1412     di << "\n ShapeExtend_DONE5 : non adjacent intersection fixed by increasing tolerance of vertex(vertices)";
1413
1414   status = ShapeExtend_DONE6;
1415   if(sfw->StatusSelfIntersection(status))
1416     di << "\n ShapeExtend_DONE6 : tolerance of edge was increased to hide intersection";
1417
1418   status = ShapeExtend_DONE7;
1419   if(sfw->StatusSelfIntersection(status))
1420     di << "\n ShapeExtend_DONE7 : range of some edges was decreased to avoid intersection";
1421
1422   TopoDS_Wire finalwire = sfw->Wire();
1423
1424   DBRep::Set(argv[1],finalwire);
1425
1426
1427   advWA->Load(TopoDS::Wire(finalwire));
1428
1429   TopoDS_Face fface =
1430     BRepBuilderAPI_MakeFace(TopoDS::Wire(finalwire),Standard_True);
1431
1432   DBRep::Set(argv[2],fface);
1433
1434   advWA->SetFace(fface);
1435
1436   TopTools_IndexedMapOfShape fM;
1437   TopExp::MapShapes(finalwire, TopAbs_EDGE, fM);
1438
1439   for(j= 1; j<=fM.Extent(); ++j) {
1440     Standard_Integer num=1;
1441     IntRes2d_SequenceOfIntersectionPoint points2d;
1442     TColgp_SequenceOfPnt points3d;
1443
1444     di << "\n j =" << j << ",  CheckSelfIntersectingEdge = " <<
1445       (Standard_Integer) advWA->CheckSelfIntersectingEdge(j, points2d, points3d);
1446
1447
1448     status = ShapeExtend_FAIL1;
1449     if(advWA->StatusSelfIntersection(status))
1450
1451       status = ShapeExtend_FAIL2;
1452     if(advWA->StatusSelfIntersection(status))
1453       di << "\n No Vertices found in the edge";
1454
1455     status = ShapeExtend_DONE1;
1456     if(advWA->StatusSelfIntersection(status))
1457       di << "\n Self-intersection found in the edge";
1458
1459     num = points2d.Length();
1460     di << "\n No. of self-intersecting edges : " << num;
1461
1462     for(int i=1; i<=num; ++i) {
1463       gp_Pnt pt = points3d(i);
1464       di << "\n Intersecting pt : (" << pt.X() << ", " << pt.Y() << ", " << pt.Z() << ")";
1465     }
1466
1467   }
1468
1469   di << "\n";
1470
1471   return 0;
1472 }
1473
1474 void QABugs::Commands_17(Draw_Interpretor& theCommands) {
1475   const char *group = "QABugs";
1476
1477   theCommands.Add ("BUC60842", "BUC60842", __FILE__, BUC60842, group);
1478   theCommands.Add ("BUC60843", "BUC60843 result_shape name_of_circle name_of_curve [ par1 [ tol ] ]", __FILE__, BUC60843, group);
1479   theCommands.Add ("BUC60970", "BUC60970 shape result", __FILE__, BUC60970, group);
1480   theCommands.Add ("BUC60915", "BUC60915", __FILE__, BUC60915_1, group);
1481   theCommands.Add ("OCC138", "OCC138", __FILE__, OCC138, group);
1482   theCommands.Add ("OCC353","OCC353",__FILE__,OCC353,group);
1483   theCommands.Add ("OCC138LC", "OCC138LC", __FILE__, OCC138LC, group);
1484   theCommands.Add ("OCC566", "OCC566 shape [ xmin ymin zmin xmax ymax zmax] ; print bounding box", __FILE__, OCC566, group);
1485   theCommands.Add ("OCC570", "OCC570 result", __FILE__, OCC570, group);
1486
1487   theCommands.Add("OCC570mkevol",
1488     "OCC570mkevol result object (then use updatevol) [R/Q/P]; mkevol",
1489     __FILE__, MKEVOL,group);
1490   theCommands.Add("OCC570updatevol",
1491     "OCC570updatevol edge u1 rad1 u2 rad2 ...; updatevol",
1492     __FILE__, UPDATEVOL,group);
1493   theCommands.Add("OCC570updatevollaw",
1494     "OCC570updatevollaw edge u1 rad1 u2 rad2 ...; updatevollaw",
1495     __FILE__, UPDATEVOL,group);
1496   theCommands.Add("OCC570buildevol",
1497     "OCC570buildevol; end of the evol fillet computation",
1498     __FILE__, BUILDEVOL,group);
1499
1500   theCommands.Add ("OCC606", "OCC606 result shape [-t]", __FILE__, OCC606, group);
1501
1502   theCommands.Add ("OCC813", "OCC813 U V", __FILE__, OCC813, group);
1503   theCommands.Add ("OCC814", "OCC814", __FILE__, OCC814, group);
1504   theCommands.Add ("OCC884", "OCC884 result shape [toler [maxtoler]]", __FILE__, OCC884, group);
1505
1506   theCommands.Add ("OCC1174_1", "OCC1174_1 shape", __FILE__, OCC1174_1, group);
1507   theCommands.Add ("OCC1174_2", "OCC1174_2 shape", __FILE__, OCC1174_2, group);
1508
1509   theCommands.Add ("OCCN1", "OCCN1 angle fuse(1 for boss / 0 for slot) length", __FILE__, OCCN1, group);
1510   theCommands.Add ("OCCN2", "OCCN2", __FILE__, OCCN2, group);
1511
1512   theCommands.Add ("OCC2569", "OCC2569 nbpoles result", __FILE__, OCC2569, group);
1513
1514   theCommands.Add ("OCC1642", "OCC1642 FinalWare FinalFace InitWare InitFace shape FixReorder FixDegenerated FixConnected FixSelfIntersection", __FILE__, OCC1642, group);
1515
1516   return;
1517 }