0030239: Visualization, Graphic3d_ArrayOfPrimitives - pass Graphic3d_ArrayFlags bitma...
[occt.git] / src / ViewerTest / ViewerTest_ObjectCommands.cxx
1 // Created on: 1998-11-12
2 // Created by: Robert COUBLANC
3 // Copyright (c) 1998-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17 #include <ViewerTest.hxx>
18
19 #include <AIS_PlaneTrihedron.hxx>
20
21 #include <Quantity_NameOfColor.hxx>
22 #include <Draw_Interpretor.hxx>
23 #include <Draw.hxx>
24 #include <Draw_Appli.hxx>
25 #include <DBRep.hxx>
26 #include <DBRep_DrawableShape.hxx>
27
28 #include <Font_BRepFont.hxx>
29 #include <Font_BRepTextBuilder.hxx>
30 #include <Font_FontMgr.hxx>
31
32 #include <NCollection_List.hxx>
33
34 #include <OSD_Chronometer.hxx>
35 #include <TCollection_AsciiString.hxx>
36 #include <V3d_Viewer.hxx>
37 #include <V3d_View.hxx>
38 #include <V3d.hxx>
39
40 #include <AIS_Shape.hxx>
41 #include <AIS_DisplayMode.hxx>
42 #include <AIS_PointCloud.hxx>
43 #include <TColStd_MapOfInteger.hxx>
44 #include <AIS_MapOfInteractive.hxx>
45 #include <ViewerTest_AutoUpdater.hxx>
46 #include <ViewerTest_DoubleMapOfInteractiveAndName.hxx>
47 #include <ViewerTest_DoubleMapIteratorOfDoubleMapOfInteractiveAndName.hxx>
48 #include <ViewerTest_EventManager.hxx>
49
50 #include <TopoDS_Solid.hxx>
51 #include <BRepTools.hxx>
52 #include <BRep_Builder.hxx>
53 #include <TopAbs_ShapeEnum.hxx>
54
55 #include <TopoDS.hxx>
56 #include <BRep_Tool.hxx>
57 #include <TopExp_Explorer.hxx>
58
59 #include <BRepAdaptor_Curve.hxx>
60 #include <BRepAdaptor_Surface.hxx>
61
62 #include <TopAbs.hxx>
63 #include <TopExp.hxx>
64 #include <TopoDS_Vertex.hxx>
65 #include <TopoDS_Shape.hxx>
66 #include <TopoDS_Face.hxx>
67
68 #include <Draw_Window.hxx>
69 #include <AIS_ListIteratorOfListOfInteractive.hxx>
70 #include <AIS_ListOfInteractive.hxx>
71 #include <AIS_ColoredShape.hxx>
72 #include <AIS_DisplayMode.hxx>
73 #include <AIS_Shape.hxx>
74
75 #include <AIS_InteractiveContext.hxx>
76 #include <Geom_Plane.hxx>
77 #include <gp_Pln.hxx>
78 #include <TCollection_ExtendedString.hxx>
79 #include <TCollection_HAsciiString.hxx>
80 #include <GC_MakePlane.hxx>
81 #include <gp_Circ.hxx>
82 #include <AIS_Axis.hxx>
83 #include <Geom_Axis2Placement.hxx>
84 #include <Geom_Axis1Placement.hxx>
85 #include <AIS_Trihedron.hxx>
86 #include <AIS_Axis.hxx>
87 #include <gp_Trsf.hxx>
88 #include <gp_Quaternion.hxx>
89 #include <TopLoc_Location.hxx>
90
91 #include <HLRAlgo_Projector.hxx>
92 #include <HLRBRep_Algo.hxx>
93 #include <HLRBRep_PolyAlgo.hxx>
94 #include <HLRBRep_PolyHLRToShape.hxx>
95 #include <HLRBRep_HLRToShape.hxx>
96 #include <Aspect_Window.hxx>
97
98 #include <Graphic3d_ArrayOfPoints.hxx>
99 #include <Graphic3d_ArrayOfSegments.hxx>
100 #include <Graphic3d_ArrayOfPolylines.hxx>
101 #include <Graphic3d_ArrayOfTriangles.hxx>
102 #include <Graphic3d_ArrayOfTriangleFans.hxx>
103 #include <Graphic3d_ArrayOfTriangleStrips.hxx>
104 #include <Graphic3d_ArrayOfQuadrangles.hxx>
105 #include <Graphic3d_ArrayOfQuadrangleStrips.hxx>
106 #include <Graphic3d_ArrayOfPolygons.hxx>
107 #include <Graphic3d_AspectMarker3d.hxx>
108 #include <Graphic3d_Group.hxx>
109 #include <Standard_Real.hxx>
110
111 #include <AIS_Circle.hxx>
112 #include <BRepBuilderAPI_MakeEdge.hxx>
113 #include <BRepBuilderAPI_MakeFace.hxx>
114 #include <BRepBuilderAPI_MakeWire.hxx>
115 #include <Geom_Circle.hxx>
116 #include <GC_MakeCircle.hxx>
117 #include <Select3D_SensitiveCircle.hxx>
118 #include <SelectMgr_EntityOwner.hxx>
119 #include <SelectMgr_Selection.hxx>
120 #include <StdFail_NotDone.hxx>
121 #include <StdPrs_ShadedShape.hxx>
122 #include <TopoDS_Wire.hxx>
123
124 #include <AIS_MultipleConnectedInteractive.hxx>
125 #include <AIS_ConnectedInteractive.hxx>
126 #include <AIS_TextLabel.hxx>
127 #include <TopLoc_Location.hxx>
128 #include <TColStd_ListOfInteger.hxx>
129 #include <TColStd_ListIteratorOfListOfInteger.hxx>
130
131 #include <Select3D_SensitiveSegment.hxx>
132 #include <Select3D_SensitivePrimitiveArray.hxx>
133 #include <Select3D_SensitiveTriangle.hxx>
134 #include <Select3D_SensitiveCurve.hxx>
135 #include <Select3D_SensitivePoint.hxx>
136 #include <BRepAdaptor_Curve.hxx>
137 #include <StdPrs_Curve.hxx>
138
139 #include <BRepExtrema_ExtPC.hxx>
140 #include <BRepExtrema_ExtPF.hxx>
141
142 #include <Prs3d_Arrow.hxx>
143 #include <Prs3d_ArrowAspect.hxx>
144 #include <Prs3d_DatumAttribute.hxx>
145 #include <Prs3d_DatumAspect.hxx>
146 #include <Prs3d_Drawer.hxx>
147 #include <Prs3d_VertexDrawMode.hxx>
148 #include <Prs3d_LineAspect.hxx>
149 #include <Prs3d_PointAspect.hxx>
150 #include <Prs3d_Presentation.hxx>
151 #include <Prs3d_TextAspect.hxx>
152
153 #include <Image_AlienPixMap.hxx>
154 #include <TColStd_HArray1OfAsciiString.hxx>
155 #include <TColStd_HSequenceOfAsciiString.hxx>
156
157 extern ViewerTest_DoubleMapOfInteractiveAndName& GetMapOfAIS();
158 extern Standard_Boolean VDisplayAISObject (const TCollection_AsciiString& theName,
159                                            const Handle(AIS_InteractiveObject)& theAISObj,
160                                            Standard_Boolean theReplaceIfExists = Standard_True);
161 extern int ViewerMainLoop(Standard_Integer argc, const char** argv);
162 extern Handle(AIS_InteractiveContext)& TheAISContext();
163
164 namespace
165 {
166   static bool convertToColor (const Handle(TColStd_HSequenceOfAsciiString)& theColorValues,
167                               Quantity_Color& theColor)
168   {
169     const char* anArgs[3] =
170     {
171       theColorValues->Size() >= 1 ? theColorValues->Value (1).ToCString() : "",
172       theColorValues->Size() >= 2 ? theColorValues->Value (2).ToCString() : "",
173       theColorValues->Size() >= 3 ? theColorValues->Value (3).ToCString() : ""
174     };
175     return ViewerTest::ParseColor (theColorValues->Size(), anArgs, theColor) != 0;
176   }
177
178   static bool convertToDatumPart (const TCollection_AsciiString& theValue,
179                                   Prs3d_DatumParts& theDatumPart)
180   {
181     TCollection_AsciiString aValue = theValue;
182     aValue.LowerCase();
183     if      (aValue == "origin")  theDatumPart = Prs3d_DP_Origin;
184     else if (aValue == "xaxis")   theDatumPart = Prs3d_DP_XAxis;
185     else if (aValue == "yaxis")   theDatumPart = Prs3d_DP_YAxis;
186     else if (aValue == "zaxis")   theDatumPart = Prs3d_DP_ZAxis;
187     else if (aValue == "xarrow")  theDatumPart = Prs3d_DP_XArrow;
188     else if (aValue == "yarrow")  theDatumPart = Prs3d_DP_YArrow;
189     else if (aValue == "zarrow")  theDatumPart = Prs3d_DP_ZArrow;
190     else if (aValue == "xoyaxis") theDatumPart = Prs3d_DP_XOYAxis;
191     else if (aValue == "yozaxis") theDatumPart = Prs3d_DP_YOZAxis;
192     else if (aValue == "xozaxis") theDatumPart = Prs3d_DP_XOZAxis;
193     else if (aValue == "whole")   theDatumPart = Prs3d_DP_None;
194     else
195     {
196       return false;
197     }
198     return true;
199   }
200
201   static void convertToDatumParts (const TCollection_AsciiString& theValue,
202                                    NCollection_List<Prs3d_DatumParts>& theParts)
203   {
204     TCollection_AsciiString aValue = theValue;
205     const Standard_Integer aSplitPos = theValue.Search ("|");
206     Prs3d_DatumParts aPart = Prs3d_DP_None;
207     if (aSplitPos > 0)
208     {
209       convertToDatumParts (theValue.SubString (aSplitPos + 1, theValue.Length()), theParts);
210       if (aSplitPos == 1) // first symbol
211       {
212         return;
213       }
214       aValue = theValue.SubString (1, aSplitPos - 1);
215     }
216     if (convertToDatumPart (aValue, aPart))
217     {
218       theParts.Append (aPart);
219     }
220   }
221
222   static bool convertToDatumAttribute (const TCollection_AsciiString& theValue,
223                                        Prs3d_DatumAttribute& theAttribute)
224   {
225     TCollection_AsciiString aValue = theValue;
226     aValue.LowerCase();
227     if      (aValue == "xaxislength")       theAttribute = Prs3d_DA_XAxisLength;
228     else if (aValue == "yaxislength")       theAttribute = Prs3d_DA_YAxisLength;
229     else if (aValue == "zaxislength")       theAttribute = Prs3d_DA_ZAxisLength;
230     else if (aValue == "tuberadiuspercent") theAttribute = Prs3d_DP_ShadingTubeRadiusPercent;
231     else if (aValue == "coneradiuspercent") theAttribute = Prs3d_DP_ShadingConeRadiusPercent;
232     else if (aValue == "conelengthpercent") theAttribute = Prs3d_DP_ShadingConeLengthPercent;
233     else if (aValue == "originradiuspercent") theAttribute = Prs3d_DP_ShadingOriginRadiusPercent;
234     else if (aValue == "shadingnumberoffacettes") theAttribute = Prs3d_DP_ShadingNumberOfFacettes;
235     else
236       return false;
237     return true;
238   }
239
240   static void convertToDatumAttributes (const TCollection_AsciiString& theValue,
241                                         NCollection_List<Prs3d_DatumAttribute>& theAttributes)
242   {
243     TCollection_AsciiString aValue = theValue;
244     const Standard_Integer aSplitPos = theValue.Search ("|");
245     Prs3d_DatumAttribute anAttribute = Prs3d_DA_XAxisLength;
246     if (aSplitPos > 0)
247     {
248       convertToDatumAttributes (theValue.SubString (aSplitPos + 1, theValue.Length()), theAttributes);
249       if (aSplitPos == 1) // first symbol
250       {
251         return;
252       }
253       aValue = theValue.SubString (1, aSplitPos - 1);
254     }
255     if (convertToDatumAttribute (aValue, anAttribute))
256     {
257       theAttributes.Append (anAttribute);
258     }
259   }
260
261   static bool convertToDatumAxes (const TCollection_AsciiString& theValue,
262                                   Prs3d_DatumAxes& theDatumAxes)
263   {
264     TCollection_AsciiString aValue = theValue;
265     aValue.LowerCase();
266     if      (aValue == "x")   theDatumAxes = Prs3d_DA_XAxis;
267     else if (aValue == "y")   theDatumAxes = Prs3d_DA_YAxis;
268     else if (aValue == "z")   theDatumAxes = Prs3d_DA_ZAxis;
269     else if (aValue == "xy")  theDatumAxes = Prs3d_DA_XYAxis;
270     else if (aValue == "zy")  theDatumAxes = Prs3d_DA_YZAxis;
271     else if (aValue == "xz")  theDatumAxes = Prs3d_DA_XZAxis;
272     else if (aValue == "xyz") theDatumAxes = Prs3d_DA_XYZAxis;
273     else
274     {
275       return false;
276     }
277     return true;
278   }
279
280   static Standard_Boolean setTrihedronParams (Standard_Integer  theArgsNb,
281                                               const char** theArgVec,
282                                               Handle(AIS_Trihedron) theTrihedron)
283   {
284     NCollection_DataMap<TCollection_AsciiString, Handle(TColStd_HSequenceOfAsciiString)> aMapOfArgs;
285     TCollection_AsciiString aParseKey;
286     for (Standard_Integer anArgIt = 1; anArgIt < theArgsNb; ++anArgIt)
287     {
288       TCollection_AsciiString anArg (theArgVec [anArgIt]);
289       if (anArg.Value (1) == '-'
290       && !anArg.IsRealValue())
291       {
292         aParseKey = anArg;
293         aParseKey.Remove (1);
294         aParseKey.LowerCase();
295         std::string aKey = aParseKey.ToCString();
296         aMapOfArgs.Bind (aParseKey, new TColStd_HSequenceOfAsciiString());
297         continue;
298       }
299
300       if (aParseKey.IsEmpty())
301       {
302         continue;
303       }
304
305       aMapOfArgs (aParseKey)->Append (anArg);
306     }
307
308     // Check parameters
309     if ((aMapOfArgs.IsBound ("xaxis") && !aMapOfArgs.IsBound ("zaxis"))
310     || (!aMapOfArgs.IsBound ("xaxis") &&  aMapOfArgs.IsBound ("zaxis")))
311     {
312       std::cout << "Syntax error: -xaxis and -zaxis parameters are to set together.\n";
313       return Standard_False;
314     }
315
316     Handle(TColStd_HSequenceOfAsciiString) aValues;
317     Handle(Geom_Axis2Placement) aComponent = theTrihedron->Component();
318     if (aMapOfArgs.Find ("origin", aValues))
319     {
320       aComponent->SetLocation (gp_Pnt (aValues->Value (1).RealValue(),
321                                        aValues->Value (2).RealValue(),
322                                        aValues->Value (3).RealValue()));
323     }
324     Handle(TColStd_HSequenceOfAsciiString) aXValues, aZValues;
325     if (aMapOfArgs.Find ("xaxis", aXValues) && aMapOfArgs.Find ("zaxis", aZValues))
326     {
327       gp_Dir aXDir (aXValues->Value (1).RealValue(),
328                     aXValues->Value (2).RealValue(),
329                     aXValues->Value (3).RealValue());
330
331       gp_Dir aZDir (aZValues->Value (1).RealValue(),
332                     aZValues->Value (2).RealValue(),
333                     aZValues->Value (3).RealValue());
334
335       if (!aZDir.IsNormal (aXDir, M_PI / 180.0))
336       {
337         std::cout << "Syntax error - parameters 'xaxis' and 'zaxis' are not applied as VectorX is not normal to VectorZ\n";
338         return Standard_False;
339       }
340
341       aComponent->SetXDirection(aXDir);
342       aComponent->SetDirection (aZDir);
343     }
344
345     if (aMapOfArgs.Find ("dispmode", aValues))
346     {
347       TCollection_AsciiString aValue (aValues->Value (1));
348       bool isWireframe = true;
349       if (aValue.IsEqual ("sh") || aValue.IsEqual ("shading"))
350         isWireframe = false;
351       theTrihedron->SetDatumDisplayMode (isWireframe ? Prs3d_DM_WireFrame
352                                                      : Prs3d_DM_Shaded);
353     }
354
355     if (aMapOfArgs.Find ("hidelabels", aValues))
356     {
357       Standard_Boolean toHideLabels = Standard_True;
358       if (aValues->Size() == 1)
359       {
360         ViewerTest::ParseOnOff (aValues->First().ToCString(), toHideLabels);
361       }
362       else if (aValues->Size() != 0)
363       {
364         std::cout << "Syntax error: -hidelabels expects parameter 'on' or 'off' after.\n";
365         return Standard_False;
366       }
367
368       if (!theTrihedron->Attributes()->HasOwnDatumAspect())
369       {
370         theTrihedron->Attributes()->SetDatumAspect (new Prs3d_DatumAspect());
371       }
372       theTrihedron->Attributes()->DatumAspect()->SetDrawLabels (!toHideLabels);
373     }
374
375     if (aMapOfArgs.Find ("hidearrows", aValues))
376     {
377       Standard_Boolean toHideArrows = Standard_True;
378       if (aValues->Size() == 1)
379       {
380         ViewerTest::ParseOnOff (aValues->First().ToCString(), toHideArrows);
381       }
382       else if (aValues->Size() != 0)
383       {
384         std::cout << "Syntax error: -hidearrows expects parameter 'on' or 'off' after.\n";
385         return Standard_False;
386       }
387
388       if (!theTrihedron->Attributes()->HasOwnDatumAspect())
389       {
390         theTrihedron->Attributes()->SetDatumAspect (new Prs3d_DatumAspect());
391       }
392       theTrihedron->Attributes()->DatumAspect()->SetDrawArrows (!toHideArrows);
393     }
394
395     if (aMapOfArgs.Find ("color", aValues))
396     {
397       NCollection_List<Prs3d_DatumParts> aParts;
398       if (aValues->Size() < 2)
399       {
400         std::cout << "Syntax error: -color wrong parameters.\n";
401         return Standard_False;
402       }
403
404       convertToDatumParts (aValues->Value(1), aParts);
405       aValues->Remove (1); // datum part is processed
406       Quantity_Color aColor;
407       if (!convertToColor (aValues, aColor))
408       {
409         std::cout << "Syntax error: -color wrong parameters.\n";
410         return Standard_False;
411       }
412
413       for (NCollection_List<Prs3d_DatumParts>::Iterator anIterator (aParts); anIterator.More(); anIterator.Next())
414       {
415         Prs3d_DatumParts aDatumPart = anIterator.Value();
416         if (aDatumPart == Prs3d_DP_None)
417         {
418           theTrihedron->SetColor (aColor);
419         }
420         else
421         {
422           theTrihedron->SetDatumPartColor (aDatumPart, aColor);
423         }
424       }
425     }
426
427     if (aMapOfArgs.Find ("textcolor", aValues))
428     {
429       Quantity_Color aColor;
430       if (!convertToColor (aValues, aColor))
431       {
432         std::cout << "Syntax error: -textcolor wrong parameters.\n";
433         return Standard_False;
434       }
435       theTrihedron->SetTextColor (aColor);
436     }
437
438     if (aMapOfArgs.Find ("arrowcolor", aValues))
439     {
440       Quantity_Color aColor;
441       if (!convertToColor (aValues, aColor))
442       {
443         std::cout << "Syntax error: -arrowcolor wrong parameters.\n";
444         return Standard_False;
445       }
446       theTrihedron->SetArrowColor (aColor);
447     }
448
449     if (aMapOfArgs.Find ("attribute", aValues))
450     {
451       NCollection_List<Prs3d_DatumAttribute> anAttributes;
452       if (aValues->Size() != 2)
453       {
454         std::cout << "Syntax error: -attribute wrong parameters.\n";
455         return Standard_False;
456       }
457
458       convertToDatumAttributes (aValues->Value (1), anAttributes);
459       if (!theTrihedron->Attributes()->HasOwnDatumAspect())
460         theTrihedron->Attributes()->SetDatumAspect(new Prs3d_DatumAspect());
461       for (NCollection_List<Prs3d_DatumAttribute>::Iterator anIterator (anAttributes); anIterator.More(); anIterator.Next())
462       {
463         theTrihedron->Attributes()->DatumAspect()->SetAttribute (anIterator.Value(), aValues->Value (2).RealValue());
464       }
465     }
466
467     if (aMapOfArgs.Find ("priority", aValues))
468     {
469       Prs3d_DatumParts aDatumPart;
470       if (aValues->Size() < 2
471       || !convertToDatumPart (aValues->Value (1), aDatumPart))
472       {
473         std::cout << "Syntax error: -priority wrong parameters.\n";
474         return Standard_False;
475       }
476       theTrihedron->SetSelectionPriority (aDatumPart, aValues->Value (2).IntegerValue());
477     }
478
479     if (aMapOfArgs.Find ("labels", aValues))
480     {
481       Prs3d_DatumParts aDatumPart = Prs3d_DP_None;
482       if (aValues->Size() > 2
483        && convertToDatumPart(aValues->Value(1), aDatumPart)
484        && aDatumPart >= Prs3d_DP_XAxis
485        && aDatumPart <= Prs3d_DP_ZAxis) // labels are set to axes only
486       {
487         theTrihedron->SetLabel (aDatumPart, aValues->Value (2));
488       }
489       else
490       {
491         std::cout << "Syntax error: -labels wrong parameters.\n";
492         return Standard_False;
493       }
494     }
495
496     if (aMapOfArgs.Find ("drawaxes", aValues))
497     {
498       Prs3d_DatumAxes aDatumAxes = Prs3d_DA_XAxis;
499       if (aValues->Size() < 1
500       || !convertToDatumAxes (aValues->Value (1), aDatumAxes))
501       {
502         std::cout << "Syntax error: -drawaxes wrong parameters.\n";
503         return Standard_False;
504       }
505       if (!theTrihedron->Attributes()->HasOwnDatumAspect())
506         theTrihedron->Attributes()->SetDatumAspect(new Prs3d_DatumAspect());
507       theTrihedron->Attributes()->DatumAspect()->SetDrawDatumAxes (aDatumAxes);
508     }
509     return Standard_True;
510   }
511 }
512
513 //==============================================================================
514 //function : Vtrihedron 2d
515 //purpose  : Create a plane with a 2D  trihedron from a faceselection
516 //Draw arg : vtri2d  name
517 //==============================================================================
518 static int VTrihedron2D (Draw_Interpretor& /*theDI*/,
519                          Standard_Integer  theArgsNum,
520                          const char**      theArgVec)
521 {
522   if (theArgsNum != 2)
523   {
524     std::cerr << theArgVec[0]<< " error.\n";
525     return 1;
526   }
527
528   TopTools_ListOfShape aShapes;
529   ViewerTest::GetSelectedShapes (aShapes);
530
531   if (aShapes.Extent() != 1)
532   {
533     std::cerr << "Error: wrong number of selected shapes.\n";
534     return 1;
535   }
536
537   const TopoDS_Shape& aShape = aShapes.First();
538
539   TopoDS_Face     aFace = TopoDS::Face (aShape);
540   TopExp_Explorer aFaceExp (aFace, TopAbs_EDGE);
541   TopoDS_Edge     anEdge0 = TopoDS::Edge (aFaceExp.Current());
542
543   gp_Pnt A,B,C;
544   if (aFaceExp.More())
545   {
546     aFaceExp.Next();
547     TopoDS_Edge anEdge1 = TopoDS::Edge (aFaceExp.Current() );
548     BRepAdaptor_Curve aCurve0 (anEdge0);
549     BRepAdaptor_Curve aCurve1 (anEdge1);
550     A = aCurve1.Value (0.1);
551     B = aCurve1.Value (0.9);
552     C = aCurve0.Value (0.5);
553   }
554   else
555   {
556     BRepAdaptor_Curve aCurve0 (anEdge0);
557     A = aCurve0.Value (0.1);
558     B = aCurve0.Value (0.9);
559     C = aCurve0.Value (0.5);
560   }
561
562   GC_MakePlane aMkPlane (A,B,C);
563
564   Handle(AIS_PlaneTrihedron) anAISPlaneTri = new AIS_PlaneTrihedron (aMkPlane.Value());
565   TCollection_AsciiString aName (theArgVec[1]);
566
567   VDisplayAISObject (aName, anAISPlaneTri);
568
569   return 0;
570 }
571
572 //=======================================================================
573 //function : VTrihedron
574 //purpose  :
575 //=======================================================================
576 static int VTrihedron (Draw_Interpretor& ,
577                        Standard_Integer theArgsNb,
578                        const char** theArgVec)
579 {
580   if (theArgsNb < 2)
581   {
582     std::cout << "Syntax error: the wrong number of input parameters.\n";
583     return 1;
584   }
585
586   TCollection_AsciiString aName (theArgVec[1]);
587   gp_Pln aWorkingPlane;
588   Standard_Boolean toUpdate = Standard_True;
589
590   NCollection_DataMap<TCollection_AsciiString, Standard_Real> aRealParams;
591   NCollection_DataMap<TCollection_AsciiString, TCollection_AsciiString> aStringParams;
592
593   Handle(AIS_Trihedron) aTrihedron;
594   Handle(AIS_InteractiveObject) anObject;
595   if (GetMapOfAIS().Find2 (aName, anObject))
596   {
597     aTrihedron = Handle(AIS_Trihedron)::DownCast (anObject);
598     if (aTrihedron.IsNull())
599     {
600       std::cout << "Syntax error: no trihedron with this name.\n";
601       return 1;
602     }
603   }
604   else
605   {
606     Handle(Geom_Axis2Placement) aPlacement = new Geom_Axis2Placement (gp_Pnt (0.0, 0.0, 0.0),
607                                                                       gp::DZ(), gp::DX());
608     aTrihedron = new AIS_Trihedron (aPlacement);
609   }
610
611   if (!setTrihedronParams (theArgsNb, theArgVec, aTrihedron))
612   {
613     return 1;
614   }
615
616   // Redisplay a dimension after parameter changing.
617   if (ViewerTest::GetAISContext()->IsDisplayed (aTrihedron))
618   {
619     ViewerTest::GetAISContext()->Redisplay (aTrihedron, toUpdate);
620   }
621   else
622   {
623     VDisplayAISObject (theArgVec[1], aTrihedron);
624   }
625
626   return 0;
627 }
628
629 //==============================================================================
630 //function : VSize
631 //author   : ege
632 //purpose  : Change the size of a named or selected trihedron
633 //           if no name : it affects the trihedrons witch are selected otherwise nothing is donne
634 //           if no value, the value is set at 100 by default
635 //Draw arg : vsize [name] [size]
636 //==============================================================================
637
638 static int VSize (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
639
640 {
641   // Declaration de booleens
642   Standard_Boolean             ThereIsName;
643   Standard_Boolean             ThereIsCurrent;
644   Standard_Real                value;
645   Standard_Boolean             hascol;
646
647   Quantity_Color col = Quantity_NOC_BLACK;
648
649   // Verification des arguments
650   if ( argc>3 ) {di<<argv[0]<<" Syntaxe error\n"; return 1;}
651
652   // Verification du nombre d'arguments
653   if (argc==1)      {ThereIsName=Standard_False;value=100;}
654   else if (argc==2) {ThereIsName=Standard_False;value=Draw::Atof(argv[1]);}
655   else              {ThereIsName=Standard_True;value=Draw::Atof(argv[2]);}
656
657   // On set le booleen ThereIsCurrent
658   if (TheAISContext() -> NbSelected() > 0) {ThereIsCurrent=Standard_True;}
659   else {ThereIsCurrent=Standard_False;}
660
661
662
663   //===============================================================
664   // Il n'y a pas de nom  mais des objets selectionnes
665   //===============================================================
666   if (!ThereIsName && ThereIsCurrent)
667   {
668
669     ViewerTest_DoubleMapIteratorOfDoubleMapOfInteractiveAndName
670       it (GetMapOfAIS());
671
672     while ( it.More() ) {
673
674       Handle(AIS_InteractiveObject) aShape = it.Key1();
675
676       if (!aShape.IsNull() &&  TheAISContext()->IsSelected(aShape) )
677       {
678
679         // On verifie que l'AIS InteraciveObject selectionne est bien
680         // un AIS_Trihedron
681         if (aShape->Type()==AIS_KOI_Datum && aShape->Signature()==3) {
682
683           if (aShape->HasColor())
684           {
685             hascol = Standard_True;
686
687             // On recupere la couleur de aShape
688             aShape->Color (col);
689           }
690           else
691           {
692             hascol = Standard_False;
693           }
694
695           // On downcast aShape  de AIS_InteractiveObject a AIS_Trihedron
696           // pour lui appliquer la methode SetSize()
697           Handle(AIS_Trihedron) aTrihedron = Handle(AIS_Trihedron)::DownCast (aShape);
698
699           // C'est bien un triedre,on chage sa valeur!
700           aTrihedron->SetSize(value);
701
702           // On donne la couleur au Trihedron
703           if(hascol)   aTrihedron->SetColor(col);
704           else         aTrihedron->UnsetColor();
705
706
707           // The trihedron hasn't be errased from the map
708           // so you just have to redisplay it
709           TheAISContext() ->Redisplay(aTrihedron,Standard_False);
710
711         }
712
713       }
714
715       it.Next();
716     }
717
718     TheAISContext() ->UpdateCurrentViewer();
719   }
720
721   //===============================================================
722   // Il n'y a pas d'arguments et aucuns objets selectionne Rien A Faire!
723   //===============================================================
724
725
726
727   //===============================================================
728   // Il y a un nom de triedre passe en argument
729   //===============================================================
730   if (ThereIsName) {
731     TCollection_AsciiString name=argv[1];
732
733     // on verifie que ce nom correspond bien a une shape
734     Handle(AIS_InteractiveObject) aShape;
735     if (GetMapOfAIS().Find2(name, aShape))
736     {
737       // On verifie que l'AIS InteraciveObject est bien
738       // un AIS_Trihedron
739       if (!aShape.IsNull() &&
740         aShape->Type()==AIS_KOI_Datum && aShape->Signature()==3)
741       {
742         if (aShape->HasColor())
743         {
744           hascol=Standard_True;
745
746           // On recupere la couleur de aShape
747           aShape->Color (col);
748         }
749         else
750         {
751           hascol = Standard_False;
752         }
753
754         // On downcast aShape de AIS_InteractiveObject a AIS_Trihedron
755         // pour lui appliquer la methode SetSize()
756         Handle(AIS_Trihedron) aTrihedron = Handle(AIS_Trihedron)::DownCast (aShape);
757
758         // C'est bien un triedre,on chage sa valeur
759         aTrihedron->SetSize(value);
760
761         // On donne la couleur au Trihedron
762         if(hascol)   aTrihedron->SetColor(col);
763         else         aTrihedron->UnsetColor();
764
765         // The trihedron hasn't be errased from the map
766         // so you just have to redisplay it
767         TheAISContext() ->Redisplay(aTrihedron,Standard_False);
768
769         TheAISContext() ->UpdateCurrentViewer();
770       }
771     }
772   }
773   return 0;
774 }
775
776
777 //==============================================================================
778
779 //==============================================================================
780 //function : VPlaneTrihedron
781 //purpose  : Create a plane from a trihedron selection. If no arguments are set, the default
782 //Draw arg : vplanetri  name
783 //==============================================================================
784 #include <AIS_Plane.hxx>
785
786
787
788 static int VPlaneTrihedron (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
789
790 {
791   // Verification des arguments
792   if ( argc!=2) {di<<argv[0]<<" error\n"; return 1;}
793
794   if (TheAISContext()->NbSelected() != 1)
795   {
796     std::cerr << "Error: Wrong number of selected shapes.\n";
797     return 1;
798   }
799
800   TheAISContext()->InitSelected();
801   Handle(AIS_InteractiveObject) aTest = TheAISContext()->SelectedInteractive();
802   Handle(AIS_Plane) aPlane = Handle(AIS_Plane)::DownCast (aTest);
803   if (aPlane.IsNull())
804   {
805     std::cerr << "Error: Selected shape is not a plane.\n";
806     return 1;
807   }
808
809   VDisplayAISObject (argv[1], aPlane);
810
811   return 0;
812 }
813
814
815
816 //==============================================================================
817 // Fonction        First click      2de click
818 //
819 // vaxis           vertex           vertex
820 //                 edge             None
821 // vaxispara       edge             vertex
822 // vaxisortho      edge             Vertex
823 // vaxisinter      Face             Face
824 //==============================================================================
825
826 //==============================================================================
827 //function : VAxisBuilder
828 //purpose  :
829 //Draw arg : vaxis AxisName Xa Ya Za Xb Yb Zb
830 //==============================================================================
831 #include <TopoDS_Edge.hxx>
832 #include <TopoDS_Vertex.hxx>
833 #include <TopExp.hxx>
834 #include <Geom_Line.hxx>
835
836 static int VAxisBuilder(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
837 {
838   // Declarations
839   Standard_Boolean HasArg;
840   TCollection_AsciiString name;
841
842   // Verification
843   if (argc<2 || argc>8 ) {di<<" Syntaxe error\n";return 1;}
844   if (argc==8) HasArg=Standard_True;
845   else HasArg=Standard_False;
846
847   name=argv[1];
848
849   TopTools_ListOfShape aShapes;
850   ViewerTest::GetSelectedShapes (aShapes);
851
852   // Cas ou il y a des arguments
853   // Purpose: Teste le constructeur AIS_Axis::AIS_Axis(x: Line from Geom)
854   if (HasArg) {
855     Standard_Real coord[6];
856     for(Standard_Integer i=0;i<=5;i++){
857       coord[i]=Draw::Atof(argv[2+i]);
858     }
859     gp_Pnt p1(coord[0],coord[1],coord[2]), p2(coord[3],coord[4],coord[5]) ;
860
861     gp_Vec myVect (p1,p2);
862     Handle(Geom_Line) myLine=new Geom_Line (p1 ,myVect );
863     Handle(AIS_Axis) TheAxis=new AIS_Axis (myLine );
864     GetMapOfAIS().Bind (TheAxis,name);
865     TheAISContext()->Display(TheAxis, Standard_True);
866   }
867
868   // Pas d'arguments
869   else {
870     // fonction vaxis
871     // Purpose: Teste le constructeur AIS_Axis::AIS_Axis (x:Axis1Placement from Geom)
872     if ( !strcasecmp(argv[0], "vaxis")) {
873       if (aShapes.Extent() != 2 && aShapes.Extent() != 1)
874       {
875         std::cerr << "Error: Wrong number of selected shapes.\n";
876         return 1;
877       }
878
879       const TopoDS_Shape& aShapeA = aShapes.First();
880       if (aShapeA.ShapeType() == TopAbs_VERTEX)
881       {
882         if (aShapes.Extent() != 2)
883         {
884           std::cerr << "Error: Wron number of selected shapes.\n";
885           return 1;
886         }
887
888         const TopoDS_Shape& aShapeB = aShapes.Last();
889         if (aShapeB.ShapeType() != TopAbs_VERTEX)
890         {
891           std::cerr << "Syntax error: You should select two vertices or one edge.\n";
892           return 1;
893         }
894
895         // Construction de l'axe
896         gp_Pnt A = BRep_Tool::Pnt (TopoDS::Vertex (aShapeA));
897         gp_Pnt B = BRep_Tool::Pnt (TopoDS::Vertex (aShapeB));
898         gp_Vec V (A,B);
899         gp_Dir D (V);
900         Handle(Geom_Axis1Placement) OrigineAndVect=new Geom_Axis1Placement (A,D);
901         Handle(AIS_Axis) TheAxis=new AIS_Axis (OrigineAndVect);
902         GetMapOfAIS().Bind (TheAxis,name);
903         TheAISContext()->Display (TheAxis, Standard_True);
904       }
905       else
906       {
907         TopoDS_Edge    ed =TopoDS::Edge (aShapeA);
908         TopoDS_Vertex  Va,Vb;
909         TopExp::Vertices(ed,Va,Vb );
910         gp_Pnt A=BRep_Tool::Pnt(Va);
911         gp_Pnt B=BRep_Tool::Pnt(Vb);
912         gp_Vec  V (A,B);
913         gp_Dir   D (V);
914         Handle(Geom_Axis1Placement) OrigineAndVect=new Geom_Axis1Placement (A,D);
915         Handle(AIS_Axis) TheAxis=new AIS_Axis (OrigineAndVect);
916         GetMapOfAIS().Bind (TheAxis,name);
917         TheAISContext()->Display (TheAxis, Standard_True);
918       }
919
920     }
921
922     // Fonction axispara
923     // Purpose: Teste le constructeur AIS_Axis::AIS_Axis(x: Axis2Placement from Geom, y: TypeOfAxis from AIS)
924     else if ( !strcasecmp(argv[0], "vaxispara"))
925     {
926       if (aShapes.Extent() != 2)
927       {
928         std::cerr << "Error: Wrong number of selected shapes.\n";
929         return 1;
930       }
931
932       const TopoDS_Shape& aShapeA = aShapes.First();
933       const TopoDS_Shape& aShapeB = aShapes.Last();
934       if (!(aShapeA.ShapeType() == TopAbs_EDGE
935          && aShapeB.ShapeType() == TopAbs_VERTEX))
936       {
937         std::cerr << "Syntax error: You should select face and then vertex.\n";
938         return 1;
939       }
940
941       TopoDS_Edge    ed=TopoDS::Edge (aShapeA);
942       gp_Pnt B=BRep_Tool::Pnt(TopoDS::Vertex(aShapeB));
943       TopoDS_Vertex  Va,Vc;
944       TopExp::Vertices(ed,Va,Vc );
945       gp_Pnt A=BRep_Tool::Pnt(Va);
946       gp_Pnt C=BRep_Tool::Pnt(Vc);
947       gp_Vec  V (A,C);
948       gp_Dir   D (V);
949       Handle(Geom_Axis1Placement) OrigineAndVect=new Geom_Axis1Placement (B,D);
950       Handle(AIS_Axis) TheAxis=new AIS_Axis (OrigineAndVect);
951       GetMapOfAIS().Bind (TheAxis,name);
952       TheAISContext()->Display (TheAxis, Standard_True);
953
954     }
955
956     // Fonction axisortho
957     else
958     {
959       if (aShapes.Extent() != 2)
960       {
961         std::cerr << "Error: Wrong number of selected shapes.\n";
962         return 1;
963       }
964
965       const TopoDS_Shape& aShapeA = aShapes.First();
966       const TopoDS_Shape& aShapeB = aShapes.Last();
967       if (!(aShapeA.ShapeType() == TopAbs_EDGE
968          && aShapeB.ShapeType() == TopAbs_VERTEX))
969       {
970         std::cerr << "Syntax error: You should select face and then vertex.\n";
971         return 1;
972       }
973
974       // Construction de l'axe
975       TopoDS_Edge    ed=TopoDS::Edge(aShapeA) ;
976       gp_Pnt B=BRep_Tool::Pnt(TopoDS::Vertex(aShapeB) );
977       TopoDS_Vertex  Va,Vc;
978       TopExp::Vertices(ed,Va,Vc );
979       gp_Pnt A=BRep_Tool::Pnt(Va);
980       gp_Pnt C=BRep_Tool::Pnt(Vc);
981       gp_Pnt E(A.Y()+A.Z()-C.Y()-C.Z()  ,C.X()-A.X() ,C.X()-A.X() );
982       gp_Vec  V (A,E);
983       gp_Dir   D (V);
984       Handle(Geom_Axis1Placement) OrigineAndVect=new Geom_Axis1Placement (B,D);
985       Handle(AIS_Axis) TheAxis=new AIS_Axis (OrigineAndVect);
986       GetMapOfAIS().Bind (TheAxis,name);
987       TheAISContext()->Display (TheAxis, Standard_True);
988
989     }
990
991   }
992   return 0;
993 }
994
995
996 //==============================================================================
997 // Fonction        First click      Result
998 //
999 // vpoint          vertex           AIS_Point=Vertex
1000 //                 edge             AIS_Point=Middle of the edge
1001 //==============================================================================
1002
1003 //==============================================================================
1004 //function : VPointBuilder
1005 //purpose  : Build an AIS_Point from coordinates or with a selected vertex or edge
1006 //Draw arg : vpoint PoinName [Xa] [Ya] [Za]
1007 //==============================================================================
1008 #include <TopoDS_Edge.hxx>
1009 #include <TopoDS_Vertex.hxx>
1010 #include <TopExp.hxx>
1011 #include <AIS_Point.hxx>
1012 #include <Geom_CartesianPoint.hxx>
1013
1014 static int VPointBuilder(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
1015 {
1016   // Declarations
1017   Standard_Boolean HasArg;
1018   TCollection_AsciiString name;
1019
1020   // Verification
1021   if (argc<2 || argc>5 ) {di<<" Syntaxe error\n";return 1;}
1022   if (argc==5) HasArg=Standard_True;
1023   else HasArg=Standard_False;
1024
1025   name=argv[1];
1026
1027   // Il y a des arguments: teste l'unique constructeur AIS_Pnt::AIS_Pnt(Point from Geom)
1028   if (HasArg) {
1029     Standard_Real thecoord[3];
1030     for(Standard_Integer i=0;i<=2;i++)
1031       thecoord[i]=Draw::Atof(argv[2+i]);
1032     Handle(Geom_CartesianPoint )  myGeomPoint= new Geom_CartesianPoint (thecoord[0],thecoord[1],thecoord[2]);
1033     Handle(AIS_Point)  myAISPoint=new AIS_Point(myGeomPoint );
1034     GetMapOfAIS().Bind (myAISPoint,name);
1035     TheAISContext()->Display (myAISPoint, Standard_True);
1036   }
1037
1038   // Il n'a pas d'arguments
1039   else
1040   {
1041     TopTools_ListOfShape aShapes;
1042     ViewerTest::GetSelectedShapes (aShapes);
1043
1044     if (aShapes.Extent() != 1)
1045     {
1046       std::cerr << "Error: Wrong number of selected shapes.\n";
1047       std::cerr << "\tYou should select one edge or vertex.\n";
1048       return 1;
1049     }
1050
1051     const TopoDS_Shape& aShapeA = aShapes.First();
1052
1053     if (aShapeA.ShapeType()==TopAbs_VERTEX )
1054     {
1055       gp_Pnt A=BRep_Tool::Pnt(TopoDS::Vertex(aShapeA ) );
1056       Handle(Geom_CartesianPoint) myGeomPoint= new Geom_CartesianPoint (A );
1057       Handle(AIS_Point)  myAISPoint = new AIS_Point  (myGeomPoint );
1058       GetMapOfAIS().Bind(myAISPoint,name);
1059       TheAISContext()->Display (myAISPoint, Standard_True);
1060     }
1061     else
1062     {
1063       TopoDS_Edge myEdge=TopoDS::Edge(aShapeA);
1064       TopoDS_Vertex myVertexA,myVertexB;
1065       TopExp::Vertices (myEdge ,myVertexA ,myVertexB );
1066       gp_Pnt A=BRep_Tool::Pnt(myVertexA );
1067       gp_Pnt B=BRep_Tool::Pnt(myVertexB );
1068       // M est le milieu de [AB]
1069       Handle(Geom_CartesianPoint) myGeomPointM= new Geom_CartesianPoint ( (A.X()+B.X())/2  , (A.Y()+B.Y())/2  , (A.Z()+B.Z())/2  );
1070       Handle(AIS_Point)  myAISPointM = new AIS_Point  (myGeomPointM );
1071       GetMapOfAIS().Bind(myAISPointM,name);
1072       TheAISContext()->Display (myAISPointM, Standard_True);
1073     }
1074
1075   }
1076   return 0;
1077
1078 }
1079
1080 //==============================================================================
1081 // Function        1st click   2de click  3de click
1082 // vplane          Vertex      Vertex     Vertex
1083 //                 Vertex      Edge
1084 //                 Edge        Vertex
1085 //                 Face
1086 // vplanepara      Face        Vertex
1087 //                 Vertex      Face
1088 // vplaneortho     Face        Edge
1089 //                 Edge        Face
1090 //==============================================================================
1091
1092 //==============================================================================
1093 //function : VPlaneBuilder
1094 //purpose  : Build an AIS_Plane from selected entities or Named AIS components
1095 //Draw arg : vplane PlaneName [AxisName]  [PointName] [TypeOfSensitivity]
1096 //                            [PointName] [PointName] [PointName] [TypeOfSensitivity]
1097 //                            [PlaneName] [PointName] [TypeOfSensitivity]
1098 //==============================================================================
1099
1100 static Standard_Integer VPlaneBuilder (Draw_Interpretor& /*di*/,
1101                                        Standard_Integer argc,
1102                                        const char** argv)
1103 {
1104   // Declarations
1105   Standard_Boolean hasArg;
1106   TCollection_AsciiString aName;
1107
1108   // Verification
1109   if (argc<2 || argc>6 )
1110   {
1111     std::cout<<" Syntax error\n";
1112     return 1;
1113   }
1114   if (argc == 6 || argc==5 || argc==4)
1115     hasArg=Standard_True;
1116   else 
1117     hasArg=Standard_False;
1118
1119   aName=argv[1];
1120
1121   // There are some arguments
1122   if (hasArg)
1123   {
1124     Handle(AIS_InteractiveObject) aShapeA;
1125     if (!GetMapOfAIS().Find2 (argv[2], aShapeA))
1126     {
1127       std::cout<<"vplane: error 1st name doesn't exist in the GetMapOfAIS()\n";
1128       return 1;
1129     }
1130
1131     // The first argument is an AIS_Point
1132     if (!aShapeA.IsNull() &&
1133         aShapeA->Type()==AIS_KOI_Datum &&
1134         aShapeA->Signature()==1)
1135     {
1136         // The second argument must also be an AIS_Point
1137         Handle(AIS_InteractiveObject) aShapeB;
1138         if (argc<5 || !GetMapOfAIS().Find2 (argv[3], aShapeB))
1139         {
1140           std::cout<<"vplane: error 2nd name doesn't exist in the GetMapOfAIS()\n";
1141           return 1;
1142         }
1143         // If B is not an AIS_Point
1144         if (aShapeB.IsNull() ||
1145           (!(aShapeB->Type()==AIS_KOI_Datum && aShapeB->Signature()==1)))
1146         {
1147           std::cout<<"vplane: error 2nd object is expected to be an AIS_Point.\n";
1148           return 1;
1149         }
1150         // The third object is an AIS_Point
1151         Handle(AIS_InteractiveObject) aShapeC;
1152         if (!GetMapOfAIS().Find2(argv[4], aShapeC)) 
1153         {
1154           std::cout<<"vplane: error 3d name doesn't exist in the GetMapOfAIS().\n";
1155           return 1; 
1156         }
1157         // If C is not an AIS_Point
1158         if (aShapeC.IsNull() ||
1159           (!(aShapeC->Type()==AIS_KOI_Datum && aShapeC->Signature()==1)))
1160         {
1161           std::cout<<"vplane: error 3d object is expected to be an AIS_Point.\n";
1162           return 1;
1163         }
1164
1165         // Treatment of objects A, B, C
1166         // Downcast an AIS_IO to AIS_Point
1167         Handle(AIS_Point) anAISPointA = Handle(AIS_Point)::DownCast( aShapeA);
1168         Handle(AIS_Point) anAISPointB = Handle(AIS_Point)::DownCast( aShapeB);
1169         Handle(AIS_Point) anAISPointC = Handle(AIS_Point)::DownCast( aShapeC);
1170
1171         Handle(Geom_CartesianPoint ) aCartPointA = 
1172           Handle(Geom_CartesianPoint)::DownCast( anAISPointA->Component());
1173
1174         Handle(Geom_CartesianPoint ) aCartPointB = 
1175           Handle(Geom_CartesianPoint)::DownCast( anAISPointB->Component());
1176
1177         Handle(Geom_CartesianPoint ) aCartPointC = 
1178           Handle(Geom_CartesianPoint)::DownCast( anAISPointC->Component());
1179
1180         // Verification that the three points are different
1181         if(Abs(aCartPointB->X()-aCartPointA->X())<=Precision::Confusion() &&
1182            Abs(aCartPointB->Y()-aCartPointA->Y())<=Precision::Confusion() &&
1183            Abs(aCartPointB->Z()-aCartPointA->Z())<=Precision::Confusion())
1184         {
1185           // B=A
1186           std::cout<<"vplane error: same points\n";return 1;
1187         }
1188         if(Abs(aCartPointC->X()-aCartPointA->X())<=Precision::Confusion() &&
1189            Abs(aCartPointC->Y()-aCartPointA->Y())<=Precision::Confusion() &&
1190            Abs(aCartPointC->Z()-aCartPointA->Z())<=Precision::Confusion())
1191         {
1192           // C=A
1193           std::cout<<"vplane error: same points\n";return 1;
1194         }
1195         if(Abs(aCartPointC->X()-aCartPointB->X())<=Precision::Confusion() &&
1196            Abs(aCartPointC->Y()-aCartPointB->Y())<=Precision::Confusion() &&
1197            Abs(aCartPointC->Z()-aCartPointB->Z())<=Precision::Confusion())
1198         {
1199           // C=B
1200           std::cout<<"vplane error: same points\n";return 1;
1201         }
1202
1203         gp_Pnt A = aCartPointA->Pnt();
1204         gp_Pnt B = aCartPointB->Pnt();
1205         gp_Pnt C = aCartPointC->Pnt();
1206
1207         // Construction of AIS_Plane
1208         GC_MakePlane MkPlane (A,B,C);
1209         Handle(Geom_Plane) aGeomPlane = MkPlane.Value();
1210         Handle(AIS_Plane)  anAISPlane = new AIS_Plane(aGeomPlane );
1211         GetMapOfAIS().Bind (anAISPlane,aName );
1212         if (argc == 6)
1213         {
1214           Standard_Integer aType = Draw::Atoi (argv[5]);
1215           if (aType != 0 && aType != 1)
1216           {
1217             std::cout << "vplane error: wrong type of sensitivity!\n"
1218                       << "Should be one of the following values:\n"
1219                       << "0 - Interior\n"
1220                       << "1 - Boundary"
1221                       << std::endl;
1222             return 1;
1223           }
1224           else
1225           {
1226             anAISPlane->SetTypeOfSensitivity (Select3D_TypeOfSensitivity (aType));
1227           }
1228         }
1229         TheAISContext()->Display (anAISPlane, Standard_True);
1230       }
1231
1232       // The first argument is an AIS_Axis
1233       // Creation of a plane orthogonal to the axis through a point
1234     else if (aShapeA->Type()==AIS_KOI_Datum && aShapeA->Signature()==2 ) {
1235       // The second argument should be an AIS_Point
1236       Handle(AIS_InteractiveObject) aShapeB;
1237       if (argc!=4 || !GetMapOfAIS().Find2 (argv[3], aShapeB))
1238       {
1239         std::cout<<"vplane: error 2d name doesn't exist in the GetMapOfAIS()\n";
1240         return 1;
1241       }
1242       // If B is not an AIS_Point
1243       if (aShapeB.IsNull() ||
1244         (!(aShapeB->Type()==AIS_KOI_Datum && aShapeB->Signature()==1)))
1245       {
1246         std::cout<<"vplane: error 2d object is expected to be an AIS_Point\n";
1247         return 1;
1248       }
1249
1250       // Treatment of objects A and B
1251       Handle(AIS_Axis) anAISAxisA = Handle(AIS_Axis)::DownCast(aShapeA);
1252       Handle(AIS_Point) anAISPointB = Handle(AIS_Point)::DownCast(aShapeB);
1253
1254       Handle(Geom_Line ) aGeomLineA = anAISAxisA ->Component();
1255       Handle(Geom_Point) aGeomPointB = anAISPointB->Component()  ;
1256
1257       gp_Ax1 anAxis = aGeomLineA->Position();
1258       Handle(Geom_CartesianPoint) aCartPointB = 
1259         Handle(Geom_CartesianPoint)::DownCast(aGeomPointB);
1260
1261       gp_Dir D =anAxis.Direction();
1262       gp_Pnt B = aCartPointB->Pnt();
1263
1264       // Construction of AIS_Plane
1265       Handle(Geom_Plane) aGeomPlane = new Geom_Plane(B,D);
1266       Handle(AIS_Plane) anAISPlane = new AIS_Plane(aGeomPlane,B );
1267       GetMapOfAIS().Bind (anAISPlane,aName );
1268       if (argc == 5)
1269       {
1270         Standard_Integer aType = Draw::Atoi (argv[4]);
1271         if (aType != 0 && aType != 1)
1272         {
1273           std::cout << "vplane error: wrong type of sensitivity!\n"
1274                     << "Should be one of the following values:\n"
1275                     << "0 - Interior\n"
1276                     << "1 - Boundary"
1277                     << std::endl;
1278           return 1;
1279         }
1280         else
1281         {
1282           anAISPlane->SetTypeOfSensitivity (Select3D_TypeOfSensitivity (aType));
1283         }
1284       }
1285       TheAISContext()->Display (anAISPlane, Standard_True);
1286
1287     }
1288     // The first argumnet is an AIS_Plane
1289     // Creation of a plane parallel to the plane passing through the point
1290     else if (aShapeA->Type()==AIS_KOI_Datum && aShapeA->Signature()==7)
1291     {
1292       // The second argument should be an AIS_Point
1293       Handle(AIS_InteractiveObject) aShapeB;
1294       if (argc!=4 || !GetMapOfAIS().Find2 (argv[3], aShapeB))
1295       {
1296         std::cout<<"vplane: error 2d name doesn't exist in the GetMapOfAIS()\n";
1297         return 1;
1298       }
1299       // B should be an AIS_Point
1300       if (aShapeB.IsNull() ||
1301          (!(aShapeB->Type()==AIS_KOI_Datum && aShapeB->Signature()==1)))
1302       {
1303         std::cout<<"vplane: error 2d object is expected to be an AIS_Point\n";
1304         return 1;
1305       }
1306
1307       // Treatment of objects A and B
1308       Handle(AIS_Plane) anAISPlaneA = Handle(AIS_Plane)::DownCast(aShapeA);
1309       Handle(AIS_Point) anAISPointB = Handle(AIS_Point)::DownCast(aShapeB);
1310
1311       Handle(Geom_Plane) aNewGeomPlane= anAISPlaneA->Component();
1312       Handle(Geom_Point) aGeomPointB = anAISPointB->Component();
1313
1314       Handle(Geom_CartesianPoint) aCartPointB = 
1315         Handle(Geom_CartesianPoint)::DownCast(aGeomPointB);
1316       gp_Pnt B= aCartPointB->Pnt();
1317
1318       // Construction of an AIS_Plane
1319       Handle(AIS_Plane) anAISPlane = new AIS_Plane(aNewGeomPlane, B);
1320       GetMapOfAIS().Bind (anAISPlane, aName);
1321       if (argc == 5)
1322       {
1323         Standard_Integer aType = Draw::Atoi (argv[4]);
1324         if (aType != 0 && aType != 1)
1325         {
1326           std::cout << "vplane error: wrong type of sensitivity!\n"
1327                     << "Should be one of the following values:\n"
1328                     << "0 - Interior\n"
1329                     << "1 - Boundary"
1330                     << std::endl;
1331           return 1;
1332         }
1333         else
1334         {
1335           anAISPlane->SetTypeOfSensitivity (Select3D_TypeOfSensitivity (aType));
1336         }
1337       }
1338       TheAISContext()->Display (anAISPlane, Standard_True);
1339     }
1340     // Error
1341     else
1342     {
1343       std::cout<<"vplane: error 1st object is not an AIS\n";
1344       return 1;
1345     }
1346   }
1347   // There are no arguments
1348   else 
1349   {
1350     TopTools_ListOfShape aShapes;
1351     ViewerTest::GetSelectedShapes (aShapes);
1352
1353     // Function vplane
1354     // Test the constructor AIS_Plane::AIS_Plane(Geom_Plane, Standard_Boolean )
1355     if (!strcasecmp(argv[0], "vplane"))
1356     {
1357       if (aShapes.Extent() < 1 || aShapes.Extent() > 3)
1358       {
1359         std::cerr << "Error: Wront number of selected shapes.\n";
1360         std::cerr << "\tYou should one of variant: face, edge and vertex or three vertices.\n";
1361         return 1;
1362       }
1363
1364       const TopoDS_Shape& aShapeA = aShapes.First();
1365       if (aShapeA.ShapeType() == TopAbs_VERTEX)
1366       {
1367         if (aShapes.Extent() == 2)
1368         {
1369           const TopoDS_Shape& aShapeB = aShapes.Last();
1370           if (aShapeB.ShapeType() != TopAbs_EDGE)
1371           {
1372             std::cerr << "Syntax error: Together with vertex should be edge.\n";
1373             return 1;
1374           }
1375
1376           // Verify that the vertex is not on the edge ShapeB
1377           TopoDS_Edge anEdgeB = TopoDS::Edge(aShapeB);
1378           TopoDS_Vertex aVertA = TopoDS::Vertex(aShapeA);
1379
1380           BRepExtrema_ExtPC OrthoProj(aVertA, anEdgeB);
1381           if (OrthoProj.SquareDistance(1)<Precision::Approximation())
1382           {
1383             // The vertex is on the edge
1384             std::cout<<" vplane: error point is on the edge\n";
1385             return 1;
1386           }
1387           else
1388           {
1389             gp_Pnt A = BRep_Tool::Pnt(aVertA);
1390             TopoDS_Vertex aVBa, aVBb;
1391             TopExp::Vertices(anEdgeB ,aVBa ,aVBb);
1392             gp_Pnt aBa = BRep_Tool::Pnt(aVBa);
1393             gp_Pnt aBb = BRep_Tool::Pnt(aVBb);
1394             GC_MakePlane MkPlane (A, aBa, aBb);
1395             Handle(Geom_Plane) aGeomPlane = MkPlane.Value();
1396             Handle(AIS_Plane) anAISPlane = new AIS_Plane (aGeomPlane);
1397             GetMapOfAIS().Bind (anAISPlane, aName);
1398             TheAISContext()->Display (anAISPlane, Standard_True);
1399           }
1400         }
1401         else if (aShapes.Extent() == 3)
1402         {
1403           TopTools_ListOfShape::Iterator anIter (aShapes);
1404
1405           anIter.Next();
1406           const TopoDS_Shape& aShapeB = anIter.Value();
1407
1408           anIter.Next();
1409           const TopoDS_Shape& aShapeC = anIter.Value();
1410
1411           if (!(aShapeB.ShapeType() == TopAbs_VERTEX
1412              && aShapeC.ShapeType() == TopAbs_VERTEX))
1413           {
1414             std::cerr << "Syntax error: You should one of variant: face, edge and vertex or three vertices.\n";
1415             return 1;
1416           }
1417
1418           gp_Pnt A = BRep_Tool::Pnt(TopoDS::Vertex(aShapeA));
1419           gp_Pnt B = BRep_Tool::Pnt(TopoDS::Vertex(aShapeB));
1420           gp_Pnt C = BRep_Tool::Pnt(TopoDS::Vertex(aShapeC));
1421           GC_MakePlane MkPlane(A, B, C);
1422           Handle(Geom_Plane) aGeomPlane = MkPlane.Value();
1423           Handle(AIS_Plane) anAISPlane = new AIS_Plane (aGeomPlane);
1424           GetMapOfAIS().Bind (anAISPlane, aName);
1425           TheAISContext()->Display (anAISPlane, Standard_True);
1426         }
1427         else
1428         {
1429           std::cerr << "Syntax error: You should one of variant: face, edge and vertex or three vertices.\n";
1430           return 1;
1431         }
1432       }
1433       else if (aShapeA.ShapeType() == TopAbs_EDGE)
1434       {
1435         if (aShapes.Extent() != 2)
1436         {
1437           std::cerr << "Error: wrong number of selected shapes.\n";
1438           return 1;
1439         }
1440
1441         const TopoDS_Shape& aShapeB = aShapes.Last();
1442         if (aShapeB.ShapeType() != TopAbs_VERTEX)
1443         {
1444           std::cerr << "Syntax error: Together with edge should be vertex.\n";
1445           return 1;
1446         }
1447
1448         // Check that the vertex aShapeB is not on the edge
1449         TopoDS_Edge anEdgeA = TopoDS::Edge(aShapeA);
1450         TopoDS_Vertex aVertB = TopoDS::Vertex(aShapeB);
1451
1452         BRepExtrema_ExtPC OrthoProj (aVertB, anEdgeA);
1453         if (OrthoProj.SquareDistance(1)<Precision::Approximation())
1454         {
1455           // The vertex is on the edge
1456           std::cout<<" vplane: error point is on the edge\n";
1457           return 1;
1458         }
1459
1460         gp_Pnt B = BRep_Tool::Pnt(aVertB);
1461         TopoDS_Vertex aVAa, aVAb;
1462         TopExp::Vertices(anEdgeA, aVAa, aVAb);
1463         gp_Pnt Aa = BRep_Tool::Pnt(aVAa);
1464         gp_Pnt Ab = BRep_Tool::Pnt(aVAb);
1465         GC_MakePlane MkPlane (B,Aa,Ab);
1466         Handle(Geom_Plane) aGeomPlane = MkPlane.Value();
1467         Handle(AIS_Plane) anAISPlane = new AIS_Plane (aGeomPlane);
1468         GetMapOfAIS().Bind (anAISPlane ,aName);
1469         TheAISContext()->Display (anAISPlane, Standard_True);
1470       }
1471       else if (aShapeA.ShapeType() == TopAbs_FACE)
1472       {
1473         TopoDS_Face aFace = TopoDS::Face(aShapeA);
1474         BRepAdaptor_Surface aSurface (aFace, Standard_False);
1475         if (aSurface.GetType()==GeomAbs_Plane)
1476         {
1477           gp_Pln aPlane = aSurface.Plane();
1478           Handle(Geom_Plane) aGeomPlane = new Geom_Plane(aPlane);
1479           Handle(AIS_Plane) anAISPlane = new AIS_Plane(aGeomPlane);
1480           GetMapOfAIS().Bind (anAISPlane, aName);
1481           TheAISContext()->Display (anAISPlane, Standard_True);
1482         }
1483         else
1484         {
1485           std::cout<<" vplane: error\n";
1486           return 1;
1487         }
1488       }
1489       else
1490       {
1491         std::cerr << "Syntax error: You should one of variant: face, edge and vertex or three vertices.\n";
1492         return 1;
1493       }
1494     }
1495
1496     // Function vPlanePara
1497     // ===================
1498     // test the constructor AIS_Plane::AIS_Plane(Geom_Plane,gp_Pnt)
1499     else if (!strcasecmp(argv[0], "vplanepara"))
1500     {
1501       if (aShapes.Extent() != 2)
1502       {
1503         std::cerr << "Error: Wrong number of selected shapes.\n";
1504         return 1;
1505       }
1506
1507       const TopoDS_Shape* aShapeA = &aShapes.First();
1508       const TopoDS_Shape* aShapeB = &aShapes.Last();
1509       if (aShapeA->ShapeType() != TopAbs_VERTEX)
1510       {
1511         std::swap (aShapeA, aShapeB);
1512       }
1513
1514       if (!(aShapeA->ShapeType() == TopAbs_VERTEX
1515          && aShapeB->ShapeType() == TopAbs_FACE))
1516       {
1517         std::cerr << "Syntax error: you should select face and vertex.\n";
1518         return 1;
1519       }
1520
1521       gp_Pnt A = BRep_Tool::Pnt(TopoDS::Vertex(*aShapeA));
1522
1523       TopoDS_Face aFace = TopoDS::Face(*aShapeB);
1524       BRepAdaptor_Surface aSurface (aFace, Standard_False);
1525       if (aSurface.GetType() == GeomAbs_Plane)
1526       {
1527         gp_Pln aPlane = aSurface.Plane();
1528         // Construct a plane parallel to aGeomPlane through A
1529         aPlane.SetLocation(A);
1530         Handle(Geom_Plane) aGeomPlane = new Geom_Plane (aPlane);
1531         Handle(AIS_Plane) aAISPlane = new AIS_Plane (aGeomPlane, A);
1532         GetMapOfAIS().Bind (aAISPlane ,aName);
1533         TheAISContext()->Display (aAISPlane, Standard_True);
1534       }
1535       else
1536       {
1537         std::cerr << "Error: Builded surface is not a plane.\n";
1538         return 1;
1539       }
1540     }
1541
1542     // Function vplaneortho
1543     // ====================
1544     // test the constructor AIS_Plane::AIS_Plane(Geom_Plane,gp_Pnt,gp_Pnt,gp_Pnt)
1545     else
1546     {
1547       if (aShapes.Extent() != 2)
1548       {
1549         std::cerr << "Error: wrong number of selected shapes.\n";
1550         return 1;
1551       }
1552
1553       const TopoDS_Shape* aShapeA = &aShapes.First();
1554       const TopoDS_Shape* aShapeB = &aShapes.Last();
1555
1556       if (aShapeA->ShapeType() != TopAbs_EDGE)
1557       {
1558         std::swap (aShapeA, aShapeB);
1559       }
1560
1561       if (!(aShapeA->ShapeType() == TopAbs_EDGE
1562          && aShapeB->ShapeType() == TopAbs_FACE))
1563       {
1564         std::cerr << "Error: you should select edge and face.\n";
1565         return 1;
1566       }
1567
1568       // Construction of plane
1569       TopoDS_Edge anEdgeA = TopoDS::Edge(*aShapeA);
1570       TopoDS_Vertex aVAa, aVAb;
1571       TopExp::Vertices(anEdgeA, aVAa, aVAb);
1572       gp_Pnt Aa = BRep_Tool::Pnt(aVAa);
1573       gp_Pnt Ab = BRep_Tool::Pnt(aVAb);
1574       gp_Vec ab (Aa,Ab);
1575
1576       gp_Dir Dab (ab);
1577       // Creation of rotation axis
1578       gp_Ax1 aRotAxis (Aa,Dab);
1579
1580       TopoDS_Face aFace = TopoDS::Face(*aShapeB);
1581       // The edge must be parallel to the face
1582       BRepExtrema_ExtPF aHeightA (aVAa, aFace);
1583       BRepExtrema_ExtPF aHeightB (aVAb, aFace);
1584       // Compare to heights
1585       if (fabs(sqrt(aHeightA.SquareDistance(1)) - sqrt(aHeightB.SquareDistance(1)))
1586           >Precision::Confusion())
1587       {
1588         // the edge is not parallel to the face
1589         std::cout<<" vplaneortho error: the edge is not parallel to the face\n";
1590         return 1;
1591       }
1592       // the edge is OK
1593       BRepAdaptor_Surface aSurface (aFace, Standard_False);
1594       if (aSurface.GetType()==GeomAbs_Plane)
1595       {
1596         gp_Pln aPlane = aSurface.Plane();
1597         // It rotates a half turn round the axis of rotation
1598         aPlane.Rotate(aRotAxis , M_PI/2);
1599
1600         Handle(Geom_Plane) aGeomPlane = new Geom_Plane (aPlane);
1601         // constructed aGeomPlane parallel to a plane containing the edge (center mid-edge)
1602         gp_Pnt aMiddle ((Aa.X()+Ab.X() )/2 ,(Aa.Y()+Ab.Y() )/2 ,(Aa.Z()+Ab.Z() )/2 );
1603         Handle(AIS_Plane) anAISPlane = new AIS_Plane (aGeomPlane, aMiddle);
1604         GetMapOfAIS().Bind (anAISPlane, aName);
1605         TheAISContext()->Display (anAISPlane, Standard_True);
1606       }
1607       else
1608       {
1609         std::cout<<" vplaneortho: error\n";
1610         return 1;
1611       }
1612     }
1613   }
1614   return 0;
1615 }
1616
1617 //===============================================================================================
1618 //function : VChangePlane
1619 //purpose  :
1620 //===============================================================================================
1621 static int VChangePlane (Draw_Interpretor& /*theDi*/, Standard_Integer theArgsNb, const char** theArgVec)
1622 {
1623   Handle(AIS_InteractiveContext) aContextAIS = ViewerTest::GetAISContext();
1624   if (aContextAIS.IsNull())
1625   {
1626     std::cout << theArgVec[0] << "AIS context is not available.\n";
1627     return 1;
1628   }
1629
1630   if (theArgsNb < 3 || theArgsNb > 11)
1631   {
1632     std::cerr << theArgVec[0] 
1633               << ": incorrect number of command arguments.\n"
1634               << "Type help for more information.\n";
1635     return 1;
1636   }
1637
1638   TCollection_AsciiString aName (theArgVec[1]);
1639
1640   Handle(AIS_Plane) aPlane = GetMapOfAIS().IsBound2(aName)
1641     ? Handle(AIS_Plane)::DownCast (GetMapOfAIS().Find2 (aName))
1642     : NULL;
1643
1644   if ( aPlane.IsNull() )
1645   {
1646     std::cout << theArgVec[0] 
1647               << ": there is no interactive plane with the given name."
1648               << "Type help for more information.\n";
1649     return 1;
1650   }
1651
1652   Standard_Real aCenterX = aPlane->Center().X();
1653   Standard_Real aCenterY = aPlane->Center().Y();
1654   Standard_Real aCenterZ = aPlane->Center().Z();
1655
1656   Standard_Real aDirX = aPlane->Component()->Axis().Direction().X();
1657   Standard_Real aDirY = aPlane->Component()->Axis().Direction().Y();
1658   Standard_Real aDirZ = aPlane->Component()->Axis().Direction().Z();
1659
1660   Standard_Real aSizeX = 0.0;
1661   Standard_Real aSizeY = 0.0;
1662   aPlane->Size (aSizeX, aSizeY);
1663   Standard_Boolean isUpdate = Standard_True;
1664
1665   TCollection_AsciiString aPName, aPValue;
1666   for (Standard_Integer anArgIt = 1; anArgIt < theArgsNb; ++anArgIt)
1667   {
1668     const TCollection_AsciiString anArg = theArgVec[anArgIt];
1669     TCollection_AsciiString anArgCase = anArg;
1670     anArgCase.UpperCase();
1671     if (ViewerTest::SplitParameter (anArg, aPName, aPValue))
1672     {
1673       aPName.UpperCase();
1674       if (aPName.IsEqual ("X"))
1675       {
1676         aCenterX = aPValue.RealValue();
1677       }
1678       else if (aPName.IsEqual ("Y"))
1679       {
1680         aCenterY = aPValue.RealValue();
1681       }
1682       else if (aPName.IsEqual ("Z"))
1683       {
1684         aCenterZ = aPValue.RealValue();
1685       }
1686       else if (aPName.IsEqual ("DX"))
1687       {
1688         aDirX = aPValue.RealValue();
1689       }
1690       else if (aPName.IsEqual ("DY"))
1691       {
1692         aDirY = aPValue.RealValue();
1693       }
1694       else if (aPName.IsEqual ("DZ"))
1695       {
1696         aDirZ = aPValue.RealValue();
1697       }
1698       else if (aPName.IsEqual ("SX"))
1699       {
1700         aSizeX = aPValue.RealValue();
1701       }
1702       else if (aPName.IsEqual ("SY"))
1703       {
1704         aSizeY = aPValue.RealValue();
1705       }
1706     }
1707     else if (anArg.IsEqual ("NOUPDATE"))
1708     {
1709       isUpdate = Standard_False;
1710     }
1711   }
1712
1713   gp_Dir aDirection (aDirX, aDirY, aDirZ);
1714   gp_Pnt aCenterPnt (aCenterX, aCenterY, aCenterZ);
1715   aPlane->SetCenter (aCenterPnt);
1716   aPlane->SetComponent (new Geom_Plane (aCenterPnt, aDirection));
1717   aPlane->SetSize (aSizeX, aSizeY);
1718
1719   aContextAIS->Update (aPlane, isUpdate);
1720
1721   return 0;
1722 }
1723
1724 //==============================================================================
1725 // Fonction  vline
1726 // ---------------  Uniquement par parametre. Pas de selection dans le viewer.
1727 //==============================================================================
1728
1729 //==============================================================================
1730 //function : VLineBuilder
1731 //purpose  : Build an AIS_Line
1732 //Draw arg : vline LineName  [AIS_PointName] [AIS_PointName]
1733 //                           [Xa] [Ya] [Za]   [Xb] [Yb] [Zb]
1734 //==============================================================================
1735 #include <Geom_CartesianPoint.hxx>
1736 #include <AIS_Line.hxx>
1737
1738
1739 static int VLineBuilder(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
1740 {
1741   // Verifications
1742   if (argc!=4 && argc!=8 && argc!=2 )  {di<<"vline error: number of arguments not correct \n";return 1; }
1743
1744   // On recupere les parametres
1745   Handle(AIS_InteractiveObject) theShapeA;
1746   Handle(AIS_InteractiveObject) theShapeB;
1747
1748   // Parametres: AIS_Point AIS_Point
1749   // ===============================
1750   if (argc==4) {
1751     GetMapOfAIS().Find2 (argv[2], theShapeA);
1752     // On verifie que c'est bien une AIS_Point
1753     if (!theShapeA.IsNull() &&
1754       theShapeA->Type()==AIS_KOI_Datum && theShapeA->Signature()==1) {
1755         // on recupere le deuxieme AIS_Point
1756         GetMapOfAIS().Find2 (argv[3], theShapeB);
1757         if (theShapeB.IsNull() ||
1758           (!(theShapeB->Type()==AIS_KOI_Datum && theShapeB->Signature()==1)))
1759         {
1760           di <<"vline error: wrong type of 2de argument.\n";
1761           return 1;
1762         }
1763       }
1764     else {di <<"vline error: wrong type of 1st argument.\n";return 1; }
1765     // Les deux parametres sont du bon type. On verifie que les points ne sont pas confondus
1766     Handle(AIS_Point) theAISPointA= Handle(AIS_Point)::DownCast (theShapeA);
1767     Handle(AIS_Point) theAISPointB= Handle(AIS_Point)::DownCast (theShapeB);
1768
1769     Handle(Geom_Point ) myGeomPointBA=  theAISPointA->Component();
1770     Handle(Geom_CartesianPoint ) myCartPointA= Handle(Geom_CartesianPoint)::DownCast (myGeomPointBA);
1771     //    Handle(Geom_CartesianPoint ) myCartPointA= *(Handle(Geom_CartesianPoint)*)& (theAISPointA->Component() ) ;
1772
1773     Handle(Geom_Point ) myGeomPointB=  theAISPointB->Component();
1774     Handle(Geom_CartesianPoint ) myCartPointB= Handle(Geom_CartesianPoint)::DownCast (myGeomPointB);
1775     //    Handle(Geom_CartesianPoint ) myCartPointB= *(Handle(Geom_CartesianPoint)*)& (theAISPointB->Component() ) ;
1776
1777     if (myCartPointB->X()==myCartPointA->X() && myCartPointB->Y()==myCartPointA->Y() && myCartPointB->Z()==myCartPointA->Z() ) {
1778       // B=A
1779       di<<"vline error: same points\n";return 1;
1780     }
1781     // Les deux points sont OK...Construction de l'AIS_Line (en faite, le segment AB)
1782     Handle(AIS_Line) theAISLine= new AIS_Line(myCartPointA,myCartPointB );
1783     GetMapOfAIS().Bind(theAISLine,argv[1] );
1784     TheAISContext()->Display (theAISLine, Standard_True);
1785
1786   }
1787
1788   // Parametres 6 Reals
1789   // ==================
1790
1791   else if (argc==8) {
1792     // On verifie que les deux points ne sont pas confondus
1793
1794     Standard_Real coord[6];
1795     for(Standard_Integer i=0;i<=2;i++){
1796       coord[i]=Draw::Atof(argv[2+i]);
1797       coord[i+3]=Draw::Atof(argv[5+i]);
1798     }
1799
1800     Handle(Geom_CartesianPoint ) myCartPointA=new Geom_CartesianPoint (coord[0],coord[1],coord[2] );
1801     Handle(Geom_CartesianPoint ) myCartPointB=new Geom_CartesianPoint (coord[3],coord[4],coord[5] );
1802
1803     Handle(AIS_Line) theAISLine= new AIS_Line(myCartPointA,myCartPointB );
1804     GetMapOfAIS().Bind(theAISLine,argv[1] );
1805     TheAISContext()->Display (theAISLine, Standard_True);
1806
1807   }
1808
1809   // Pas de parametres: Selection dans le viewer.
1810   // ============================================
1811
1812   else
1813   {
1814     TopTools_ListOfShape aShapes;
1815     ViewerTest::GetSelectedShapes (aShapes);
1816     if (aShapes.Extent() != 2)
1817     {
1818       std::cerr << "Error: wrong number of selected shapes.\n";
1819       return 1;
1820     }
1821
1822     const TopoDS_Shape& aShapeA = aShapes.First();
1823     const TopoDS_Shape& aShapeB = aShapes.Last();
1824
1825     if (!(aShapeA.ShapeType() == TopAbs_VERTEX
1826        && aShapeB.ShapeType() == TopAbs_VERTEX))
1827     {
1828       std::cerr << "Error: you should select two different vertex.\n";
1829       return 1;
1830     }
1831
1832     // Construction de la line
1833     gp_Pnt A = BRep_Tool::Pnt (TopoDS::Vertex (aShapeA));
1834     gp_Pnt B = BRep_Tool::Pnt (TopoDS::Vertex (aShapeB));
1835
1836     Handle(Geom_CartesianPoint ) myCartPointA=new Geom_CartesianPoint(A);
1837     Handle(Geom_CartesianPoint ) myCartPointB=new Geom_CartesianPoint(B);
1838
1839     Handle(AIS_Line) theAISLine= new AIS_Line(myCartPointA,myCartPointB );
1840     GetMapOfAIS().Bind(theAISLine,argv[1] );
1841     TheAISContext()->Display (theAISLine, Standard_True);
1842   }
1843
1844   return 0;
1845 }
1846
1847 //==============================================================================
1848 // class   : FilledCircle
1849 // purpose : creates filled circle based on AIS_InteractiveObject 
1850 //           and Geom_Circle.
1851 //           This class is used to check method Matches() of class 
1852 //           Select3D_SensitiveCircle with member myFillStatus = Standard_True, 
1853 //           because none of AIS classes provides creation of 
1854 //           Select3D_SensitiveCircle with member myFillStatus = Standard_True 
1855 //           (look method ComputeSelection() )
1856 //============================================================================== 
1857
1858 Handle(Geom_Circle) CreateCircle(gp_Pnt theCenter, Standard_Real theRadius) 
1859 {
1860   gp_Ax2 anAxes(theCenter, gp_Dir(gp_Vec(0., 0., 1.))); 
1861   gp_Circ aCirc(anAxes, theRadius);
1862   Handle(Geom_Circle) aCircle = new Geom_Circle(aCirc);
1863   return aCircle;
1864 }
1865
1866 class FilledCircle : public AIS_InteractiveObject 
1867 {
1868 public:
1869     // CASCADE RTTI
1870     DEFINE_STANDARD_RTTI_INLINE(FilledCircle,AIS_InteractiveObject); 
1871
1872     FilledCircle(gp_Pnt theCenter, Standard_Real theRadius);
1873     FilledCircle(Handle(Geom_Circle) theCircle);
1874
1875 private:
1876     TopoDS_Face ComputeFace();
1877
1878     // Virtual methods implementation
1879     void Compute (  const Handle(PrsMgr_PresentationManager3d)& thePresentationManager,
1880                   const Handle(Prs3d_Presentation)& thePresentation,
1881                   const Standard_Integer theMode) Standard_OVERRIDE;
1882
1883     void ComputeSelection (  const Handle(SelectMgr_Selection)& theSelection, 
1884                            const Standard_Integer theMode) Standard_OVERRIDE;
1885
1886 protected:
1887     Handle(Geom_Circle) myCircle;
1888     Standard_Boolean myFilledStatus;
1889
1890 }; 
1891
1892
1893 FilledCircle::FilledCircle(gp_Pnt theCenter, Standard_Real theRadius) 
1894 {
1895   myCircle = CreateCircle(theCenter, theRadius);
1896   myFilledStatus = Standard_True;
1897 }
1898
1899 FilledCircle::FilledCircle(Handle(Geom_Circle) theCircle) 
1900 {
1901   myCircle = theCircle;
1902   myFilledStatus = Standard_True;
1903 }
1904
1905 TopoDS_Face FilledCircle::ComputeFace() 
1906 {
1907   // Create edge from myCircle 
1908   BRepBuilderAPI_MakeEdge anEdgeMaker(myCircle->Circ());
1909   TopoDS_Edge anEdge = anEdgeMaker.Edge(); 
1910
1911   // Create wire from anEdge 
1912   BRepBuilderAPI_MakeWire aWireMaker(anEdge);
1913   TopoDS_Wire aWire = aWireMaker.Wire();
1914
1915   // Create face from aWire
1916   BRepBuilderAPI_MakeFace aFaceMaker(aWire);
1917   TopoDS_Face aFace = aFaceMaker.Face();
1918
1919   return aFace;
1920 }
1921
1922 void FilledCircle::Compute(const Handle(PrsMgr_PresentationManager3d) &/*thePresentationManager*/, 
1923                            const Handle(Prs3d_Presentation) &thePresentation, 
1924                            const Standard_Integer theMode) 
1925 {
1926   thePresentation->Clear();
1927
1928   TopoDS_Face aFace = ComputeFace();
1929
1930   if (aFace.IsNull()) return;
1931   if (theMode != 0) return;
1932
1933   StdPrs_ShadedShape::Add(thePresentation, aFace, myDrawer);
1934 }
1935
1936 void FilledCircle::ComputeSelection(const Handle(SelectMgr_Selection) &theSelection, 
1937                                     const Standard_Integer /*theMode*/)
1938 {
1939   Handle(SelectMgr_EntityOwner) anEntityOwner = new SelectMgr_EntityOwner(this);
1940   Handle(Select3D_SensitiveCircle) aSensitiveCircle = new Select3D_SensitiveCircle(anEntityOwner, 
1941       myCircle, myFilledStatus);
1942   theSelection->Add(aSensitiveCircle);
1943 }
1944
1945 //==============================================================================
1946 // Fonction  vcircle
1947 // -----------------  Uniquement par parametre. Pas de selection dans le viewer.
1948 //==============================================================================
1949
1950 //==============================================================================
1951 //function : VCircleBuilder
1952 //purpose  : Build an AIS_Circle
1953 //Draw arg : vcircle CircleName PlaneName PointName Radius IsFilled
1954 //                              PointName PointName PointName IsFilled
1955 //==============================================================================
1956
1957 void DisplayCircle (Handle (Geom_Circle) theGeomCircle,
1958                     TCollection_AsciiString theName, 
1959                     Standard_Boolean isFilled) 
1960 {
1961   Handle(AIS_InteractiveObject) aCircle;
1962   if (isFilled) 
1963   {
1964     aCircle = new FilledCircle(theGeomCircle);
1965   }
1966   else
1967   {
1968     aCircle = new AIS_Circle(theGeomCircle);
1969     Handle(AIS_Circle)::DownCast (aCircle)->SetFilledCircleSens (Standard_False);
1970   }
1971
1972   // Check if there is an object with given name
1973   // and remove it from context
1974   if (GetMapOfAIS().IsBound2(theName)) 
1975   {
1976     Handle(AIS_InteractiveObject) anInterObj = GetMapOfAIS().Find2(theName);
1977     TheAISContext()->Remove(anInterObj, Standard_False);
1978     GetMapOfAIS().UnBind2(theName);
1979    }
1980
1981    // Bind the circle to its name
1982    GetMapOfAIS().Bind(aCircle, theName);
1983
1984    // Display the circle
1985    TheAISContext()->Display (aCircle, Standard_True);
1986   
1987 }
1988
1989 static int VCircleBuilder(Draw_Interpretor& /*di*/, Standard_Integer argc, const char** argv)
1990 {
1991   // Verification of the arguments
1992   if (argc>6 || argc<2) 
1993   { 
1994     std::cout << "vcircle error: expect 4 arguments.\n"; 
1995     return 1; // TCL_ERROR 
1996   }
1997
1998   // There are all arguments
1999   if (argc == 6) 
2000   {
2001     // Get arguments
2002     TCollection_AsciiString aName(argv[1]);
2003     Standard_Boolean isFilled = Draw::Atoi(argv[5]) != 0;
2004
2005     Handle(AIS_InteractiveObject) theShapeA, theShapeB;
2006     GetMapOfAIS().Find2 (argv[2], theShapeA);
2007     GetMapOfAIS().Find2 (argv[3], theShapeB);
2008
2009     // Arguments: AIS_Point AIS_Point AIS_Point
2010     // ========================================
2011     if (!theShapeA.IsNull() && !theShapeB.IsNull() &&
2012       theShapeA->Type()==AIS_KOI_Datum && theShapeA->Signature()==1)
2013     {
2014       if (theShapeB->Type()!=AIS_KOI_Datum || theShapeB->Signature()!=1 ) 
2015       {
2016         std::cout << "vcircle error: 2d argument is unexpected to be a point.\n";
2017         return 1; // TCL_ERROR 
2018       }
2019       // The third object must be a point
2020       Handle(AIS_InteractiveObject) theShapeC;
2021       GetMapOfAIS().Find2 (argv[4], theShapeC);
2022       if (theShapeC.IsNull() ||
2023         theShapeC->Type()!=AIS_KOI_Datum || theShapeC->Signature()!=1 ) 
2024       {
2025         std::cout << "vcircle error: 3d argument is unexpected to be a point.\n";
2026         return 1; // TCL_ERROR 
2027       }
2028         // tag
2029         // Verify that the three points are different
2030         Handle(AIS_Point) theAISPointA = Handle(AIS_Point)::DownCast(theShapeA);
2031         Handle(AIS_Point) theAISPointB = Handle(AIS_Point)::DownCast(theShapeB);
2032         Handle(AIS_Point) theAISPointC = Handle(AIS_Point)::DownCast(theShapeC);
2033         
2034         Handle(Geom_Point) myGeomPointA = theAISPointA->Component();
2035         Handle(Geom_CartesianPoint) myCartPointA = 
2036           Handle(Geom_CartesianPoint)::DownCast(myGeomPointA);
2037
2038         Handle(Geom_Point) myGeomPointB = theAISPointB->Component();
2039         Handle(Geom_CartesianPoint) myCartPointB =
2040           Handle(Geom_CartesianPoint)::DownCast(myGeomPointB);
2041
2042         Handle(Geom_Point) myGeomPointC = theAISPointC->Component();
2043         Handle(Geom_CartesianPoint) myCartPointC =
2044           Handle(Geom_CartesianPoint)::DownCast(myGeomPointC);
2045
2046         // Test A=B
2047         if (Abs(myCartPointA->X()-myCartPointB->X()) <= Precision::Confusion() && 
2048             Abs(myCartPointA->Y()-myCartPointB->Y()) <= Precision::Confusion() && 
2049             Abs(myCartPointA->Z()-myCartPointB->Z()) <= Precision::Confusion() ) 
2050         {
2051           std::cout << "vcircle error: Same points.\n"; 
2052           return 1; // TCL_ERROR 
2053         }
2054         // Test A=C
2055         if (Abs(myCartPointA->X()-myCartPointC->X()) <= Precision::Confusion() &&
2056             Abs(myCartPointA->Y()-myCartPointC->Y()) <= Precision::Confusion() && 
2057             Abs(myCartPointA->Z()-myCartPointC->Z()) <= Precision::Confusion() ) 
2058         {
2059           std::cout << "vcircle error: Same points.\n"; 
2060           return 1; // TCL_ERROR 
2061         }
2062         // Test B=C
2063         if (Abs(myCartPointB->X()-myCartPointC->X()) <= Precision::Confusion() && 
2064             Abs(myCartPointB->Y()-myCartPointC->Y()) <= Precision::Confusion() && 
2065             Abs(myCartPointB->Z()-myCartPointC->Z()) <= Precision::Confusion() ) 
2066         {
2067           std::cout << "vcircle error: Same points.\n"; 
2068           return 1;// TCL_ERROR 
2069         }
2070         // Construction of the circle
2071         GC_MakeCircle Cir = GC_MakeCircle (myCartPointA->Pnt(), 
2072           myCartPointB->Pnt(), myCartPointC->Pnt() );
2073         Handle (Geom_Circle) theGeomCircle;
2074         try 
2075         {
2076           theGeomCircle = Cir.Value();
2077         }
2078         catch (StdFail_NotDone)
2079         {
2080           std::cout << "vcircle error: can't create circle\n";
2081           return -1; // TCL_ERROR
2082         }
2083         
2084         DisplayCircle(theGeomCircle, aName, isFilled);
2085     }
2086
2087     // Arguments: AIS_Plane AIS_Point Real
2088     // ===================================
2089     else if (theShapeA->Type() == AIS_KOI_Datum && 
2090       theShapeA->Signature() == 7 ) 
2091     {
2092       if (theShapeB->Type() != AIS_KOI_Datum || 
2093         theShapeB->Signature() != 1 ) 
2094       {
2095         std::cout << "vcircle error: 2d element is a unexpected to be a point.\n"; 
2096         return 1; // TCL_ERROR 
2097       }
2098       // Check that the radius is >= 0
2099       if (Draw::Atof(argv[4]) <= 0 ) 
2100       {
2101         std::cout << "vcircle error: the radius must be >=0.\n"; 
2102         return 1; // TCL_ERROR 
2103       }
2104
2105       // Recover the normal to the plane
2106       Handle(AIS_Plane) theAISPlane = Handle(AIS_Plane)::DownCast(theShapeA);
2107       Handle(AIS_Point) theAISPointB = Handle(AIS_Point)::DownCast(theShapeB); 
2108
2109       Handle(Geom_Plane) myGeomPlane = theAISPlane->Component();
2110       Handle(Geom_Point) myGeomPointB = theAISPointB->Component();
2111       Handle(Geom_CartesianPoint) myCartPointB = 
2112         Handle(Geom_CartesianPoint)::DownCast(myGeomPointB);
2113
2114       gp_Pln mygpPlane = myGeomPlane->Pln();
2115       gp_Ax1 thegpAxe = mygpPlane.Axis();
2116       gp_Dir theDir = thegpAxe.Direction();
2117       gp_Pnt theCenter = myCartPointB->Pnt();
2118       Standard_Real TheR = Draw::Atof(argv[4]);
2119       GC_MakeCircle Cir = GC_MakeCircle (theCenter, theDir ,TheR);
2120       Handle (Geom_Circle) theGeomCircle;
2121       try 
2122       {
2123         theGeomCircle = Cir.Value();
2124       }
2125       catch (StdFail_NotDone)
2126       {
2127         std::cout << "vcircle error: can't create circle\n";
2128         return -1; // TCL_ERROR
2129       }
2130
2131       DisplayCircle(theGeomCircle, aName, isFilled);
2132
2133     }
2134
2135     // Error
2136     else
2137     {
2138       std::cout << "vcircle error: 1st argument is a unexpected type.\n"; 
2139       return 1; // TCL_ERROR 
2140     }
2141
2142   }
2143   // No arguments: selection in the viewer
2144   // =========================================
2145   else 
2146   {
2147     // Get the name of the circle 
2148     TCollection_AsciiString aName(argv[1]);
2149
2150     TopTools_ListOfShape aShapes;
2151     ViewerTest::GetSelectedShapes (aShapes);
2152     if (aShapes.Extent() != 3 && aShapes.Extent() != 2)
2153     {
2154       std::cerr << "Error: Wrong number of selected shapes.\n";
2155       return 1;
2156     }
2157
2158     const TopoDS_Shape& aShapeA = aShapes.First();
2159     if (aShapeA.ShapeType() == TopAbs_VERTEX ) 
2160     {
2161       if (aShapes.Extent() != 3)
2162       {
2163         std::cerr << "Error: wrong number of selected shapes.\n";
2164         return 1;
2165       }
2166
2167       TopTools_ListOfShape::Iterator anIter (aShapes);
2168
2169       anIter.Next();
2170       const TopoDS_Shape& aShapeB = anIter.Value();
2171
2172       anIter.Next();
2173       const TopoDS_Shape& aShapeC = anIter.Value();
2174       
2175       // Get isFilled
2176       Standard_Boolean isFilled;
2177       std::cout << "Enter filled status (0 or 1)\n";
2178       cin >> isFilled;
2179
2180       // Construction of the circle
2181       gp_Pnt A = BRep_Tool::Pnt (TopoDS::Vertex (aShapeA));
2182       gp_Pnt B = BRep_Tool::Pnt (TopoDS::Vertex (aShapeB));
2183       gp_Pnt C = BRep_Tool::Pnt (TopoDS::Vertex (aShapeC));
2184
2185       GC_MakeCircle Cir = GC_MakeCircle (A, B, C);
2186       Handle (Geom_Circle) theGeomCircle;
2187       try 
2188       {
2189         theGeomCircle = Cir.Value();
2190       }
2191       catch (StdFail_NotDone)
2192       {
2193         std::cout << "vcircle error: can't create circle\n";
2194         return -1; // TCL_ERROR
2195       }
2196
2197       DisplayCircle(theGeomCircle, aName, isFilled);
2198
2199     }
2200     else if (aShapeA.ShapeType() == TopAbs_FACE)
2201     {
2202       const TopoDS_Shape& aShapeB = aShapes.Last();
2203
2204       // Recover the radius 
2205       Standard_Real theRad;
2206       do 
2207       {
2208         std::cout << " Enter the value of the radius:\n";
2209         cin >> theRad;
2210       } while (theRad <= 0);
2211       
2212       // Get filled status
2213       Standard_Boolean isFilled;
2214       std::cout << "Enter filled status (0 or 1)\n";
2215       cin >> isFilled;
2216
2217       // Recover the normal to the plane. tag
2218       TopoDS_Face myFace = TopoDS::Face(aShapeA);
2219       BRepAdaptor_Surface mySurface (myFace, Standard_False);
2220       gp_Pln myPlane = mySurface.Plane();
2221       Handle(Geom_Plane) theGeomPlane = new Geom_Plane (myPlane);
2222       gp_Pln mygpPlane = theGeomPlane->Pln();
2223       gp_Ax1 thegpAxe = mygpPlane.Axis();
2224       gp_Dir theDir = thegpAxe.Direction();
2225
2226       // Recover the center
2227       gp_Pnt theCenter = BRep_Tool::Pnt (TopoDS::Vertex (aShapeB));
2228
2229       // Construct the circle
2230       GC_MakeCircle Cir = GC_MakeCircle (theCenter, theDir ,theRad);
2231       Handle (Geom_Circle) theGeomCircle;
2232       try 
2233       {
2234         theGeomCircle = Cir.Value();
2235       }
2236       catch (StdFail_NotDone)
2237       {
2238         std::cout << "vcircle error: can't create circle\n";
2239         return -1; // TCL_ERROR
2240       }
2241
2242       DisplayCircle(theGeomCircle, aName, isFilled);
2243     }
2244     else
2245     {
2246       std::cerr << "Error: You should select face and vertex or three vertices.\n";
2247       return 1;
2248     }
2249   }
2250
2251   return 0;
2252 }
2253
2254 //=======================================================================
2255 //function : VDrawText
2256 //purpose  :
2257 //=======================================================================
2258 static int VDrawText (Draw_Interpretor& theDI,
2259                       Standard_Integer  theArgsNb,
2260                       const char**      theArgVec)
2261 {
2262   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
2263   if (theArgsNb < 3)
2264   {
2265     std::cout << "Error: wrong number of arguments! See usage:\n";
2266     theDI.PrintHelp (theArgVec[0]);
2267     return 1;
2268   }
2269   else if (aContext.IsNull())
2270   {
2271     std::cout << "Error: no active view!\n";
2272     return 1;
2273   }
2274
2275   Standard_Integer           anArgIt = 1;
2276   TCollection_ExtendedString aName (theArgVec[anArgIt++], Standard_True);
2277   TCollection_ExtendedString aText (theArgVec[anArgIt++], Standard_True);
2278   Handle(AIS_TextLabel)      aTextPrs;
2279   ViewerTest_AutoUpdater     anAutoUpdater (aContext, ViewerTest::CurrentView());
2280
2281   Standard_Boolean isNewPrs = Standard_False;
2282   if (GetMapOfAIS().IsBound2 (aName))
2283   {
2284     aTextPrs = Handle(AIS_TextLabel)::DownCast (GetMapOfAIS().Find2 (aName));
2285   }
2286
2287   if (aTextPrs.IsNull())
2288   {
2289     isNewPrs = Standard_True;
2290     aTextPrs = new AIS_TextLabel();
2291     aTextPrs->SetFont ("Courier");
2292   }
2293
2294   aTextPrs->SetText (aText);
2295
2296   Handle(Graphic3d_TransformPers) aTrsfPers;
2297   Aspect_TypeOfDisplayText aDisplayType = Aspect_TODT_NORMAL;
2298
2299   Standard_Boolean aHasPlane = Standard_False;
2300   gp_Dir           aNormal;
2301   gp_Dir           aDirection;
2302   gp_Pnt           aPos;
2303
2304   for (; anArgIt < theArgsNb; ++anArgIt)
2305   {
2306     TCollection_AsciiString aParam (theArgVec[anArgIt]);
2307     aParam.LowerCase();
2308
2309     if (anAutoUpdater.parseRedrawMode (aParam))
2310     {
2311       continue;
2312     }
2313     else if (aParam == "-pos"
2314           || aParam == "-position")
2315     {
2316       if (anArgIt + 3 >= theArgsNb)
2317       {
2318         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2319         return 1;
2320       }
2321
2322       aPos.SetX (Draw::Atof (theArgVec[++anArgIt]));
2323       aPos.SetY (Draw::Atof (theArgVec[++anArgIt]));
2324       aPos.SetZ (Draw::Atof (theArgVec[++anArgIt]));
2325       aTextPrs->SetPosition (aPos);
2326     }
2327     else if (aParam == "-color")
2328     {
2329       if (anArgIt + 1 >= theArgsNb)
2330       {
2331         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2332         return 1;
2333       }
2334
2335       TCollection_AsciiString aColor (theArgVec[anArgIt + 1]);
2336       Quantity_NameOfColor aNameOfColor = Quantity_NOC_BLACK;
2337       if (Quantity_Color::ColorFromName (aColor.ToCString(), aNameOfColor))
2338       {
2339         anArgIt += 1;
2340         aTextPrs->SetColor (aNameOfColor);
2341         continue;
2342       }
2343       else if (anArgIt + 3 >= theArgsNb)
2344       {
2345         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2346         return 1;
2347       }
2348
2349       TCollection_AsciiString aGreen (theArgVec[anArgIt + 2]);
2350       TCollection_AsciiString aBlue  (theArgVec[anArgIt + 3]);
2351       if (!aColor.IsRealValue()
2352        || !aGreen.IsRealValue()
2353        || !aBlue.IsRealValue())
2354       {
2355         std::cout << "Error: wrong syntax at '" << aParam.ToCString() << "'.\n";
2356         return 1;
2357       }
2358
2359       const Graphic3d_Vec3d anRGB (aColor.RealValue(),
2360                                    aGreen.RealValue(),
2361                                    aBlue.RealValue());
2362
2363       aTextPrs->SetColor (Quantity_Color (anRGB.r(), anRGB.g(), anRGB.b(), Quantity_TOC_RGB));
2364       anArgIt += 3;
2365     }
2366     else if (aParam == "-halign")
2367     {
2368       if (++anArgIt >= theArgsNb)
2369       {
2370         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2371         return 1;
2372       }
2373
2374       TCollection_AsciiString aType (theArgVec[anArgIt]);
2375       aType.LowerCase();
2376       if (aType == "left")
2377       {
2378         aTextPrs->SetHJustification (Graphic3d_HTA_LEFT);
2379       }
2380       else if (aType == "center")
2381       {
2382         aTextPrs->SetHJustification (Graphic3d_HTA_CENTER);
2383       }
2384       else if (aType == "right")
2385       {
2386         aTextPrs->SetHJustification (Graphic3d_HTA_RIGHT);
2387       }
2388       else
2389       {
2390         std::cout << "Error: wrong syntax at '" << aParam.ToCString() << "'.\n";
2391         return 1;
2392       }
2393     }
2394     else if (aParam == "-valign")
2395     {
2396       if (++anArgIt >= theArgsNb)
2397       {
2398         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2399         return 1;
2400       }
2401
2402       TCollection_AsciiString aType (theArgVec[anArgIt]);
2403       aType.LowerCase();
2404       if (aType == "top")
2405       {
2406         aTextPrs->SetVJustification (Graphic3d_VTA_TOP);
2407       }
2408       else if (aType == "center")
2409       {
2410         aTextPrs->SetVJustification (Graphic3d_VTA_CENTER);
2411       }
2412       else if (aType == "bottom")
2413       {
2414         aTextPrs->SetVJustification (Graphic3d_VTA_BOTTOM);
2415       }
2416       else if (aType == "topfirstline")
2417       {
2418         aTextPrs->SetVJustification (Graphic3d_VTA_TOPFIRSTLINE);
2419       }
2420       else
2421       {
2422         std::cout << "Error: wrong syntax at '" << aParam.ToCString() << "'.\n";
2423         return 1;
2424       }
2425     }
2426     else if (aParam == "-angle")
2427     {
2428       if (++anArgIt >= theArgsNb)
2429       {
2430         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2431         return 1;
2432       }
2433
2434       aTextPrs->SetAngle (Draw::Atof (theArgVec[anArgIt]) * (M_PI / 180.0));
2435     }
2436     else if (aParam == "-zoom")
2437     {
2438       if (++anArgIt >= theArgsNb)
2439       {
2440         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2441         return 1;
2442       }
2443
2444       aTextPrs->SetZoomable (Draw::Atoi (theArgVec[anArgIt]) == 1);
2445     }
2446     else if (aParam == "-height")
2447     {
2448       if (++anArgIt >= theArgsNb)
2449       {
2450         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2451         return 1;
2452       }
2453
2454       aTextPrs->SetHeight (Draw::Atof(theArgVec[anArgIt]));
2455     }
2456     else if (aParam == "-aspect")
2457     {
2458       if (++anArgIt >= theArgsNb)
2459       {
2460         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2461         return 1;
2462       }
2463
2464       TCollection_AsciiString anOption (theArgVec[anArgIt]);
2465       anOption.LowerCase();
2466       if (anOption.IsEqual ("regular"))
2467       {
2468         aTextPrs->SetFontAspect (Font_FA_Regular);
2469       }
2470       else if (anOption.IsEqual ("bold"))
2471       {
2472         aTextPrs->SetFontAspect (Font_FA_Bold);
2473       }
2474       else if (anOption.IsEqual ("italic"))
2475       {
2476         aTextPrs->SetFontAspect (Font_FA_Italic);
2477       }
2478       else if (anOption.IsEqual ("bolditalic"))
2479       {
2480         aTextPrs->SetFontAspect (Font_FA_BoldItalic);
2481       }
2482     }
2483     else if (aParam == "-font")
2484     {
2485       if (++anArgIt >= theArgsNb)
2486       {
2487         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2488         return 1;
2489       }
2490
2491       aTextPrs->SetFont (theArgVec[anArgIt]);
2492     }
2493     else if (aParam == "-plane")
2494     {
2495       if (anArgIt + 6 >= theArgsNb)
2496       {
2497         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2498         return 1;
2499       }
2500
2501       Standard_Real aX = Draw::Atof (theArgVec[++anArgIt]);
2502       Standard_Real aY = Draw::Atof (theArgVec[++anArgIt]);
2503       Standard_Real aZ = Draw::Atof (theArgVec[++anArgIt]);
2504       aNormal.SetCoord (aX, aY, aZ);
2505
2506       aX = Draw::Atof (theArgVec[++anArgIt]);
2507       aY = Draw::Atof (theArgVec[++anArgIt]);
2508       aZ = Draw::Atof (theArgVec[++anArgIt]);
2509       aDirection.SetCoord (aX, aY, aZ);
2510
2511       aHasPlane = Standard_True;
2512     }
2513     else if (aParam == "-flipping")
2514     {
2515       aTextPrs->SetFlipping (Standard_True);
2516     }
2517     else if (aParam == "-disptype"
2518           || aParam == "-displaytype")
2519     {
2520       if (++anArgIt >= theArgsNb)
2521       {
2522         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2523         return 1;
2524       }
2525       TCollection_AsciiString aType (theArgVec[anArgIt]);
2526       aType.LowerCase();
2527       if (aType == "subtitle")
2528         aDisplayType = Aspect_TODT_SUBTITLE;
2529       else if (aType == "decal")
2530         aDisplayType = Aspect_TODT_DEKALE;
2531       else if (aType == "blend")
2532         aDisplayType = Aspect_TODT_BLEND;
2533       else if (aType == "dimension")
2534         aDisplayType = Aspect_TODT_DIMENSION;
2535       else if (aType == "normal")
2536         aDisplayType = Aspect_TODT_NORMAL;
2537       else if (aType == "shadow")
2538         aDisplayType = Aspect_TODT_SHADOW;
2539       else
2540       {
2541         std::cout << "Error: wrong display type '" << aType << "'.\n";
2542         return 1;
2543       }
2544     }
2545     else if (aParam == "-subcolor"
2546           || aParam == "-subtitlecolor")
2547     {
2548       if (anArgIt + 1 >= theArgsNb)
2549       {
2550         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2551         return 1;
2552       }
2553
2554       TCollection_AsciiString aColor (theArgVec[anArgIt + 1]);
2555       Quantity_NameOfColor aNameOfColor = Quantity_NOC_BLACK;
2556       if (Quantity_Color::ColorFromName (aColor.ToCString(), aNameOfColor))
2557       {
2558         anArgIt += 1;
2559         aTextPrs->SetColorSubTitle (aNameOfColor);
2560         continue;
2561       }
2562       else if (anArgIt + 3 >= theArgsNb)
2563       {
2564         std::cout << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2565         return 1;
2566       }
2567
2568       TCollection_AsciiString aGreen (theArgVec[anArgIt + 2]);
2569       TCollection_AsciiString aBlue  (theArgVec[anArgIt + 3]);
2570       if (!aColor.IsRealValue()
2571        || !aGreen.IsRealValue()
2572        || !aBlue.IsRealValue())
2573       {
2574         std::cout << "Error: wrong syntax at '" << aParam.ToCString() << "'.\n";
2575         return 1;
2576       }
2577
2578       const Graphic3d_Vec3d anRGB (aColor.RealValue(),
2579                                    aGreen.RealValue(),
2580                                    aBlue.RealValue());
2581
2582       aTextPrs->SetColorSubTitle (Quantity_Color (anRGB.r(), anRGB.g(), anRGB.b(), Quantity_TOC_RGB));
2583       anArgIt += 3;
2584     }
2585     else if (aParam == "-2d")
2586     {
2587       aTrsfPers = new Graphic3d_TransformPers (Graphic3d_TMF_2d);
2588     }
2589     else if (aParam == "-trsfperspos"
2590           || aParam == "-perspos")
2591     {
2592       if (anArgIt + 2 >= theArgsNb)
2593       {
2594         std::cerr << "Error: wrong number of values for parameter '" << aParam.ToCString() << "'.\n";
2595         return 1;
2596       }
2597
2598       TCollection_AsciiString aX (theArgVec[++anArgIt]);
2599       TCollection_AsciiString aY (theArgVec[++anArgIt]);
2600       TCollection_AsciiString aZ = "0";
2601       if (!aX.IsIntegerValue()
2602        || !aY.IsIntegerValue())
2603       {
2604         std::cerr << "Error: wrong syntax at '" << aParam << "'.\n";
2605         return 1;
2606       }
2607       if (anArgIt + 1 < theArgsNb)
2608       {
2609         TCollection_AsciiString aTemp = theArgVec[anArgIt + 1];
2610         if (aTemp.IsIntegerValue())
2611         {
2612           aZ = aTemp;
2613           ++anArgIt;
2614         }
2615       }
2616
2617       aTrsfPers = Graphic3d_TransformPers::FromDeprecatedParams (Graphic3d_TMF_2d, gp_Pnt (aX.IntegerValue(), aY.IntegerValue(), aZ.IntegerValue()));
2618     }
2619     else
2620     {
2621       std::cout << "Error: unknown argument '" << aParam << "'\n";
2622       return 1;
2623     }
2624   }
2625
2626   if (aHasPlane)
2627   {
2628     aTextPrs->SetOrientation3D (gp_Ax2 (aPos, aNormal, aDirection));
2629   }
2630
2631   aTextPrs->SetDisplayType (aDisplayType);
2632
2633   if (!aTrsfPers.IsNull())
2634   {
2635     aContext->SetTransformPersistence (aTextPrs, aTrsfPers);
2636     aTextPrs->SetZLayer(Graphic3d_ZLayerId_TopOSD);
2637     if (aTextPrs->Position().Z() != 0)
2638     {
2639       aTextPrs->SetPosition (gp_Pnt(aTextPrs->Position().X(), aTextPrs->Position().Y(), 0));
2640     }
2641   }
2642   else if (!aTextPrs->TransformPersistence().IsNull())
2643   {
2644     aContext->SetTransformPersistence (aTextPrs, Handle(Graphic3d_TransformPers)());
2645   }
2646
2647   if (isNewPrs)
2648   {
2649     ViewerTest::Display (aName, aTextPrs, Standard_False);
2650   }
2651   else
2652   {
2653     aContext->Redisplay (aTextPrs, Standard_False, Standard_True);
2654   }
2655   return 0;
2656 }
2657
2658 #include <math.h>
2659 #include <gp_Pnt.hxx>
2660 #include <Graphic3d_ArrayOfPoints.hxx>
2661 #include <Graphic3d_ArrayOfPrimitives.hxx>
2662 #include <Graphic3d_ArrayOfTriangles.hxx>
2663 #include <Poly_Array1OfTriangle.hxx>
2664 #include <Poly_Triangle.hxx>
2665 #include <Poly_Triangulation.hxx>
2666 #include <TColgp_Array1OfPnt.hxx>
2667 #include <TShort_Array1OfShortReal.hxx>
2668 #include <TShort_HArray1OfShortReal.hxx>
2669
2670 #include <AIS_Triangulation.hxx>
2671 #include <StdPrs_ToolTriangulatedShape.hxx>
2672 #include <Poly_Connect.hxx>
2673 #include <TColgp_Array1OfDir.hxx>
2674 #include <Graphic3d_GraphicDriver.hxx>
2675
2676 #include <TColStd_Array1OfInteger.hxx>
2677 #include <TColStd_HArray1OfInteger.hxx>
2678 #include <Prs3d_ShadingAspect.hxx>
2679 #include <Graphic3d_MaterialAspect.hxx>
2680 #include <Graphic3d_AspectFillArea3d.hxx>
2681
2682 #include <BRepPrimAPI_MakeCylinder.hxx>
2683 #include <TopoDS_Shape.hxx>
2684 #include <TopExp_Explorer.hxx>
2685 #include <TopAbs.hxx>
2686 #include <StdSelect_ShapeTypeFilter.hxx>
2687 #include <AIS_InteractiveObject.hxx>
2688
2689
2690 //===============================================================================================
2691 //function : CalculationOfSphere
2692 //author   : psn
2693 //purpose  : Create a Sphere
2694 //===============================================================================================
2695
2696 Handle( Poly_Triangulation ) CalculationOfSphere( double X , double Y , double Z ,
2697                                                   int res ,
2698                                                   double Radius ){
2699   double mRadius = Radius;
2700   double mCenter[3] = {X,Y,Z};
2701   int mThetaResolution;
2702   int mPhiResolution;
2703   double mStartTheta = 0;//StartTheta;
2704   double mEndTheta = 360;//EndTheta;
2705   double mStartPhi = 0;//StartPhi;
2706   double mEndPhi = 180;//EndPhi;
2707   res = res < 4 ? 4 : res;
2708
2709   mThetaResolution = res;
2710   mPhiResolution = res;
2711
2712   int i, j;
2713   int jStart, jEnd, numOffset;
2714   double x[3], n[3], deltaPhi, deltaTheta, phi, theta, radius;
2715   double startTheta, endTheta, startPhi, endPhi;
2716   int base, numPoles=0, thetaResolution, phiResolution;
2717
2718   int pts[3];
2719   int piece = -1;
2720   int numPieces = 1;
2721   if ( numPieces > mThetaResolution ) {
2722     numPieces = mThetaResolution;
2723   }
2724
2725   int localThetaResolution =  mThetaResolution;
2726   double localStartTheta =  mStartTheta;
2727   double localEndTheta =  mEndTheta;
2728
2729   while ( localEndTheta < localStartTheta ) {
2730     localEndTheta += 360.0;
2731   }
2732
2733   deltaTheta = (localEndTheta - localStartTheta) / localThetaResolution;
2734
2735   // Change the ivars based on pieces.
2736   int start, end;
2737   start = piece * localThetaResolution / numPieces;
2738   end = (piece+1) * localThetaResolution / numPieces;
2739   localEndTheta = localStartTheta + (double)(end) * deltaTheta;
2740   localStartTheta = localStartTheta + (double)(start) * deltaTheta;
2741   localThetaResolution = end - start;
2742
2743   // Create north pole if needed
2744   int number_point = 0;
2745   int number_pointArray = 0;
2746
2747   if ( mStartPhi <= 0.0 ) {
2748     number_pointArray++;
2749     numPoles++;
2750   }
2751   if ( mEndPhi >= 180.0 ) {
2752     number_pointArray++;
2753     numPoles++;
2754   }
2755
2756   // Check data, determine increments, and convert to radians
2757   startTheta = (localStartTheta < localEndTheta ? localStartTheta : localEndTheta);
2758   startTheta *= M_PI  / 180.0;
2759   endTheta = (localEndTheta > localStartTheta ? localEndTheta : localStartTheta);
2760   endTheta *= M_PI  / 180.0;
2761
2762
2763   startPhi = ( mStartPhi <  mEndPhi ?  mStartPhi :  mEndPhi);
2764   startPhi *= M_PI  / 180.0;
2765   endPhi = ( mEndPhi >  mStartPhi ?  mEndPhi :  mStartPhi);
2766   endPhi *= M_PI  / 180.0;
2767
2768   phiResolution =  mPhiResolution - numPoles;
2769   deltaPhi = (endPhi - startPhi) / ( mPhiResolution - 1);
2770   thetaResolution = localThetaResolution;
2771   if ( fabs(localStartTheta - localEndTheta) < 360.0 ) {
2772     ++localThetaResolution;
2773   }
2774   deltaTheta = (endTheta - startTheta) / thetaResolution;
2775
2776   jStart = ( mStartPhi <= 0.0 ? 1 : 0);
2777   jEnd = ( mEndPhi >= 180.0 ?  mPhiResolution - 1 :  mPhiResolution);
2778
2779   // Create intermediate points
2780   for ( i = 0; i < localThetaResolution; i++ ) {
2781     for ( j = jStart; j < jEnd; j++ ) {
2782         number_pointArray++;
2783     }
2784   }
2785
2786   //Generate mesh connectivity
2787   base = phiResolution * localThetaResolution;
2788
2789   int number_triangle = 0 ;
2790   if ( mStartPhi <= 0.0 ) { // around north pole
2791     number_triangle += localThetaResolution;
2792   }
2793
2794   if ( mEndPhi >= 180.0 ) { // around south pole
2795     number_triangle += localThetaResolution;
2796   }
2797
2798   // bands in-between poles
2799   for ( i=0; i < localThetaResolution; i++){
2800     for ( j=0; j < (phiResolution-1); j++){
2801        number_triangle +=2;
2802     }
2803   }
2804
2805   Handle( Poly_Triangulation ) polyTriangulation = new Poly_Triangulation(number_pointArray, number_triangle, false);
2806   TColgp_Array1OfPnt& PointsOfArray = polyTriangulation->ChangeNodes();
2807   Poly_Array1OfTriangle& pArrayTriangle = polyTriangulation->ChangeTriangles();
2808
2809   if (  mStartPhi <= 0.0 ){
2810       x[0] =  mCenter[0];
2811       x[1] =  mCenter[1];
2812       x[2] =  mCenter[2] +  mRadius;
2813       PointsOfArray.SetValue(1,gp_Pnt(x[0],x[1],x[2]));
2814   }
2815
2816   // Create south pole if needed
2817   if (  mEndPhi >= 180.0 ){
2818       x[0] =  mCenter[0];
2819       x[1] =  mCenter[1];
2820       x[2] =  mCenter[2] -  mRadius;
2821       PointsOfArray.SetValue(2,gp_Pnt(x[0],x[1],x[2]));
2822   }
2823
2824   number_point = 3;
2825   for ( i=0; i < localThetaResolution; i++){
2826     theta = localStartTheta * M_PI / 180.0 + i*deltaTheta;
2827     for ( j = jStart; j < jEnd; j++){
2828         phi = startPhi + j*deltaPhi;
2829         radius =  mRadius * sin((double)phi);
2830         n[0] = radius * cos((double)theta);
2831         n[1] = radius * sin((double)theta);
2832         n[2] =  mRadius * cos((double)phi);
2833         x[0] = n[0] +  mCenter[0];
2834         x[1] = n[1] +  mCenter[1];
2835         x[2] = n[2] +  mCenter[2];
2836         PointsOfArray.SetValue(number_point,gp_Pnt(x[0],x[1],x[2]));
2837         number_point++;
2838       }
2839     }
2840
2841   numPoles = 3;
2842   number_triangle = 1;
2843   if ( mStartPhi <= 0.0 ){// around north pole
2844     for (i=0; i < localThetaResolution; i++){
2845         pts[0] = phiResolution*i + numPoles;
2846         pts[1] = (phiResolution*(i+1) % base) + numPoles;
2847         pts[2] = 1;
2848         pArrayTriangle.SetValue(number_triangle,Poly_Triangle(pts[0],pts[1],pts[2]));
2849         number_triangle++;
2850       }
2851     }
2852
2853   if (  mEndPhi >= 180.0 ){ // around south pole
2854     numOffset = phiResolution - 1 + numPoles;
2855     for (i=0; i < localThetaResolution; i++){
2856         pts[0] = phiResolution*i + numOffset;
2857         pts[2] = ((phiResolution*(i+1)) % base) + numOffset;
2858         pts[1] = numPoles - 1;
2859         pArrayTriangle.SetValue(number_triangle,Poly_Triangle(pts[0],pts[1],pts[2]));
2860         number_triangle++;
2861       }
2862     }
2863
2864   // bands in-between poles
2865
2866   for (i=0; i < localThetaResolution; i++){
2867     for (j=0; j < (phiResolution-1); j++){
2868         pts[0] = phiResolution*i + j + numPoles;
2869         pts[1] = pts[0] + 1;
2870         pts[2] = ((phiResolution*(i+1)+j) % base) + numPoles + 1;
2871         pArrayTriangle.SetValue(number_triangle,Poly_Triangle(pts[0],pts[1],pts[2]));
2872         number_triangle++;
2873         pts[1] = pts[2];
2874         pts[2] = pts[1] - 1;
2875         pArrayTriangle.SetValue(number_triangle,Poly_Triangle(pts[0],pts[1],pts[2]));
2876         number_triangle++;
2877       }
2878     }
2879
2880   Poly_Connect* pc = new Poly_Connect(polyTriangulation);
2881
2882   Handle(TShort_HArray1OfShortReal) Normals = new TShort_HArray1OfShortReal(1, polyTriangulation->NbNodes() * 3);
2883
2884   Standard_Integer index[3];
2885   Standard_Real Tol = Precision::Confusion();
2886
2887   gp_Dir Nor;
2888   for (i = PointsOfArray.Lower(); i <= PointsOfArray.Upper(); i++) {
2889       gp_XYZ eqPlan(0, 0, 0);
2890       for ( pc->Initialize(i); pc->More(); pc->Next()) {
2891         pArrayTriangle(pc->Value()).Get(index[0], index[1], index[2]);
2892         gp_XYZ v1(PointsOfArray(index[1]).Coord()-PointsOfArray(index[0]).Coord());
2893         gp_XYZ v2(PointsOfArray(index[2]).Coord()-PointsOfArray(index[1]).Coord());
2894         gp_XYZ vv = v1^v2;
2895         Standard_Real mod = vv.Modulus();
2896         if(mod < Tol) continue;
2897         eqPlan += vv/mod;
2898       }
2899
2900       Standard_Real modmax = eqPlan.Modulus();
2901
2902       if(modmax > Tol)
2903         Nor = gp_Dir(eqPlan);
2904       else
2905         Nor = gp_Dir(0., 0., 1.);
2906       
2907       Standard_Integer k = (i - PointsOfArray.Lower()) * 3;
2908       Normals->SetValue(k + 1, (Standard_ShortReal)Nor.X());
2909       Normals->SetValue(k + 2, (Standard_ShortReal)Nor.Y());
2910       Normals->SetValue(k + 3, (Standard_ShortReal)Nor.Z());
2911   }
2912
2913   delete pc;
2914   polyTriangulation->SetNormals(Normals);
2915
2916   return polyTriangulation;
2917 }
2918
2919 //===============================================================================================
2920 //function : VDrawSphere
2921 //author   : psn
2922 //purpose  : Create an AIS shape.
2923 //===============================================================================================
2924 static int VDrawSphere (Draw_Interpretor& /*di*/, Standard_Integer argc, const char** argv)
2925 {
2926   // check for errors
2927   Handle(AIS_InteractiveContext) aContextAIS = ViewerTest::GetAISContext();
2928   if (aContextAIS.IsNull())
2929   {
2930     std::cout << "Call vinit before!\n";
2931     return 1;
2932   }
2933   else if (argc < 3)
2934   {
2935     std::cout << "Use: " << argv[0]
2936               << " shapeName Fineness [X=0.0 Y=0.0 Z=0.0] [Radius=100.0] [ToShowEdges=0]\n";
2937     return 1;
2938   }
2939
2940   // read the arguments
2941   TCollection_AsciiString aShapeName (argv[1]);
2942   Standard_Integer aResolution = Draw::Atoi (argv[2]);
2943   Standard_Real aCenterX = (argc > 5) ? Draw::Atof (argv[3]) : 0.0;
2944   Standard_Real aCenterY = (argc > 5) ? Draw::Atof (argv[4]) : 0.0;
2945   Standard_Real aCenterZ = (argc > 5) ? Draw::Atof (argv[5]) : 0.0;
2946   Standard_Real aRadius =  (argc > 6) ? Draw::Atof (argv[6]) : 100.0;
2947   Standard_Boolean toShowEdges = (argc > 7) ? Draw::Atoi (argv[7]) == 1 : Standard_False;
2948   Standard_Boolean toPrintInfo = (argc > 8) ? Draw::Atoi (argv[8]) == 1 : Standard_True;
2949
2950   // remove AIS object with given name from map
2951   VDisplayAISObject (aShapeName, Handle(AIS_InteractiveObject)());
2952
2953   if (toPrintInfo)
2954     std::cout << "Compute Triangulation...\n";
2955   Handle(AIS_Triangulation) aShape
2956     = new AIS_Triangulation (CalculationOfSphere (aCenterX, aCenterY, aCenterZ,
2957                                                   aResolution,
2958                                                   aRadius));
2959   Standard_Integer aNumberPoints    = aShape->GetTriangulation()->Nodes().Length();
2960   Standard_Integer aNumberTriangles = aShape->GetTriangulation()->Triangles().Length();
2961
2962   // stupid initialization of Green color in RGBA space as integer
2963   // probably wrong for big-endian CPUs
2964   const Graphic3d_Vec4ub aColor (0, 255, 0, 0);
2965
2966   // setup colors array per vertex
2967   Handle(TColStd_HArray1OfInteger) aColorArray = new TColStd_HArray1OfInteger (1, aNumberPoints);
2968   for (Standard_Integer aNodeId = 1; aNodeId <= aNumberPoints; ++aNodeId)
2969   {
2970     aColorArray->SetValue (aNodeId, *reinterpret_cast<const Standard_Integer*> (aColor.GetData()));
2971   }
2972   aShape->SetColors (aColorArray);
2973
2974   // show statistics
2975   Standard_Integer aPointsSize      = aNumberPoints * 3 * sizeof(float);  // 3x GLfloat
2976   Standard_Integer aNormalsSize     = aNumberPoints * 3 * sizeof(float);  // 3x GLfloat
2977   Standard_Integer aColorsSize      = aNumberPoints * 3 * sizeof(float);  // 3x GLfloat without alpha
2978   Standard_Integer aTrianglesSize   = aNumberTriangles * 3 * sizeof(int); // 3x GLint
2979   Standard_Integer aPolyConnectSize = aNumberPoints * 4 + aNumberTriangles * 6 * 4;
2980   Standard_Integer aTotalSize       = aPointsSize + aNormalsSize + aColorsSize + aTrianglesSize;
2981   aTotalSize >>= 20; //MB
2982   aNormalsSize >>= 20;
2983   aColorsSize >>= 20;
2984   aTrianglesSize >>= 20;
2985   aPolyConnectSize >>= 20;
2986   if (toPrintInfo)
2987   {
2988     std::cout << "NumberOfPoints:    " << aNumberPoints << "\n"
2989       << "NumberOfTriangles: " << aNumberTriangles << "\n"
2990       << "Amount of memory required for PolyTriangulation without Normals: " << (aTotalSize - aNormalsSize) << " Mb\n"
2991       << "Amount of memory for colors: " << aColorsSize << " Mb\n"
2992       << "Amount of memory for PolyConnect: " << aPolyConnectSize << " Mb\n"
2993       << "Amount of graphic card memory required: " << aTotalSize << " Mb\n";
2994   }
2995
2996   // Setting material properties, very important for desirable visual result!
2997   Graphic3d_MaterialAspect aMat (Graphic3d_NOM_PLASTIC);
2998   aMat.SetAmbient  (0.2f);
2999   aMat.SetSpecular (0.5f);
3000   Handle(Graphic3d_AspectFillArea3d) anAspect
3001     = new Graphic3d_AspectFillArea3d (Aspect_IS_SOLID,
3002                                       Quantity_NOC_RED,
3003                                       Quantity_NOC_YELLOW,
3004                                       Aspect_TOL_SOLID,
3005                                       1.0,
3006                                       aMat,
3007                                       aMat);
3008   Handle(Prs3d_ShadingAspect) aShAsp = new Prs3d_ShadingAspect();
3009   if (toShowEdges)
3010   {
3011     anAspect->SetEdgeOn();
3012   }
3013   else
3014   {
3015     anAspect->SetEdgeOff();
3016   }
3017   aShAsp->SetAspect (anAspect);
3018   aShape->Attributes()->SetShadingAspect (aShAsp);
3019
3020   VDisplayAISObject (aShapeName, aShape);
3021   return 0;
3022 }
3023
3024 //=============================================================================
3025 //function : VComputeHLR
3026 //purpose  :
3027 //=============================================================================
3028
3029 static int VComputeHLR (Draw_Interpretor& ,
3030                         Standard_Integer theArgNb,
3031                         const char** theArgVec)
3032 {
3033   TCollection_AsciiString aShapeName, aHlrName;
3034   TopoDS_Shape aSh;
3035   gp_Pnt anEye;
3036   gp_Dir aDir;
3037   gp_Ax2 aProjAx;
3038   bool hasViewDirArg = false;
3039   Prs3d_TypeOfHLR anAlgoType = Prs3d_TOH_PolyAlgo;
3040   bool toShowTangentEdges = false, toShowHiddenEdges = false;
3041   int aNbIsolines = 0;
3042   if (Handle(V3d_Viewer) aViewer = ViewerTest::GetViewerFromContext())
3043   {
3044     gp_Dir aRight;
3045     Handle(V3d_View) aView = ViewerTest::CurrentView();
3046     Standard_Integer aWidth, aHeight;
3047     Standard_Real aCentX, aCentY, aCentZ, aDirX, aDirY, aDirZ;
3048     Standard_Real aRightX, aRightY, aRightZ;
3049     aView->Window()->Size (aWidth, aHeight);
3050
3051     aView->ConvertWithProj (aWidth, aHeight/2, 
3052                             aRightX, aRightY, aRightZ,
3053                             aDirX, aDirY, aDirZ);
3054     aView->ConvertWithProj (aWidth/2, aHeight/2, 
3055                             aCentX, aCentY, aCentZ,
3056                             aDirX, aDirY, aDirZ);
3057
3058     anEye.SetCoord (-aCentX, -aCentY, -aCentZ);
3059     aDir.SetCoord (-aDirX, -aDirY, -aDirZ);
3060     aRight.SetCoord (aRightX - aCentX, aRightY - aCentY, aRightZ - aCentZ);
3061     aProjAx.SetLocation (anEye);
3062     aProjAx.SetDirection (aDir);
3063     aProjAx.SetXDirection (aRight);
3064   }
3065   for (Standard_Integer anArgIter = 1; anArgIter < theArgNb; ++anArgIter)
3066   {
3067     TCollection_AsciiString anArgCase (theArgVec[anArgIter]);
3068     anArgCase.LowerCase();
3069     if (anArgIter + 1 < theArgNb
3070      && (anArgCase == "-algotype"
3071       || anArgCase == "-algo"
3072       || anArgCase == "-type"))
3073     {
3074       TCollection_AsciiString anArgNext (theArgVec[++anArgIter]);
3075       anArgNext.LowerCase();
3076       if (anArgNext == "polyalgo")
3077       {
3078         anAlgoType = Prs3d_TOH_PolyAlgo;
3079       }
3080       else if (anArgNext == "algo")
3081       {
3082         anAlgoType = Prs3d_TOH_Algo;
3083       }
3084       else
3085       {
3086         std::cout << "Syntax error: unknown algo type '" << anArgNext << "'\n";
3087         return 1;
3088       }
3089     }
3090     else if (anArgCase == "-showhiddenedges"
3091           || anArgCase == "-hiddenedges"
3092           || anArgCase == "-hidden")
3093     {
3094       toShowHiddenEdges = true;
3095       if (anArgIter + 1 < theArgNb
3096        && ViewerTest::ParseOnOff (theArgVec[anArgIter + 1], toShowHiddenEdges))
3097       {
3098         ++anArgIter;
3099       }
3100     }
3101     else if (anArgCase == "-showtangentedges"
3102           || anArgCase == "-tangentedges"
3103           || anArgCase == "-tangent")
3104     {
3105       toShowTangentEdges = true;
3106       if (anArgIter + 1 < theArgNb
3107        && ViewerTest::ParseOnOff (theArgVec[anArgIter + 1], toShowTangentEdges))
3108       {
3109         ++anArgIter;
3110       }
3111     }
3112     else if (anArgIter + 1 < theArgNb
3113           && (anArgCase == "-nbiso"
3114            || anArgCase == "-nbisolines"))
3115     {
3116       aNbIsolines = Draw::Atoi (theArgVec[++anArgIter]);
3117     }
3118     else if (aSh.IsNull())
3119     {
3120       aSh = DBRep::Get (theArgVec[anArgIter]);
3121       aShapeName = theArgVec[anArgIter];
3122       if (aSh.IsNull())
3123       {
3124         BRep_Builder aBrepBuilder;
3125         BRepTools::Read (aSh, theArgVec[anArgIter], aBrepBuilder);
3126         if (aSh.IsNull())
3127         {
3128           std::cout << "Syntax error: no shape with name " << theArgVec[anArgIter] << " found\n";
3129           return 1;
3130         }
3131       }
3132     }
3133     else if (aHlrName.IsEmpty())
3134     {
3135       aHlrName = theArgVec[anArgIter];
3136     }
3137     else if (!hasViewDirArg
3138           && anArgIter + 8 < theArgNb)
3139     {
3140       hasViewDirArg = true;
3141       gp_Dir anUp;
3142       anEye.SetCoord (Draw::Atof (theArgVec[anArgIter + 0]), Draw::Atof (theArgVec[anArgIter + 1]), Draw::Atof (theArgVec[anArgIter + 2]));
3143       aDir .SetCoord (Draw::Atof (theArgVec[anArgIter + 3]), Draw::Atof (theArgVec[anArgIter + 4]), Draw::Atof (theArgVec[anArgIter + 5]));
3144       anUp .SetCoord (Draw::Atof (theArgVec[anArgIter + 6]), Draw::Atof (theArgVec[anArgIter + 7]), Draw::Atof (theArgVec[anArgIter + 8]));
3145       aProjAx.SetLocation (anEye);
3146       aProjAx.SetDirection (aDir);
3147       aProjAx.SetYDirection (anUp);
3148       anArgIter += 8;
3149     }
3150     else
3151     {
3152       std::cout << "Syntax error: unknown argument '" << theArgVec[anArgIter] << "'\n";
3153       return 1;
3154     }
3155   }
3156
3157   if (aHlrName.IsEmpty() || aSh.IsNull()
3158    || (ViewerTest::GetAISContext().IsNull() && hasViewDirArg))
3159   {
3160     std::cout << "Syntax error: wrong number of arguments\n";
3161     return 1;
3162   }
3163
3164   HLRAlgo_Projector aProjector (aProjAx);
3165   TopoDS_Shape aVisible[6];
3166   TopoDS_Shape aHidden[6];
3167   if (anAlgoType == Prs3d_TOH_PolyAlgo)
3168   {
3169     Handle(HLRBRep_PolyAlgo) aPolyAlgo = new HLRBRep_PolyAlgo();
3170     aPolyAlgo->Projector (aProjector);
3171     aPolyAlgo->Load (aSh);
3172     aPolyAlgo->Update();
3173
3174     HLRBRep_PolyHLRToShape aHLRToShape;
3175     aHLRToShape.Update (aPolyAlgo);
3176
3177     aVisible[HLRBRep_Sharp]   = aHLRToShape.VCompound();
3178     aVisible[HLRBRep_OutLine] = aHLRToShape.OutLineVCompound(); // extract visible outlines
3179     aVisible[HLRBRep_RgNLine] = aHLRToShape.RgNLineVCompound();
3180     if (toShowTangentEdges)
3181     {
3182       aVisible[HLRBRep_Rg1Line] = aHLRToShape.Rg1LineVCompound();
3183     }
3184     if (toShowHiddenEdges)
3185     {
3186       aHidden[HLRBRep_Sharp]   = aHLRToShape.HCompound();
3187       aHidden[HLRBRep_OutLine] = aHLRToShape.OutLineHCompound();
3188       aHidden[HLRBRep_RgNLine] = aHLRToShape.RgNLineHCompound();
3189       if (toShowTangentEdges)
3190       {
3191         aHidden[HLRBRep_Rg1Line] = aHLRToShape.Rg1LineHCompound();
3192       }
3193     }
3194   }
3195   else
3196   {
3197     Handle(HLRBRep_Algo) aHlrAlgo = new HLRBRep_Algo();
3198     aHlrAlgo->Add (aSh, aNbIsolines);
3199     aHlrAlgo->Projector (aProjector);
3200     aHlrAlgo->Update();
3201     aHlrAlgo->Hide();
3202
3203     HLRBRep_HLRToShape aHLRToShape (aHlrAlgo);
3204     aVisible[HLRBRep_Sharp]   = aHLRToShape.VCompound();
3205     aVisible[HLRBRep_OutLine] = aHLRToShape.OutLineVCompound();
3206     aVisible[HLRBRep_RgNLine] = aHLRToShape.RgNLineVCompound();
3207     if (toShowTangentEdges)
3208     {
3209       aVisible[HLRBRep_Rg1Line] = aHLRToShape.Rg1LineVCompound();
3210     }
3211     aVisible[HLRBRep_IsoLine] = aHLRToShape.IsoLineVCompound();
3212
3213     if (toShowHiddenEdges)
3214     {
3215       aHidden[HLRBRep_Sharp]   = aHLRToShape.HCompound();
3216       aHidden[HLRBRep_OutLine] = aHLRToShape.OutLineHCompound();
3217       aHidden[HLRBRep_RgNLine] = aHLRToShape.RgNLineHCompound();
3218       if (toShowTangentEdges)
3219       {
3220         aHidden[HLRBRep_Rg1Line] = aHLRToShape.Rg1LineHCompound();
3221       }
3222       aHidden[HLRBRep_IsoLine] = aHLRToShape.IsoLineHCompound();
3223     }
3224     // extract 3d
3225     //aVisible[HLRBRep_Sharp]   = aHLRToShape.CompoundOfEdges (HLRBRep_Sharp, Standard_True, Standard_True);
3226     //aVisible[HLRBRep_OutLine] = aHLRToShape.OutLineVCompound3d();
3227   }
3228
3229   TopoDS_Compound aCompRes, aCompVis, aCompHid;
3230   BRep_Builder aBuilder;
3231   aBuilder.MakeCompound (aCompVis);
3232   aBuilder.MakeCompound (aCompHid);
3233   aBuilder.MakeCompound (aCompRes);
3234   for (int aTypeIter = 0; aTypeIter < 6; ++aTypeIter)
3235   {
3236     if (!aVisible[aTypeIter].IsNull())
3237     {
3238       aBuilder.Add (aCompVis, aVisible[aTypeIter]);
3239     }
3240     if (!aHidden[aTypeIter].IsNull())
3241     {
3242       aBuilder.Add (aCompHid, aHidden[aTypeIter]);
3243     }
3244   }
3245   aBuilder.Add (aCompRes, aCompVis);
3246   aBuilder.Add (aCompRes, aCompHid);
3247
3248   // create an AIS shape and display it
3249   if (!ViewerTest::GetAISContext().IsNull())
3250   {
3251     Handle(AIS_ColoredShape) anObject = new AIS_ColoredShape (aCompRes);
3252     if (toShowHiddenEdges)
3253     {
3254       Handle(Prs3d_LineAspect) aLineAspect = new Prs3d_LineAspect (Quantity_Color (Quantity_NOC_RED), Aspect_TOL_DASH, 1.0f);
3255       for (int aTypeIter = 0; aTypeIter < 6; ++aTypeIter)
3256       {
3257         if (!aHidden[aTypeIter].IsNull())
3258         {
3259           Handle(AIS_ColoredDrawer) aDrawer = anObject->CustomAspects (aHidden[aTypeIter]);
3260           aDrawer->SetLineAspect (aLineAspect);
3261           aDrawer->SetWireAspect (aLineAspect);
3262           aDrawer->SetFreeBoundaryAspect (aLineAspect);
3263           aDrawer->SetUnFreeBoundaryAspect (aLineAspect);
3264         }
3265       }
3266     }
3267     ViewerTest::Display (aHlrName, anObject, true);
3268   }
3269
3270   DBRep::Set (aHlrName.ToCString(), aCompRes);
3271   return 0;
3272 }
3273
3274 // This class is a wrap for Graphic3d_ArrayOfPrimitives; it is used for
3275 // manipulating and displaying such an array with AIS context
3276
3277 class MyPArrayObject : public AIS_InteractiveObject
3278 {
3279
3280 public:
3281
3282   MyPArrayObject (const Handle(Graphic3d_ArrayOfPrimitives)& thePArray) : myPArray (thePArray) {}
3283
3284   MyPArrayObject (Handle(TColStd_HArray1OfAsciiString) theArrayDescription,
3285                   Handle(Graphic3d_AspectMarker3d) theMarkerAspect = NULL);
3286
3287   DEFINE_STANDARD_RTTI_INLINE(MyPArrayObject,AIS_InteractiveObject);
3288
3289   virtual Standard_Boolean AcceptDisplayMode (const Standard_Integer theMode) const Standard_OVERRIDE { return theMode == 0; }
3290
3291 private:
3292
3293   void Compute (const Handle(PrsMgr_PresentationManager3d)& aPresentationManager,
3294                 const Handle(Prs3d_Presentation)& aPresentation,
3295                 const Standard_Integer aMode) Standard_OVERRIDE;
3296
3297   void ComputeSelection (const Handle(SelectMgr_Selection)& theSelection,
3298                          const Standard_Integer /*theMode*/) Standard_OVERRIDE;
3299
3300   bool CheckInputCommand (const TCollection_AsciiString theCommand,
3301                           const Handle(TColStd_HArray1OfAsciiString) theArgsArray,
3302                           Standard_Integer &theArgIndex,
3303                           Standard_Integer theArgCount,
3304                           Standard_Integer theMaxArgs);
3305
3306 protected:
3307
3308   Handle(TColStd_HArray1OfAsciiString) myArrayDescription;
3309   Handle(Graphic3d_AspectMarker3d) myMarkerAspect;
3310   Handle(Graphic3d_ArrayOfPrimitives) myPArray;
3311
3312 };
3313
3314 void MyPArrayObject::Compute (const Handle(PrsMgr_PresentationManager3d)& /*aPresentationManager*/,
3315                               const Handle(Prs3d_Presentation)& thePrs,
3316                               const Standard_Integer theMode)
3317 {
3318   if (myPArray.IsNull() || theMode != 0)
3319   {
3320     return;
3321   }
3322
3323   Handle(Graphic3d_Group) aGroup = thePrs->NewGroup();
3324   if (!myMarkerAspect.IsNull())
3325   {
3326     aGroup->SetGroupPrimitivesAspect (myMarkerAspect);
3327   }
3328   aGroup->SetGroupPrimitivesAspect (myDrawer->LineAspect()->Aspect());
3329   aGroup->SetGroupPrimitivesAspect (myDrawer->ShadingAspect()->Aspect());
3330   aGroup->AddPrimitiveArray (myPArray);
3331 }
3332
3333 MyPArrayObject::MyPArrayObject (Handle(TColStd_HArray1OfAsciiString) theArrayDescription,
3334                                 Handle(Graphic3d_AspectMarker3d) theMarkerAspect)
3335 {
3336   myArrayDescription = theArrayDescription;
3337   myMarkerAspect = theMarkerAspect;
3338
3339   // Parsing array description
3340   Standard_Integer aVertexNum = 0, aBoundNum = 0, aEdgeNum = 0;
3341   Graphic3d_ArrayFlags anArrayFlags = Graphic3d_ArrayFlags_None;
3342
3343   Standard_Integer anArgIndex = 0;
3344   Standard_Integer anArgsCount = myArrayDescription->Length();
3345   TCollection_AsciiString anArrayType = myArrayDescription->Value (anArgIndex++);
3346
3347   TCollection_AsciiString aCommand;
3348   while (anArgIndex < anArgsCount)
3349   {
3350     aCommand = myArrayDescription->Value (anArgIndex);
3351     aCommand.LowerCase();
3352
3353     // vertex command
3354     if (CheckInputCommand ("v", myArrayDescription, anArgIndex, 3, anArgsCount))
3355     {
3356       // vertex has a normal or normal with color or texel
3357       if (CheckInputCommand ("n", myArrayDescription, anArgIndex, 3, anArgsCount))
3358       {
3359         anArrayFlags = anArrayFlags | Graphic3d_ArrayFlags_VertexNormal;
3360       }
3361
3362       // vertex has a color
3363       if (CheckInputCommand ("c", myArrayDescription, anArgIndex, 3, anArgsCount))
3364       {
3365         anArrayFlags = anArrayFlags | Graphic3d_ArrayFlags_VertexColor;
3366       }
3367
3368       // vertex has a texel
3369       if (CheckInputCommand ("t", myArrayDescription, anArgIndex, 2, anArgsCount))
3370       {
3371         anArrayFlags = anArrayFlags | Graphic3d_ArrayFlags_VertexTexel;
3372       }
3373
3374       aVertexNum++;
3375     }
3376     // bound command
3377     else if (CheckInputCommand ("b", myArrayDescription, anArgIndex, 1, anArgsCount))
3378     {
3379       // bound has color
3380       if (CheckInputCommand ("c", myArrayDescription, anArgIndex, 3, anArgsCount))
3381       {
3382         anArrayFlags = anArrayFlags | Graphic3d_ArrayFlags_BoundColor;
3383       }
3384
3385       aBoundNum++;
3386     }
3387     // edge command
3388     else if (CheckInputCommand ("e", myArrayDescription, anArgIndex, 1, anArgsCount))
3389     {
3390       aEdgeNum++;
3391     }
3392     // unknown command
3393     else
3394       anArgIndex++;
3395   }
3396
3397   Handle(Graphic3d_ArrayOfPrimitives) anArray;
3398   if (anArrayType == "points")
3399   {
3400     anArray = new Graphic3d_ArrayOfPoints (aVertexNum);
3401   }
3402   else if (anArrayType == "segments")
3403     anArray = new Graphic3d_ArrayOfSegments (aVertexNum, aEdgeNum, anArrayFlags);
3404   else if (anArrayType == "polylines")
3405     anArray = new Graphic3d_ArrayOfPolylines (aVertexNum, aBoundNum, aEdgeNum, anArrayFlags);
3406   else if (anArrayType == "triangles")
3407     anArray = new Graphic3d_ArrayOfTriangles (aVertexNum, aEdgeNum, anArrayFlags);
3408   else if (anArrayType == "trianglefans")
3409     anArray = new Graphic3d_ArrayOfTriangleFans (aVertexNum, aBoundNum, anArrayFlags);
3410   else if (anArrayType == "trianglestrips")
3411     anArray = new Graphic3d_ArrayOfTriangleStrips (aVertexNum, aBoundNum, anArrayFlags);
3412   else if (anArrayType == "quads")
3413     anArray = new Graphic3d_ArrayOfQuadrangles (aVertexNum, aEdgeNum, anArrayFlags);
3414   else if (anArrayType == "quadstrips")
3415     anArray = new Graphic3d_ArrayOfQuadrangleStrips (aVertexNum, aBoundNum, anArrayFlags);
3416   else if (anArrayType == "polygons")
3417     anArray = new Graphic3d_ArrayOfPolygons (aVertexNum, aBoundNum, aEdgeNum, anArrayFlags);
3418
3419   anArgIndex = 1;