Salome HOME
Fix for issue #411: remove result shape to avoid the naming selection problem (select...
[modules/shaper.git] / src / Model / Model_AttributeSelection.cpp
1 // Copyright (C) 2014-20xx CEA/DEN, EDF R&D
2
3 // File:        Model_AttributeSelection.h
4 // Created:     2 Oct 2014
5 // Author:      Mikhail PONIKAROV
6
7 #include "Model_AttributeSelection.h"
8 #include "Model_Application.h"
9 #include "Model_Events.h"
10 #include "Model_Data.h"
11 #include "Model_Document.h"
12 #include <ModelAPI_Feature.h>
13 #include <ModelAPI_ResultBody.h>
14 #include <ModelAPI_ResultConstruction.h>
15 #include <ModelAPI_CompositeFeature.h>
16 #include <GeomAPI_Shape.h>
17 #include <GeomAPI_PlanarEdges.h>
18 #include <GeomAlgoAPI_SketchBuilder.h>
19 #include <Events_Error.h>
20
21 #include <TNaming_Selector.hxx>
22 #include <TNaming_NamedShape.hxx>
23 #include <TNaming_Tool.hxx>
24 #include <TNaming_Builder.hxx>
25 #include <TNaming_Localizer.hxx>
26 #include <TopoDS_Shape.hxx>
27 #include <TopoDS_Compound.hxx>
28 #include <TDataStd_IntPackedMap.hxx>
29 #include <TDataStd_Integer.hxx>
30 #include <TDataStd_UAttribute.hxx>
31 #include <TDataStd_Name.hxx>
32 #include <TopTools_MapOfShape.hxx>
33 #include <TopTools_IndexedMapOfShape.hxx>
34 #include <TopTools_MapIteratorOfMapOfShape.hxx>
35 #include <TopTools_ListOfShape.hxx>
36 #include <TopExp_Explorer.hxx>
37 #include <TDF_LabelMap.hxx>
38 #include <BRep_Tool.hxx>
39 #include <BRep_Builder.hxx>
40 #include <TopoDS_Edge.hxx>
41 #include <TopoDS.hxx>
42 #include <TopExp.hxx>
43 #include <TColStd_MapOfTransient.hxx>
44 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
45 #include <TopTools_ListIteratorOfListOfShape.hxx>
46 #include <TColStd_MapIteratorOfPackedMapOfInteger.hxx>
47 #include <gp_Pnt.hxx>
48 #include <Precision.hxx>
49 #include <TDF_ChildIterator.hxx>
50 #include <TDF_ChildIDIterator.hxx>
51 #include <TDataStd_Name.hxx>
52 #include <TopAbs_ShapeEnum.hxx>
53 #include <TopoDS_Iterator.hxx>
54 using namespace std;
55 //#define DEB_NAMING 1
56 #ifdef DEB_NAMING
57 #include <BRepTools.hxx>
58 #endif
59 /// adeed to the index in the packed map to signalize that the vertex of edge is seleted
60 /// (multiplied by the index of the edge)
61 static const int kSTART_VERTEX_DELTA = 1000000;
62 // identifier that there is simple reference: selection equals to context
63 Standard_GUID kSIMPLE_REF_ID("635eacb2-a1d6-4dec-8348-471fae17cb29");
64
65 // on this label is stored:
66 // TNaming_NamedShape - selected shape
67 // TNaming_Naming - topological selection information (for the body)
68 // TDataStd_IntPackedMap - indexes of edges in composite element (for construction)
69 // TDataStd_Integer - type of the selected shape (for construction)
70 // TDF_Reference - from ReferenceAttribute, the context
71 #define DDDD 1
72 void Model_AttributeSelection::setValue(const ResultPtr& theContext,
73   const std::shared_ptr<GeomAPI_Shape>& theSubShape)
74 {
75   const std::shared_ptr<GeomAPI_Shape>& anOldShape = value();
76   bool isOldShape = 
77     (theSubShape == anOldShape || (theSubShape && anOldShape && theSubShape->isEqual(anOldShape)));
78   if (isOldShape) return; // shape is the same, so context is also unchanged
79   // update the referenced object if needed
80   bool isOldContext = theContext == myRef.value();
81
82
83   if (!isOldContext)
84     myRef.setValue(theContext);
85
86   // do noth use naming if selected shape is result shape itself, but not sub-shape
87   TDF_Label aSelLab = selectionLabel();
88   aSelLab.ForgetAttribute(kSIMPLE_REF_ID);
89   if (theContext->groupName() == ModelAPI_ResultBody::group()) {
90     // do not select the whole shape for body:it is already must be in the data framework
91     if (theContext->shape().get() && theContext->shape()->isEqual(theSubShape)) {
92       aSelLab.ForgetAllAttributes(true);
93       TDataStd_UAttribute::Set(aSelLab, kSIMPLE_REF_ID);
94     } else {
95       selectBody(theContext, theSubShape);
96     }
97   } else if (theContext->groupName() == ModelAPI_ResultConstruction::group()) {
98     selectConstruction(theContext, theSubShape);
99   }
100   myIsInitialized = true;
101
102   std::string aSelName = namingName();
103   if(!aSelName.empty())
104           TDataStd_Name::Set(selectionLabel(), aSelName.c_str()); //set name
105 #ifdef DDDD
106   //####
107   //selectSubShape("FACE",  "Extrusion_1/LateralFace_3");
108   //selectSubShape("FACE",  "Extrusion_1/TopFace");
109   //selectSubShape("EDGE", "Extrusion_1/TopFace|Extrusion_1/LateralFace_1");
110   //selectSubShape("EDGE", "Sketch_1/Edge_6");
111 #endif
112   owner()->data()->sendAttributeUpdated(this);
113 }
114
115 std::shared_ptr<GeomAPI_Shape> Model_AttributeSelection::value()
116 {
117   std::shared_ptr<GeomAPI_Shape> aResult;
118   if (myIsInitialized) {
119     TDF_Label aSelLab = selectionLabel();
120     if (aSelLab.IsAttribute(kSIMPLE_REF_ID)) { // it is just reference to shape, not sub-shape
121       ResultPtr aContext = context();
122       if (!aContext.get()) 
123         return aResult; // empty result
124       return aContext->shape();
125     }
126
127     Handle(TNaming_NamedShape) aSelection;
128     if (selectionLabel().FindAttribute(TNaming_NamedShape::GetID(), aSelection)) {
129       TopoDS_Shape aSelShape = aSelection->Get();
130       aResult = std::shared_ptr<GeomAPI_Shape>(new GeomAPI_Shape);
131       aResult->setImpl(new TopoDS_Shape(aSelShape));
132     } else { // for simple construction element: just shape of this construction element
133       ResultConstructionPtr aConstr = 
134         std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(context());
135       if (aConstr) {
136         return aConstr->shape();
137       }
138     }
139   }
140   return aResult;
141 }
142
143 Model_AttributeSelection::Model_AttributeSelection(TDF_Label& theLabel)
144   : myRef(theLabel)
145 {
146   myIsInitialized = myRef.isInitialized();
147 }
148
149 ResultPtr Model_AttributeSelection::context() {
150   return std::dynamic_pointer_cast<ModelAPI_Result>(myRef.value());
151 }
152
153
154 void Model_AttributeSelection::setObject(const std::shared_ptr<ModelAPI_Object>& theObject)
155 {
156   ModelAPI_AttributeSelection::setObject(theObject);
157   myRef.setObject(theObject);
158 }
159
160 TDF_LabelMap& Model_AttributeSelection::scope()
161 {
162   if (myScope.IsEmpty()) { // create a new scope if not yet done
163     // gets all labels with named shapes that are bofore this feature label (before in history)
164     TDF_Label aFeatureLab = std::dynamic_pointer_cast<Model_Data>(
165       owner()->data())->label().Father();
166     int aFeatureID = aFeatureLab.Tag();
167     TDF_ChildIterator aFeaturesIter(aFeatureLab.Father());
168     for(; aFeaturesIter.More(); aFeaturesIter.Next()) {
169       if (aFeaturesIter.Value().Tag() >= aFeatureID) // the left labels are created later
170         break;
171       TDF_ChildIDIterator aNSIter(aFeaturesIter.Value(), TNaming_NamedShape::GetID(), 1);
172       for(; aNSIter.More(); aNSIter.Next()) {
173         Handle(TNaming_NamedShape) aNS = Handle(TNaming_NamedShape)::DownCast(aNSIter.Value());
174         if (!aNS.IsNull() && aNS->Evolution() != TNaming_SELECTED) {
175           myScope.Add(aNS->Label());
176         }
177       }
178     }
179   }
180   return myScope;
181 }
182
183 bool Model_AttributeSelection::update()
184 {
185   ResultPtr aContext = context();
186   if (!aContext.get()) return false;
187   TDF_Label aSelLab = selectionLabel();
188   if (aSelLab.IsAttribute(kSIMPLE_REF_ID)) { // it is just reference to shape, not sub-shape
189     return aContext->shape() && !aContext->shape()->isNull();
190   }
191
192   if (aContext->groupName() == ModelAPI_ResultBody::group()) {
193     // body: just a named shape, use selection mechanism from OCCT
194     TNaming_Selector aSelector(aSelLab);
195     bool aResult = aSelector.Solve(scope()) == Standard_True;
196     owner()->data()->sendAttributeUpdated(this);
197     return aResult;
198   } else if (aContext->groupName() == ModelAPI_ResultConstruction::group()) {
199     // construction: identification by the results indexes, recompute faces and
200     // take the face that more close by the indexes
201     std::shared_ptr<GeomAPI_PlanarEdges> aWirePtr = 
202       std::dynamic_pointer_cast<GeomAPI_PlanarEdges>(
203       std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(aContext)->shape());
204     if (aWirePtr && aWirePtr->hasPlane()) { // sketch sub-element
205       TDF_Label aLab = myRef.myRef->Label();
206       // getting a type of selected shape
207       Handle(TDataStd_Integer) aTypeAttr;
208       if (!aLab.FindAttribute(TDataStd_Integer::GetID(), aTypeAttr)) {
209         return false;
210       }
211       TopAbs_ShapeEnum aShapeType = (TopAbs_ShapeEnum)(aTypeAttr->Get());
212       // selected indexes will be needed in each "if"
213       Handle(TDataStd_IntPackedMap) aSubIds;
214       std::shared_ptr<GeomAPI_Shape> aNewSelected;
215       bool aNoIndexes = 
216         !aLab.FindAttribute(TDataStd_IntPackedMap::GetID(), aSubIds) || aSubIds->Extent() == 0;
217       // for now working only with composite features
218       FeaturePtr aContextFeature = aContext->document()->feature(aContext);
219       CompositeFeaturePtr aComposite = 
220         std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aContextFeature);
221       if (!aComposite || aComposite->numberOfSubs() == 0) {
222         return false;
223       }
224
225       if (aShapeType == TopAbs_FACE) {
226         // If this is a wire with plane defined thin it is a sketch-like object
227         std::list<std::shared_ptr<GeomAPI_Shape> > aFaces;
228         GeomAlgoAPI_SketchBuilder::createFaces(aWirePtr->origin(), aWirePtr->dirX(),
229           aWirePtr->dirY(), aWirePtr->norm(), aWirePtr, aFaces);
230         if (aFaces.empty()) // no faces, update can not work correctly
231           return false;
232         // if there is no edges indexes, any face can be used: take the first
233         std::shared_ptr<GeomAPI_Shape> aNewSelected;
234         if (aNoIndexes) {
235           aNewSelected = *(aFaces.begin());
236         } else { // searching for most looks-like initial face by the indexes
237           // prepare edges of the current resut for the fast searching
238           TColStd_MapOfTransient allCurves;
239           const int aSubNum = aComposite->numberOfSubs();
240           for(int a = 0; a < aSubNum; a++) {
241             if (aSubIds->Contains(aComposite->subFeatureId(a))) {
242               FeaturePtr aSub = aComposite->subFeature(a);
243               const std::list<std::shared_ptr<ModelAPI_Result> >& aResults = aSub->results();
244               std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aRes;
245               for(aRes = aResults.cbegin(); aRes != aResults.cend(); aRes++) {
246                 ResultConstructionPtr aConstr = 
247                   std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(*aRes);
248                 if (aConstr->shape() && aConstr->shape()->isEdge()) {
249                   const TopoDS_Shape& aResShape = aConstr->shape()->impl<TopoDS_Shape>();
250                   TopoDS_Edge anEdge = TopoDS::Edge(aResShape);
251                   if (!anEdge.IsNull()) {
252                     Standard_Real aFirst, aLast;
253                     Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast);
254                     allCurves.Add(aCurve);
255                   }
256                 }
257               }
258             }
259           }
260           // iterate new result faces and searching for these edges
261           std::list<std::shared_ptr<GeomAPI_Shape> >::iterator aFacesIter = aFaces.begin();
262           double aBestFound = 0; // best percentage of found edges
263           for(; aFacesIter != aFaces.end(); aFacesIter++) {
264             int aFound = 0, aNotFound = 0;
265             TopExp_Explorer anEdgesExp((*aFacesIter)->impl<TopoDS_Shape>(), TopAbs_EDGE);
266             for(; anEdgesExp.More(); anEdgesExp.Next()) {
267               TopoDS_Edge anEdge = TopoDS::Edge(anEdgesExp.Current());
268               if (!anEdge.IsNull()) {
269                 Standard_Real aFirst, aLast;
270                 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast);
271                 if (allCurves.Contains(aCurve)) {
272                   aFound++;
273                 } else {
274                   aNotFound++;
275                 }
276               }
277             }
278             if (aFound + aNotFound != 0) {
279               double aPercentage = double(aFound) / double(aFound + aNotFound);
280               if (aPercentage > aBestFound) {
281                 aBestFound = aPercentage;
282                 aNewSelected = *aFacesIter;
283               }
284             }
285           }
286         }
287         if (aNewSelected) { // store this new selection
288           selectConstruction(aContext, aNewSelected);
289           owner()->data()->sendAttributeUpdated(this);
290           return true;
291         }
292       } else if (aShapeType == TopAbs_EDGE) {
293         // just reselect the edge by the id
294         const int aSubNum = aComposite->numberOfSubs();
295         for(int a = 0; a < aSubNum; a++) {
296           // if aSubIds take any, the first appropriate
297           if (aSubIds->IsEmpty() || aSubIds->Contains(aComposite->subFeatureId(a))) {
298             // found the appropriate feature
299             FeaturePtr aFeature = aComposite->subFeature(a);
300             std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aResIter =
301               aFeature->results().cbegin();
302             for(;aResIter != aFeature->results().cend(); aResIter++) {
303               ResultConstructionPtr aRes = 
304                 std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(*aResIter);
305               if (aRes && aRes->shape() && aRes->shape()->isEdge()) { // found!
306                 selectConstruction(aContext, aRes->shape());
307                 owner()->data()->sendAttributeUpdated(this);
308                 return true;
309               }
310             }
311           }
312         }
313       } else if (aShapeType == TopAbs_VERTEX) {
314         // just reselect the vertex by the id of edge
315         const int aSubNum = aComposite->numberOfSubs();
316         for(int a = 0; a < aSubNum; a++) {
317           // if aSubIds take any, the first appropriate
318           int aFeatureID = aComposite->subFeatureId(a);
319           if (aSubIds->IsEmpty() || aSubIds->Contains(aFeatureID) ||
320               aSubIds->Contains(aFeatureID + kSTART_VERTEX_DELTA) ||
321               aSubIds->Contains(aFeatureID + kSTART_VERTEX_DELTA * 2)) {
322             // searching for deltas
323             int aVertexNum = 0;
324             if (aSubIds->Contains(aFeatureID + kSTART_VERTEX_DELTA)) aVertexNum = 1;
325             else if (aSubIds->Contains(aFeatureID + kSTART_VERTEX_DELTA * 2)) aVertexNum = 2;
326             // found the feature with appropriate edge
327             FeaturePtr aFeature = aComposite->subFeature(a);
328             std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aResIter =
329               aFeature->results().cbegin();
330             for(;aResIter != aFeature->results().cend(); aResIter++) {
331               ResultConstructionPtr aRes = 
332                 std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(*aResIter);
333               if (aRes && aRes->shape()) {
334                 if (aRes->shape()->isVertex() && aVertexNum == 0) { // found!
335                   selectConstruction(aContext, aRes->shape());
336                   owner()->data()->sendAttributeUpdated(this);
337                   return true;
338                 } else if (aRes->shape()->isEdge() && aVertexNum > 0) {
339                   const TopoDS_Shape& anEdge = aRes->shape()->impl<TopoDS_Shape>();
340                   int aVIndex = 1;
341                   for(TopExp_Explorer aVExp(anEdge, TopAbs_VERTEX); aVExp.More(); aVExp.Next()) {
342                     if (aVIndex == aVertexNum) { // found!
343                       std::shared_ptr<GeomAPI_Shape> aVertex(new GeomAPI_Shape);
344                       aVertex->setImpl(new TopoDS_Shape(aVExp.Current()));
345                       selectConstruction(aContext, aVertex);
346                       owner()->data()->sendAttributeUpdated(this);
347                       return true;
348                     }
349                     aVIndex++;
350                   }
351                 }
352               }
353             }
354           }
355         }
356       }
357     } else { // simple construction element: the selected is that needed
358       selectConstruction(aContext, aContext->shape());
359       owner()->data()->sendAttributeUpdated(this);
360       return true;
361     }
362   }
363   return false; // unknown case
364 }
365
366
367 void Model_AttributeSelection::selectBody(
368     const ResultPtr& theContext, const std::shared_ptr<GeomAPI_Shape>& theSubShape)
369 {
370   // perform the selection
371   TNaming_Selector aSel(selectionLabel());
372   TopoDS_Shape aNewShape = theSubShape ? theSubShape->impl<TopoDS_Shape>() : TopoDS_Shape();
373   TopoDS_Shape aContext;
374
375   ResultBodyPtr aBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(myRef.value());
376   if (aBody) {
377     aContext = aBody->shape()->impl<TopoDS_Shape>();
378   } else {
379     ResultConstructionPtr aConstr = std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(myRef.value());
380     if (aConstr) {
381       aContext = aConstr->shape()->impl<TopoDS_Shape>();
382     } else {
383       Events_Error::send("A result with shape is expected");
384       return;
385     }
386   }
387   /// fix for issue 411: result modified shapes must not participate in this selection mechanism
388   FeaturePtr aFeatureOwner = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
389   if (aFeatureOwner.get())
390     aFeatureOwner->eraseResults();
391
392   aSel.Select(aNewShape, aContext); 
393 }
394
395 /// registers the name of the shape in the label (theID == 0) of sub label (theID is a tag)
396 /// if theID is zero, 
397 static void registerSubShape(TDF_Label theMainLabel, TopoDS_Shape theShape,
398   const int theID, const FeaturePtr& theContextFeature, std::shared_ptr<Model_Document> theDoc,
399   std::string theAdditionalName,
400   Handle(TDataStd_IntPackedMap) theRefs = Handle(TDataStd_IntPackedMap)())
401 {
402   TDF_Label aLab = theID == 0 ? theMainLabel : theMainLabel.FindChild(theID);
403   TNaming_Builder aBuilder(aLab);
404   aBuilder.Generated(theShape);
405   std::stringstream aName;
406   aName<<theContextFeature->name()<<"/";
407   if (!theAdditionalName.empty())
408     aName<<theAdditionalName<<"/";
409   if (theShape.ShapeType() == TopAbs_FACE) aName<<"Face";
410   else if (theShape.ShapeType() == TopAbs_EDGE) aName<<"Edge";
411   else if (theShape.ShapeType() == TopAbs_VERTEX) aName<<"Vertex";
412
413   if (theRefs.IsNull()) {
414     aName<<theID;
415   } else { // make a compisite name from all sub-elements indexes: "1_2_3_4"
416     TColStd_MapIteratorOfPackedMapOfInteger aRef(theRefs->GetMap());
417     for(; aRef.More(); aRef.Next()) {
418       aName<<"-"<<aRef.Key();
419     }
420   }
421
422   theDoc->addNamingName(aLab, aName.str());
423   TDataStd_Name::Set(aLab, aName.str().c_str());
424 }
425
426 void Model_AttributeSelection::selectConstruction(
427     const ResultPtr& theContext, const std::shared_ptr<GeomAPI_Shape>& theSubShape)
428 {
429   std::shared_ptr<Model_Document> aMyDoc = 
430     std::dynamic_pointer_cast<Model_Document>(owner()->document());
431   FeaturePtr aContextFeature = theContext->document()->feature(theContext);
432   CompositeFeaturePtr aComposite = 
433     std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aContextFeature);
434   const TopoDS_Shape& aSubShape = theSubShape->impl<TopoDS_Shape>();
435   if (!aComposite || aComposite->numberOfSubs() == 0) {
436     // saving of context is enough: result construction contains exactly the needed shape
437     TNaming_Builder aBuilder(selectionLabel());
438     aBuilder.Generated(aSubShape);
439     aMyDoc->addNamingName(selectionLabel(), theContext->data()->name());
440     TDataStd_Name::Set(selectionLabel(), theContext->data()->name().c_str());
441     return;
442   }
443   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(owner()->data());
444   TDF_Label aLab = myRef.myRef->Label();
445   // identify the reuslts of sub-object of the composite by edges
446   // save type of the selected shape in integer attribute
447   TopAbs_ShapeEnum aShapeType = aSubShape.ShapeType();
448   TDataStd_Integer::Set(aLab, (int)aShapeType);
449   gp_Pnt aVertexPos;
450   TColStd_MapOfTransient allCurves;
451   if (aShapeType == TopAbs_VERTEX) { // compare positions
452     aVertexPos = BRep_Tool::Pnt(TopoDS::Vertex(aSubShape));
453   } else { 
454     for(TopExp_Explorer anEdgeExp(aSubShape, TopAbs_EDGE); anEdgeExp.More(); anEdgeExp.Next()) {
455       TopoDS_Edge anEdge = TopoDS::Edge(anEdgeExp.Current());
456       Standard_Real aFirst, aLast;
457       Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast);
458       allCurves.Add(aCurve);
459     }
460   }
461   // iterate and store the result ids of sub-elements and sub-elements to sub-labels
462   Handle(TDataStd_IntPackedMap) aRefs = TDataStd_IntPackedMap::Set(aLab);
463   aRefs->Clear();
464   const int aSubNum = aComposite->numberOfSubs();
465   for(int a = 0; a < aSubNum; a++) {
466     FeaturePtr aSub = aComposite->subFeature(a);
467     const std::list<std::shared_ptr<ModelAPI_Result> >& aResults = aSub->results();
468     std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aRes = aResults.cbegin();
469     // there may be many shapes (circle and center): register if at least one is in selection
470     for(; aRes != aResults.cend(); aRes++) {
471       ResultConstructionPtr aConstr = 
472         std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(*aRes);
473       if (!aConstr->shape()) {
474         continue;
475       }
476       if (aShapeType == TopAbs_VERTEX) {
477         if (aConstr->shape()->isVertex()) { // compare vertices positions
478           const TopoDS_Shape& aVertex = aConstr->shape()->impl<TopoDS_Shape>();
479           gp_Pnt aPnt = BRep_Tool::Pnt(TopoDS::Vertex(aVertex));
480           if (aPnt.IsEqual(aVertexPos, Precision::Confusion())) {
481             aRefs->Add(aComposite->subFeatureId(a));
482           }
483         } else { // get first or last vertex of the edge: last is stored with negative sign
484           const TopoDS_Shape& anEdge = aConstr->shape()->impl<TopoDS_Shape>();
485           int aDelta = kSTART_VERTEX_DELTA;
486           for(TopExp_Explorer aVExp(anEdge, TopAbs_VERTEX); aVExp.More(); aVExp.Next()) {
487             gp_Pnt aPnt = BRep_Tool::Pnt(TopoDS::Vertex(aVExp.Current()));
488             if (aPnt.IsEqual(aVertexPos, Precision::Confusion())) {
489               aRefs->Add(aDelta + aComposite->subFeatureId(a));
490               break;
491             }
492             aDelta += kSTART_VERTEX_DELTA;
493           }
494         }
495       } else {
496         if (aConstr->shape()->isEdge()) {
497           const TopoDS_Shape& aResShape = aConstr->shape()->impl<TopoDS_Shape>();
498           TopoDS_Edge anEdge = TopoDS::Edge(aResShape);
499           if (!anEdge.IsNull()) {
500             Standard_Real aFirst, aLast;
501             Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast);
502             if (allCurves.Contains(aCurve)) {
503               int anID = aComposite->subFeatureId(a);
504               aRefs->Add(anID);
505               if (aShapeType != TopAbs_EDGE) { // edge do not need the sub-edges on sub-labels
506                 // add edges to sub-label to support naming for edges selection
507                 TopExp_Explorer anEdgeExp(aSubShape, TopAbs_EDGE);
508                 for(; anEdgeExp.More(); anEdgeExp.Next()) {
509                   TopoDS_Edge anEdge = TopoDS::Edge(anEdgeExp.Current());
510                   Standard_Real aFirst, aLast;
511                   Handle(Geom_Curve) aFaceCurve = BRep_Tool::Curve(anEdge, aFirst, aLast);
512                   if (aFaceCurve == aCurve) {
513                     registerSubShape(selectionLabel(), anEdge, anID, aContextFeature, aMyDoc, "");
514                   }
515                 }
516               } else { // put vertices of the selected edge to sub-labels
517                 // add edges to sub-label to support naming for edges selection
518                 TopExp_Explorer anEdgeExp(aSubShape, TopAbs_VERTEX);
519                 int aTagIndex = anID + kSTART_VERTEX_DELTA;
520                 for(; anEdgeExp.More(); anEdgeExp.Next(), aTagIndex += kSTART_VERTEX_DELTA) {
521                   TopoDS_Vertex aV = TopoDS::Vertex(anEdgeExp.Current());
522
523                   std::stringstream anAdditionalName; 
524                   registerSubShape(selectionLabel(), aV, aTagIndex, aContextFeature, aMyDoc,
525                     "");
526                 }
527               }
528             }
529           }
530         }
531       }
532     }
533   }
534   // store the selected as primitive
535   TNaming_Builder aBuilder(selectionLabel());
536   aBuilder.Generated(aSubShape);
537   registerSubShape(selectionLabel(), aSubShape, 0, aContextFeature, aMyDoc, "", aRefs);
538 }
539
540 TDF_Label Model_AttributeSelection::selectionLabel()
541 {
542   return myRef.myRef->Label().FindChild(1);
543 }
544
545 #define FIX_BUG1 1
546 std::string GetShapeName(std::shared_ptr<Model_Document> theDoc, const TopoDS_Shape& theShape, 
547                                              const TDF_Label& theLabel)
548 {
549   std::string aName;
550   // check if the subShape is already in DF
551   Handle(TNaming_NamedShape) aNS = TNaming_Tool::NamedShape(theShape, theLabel);
552   Handle(TDataStd_Name) anAttr;
553   if(!aNS.IsNull() && !aNS->IsEmpty()) { // in the document    
554         if(aNS->Label().FindAttribute(TDataStd_Name::GetID(), anAttr)) {
555           aName = TCollection_AsciiString(anAttr->Get()).ToCString();
556           if(!aName.empty()) {      
557             const TDF_Label& aLabel = theDoc->findNamingName(aName);
558     /* MPV: the same shape with the same name may be duplicated on different labels (selection of the same construction object)
559                 if(!aLabel.IsEqual(aNS->Label())) {
560                   //aName.erase(); //something is wrong, to be checked!!!
561                   aName += "_SomethingWrong";
562                   return aName;
563                 }*/
564                 const TopoDS_Shape& aShape = aNS->Get();
565                 if(aShape.ShapeType() == TopAbs_COMPOUND) {
566                   std::string aPostFix("_");
567                   TopoDS_Iterator it(aShape);                   
568                   for (int i = 1;it.More();it.Next(), i++) {
569                     if(it.Value() == theShape) {
570                           aPostFix += TCollection_AsciiString (i).ToCString();
571                           aName    += aPostFix;
572                           break;
573                         }
574                         else continue;                                          
575                   }
576                 }
577           }     
578         }
579   }
580   return aName;
581 }
582
583 bool isTrivial (const TopTools_ListOfShape& theAncestors, TopTools_IndexedMapOfShape& theSMap)
584 {
585 // a trivial case: F1 & F2,  aNumber = 1, i.e. intersection gives 1 edge.
586   TopoDS_Compound aCmp;
587   BRep_Builder BB;
588   BB.MakeCompound(aCmp);
589   TopTools_ListIteratorOfListOfShape it(theAncestors);
590   for(;it.More();it.Next()) {
591         BB.Add(aCmp, it.Value());
592         theSMap.Add(it.Value());
593   }
594   int aNumber(0);
595   TopTools_IndexedDataMapOfShapeListOfShape aMap2;
596   TopExp::MapShapesAndAncestors(aCmp, TopAbs_EDGE, TopAbs_FACE, aMap2);
597   for (int i = 1; i <= aMap2.Extent(); i++) {
598         const TopoDS_Shape& aKey = aMap2.FindKey(i);
599         const TopTools_ListOfShape& anAncestors = aMap2.FindFromIndex(i);
600         if(anAncestors.Extent() > 1) aNumber++;
601   }
602   if(aNumber > 1) return false;
603   return true;
604 }
605 std::string Model_AttributeSelection::namingName()
606 {
607   std::string aName("");
608   if(!this->isInitialized()) return aName;
609   Handle(TDataStd_Name) anAtt;
610   if(selectionLabel().FindAttribute(TDataStd_Name::GetID(), anAtt)) {
611         aName = TCollection_AsciiString(anAtt->Get()).ToCString();
612         return aName;
613   }
614
615   std::shared_ptr<GeomAPI_Shape> aSubSh = value();
616   ResultPtr aCont = context();
617   aName = "Undefined name";
618   if(!aSubSh.get() || aSubSh->isNull() || !aCont.get() || aCont->shape()->isNull()) 
619     return aName;
620   TopoDS_Shape aSubShape = aSubSh->impl<TopoDS_Shape>();
621   TopoDS_Shape aContext  = aCont->shape()->impl<TopoDS_Shape>();
622 #ifdef DEB_NAMING
623   BRepTools::Write(aSubShape, "Selection.brep");
624   BRepTools::Write(aContext, "Context.brep");
625 #endif
626   std::shared_ptr<Model_Document> aDoc = 
627     std::dynamic_pointer_cast<Model_Document>(aCont->document());
628
629   // check if the subShape is already in DF
630   aName = GetShapeName(aDoc, aSubShape, selectionLabel());
631   if(aName.empty() ) { // not in the document!
632     TopAbs_ShapeEnum aType = aSubShape.ShapeType();
633         switch (aType) {
634           case TopAbs_FACE:
635       // the Face should be in DF. If it is not the case, it is an error ==> to be debugged             
636                 break;
637           case TopAbs_EDGE:
638                   {
639                   // name structure: F1 | F2 [| F3 | F4], where F1 & F2 the faces which gives the Edge in trivial case
640                   // if it is not atrivial case we use localization by neighbours. F3 & F4 - neighbour faces             
641                   TopTools_IndexedDataMapOfShapeListOfShape aMap;
642                   TopExp::MapShapesAndAncestors(aContext, TopAbs_EDGE, TopAbs_FACE, aMap);
643                   TopTools_IndexedMapOfShape aSMap; // map for ancestors of the sub-shape
644                   bool isTrivialCase(true);
645 /*                for (int i = 1; i <= aMap.Extent(); i++) {
646                         const TopoDS_Shape& aKey = aMap.FindKey(i);
647                         //if (aKey.IsNotEqual(aSubShape)) continue; // find exactly the selected key
648                         if (aKey.IsSame(aSubShape)) continue;
649             const TopTools_ListOfShape& anAncestors = aMap.FindFromIndex(i);
650                         // check that it is not a trivial case (F1 & F2: aNumber = 1)
651                         isTrivialCase = isTrivial(anAncestors, aSMap);                  
652                         break;
653                   }*/
654                   if(aMap.Contains(aSubShape)) {
655                         const TopTools_ListOfShape& anAncestors = aMap.FindFromKey(aSubShape);
656                         // check that it is not a trivial case (F1 & F2: aNumber = 1)
657                         isTrivialCase = isTrivial(anAncestors, aSMap);          
658                   } else 
659                           break;
660                   TopTools_ListOfShape aListOfNbs;
661                   if(!isTrivialCase) { // find Neighbors
662                         TNaming_Localizer aLocalizer;
663                         TopTools_MapOfShape aMap3;
664                         aLocalizer.FindNeighbourg(aContext, aSubShape, aMap3);
665                         //int n = aMap3.Extent();
666                         TopTools_MapIteratorOfMapOfShape it(aMap3);
667                         for(;it.More();it.Next()) {
668                           const TopoDS_Shape& aNbShape = it.Key(); // neighbor edge
669                           //TopAbs_ShapeEnum aType = aNbShape.ShapeType();
670                           const TopTools_ListOfShape& aList  = aMap.FindFromKey(aNbShape);
671                           TopTools_ListIteratorOfListOfShape it2(aList);
672                           for(;it2.More();it2.Next()) {
673                                 if(aSMap.Contains(it2.Value())) continue; // skip this Face
674                                 aListOfNbs.Append(it2.Value());
675                           }
676                         }
677                   }  // else a trivial case
678                   
679                   // build name of the sub-shape Edge
680                   for(int i=1; i <= aSMap.Extent(); i++) {
681                         const TopoDS_Shape& aFace = aSMap.FindKey(i);
682                         std::string aFaceName = GetShapeName(aDoc, aFace, selectionLabel());
683                         if(i == 1)
684                           aName = aFaceName;
685                         else 
686                           aName += "|" + aFaceName;
687                   }
688                   TopTools_ListIteratorOfListOfShape itl(aListOfNbs);
689                   for (;itl.More();itl.Next()) {
690                         std::string aFaceName = GetShapeName(aDoc, itl.Value(), selectionLabel());
691                         aName += "|" + aFaceName;
692                   }               
693                   }
694                   break;
695
696           case TopAbs_VERTEX:
697                   // name structure (Monifold Topology): 
698                   // 1) F1 | F2 | F3 - intersection of 3 faces defines a vertex - trivial case.
699                   // 2) F1 | F2 | F3 [|F4 [|Fn]] - redundant definition, but it should be kept as is to obtain safe recomputation
700                   // 2) F1 | F2      - intersection of 2 faces definses a vertex - applicable for case
701                   //                   when 1 faces is cylindrical, conical, spherical or revolution and etc.
702                   // 3) E1 | E2 | E3 - intersection of 3 edges defines a vertex - when we have case of a shell
703                   //                   or compound of 2 open faces.
704                   // 4) E1 | E2      - intesection of 2 edges defines a vertex - when we have a case of 
705                   //                   two independent edges (wire or compound)
706                   // implemented 2 first cases
707                   {
708                         TopTools_IndexedDataMapOfShapeListOfShape aMap;
709                     TopExp::MapShapesAndAncestors(aContext, TopAbs_VERTEX, TopAbs_FACE, aMap);
710                         const TopTools_ListOfShape& aList2  = aMap.FindFromKey(aSubShape);
711                         TopTools_ListOfShape aList;
712                         TopTools_MapOfShape aFMap;
713 #ifdef FIX_BUG1
714                         //int n = aList2.Extent(); //bug! duplication
715                         // fix is below
716                         TopTools_ListIteratorOfListOfShape itl2(aList2);
717                     for (int i = 1;itl2.More();itl2.Next(),i++) {
718                           if(aFMap.Add(itl2.Value()))
719                                 aList.Append(itl2.Value());
720                         }
721                         //n = aList.Extent();
722 #endif
723                         int n = aList.Extent();
724                         if(n < 3) { // open topology case => via edges
725                           TopTools_IndexedDataMapOfShapeListOfShape aMap;
726                       TopExp::MapShapesAndAncestors(aContext, TopAbs_VERTEX, TopAbs_EDGE, aMap);
727                           const TopTools_ListOfShape& aList2  = aMap.FindFromKey(aSubShape);
728                           if(aList2.Extent() >= 2)  { // regular solution
729                                 TopTools_ListIteratorOfListOfShape itl(aList2);
730                         for (int i = 1;itl.More();itl.Next(),i++) {
731                               const TopoDS_Shape& anEdge = itl.Value();
732                               std::string anEdgeName = GetShapeName(aDoc, anEdge, selectionLabel());
733                               if(i == 1)
734                                 aName = anEdgeName;
735                               else 
736                                 aName += "|" + anEdgeName;
737                                 }
738                           }
739                         } 
740                         else {
741                           TopTools_ListIteratorOfListOfShape itl(aList);
742                       for (int i = 1;itl.More();itl.Next(),i++) {
743                             const TopoDS_Shape& aFace = itl.Value();
744                             std::string aFaceName = GetShapeName(aDoc, aFace, selectionLabel());
745                             if(i == 1)
746                               aName = aFaceName;
747                             else 
748                               aName += "|" + aFaceName;
749                           }
750                         }
751                   }
752                   break;
753         }
754     // register name                    
755     // aDoc->addNamingName(selectionLabel(), aName);
756         // the selected sub-shape will not be shared and as result it will not require registration
757   }
758   return aName;
759 }
760
761 TopAbs_ShapeEnum translateType (const std::string& theType)
762 {
763   TCollection_AsciiString aStr(theType.c_str());
764   aStr.UpperCase();
765   if(aStr.IsEqual("COMP"))
766         return TopAbs_COMPOUND;
767   else if(aStr.IsEqual("COMS"))
768         return TopAbs_COMPSOLID;
769   else if(aStr.IsEqual("SOLD"))
770         return TopAbs_SOLID;
771   else if(aStr.IsEqual("SHEL"))
772         return TopAbs_SHELL;
773   else if(aStr.IsEqual("FACE"))
774         return TopAbs_FACE;
775   else if(aStr.IsEqual("WIRE"))
776         return TopAbs_WIRE;
777   else if(aStr.IsEqual("EDGE"))
778         return TopAbs_EDGE;
779   else if(aStr.IsEqual("VERT"))
780         return TopAbs_VERTEX;  
781
782   return TopAbs_SHAPE;
783 }
784
785 const TopoDS_Shape getShapeFromCompound(const std::string& theSubShapeName, const TopoDS_Shape& theCompound)
786 {
787   TopoDS_Shape aSelection;
788   std::string::size_type n = theSubShapeName.rfind('/');                        
789   if (n == std::string::npos) n = 0;
790         std::string aSubString = theSubShapeName.substr(n + 1);
791         n = aSubString.rfind('_');
792         if (n == std::string::npos) return aSelection;
793         aSubString = aSubString.substr(n+1);
794         int indx = atoi(aSubString.c_str());
795         TopoDS_Iterator it(theCompound);                        
796         for (int i = 1;it.More();it.Next(), i++) {
797           if(i == indx) {                         
798             aSelection = it.Value();                      
799             break;
800           }
801           else continue;                                                
802         }
803   return aSelection;    
804 }
805
806 const TopoDS_Shape findFaceByName(const std::string& theSubShapeName, std::shared_ptr<Model_Document> theDoc)
807 {
808   TopoDS_Shape aFace;
809   std::string::size_type n, nb = theSubShapeName.rfind('/');                    
810   if (nb == std::string::npos) nb = 0;
811   std::string aSubString = theSubShapeName.substr(nb + 1);
812   n = aSubString.rfind('_');
813   if (n != std::string::npos) {
814           std::string aSubStr2 = aSubString.substr(0, n);
815         aSubString  = theSubShapeName.substr(0, nb + 1);
816         aSubString = aSubString + aSubStr2;     
817   } else
818         aSubString = theSubShapeName;
819                                 
820   const TDF_Label& aLabel = theDoc->findNamingName(aSubString);
821   if(aLabel.IsNull()) return aFace;
822   Handle(TNaming_NamedShape) aNS;
823   if(aLabel.FindAttribute(TNaming_NamedShape::GetID(), aNS)) {
824         const TopoDS_Shape& aShape = aNS->Get();
825         if(!aShape.IsNull()) {
826           if(aShape.ShapeType() == TopAbs_COMPOUND) 
827             aFace = getShapeFromCompound(theSubShapeName, aShape);
828           else 
829                 aFace = aShape;   
830         }
831   }
832   return aFace;
833 }
834
835 int ParseName(const std::string& theSubShapeName,   std::list<std::string>& theList)
836 {
837   std::string aName = theSubShapeName;
838   std::string aLastName;
839   int n1(0), n2(0); // n1 - start position, n2 - position of the delimiter
840   while ((n2 = aName.find('|', n1)) != std::string::npos) {
841     const std::string aName1 = aName.substr(n1, n2 - n1); //name of face
842         theList.push_back(aName1);      
843         n1 = n2 + 1;
844         aLastName = aName.substr(n1);
845   }
846   if(!aLastName.empty())
847         theList.push_back(aLastName);
848   return theList.size();
849 }
850
851 const TopoDS_Shape findCommonShape(const TopAbs_ShapeEnum theType, const TopTools_ListOfShape& theList)
852 {
853   TopoDS_Shape aShape;
854   std::vector<TopTools_MapOfShape> aVec;
855   TopTools_MapOfShape aMap1, aMap2, aMap3, aMap4;
856   if(theList.Extent() > 1) {
857         aVec.push_back(aMap1);
858         aVec.push_back(aMap2);
859   }
860   if(theList.Extent() > 2)
861     aVec.push_back(aMap3);
862   if(theList.Extent() == 4)
863     aVec.push_back(aMap4);
864
865   //fill maps
866   TopTools_ListIteratorOfListOfShape it(theList);
867   for(int i = 0;it.More();it.Next(),i++) {
868         const TopoDS_Shape& aFace = it.Value();         
869         if(i < 2) {
870           TopExp_Explorer anExp (aFace, theType);
871           for(;anExp.More();anExp.Next()) {
872             const TopoDS_Shape& anEdge = anExp.Current();
873             if (!anEdge.IsNull())
874                   aVec[i].Add(anExp.Current());
875           }
876         } else {
877           TopExp_Explorer anExp (aFace, TopAbs_VERTEX);
878           for(;anExp.More();anExp.Next()) {
879             const TopoDS_Shape& aVertex = anExp.Current();
880             if (!aVertex.IsNull())
881                   aVec[i].Add(anExp.Current());
882           }
883         }
884   }
885   //trivial case: 2 faces
886   TopTools_ListOfShape aList;
887   TopTools_MapIteratorOfMapOfShape it2(aVec[0]);
888   for(;it2.More();it2.Next()) {
889         if(aVec[1].Contains(it2.Key())) {
890           aShape = it2.Key();
891           if(theList.Extent() == 2)
892             break;
893           else 
894                 aList.Append(it2.Key());
895         }
896   }
897   if(aList.Extent() && aVec.size() > 3) {// list of common edges ==> search ny neighbors 
898         if(aVec[2].Extent() && aVec[3].Extent()) {
899           TopTools_ListIteratorOfListOfShape it(aList);
900           for(;it.More();it.Next()) {
901                 const TopoDS_Shape& aCand = it.Value();
902                 // not yet completelly implemented, to be rechecked
903                 TopoDS_Vertex aV1, aV2;
904                 TopExp::Vertices(TopoDS::Edge(aCand), aV1, aV2);
905                 int aNum(0);
906                 if(aVec[2].Contains(aV1)) aNum++;
907                 else if(aVec[2].Contains(aV2)) aNum++;
908                 if(aVec[3].Contains(aV1)) aNum++;
909                 else if(aVec[3].Contains(aV2)) aNum++;
910                 if(aNum == 2) {
911                   aShape = aCand;
912                   break;
913                 }
914           }
915         }
916   }
917
918   if(aList.Extent() && aVec.size() == 3) {
919
920     TopTools_ListIteratorOfListOfShape it(aList);
921           for(;it.More();it.Next()) {
922                 const TopoDS_Shape& aCand = it.Value();
923                 if(aVec[2].Contains(aCand)) {
924                   aShape = aCand;
925                   break;
926                 }
927           }
928   }
929   return aShape;
930 }
931
932 std::string getContextName(const std::string& theSubShapeName)
933 {
934   std::string aName;
935   std::string::size_type n = theSubShapeName.find('/');                 
936   if (n == std::string::npos) return aName;
937   aName = theSubShapeName.substr(0, n);
938   return aName;
939 }
940 // type ::= COMP | COMS | SOLD | SHEL | FACE | WIRE | EDGE | VERT
941 void Model_AttributeSelection::selectSubShape(const std::string& theType, const std::string& theSubShapeName)
942 {
943   if(theSubShapeName.empty() || theType.empty()) return;
944   TopAbs_ShapeEnum aType = translateType(theType);
945   std::shared_ptr<Model_Document> aDoc =  std::dynamic_pointer_cast<Model_Document>(owner()->document());//std::dynamic_pointer_cast<Model_Document>(aCont->document());
946   std::string aContName = getContextName(theSubShapeName);
947   if(aContName.empty()) return;
948   //ResultPtr aCont = context();
949   ResultPtr aCont = aDoc->findByName(aContName);
950   if(!aCont.get() || aCont->shape()->isNull()) return;
951   TopoDS_Shape aContext  = aCont->shape()->impl<TopoDS_Shape>();
952   TopAbs_ShapeEnum aContType = aContext.ShapeType();
953   if(aType <= aContType) return; // not applicable
954  
955   
956   TopoDS_Shape aSelection;
957   switch (aType) 
958   {
959   case TopAbs_COMPOUND:
960     break;
961   case TopAbs_COMPSOLID:
962     break;
963   case TopAbs_SOLID:
964     break;
965   case TopAbs_SHELL:
966     break;
967   case TopAbs_FACE:
968     {
969           const TopoDS_Shape aSelection = findFaceByName(theSubShapeName, aDoc);
970           if(!aSelection.IsNull()) {// Select it
971                 std::shared_ptr<GeomAPI_Shape> aShapeToBeSelected(new GeomAPI_Shape());
972         aShapeToBeSelected->setImpl(new TopoDS_Shape(aSelection));
973                 setValue(aCont, aShapeToBeSelected);
974           }
975         }
976     break;
977   case TopAbs_WIRE:
978         break;
979   case TopAbs_EDGE:
980         {  
981           TopoDS_Shape aSelection;// = findFaceByName(theSubShapeName, aDoc);
982           const TDF_Label& aLabel = aDoc->findNamingName(theSubShapeName);
983       if(!aLabel.IsNull()) {
984         Handle(TNaming_NamedShape) aNS;
985         if(aLabel.FindAttribute(TNaming_NamedShape::GetID(), aNS)) {
986               const TopoDS_Shape& aShape = aNS->Get();
987               if(!aShape.IsNull()) {
988                 if(aShape.ShapeType() == TopAbs_COMPOUND) 
989                   aSelection = getShapeFromCompound(theSubShapeName, aShape);
990                         else 
991                     aSelection = aShape;          
992                   }
993                 }
994           }
995           if(aSelection.IsNull()) {
996           std::list<std::string> aListofNames;
997           int n = ParseName(theSubShapeName, aListofNames);
998           if(n > 1 && n < 5) {
999                 TopTools_ListOfShape aList;
1000         for(std::list<std::string>::iterator it =aListofNames.begin();it != aListofNames.end();it++){
1001                   const TopoDS_Shape aFace = findFaceByName(*it, aDoc);
1002                   aList.Append(aFace);          
1003                 }
1004                 aSelection = findCommonShape(TopAbs_EDGE, aList);
1005           }
1006           }
1007           if(!aSelection.IsNull()) {// Select it
1008                 std::shared_ptr<GeomAPI_Shape> aShapeToBeSelected(new GeomAPI_Shape());
1009         aShapeToBeSelected->setImpl(new TopoDS_Shape(aSelection));
1010                 setValue(aCont, aShapeToBeSelected);
1011           }
1012         }
1013     break;
1014   case TopAbs_VERTEX:
1015           {
1016           TopoDS_Shape aSelection;
1017           const TDF_Label& aLabel = aDoc->findNamingName(theSubShapeName);
1018       if(!aLabel.IsNull()) {
1019         Handle(TNaming_NamedShape) aNS;
1020         if(aLabel.FindAttribute(TNaming_NamedShape::GetID(), aNS)) {
1021               const TopoDS_Shape& aShape = aNS->Get();
1022               if(!aShape.IsNull()) {
1023                 if(aShape.ShapeType() == TopAbs_COMPOUND) 
1024                   aSelection = getShapeFromCompound(theSubShapeName, aShape);
1025                         else 
1026                     aSelection = aShape;          
1027                   }
1028                 }
1029           }
1030           if(aSelection.IsNull()) {
1031             std::list<std::string> aListofNames;
1032                 int n = ParseName(theSubShapeName, aListofNames);
1033             if(n > 1 && n < 4) { // 2 || 3
1034                 TopTools_ListOfShape aList;
1035           for(std::list<std::string>::iterator it =aListofNames.begin();it != aListofNames.end();it++){
1036                     const TopoDS_Shape aFace = findFaceByName(*it, aDoc);
1037                         if(!aFace.IsNull())
1038                       aList.Append(aFace);              
1039                   }
1040                  aSelection = findCommonShape(TopAbs_VERTEX, aList);
1041                 }
1042           }
1043           if(!aSelection.IsNull()) {// Select it
1044                 std::shared_ptr<GeomAPI_Shape> aShapeToBeSelected(new GeomAPI_Shape());
1045         aShapeToBeSelected->setImpl(new TopoDS_Shape(aSelection));
1046                 setValue(aCont, aShapeToBeSelected);
1047           }
1048           }
1049         break;
1050   default: //TopAbs_SHAPE
1051         return;
1052   }
1053
1054 }
1055
1056 int Model_AttributeSelection::Id()
1057 {
1058   std::shared_ptr<GeomAPI_Shape> aSelection = value();
1059   std::shared_ptr<GeomAPI_Shape> aContext = context()->shape();
1060   const TopoDS_Shape& aMainShape = aContext->impl<TopoDS_Shape>();
1061   const TopoDS_Shape& aSubShape = aSelection->impl<TopoDS_Shape>();
1062   int anID = 0;
1063   if (aSelection && !aSelection->isNull() &&
1064       aContext   && !aContext->isNull())
1065   {
1066     TopTools_IndexedMapOfShape aSubShapesMap;
1067     TopExp::MapShapes(aMainShape, aSubShapesMap);
1068     anID = aSubShapesMap.FindIndex(aSubShape);
1069   }
1070   return anID;
1071 }