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