src/Model/Model_AttributeSelection.cpp

   1 // Copyright (C) 2014-2019  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19
  20 #include "Model_AttributeSelection.h"
  21 #include "Model_Application.h"
  22 #include "Model_Events.h"
  23 #include "Model_Data.h"
  24 #include "Model_Document.h"
  25 #include <Model_Objects.h>
  26 #include <Model_AttributeSelectionList.h>
  27 #include <Model_ResultConstruction.h>
  28 #include <ModelAPI_Feature.h>
  29 #include <ModelAPI_ResultBody.h>
  30 #include <ModelAPI_ResultBody.h>
  31 #include <ModelAPI_ResultConstruction.h>
  32 #include <ModelAPI_ResultGroup.h>
  33 #include <ModelAPI_ResultPart.h>
  34 #include <ModelAPI_CompositeFeature.h>
  35 #include <ModelAPI_Tools.h>
  36 #include <ModelAPI_Session.h>
  37 #include <ModelAPI_Validator.h>
  38 #include <ModelGeomAlgo_Shape.h>
  39 #include <Events_InfoMessage.h>
  40 #include <GeomAPI_Edge.h>
  41 #include <GeomAPI_Pnt.h>
  42 #include <GeomAPI_ShapeIterator.h>
  43 #include <GeomAPI_Vertex.h>
  44 #include <GeomAPI_ShapeExplorer.h>
  45 #include <GeomAlgoAPI_CompoundBuilder.h>
  46 #include <GeomAlgoAPI_NExplode.h>
  47 #include <Selector_Selector.h>
  48
  49 #include <TNaming_NamedShape.hxx>
  50 #include <TNaming_Tool.hxx>
  51 #include <TNaming_Builder.hxx>
  52 #include <TNaming_SameShapeIterator.hxx>
  53 #include <TNaming_NewShapeIterator.hxx>
  54 #include <TNaming_Iterator.hxx>
  55 #include <TDataStd_Integer.hxx>
  56 #include <TDataStd_UAttribute.hxx>
  57 #include <TDataStd_Name.hxx>
  58 #include <TopTools_ListOfShape.hxx>
  59 #include <TopTools_DataMapOfShapeShape.hxx>
  60 #include <TopTools_MapOfShape.hxx>
  61 #include <TopExp_Explorer.hxx>
  62 #include <BRep_Tool.hxx>
  63 #include <TopoDS.hxx>
  64 #include <TopoDS_Edge.hxx>
  65 #include <TopExp.hxx>
  66 #include <TDF_ChildIterator.hxx>
  67 #include <TDF_ChildIDIterator.hxx>
  68 #include <TopoDS_Iterator.hxx>
  69 #include <TDF_ChildIDIterator.hxx>
  70 #include <Geom_Circle.hxx>
  71 #include <Geom_Ellipse.hxx>
  72 #include <Geom_TrimmedCurve.hxx>
  73 #include <BRep_Builder.hxx>
  74
  75 //#define DEB_NAMING 1
  76 #ifdef DEB_NAMING
  77 #include <BRepTools.hxx>
  78 #endif
  79 /// added to the index in the packed map to signalize that the vertex of edge is selected
  80 /// (multiplied by the index of the edge)
  81 static const int kSTART_VERTEX_DELTA = 1000000;
  82 // identifier that there is simple reference: selection equals to context
  83 Standard_GUID kSIMPLE_REF_ID("635eacb2-a1d6-4dec-8348-471fae17cb29");
  84 // reference to Part sub-object
  85 Standard_GUID kPART_REF_ID("635eacb2-a1d6-4dec-8348-471fae17cb27");
  86 // selection is invalid after recomputation
  87 Standard_GUID kINVALID_SELECTION("bce47fd7-80fa-4462-9d63-2f58acddd49d");
  88
  89 // identifier of the selection of the center of circle on edge
  90 Standard_GUID kCIRCLE_CENTER("d0d0e0f1-217a-4b95-8fbb-0c4132f23718");
  91 // identifier of the selection of the first focus point of ellipse on edge
  92 Standard_GUID kELLIPSE_CENTER1("f70df04c-3168-4dc9-87a4-f1f840c1275d");
  93 // identifier of the selection of the second focus point of ellipse on edge
  94 Standard_GUID kELLIPSE_CENTER2("1395ae73-8e02-4cf8-b204-06ff35873a32");
  95
  96 // prefix for the whole feature context identification
  97 const static std::string kWHOLE_FEATURE = "all-in-";
  98
  99 // on this label is stored:
 100 // TNaming_NamedShape - selected shape
 101 // TNaming_Naming - topological selection information (for the body)
 102 // TDataStd_IntPackedMap - indexes of edges in composite element (for construction)
 103 // TDataStd_Integer - type of the selected shape (for construction)
 104 // TDF_Reference - from ReferenceAttribute, the context
 105 bool Model_AttributeSelection::setValue(const ObjectPtr& theContext,
 106   const std::shared_ptr<GeomAPI_Shape>& theSubShape, const bool theTemporarily)
 107 {
 108   if (theTemporarily) {
 109     // just keep the stored without DF update
 110     myTmpContext = theContext;
 111     myTmpSubShape = theSubShape;
 112     owner()->data()->sendAttributeUpdated(this);
 113     return true;
 114   } else {
 115     myTmpContext.reset();
 116     myTmpSubShape.reset();
 117     myTmpCenterType = NOT_CENTER;
 118   }
 119
 120   CenterType aType;
 121   const std::shared_ptr<GeomAPI_Shape>& anOldShape = internalValue(aType);
 122   bool isOldContext = theContext == myRef.value();
 123   bool isOldShape = isOldContext &&
 124     (theSubShape == anOldShape || (theSubShape && anOldShape && theSubShape->isEqual(anOldShape)));
 125   if (isOldShape) return false; // shape is the same, so context is also unchanged
 126   bool aToUnblock = false;
 127   // update the referenced object if needed
 128   if (!isOldContext) {
 129     aToUnblock = !owner()->data()->blockSendAttributeUpdated(true);
 130     myRef.setValue(theContext);
 131   }
 132
 133   // do not use naming if selected shape is result shape itself, but not sub-shape
 134   TDF_Label aSelLab = selectionLabel();
 135   aSelLab.ForgetAttribute(kSIMPLE_REF_ID);
 136   aSelLab.ForgetAttribute(kINVALID_SELECTION);
 137   aSelLab.ForgetAttribute(kCIRCLE_CENTER);
 138   aSelLab.ForgetAttribute(kELLIPSE_CENTER1);
 139   aSelLab.ForgetAttribute(kELLIPSE_CENTER2);
 140
 141   bool isDegeneratedEdge = false;
 142   // do not use the degenerated edge as a shape, a null context and shape is used in the case
 143   if (theSubShape.get() && !theSubShape->isNull() && theSubShape->isEdge()) {
 144     const TopoDS_Shape& aSubShape = theSubShape->impl<TopoDS_Shape>();
 145     if (aSubShape.ShapeType() == TopAbs_EDGE)
 146       isDegeneratedEdge = BRep_Tool::Degenerated(TopoDS::Edge(aSubShape)) == Standard_True;
 147   }
 148   if (!theContext.get() || isDegeneratedEdge) {
 149     // to keep the reference attribute label
 150     TDF_Label aRefLab = myRef.myRef->Label();
 151     aSelLab.ForgetAllAttributes(true);
 152     myRef.myRef = TDF_Reference::Set(aSelLab.Father(), aSelLab.Father());
 153     if (aToUnblock)
 154       owner()->data()->blockSendAttributeUpdated(false);
 155     return false;
 156   }
 157   bool isSelectBody = theContext->groupName() == ModelAPI_ResultBody::group();
 158   if (!isSelectBody) {
 159     ResultConstructionPtr aContextConstruction =
 160       std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(theContext);
 161     isSelectBody = aContextConstruction.get() && !aContextConstruction->isInfinite();
 162   }
 163   if (isSelectBody) {
 164     ResultPtr aContextResult = std::dynamic_pointer_cast<ModelAPI_Result>(theContext);
 165     GeomShapePtr aContextShape = aContextResult->shape();
 166     // do not select the whole shape for body:it is already must be in the data framework
 167     // equal and null selected objects mean the same: object is equal to context,
 168     if (aContextShape.get() && (aContextShape->isEqual(theSubShape) || !theSubShape.get())) {
 169       aSelLab.ForgetAllAttributes(true);
 170       TDataStd_UAttribute::Set(aSelLab, kSIMPLE_REF_ID);
 171     } else {
 172       selectBody(aContextResult, theSubShape);
 173     }
 174   } else if (theContext->groupName() == ModelAPI_ResultConstruction::group()) {
 175     aSelLab.ForgetAllAttributes(true); // to remove old selection data
 176     std::shared_ptr<Model_ResultConstruction> aConstruction =
 177       std::dynamic_pointer_cast<Model_ResultConstruction>(theContext);
 178     std::shared_ptr<GeomAPI_Shape> aSubShape;
 179     if (theSubShape.get() && !aConstruction->shape()->isEqual(theSubShape))
 180       aSubShape = theSubShape; // the whole context
 181   } else if (theContext->groupName() == ModelAPI_ResultPart::group()) {
 182     aSelLab.ForgetAllAttributes(true);
 183     TDataStd_UAttribute::Set(aSelLab, kPART_REF_ID);
 184     selectPart(std::dynamic_pointer_cast<ModelAPI_Result>(theContext), theSubShape);
 185   } else if (theContext->groupName() == ModelAPI_ResultGroup::group()) {
 186     aSelLab.ForgetAllAttributes(true);
 187     TDataStd_UAttribute::Set(aSelLab, kSIMPLE_REF_ID);
 188   } else { // check the feature context: parent-Part of this feature should not be used
 189     FeaturePtr aFeatureContext = std::dynamic_pointer_cast<ModelAPI_Feature>(theContext);
 190     if (aFeatureContext.get()) {
 191       if (owner()->document() != aFeatureContext->document()) {
 192         aSelLab.ForgetAllAttributes(true);
 193         myRef.setValue(ObjectPtr());
 194         if (aToUnblock)
 195           owner()->data()->blockSendAttributeUpdated(false);
 196         return false;
 197       }
 198     }
 199   }
 200
 201   owner()->data()->sendAttributeUpdated(this);
 202
 203   if (aToUnblock)
 204     owner()->data()->blockSendAttributeUpdated(false);
 205
 206   return true;
 207 }
 208
 209 void Model_AttributeSelection::setValueCenter(
 210     const ObjectPtr& theContext, const std::shared_ptr<GeomAPI_Edge>& theEdge,
 211     const CenterType theCenterType, const bool theTemporarily)
 212 {
 213   bool anUpdated = setValue(theContext, theEdge, theTemporarily);
 214   if (theTemporarily) {
 215     myTmpCenterType = theCenterType;
 216   } else { // store in the data structure
 217     TDF_Label aSelLab = selectionLabel();
 218     switch(theCenterType) {
 219     case CIRCLE_CENTER:
 220       if (!anUpdated)
 221         anUpdated = !aSelLab.IsAttribute(kCIRCLE_CENTER);
 222       TDataStd_UAttribute::Set(aSelLab, kCIRCLE_CENTER);
 223       break;
 224     case ELLIPSE_FIRST_FOCUS:
 225       if (!anUpdated)
 226         anUpdated = !aSelLab.IsAttribute(kELLIPSE_CENTER1);
 227       TDataStd_UAttribute::Set(aSelLab, kELLIPSE_CENTER1);
 228       break;
 229     case ELLIPSE_SECOND_FOCUS:
 230       if (!anUpdated)
 231         anUpdated = !aSelLab.IsAttribute(kELLIPSE_CENTER2);
 232       TDataStd_UAttribute::Set(aSelLab, kELLIPSE_CENTER2);
 233       break;
 234     }
 235     if (anUpdated)
 236       owner()->data()->sendAttributeUpdated(this);
 237   }
 238 }
 239
 240 void Model_AttributeSelection::selectValue(
 241     const std::shared_ptr<ModelAPI_AttributeSelection>& theSource)
 242 {
 243   CenterType aType;
 244   std::shared_ptr<GeomAPI_Shape> aValue =
 245     std::dynamic_pointer_cast<Model_AttributeSelection>(theSource)->internalValue(aType);
 246   if (!aValue.get() || aType == NOT_CENTER) {
 247     setValue(theSource->context(), aValue);
 248   } else {
 249     std::shared_ptr<GeomAPI_Edge> anEdge(new GeomAPI_Edge);
 250     anEdge->setImpl(new TopoDS_Shape(aValue->impl<TopoDS_Shape>()));
 251     setValueCenter(theSource->context(), anEdge, aType);
 252   }
 253 }
 254
 255 void Model_AttributeSelection::removeTemporaryValues()
 256 {
 257   if (myTmpContext.get() || myTmpSubShape.get()) {
 258     myTmpContext.reset();
 259     myTmpSubShape.reset();
 260   }
 261 }
 262
 263 // returns the center of the edge: circular or elliptic
 264 GeomShapePtr centerByEdge(GeomShapePtr theEdge, ModelAPI_AttributeSelection::CenterType theType)
 265 {
 266   if (theType != ModelAPI_AttributeSelection::NOT_CENTER && theEdge.get() != NULL) {
 267     TopoDS_Shape aShape = theEdge->impl<TopoDS_Shape>();
 268     if (!aShape.IsNull() && aShape.ShapeType() == TopAbs_EDGE) {
 269       TopoDS_Edge anEdge = TopoDS::Edge(aShape);
 270       double aFirst, aLast;
 271       Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast);
 272       if (!aCurve.IsNull()) {
 273         TopoDS_Vertex aVertex;
 274         BRep_Builder aBuilder;
 275         if (theType == ModelAPI_AttributeSelection::CIRCLE_CENTER) {
 276           while(!aCurve.IsNull() && aCurve->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve))
 277             aCurve = Handle(Geom_TrimmedCurve)::DownCast(aCurve)->BasisCurve();
 278           Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(aCurve);
 279           if (!aCirc.IsNull()) {
 280             aBuilder.MakeVertex(aVertex, aCirc->Location(), Precision::Confusion());
 281           }
 282         } else { // ellipse
 283           while(!aCurve.IsNull() && aCurve->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve))
 284             aCurve = Handle(Geom_TrimmedCurve)::DownCast(aCurve)->BasisCurve();
 285           Handle(Geom_Ellipse) anEll = Handle(Geom_Ellipse)::DownCast(aCurve);
 286           if (!anEll.IsNull()) {
 287             aBuilder.MakeVertex(aVertex,
 288               theType == ModelAPI_AttributeSelection::ELLIPSE_FIRST_FOCUS ?
 289               anEll->Focus1() : anEll->Focus2(), Precision::Confusion());
 290           }
 291         }
 292         if (!aVertex.IsNull()) {
 293           std::shared_ptr<GeomAPI_Vertex> aResult(new GeomAPI_Vertex);
 294           aResult->setImpl(new TopoDS_Vertex(aVertex));
 295           return aResult;
 296         }
 297       }
 298     }
 299   }
 300   return theEdge; // no vertex, so, return the initial edge
 301 }
 302
 303 std::shared_ptr<GeomAPI_Shape> Model_AttributeSelection::value()
 304 {
 305   if (!ModelAPI_AttributeSelection::isInitialized() && !myTmpContext.get() && !myTmpSubShape.get())
 306     return std::shared_ptr<GeomAPI_Shape>();
 307   CenterType aType = NOT_CENTER;
 308   std::shared_ptr<GeomAPI_Shape> aResult = internalValue(aType);
 309   return centerByEdge(aResult, aType);
 310 }
 311
 312 std::shared_ptr<GeomAPI_Shape> Model_AttributeSelection::internalValue(CenterType& theType)
 313 {
 314   theType = NOT_CENTER;
 315   GeomShapePtr aResult;
 316   if (myTmpContext.get() || myTmpSubShape.get()) {
 317     theType = myTmpCenterType;
 318     ResultConstructionPtr aResulConstruction =
 319       std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(myTmpContext);
 320     if(aResulConstruction.get()) {
 321       // it is just reference to construction.
 322       return myTmpSubShape;
 323     }
 324     FeaturePtr aFeature =
 325       std::dynamic_pointer_cast<ModelAPI_Feature>(myTmpContext);
 326     if (aFeature.get()) {
 327       // it is just reference to construction.
 328       return myTmpSubShape;
 329     }
 330     return myTmpSubShape.get() ? myTmpSubShape :
 331       std::dynamic_pointer_cast<ModelAPI_Result>(myTmpContext)->shape();
 332   }
 333
 334   TDF_Label aSelLab = selectionLabel();
 335   if (aSelLab.IsAttribute(kINVALID_SELECTION))
 336     return aResult;
 337
 338   if (aSelLab.IsAttribute(kCIRCLE_CENTER))
 339     theType = CIRCLE_CENTER;
 340   else if (aSelLab.IsAttribute(kELLIPSE_CENTER1))
 341     theType = ELLIPSE_FIRST_FOCUS;
 342   else if (aSelLab.IsAttribute(kELLIPSE_CENTER2))
 343     theType = ELLIPSE_SECOND_FOCUS;
 344
 345
 346   if (myRef.isInitialized()) {
 347     if (aSelLab.IsAttribute(kSIMPLE_REF_ID)) { // it is just reference to shape, not sub-shape
 348       ResultPtr aContext = context();
 349       if (!aContext.get() || aContext->groupName() == ModelAPI_ResultConstruction::group())
 350         return aResult; // empty result, for whole construction selection also
 351       return aContext->shape();
 352     }
 353     if (aSelLab.IsAttribute(kPART_REF_ID)) {
 354       ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(context());
 355       if (!aPart.get() || !aPart->isActivated())
 356         return std::shared_ptr<GeomAPI_Shape>(); // postponed naming needed
 357       Handle(TDataStd_Integer) anIndex;
 358       if (aSelLab.FindAttribute(TDataStd_Integer::GetID(), anIndex)) {
 359         if (anIndex->Get()) { // special selection attribute was created, use it
 360           return aPart->selectionValue(anIndex->Get());
 361         } else { // face with name is already in the data model, so try to take it by name
 362           Handle(TDataStd_Name) aName;
 363           if (aSelLab.FindAttribute(TDataStd_Name::GetID(), aName)) {
 364             std::string aSubShapeName(TCollection_AsciiString(aName->Get()).ToCString());
 365             std::size_t aPartEnd = aSubShapeName.find('/');
 366             if (aPartEnd != std::string::npos && aPartEnd != aSubShapeName.rfind('/')) {
 367               std::string aNameInPart = aSubShapeName.substr(aPartEnd + 1);
 368               int anIndex;
 369               std::string aType; // to reuse already existing selection the type is not needed
 370               return aPart->shapeInPart(aNameInPart, aType, anIndex);
 371             }
 372           }
 373         }
 374       }
 375     }
 376
 377     std::shared_ptr<Model_ResultConstruction> aConstr =
 378       std::dynamic_pointer_cast<Model_ResultConstruction>(context());
 379     if (aConstr) {
 380       if (aConstr->isInfinite())
 381         return aResult; // empty result
 382     }
 383     if (!aConstr.get()) { // for construction context, return empty result as usual even
 384       // the whole feature is selected
 385       FeaturePtr aFeature = contextFeature();
 386       if (aFeature.get()) {
 387         std::list<GeomShapePtr> allShapes;
 388         std::list<ResultPtr>::const_iterator aRes = aFeature->results().cbegin();
 389         for (; aRes != aFeature->results().cend(); aRes++) {
 390           if (aRes->get() && !(*aRes)->isDisabled()) {
 391             GeomShapePtr aShape = (*aRes)->shape();
 392             if (aShape.get() && !aShape->isNull()) {
 393               allShapes.push_back(aShape);
 394             }
 395           }
 396         }
 397         return GeomAlgoAPI_CompoundBuilder::compound(allShapes);
 398       }
 399     } else {
 400       if (contextFeature().get()) {
 401         return aResult; // for the whole sketch feature selected return null => all faces
 402       }
 403     }
 404
 405     Handle(TNaming_NamedShape) aSelection;
 406     if (aSelLab.FindAttribute(TNaming_NamedShape::GetID(), aSelection)) {
 407       TopoDS_Shape aSelShape = aSelection->Get();
 408       aResult = std::shared_ptr<GeomAPI_Shape>(new GeomAPI_Shape);
 409       aResult->setImpl(new TopoDS_Shape(aSelShape));
 410     } else if (aConstr) { // simple construction element: just shape of this construction element
 411       aResult = aConstr->shape();
 412     }
 413   }
 414   return aResult;
 415 }
 416
 417 bool Model_AttributeSelection::isInvalid()
 418 {
 419   return selectionLabel().IsAttribute(kINVALID_SELECTION) == Standard_True;
 420 }
 421
 422 bool Model_AttributeSelection::isInitialized()
 423 {
 424   if (ModelAPI_AttributeSelection::isInitialized()) { // additional checks if it is initialized
 425     std::shared_ptr<GeomAPI_Shape> aResult;
 426     if (myRef.isInitialized()) {
 427       TDF_Label aSelLab = selectionLabel();
 428       // it is just reference to shape, not sub-shape
 429       if (aSelLab.IsAttribute(kSIMPLE_REF_ID) || aSelLab.IsAttribute(kPART_REF_ID)) {
 430         ResultPtr aContext = context();
 431         return aContext.get() != NULL;
 432       }
 433       Handle(TNaming_NamedShape) aSelection;
 434       if (selectionLabel().FindAttribute(TNaming_NamedShape::GetID(), aSelection)) {
 435         return !aSelection->Get().IsNull();
 436       } else { // for simple construction element: just shape of this construction element
 437         if (myRef.value().get())
 438           return true;
 439         // check that this is on open of document, so, results are not initialized yet
 440         TDF_Label aRefLab = myRef.myRef->Get();
 441         if (aRefLab.IsNull() || !owner().get())
 442           return false;
 443         std::shared_ptr<Model_Document> aMyDoc =
 444           std::dynamic_pointer_cast<Model_Document>(owner()->document());
 445         if (!aMyDoc.get())
 446           return false;
 447         // check at least the feature exists
 448         return aMyDoc->featureByLab(aRefLab).get() != NULL;
 449       }
 450     }
 451   }
 452   return false;
 453 }
 454
 455 Model_AttributeSelection::Model_AttributeSelection(TDF_Label& theLabel)
 456 : myRef(theLabel),
 457   myIsGeometricalSelection(false)
 458 {
 459   myIsInitialized = myRef.isInitialized();
 460   myParent = NULL;
 461 }
 462
 463 void Model_AttributeSelection::setID(const std::string theID)
 464 {
 465   myRef.setID(theID);
 466   ModelAPI_AttributeSelection::setID(theID);
 467   FeaturePtr aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
 468   if (myParent) {
 469     // to be able to select as geometrical selection and then - split
 470     myIsGeometricalSelection = true; // myParent->isGeometricalSelection();
 471   } else {
 472     myIsGeometricalSelection =
 473       ModelAPI_Session::get()->validators()->isGeometricalSelection(aFeature->getKind(), id());
 474   }
 475 }
 476
 477 ResultPtr Model_AttributeSelection::context()
 478 {
 479   if (!ModelAPI_AttributeSelection::isInitialized() && !myTmpContext.get() && !myTmpSubShape.get())
 480     return ResultPtr();
 481
 482   if (myTmpContext.get() || myTmpSubShape.get()) {
 483     return std::dynamic_pointer_cast<ModelAPI_Result>(myTmpContext);
 484   }
 485
 486   ResultPtr aResult = std::dynamic_pointer_cast<ModelAPI_Result>(myRef.value());
 487   // for parts there could be same-data result, so take the last enabled
 488   if (aResult.get()) {
 489     if(aResult->groupName() == ModelAPI_ResultPart::group()) {
 490       int aSize = aResult->document()->size(ModelAPI_ResultPart::group());
 491       for(int a = aSize - 1; a >= 0; a--) {
 492         ObjectPtr aPart = aResult->document()->object(ModelAPI_ResultPart::group(), a);
 493         if(aPart.get() && aPart->data() == aResult->data()) {
 494           ResultPtr aPartResult = std::dynamic_pointer_cast<ModelAPI_Result>(aPart);
 495           FeaturePtr anOwnerFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
 496           // check that this result is not this-feature result (it is forbidden to select itself)
 497           if(anOwnerFeature.get() && anOwnerFeature->firstResult() != aPartResult) {
 498             return aPartResult;
 499           }
 500         }
 501       }
 502     }
 503   } else { // if feature - construction is selected, it has only one result, return this result
 504     FeaturePtr aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(myRef.value());
 505     if (aFeature.get() && aFeature->results().size() == 1 &&
 506         aFeature->firstResult()->groupName() == ModelAPI_ResultConstruction::group())
 507       return aFeature->firstResult();
 508   }
 509   return aResult;
 510 }
 511
 512 FeaturePtr Model_AttributeSelection::contextFeature() {
 513   if (myTmpContext.get()) {
 514     return std::dynamic_pointer_cast<ModelAPI_Feature>(myTmpContext);
 515   }
 516   return std::dynamic_pointer_cast<ModelAPI_Feature>(myRef.value());
 517 }
 518 ObjectPtr Model_AttributeSelection::contextObject() {
 519   FeaturePtr aRes = contextFeature();
 520   if (aRes.get())
 521     return aRes;
 522   return context();
 523 }
 524
 525
 526 void Model_AttributeSelection::setObject(const std::shared_ptr<ModelAPI_Object>& theObject)
 527 {
 528   ModelAPI_AttributeSelection::setObject(theObject);
 529   myRef.setObject(theObject);
 530 }
 531
 532 /// Sets the invalid flag if flag is false, or removes it if "true"
 533 /// Returns theFlag
 534 static bool setInvalidIfFalse(TDF_Label& theLab, const bool theFlag) {
 535   if (theFlag) {
 536     theLab.ForgetAttribute(kINVALID_SELECTION);
 537   } else {
 538     TDataStd_UAttribute::Set(theLab, kINVALID_SELECTION);
 539   }
 540   return theFlag;
 541 }
 542
 543 void Model_AttributeSelection::split(
 544   ResultPtr theContext, TopoDS_Shape theNewShape, TopAbs_ShapeEnum theType)
 545 {
 546   TopTools_ListOfShape aSubs;
 547   for(TopoDS_Iterator anExplorer(theNewShape); anExplorer.More(); anExplorer.Next()) {
 548     if (!anExplorer.Value().IsNull() &&
 549       anExplorer.Value().ShapeType() == theType) {
 550         aSubs.Append(anExplorer.Value());
 551     } else { // invalid case; bad result shape, so, impossible to split easily
 552       aSubs.Clear();
 553       break;
 554     }
 555   }
 556   if (aSubs.Extent() > 1) { // ok to split
 557     TopTools_ListIteratorOfListOfShape aSub(aSubs);
 558     GeomShapePtr aSubSh(new GeomAPI_Shape);
 559     aSubSh->setImpl(new TopoDS_Shape(aSub.Value()));
 560     setValue(theContext, aSubSh);
 561     for(aSub.Next(); aSub.More(); aSub.Next()) {
 562       GeomShapePtr aSubSh(new GeomAPI_Shape);
 563       aSubSh->setImpl(new TopoDS_Shape(aSub.Value()));
 564       myParent->append(theContext, aSubSh);
 565     }
 566   }
 567 }
 568
 569 bool Model_AttributeSelection::update()
 570 {
 571   FeaturePtr aContextFeature = contextFeature();
 572   if (aContextFeature.get()) {
 573     owner()->data()->sendAttributeUpdated(this);  // send updated if "update" called in any way
 574     return true;
 575   }
 576   TDF_Label aSelLab = selectionLabel();
 577   ResultPtr aContext = context();
 578   if (!aContext.get())
 579     return setInvalidIfFalse(aSelLab, false);
 580   if (aSelLab.IsAttribute(kSIMPLE_REF_ID)) { // it is just reference to shape, not sub-shape
 581     return setInvalidIfFalse(aSelLab, aContext->shape() && !aContext->shape()->isNull());
 582   }
 583
 584   if (aSelLab.IsAttribute(kPART_REF_ID)) { // it is reference to the part object
 585     std::shared_ptr<GeomAPI_Shape> aNoSelection;
 586     bool aResult = selectPart(aContext, aNoSelection, true);
 587     aResult = setInvalidIfFalse(aSelLab, aResult);
 588     if (aResult) {
 589       owner()->data()->sendAttributeUpdated(this);
 590     }
 591     return aResult;
 592   }
 593
 594   if (aContext->groupName() == ModelAPI_ResultBody::group()) {
 595     // body: just a named shape, use topological selection mechanism
 596     bool aResult = false;
 597     TopoDS_Shape anOldShape;
 598     Handle(TNaming_NamedShape) aNS;
 599     if (aSelLab.FindAttribute(TNaming_NamedShape::GetID(), aNS))
 600       anOldShape = aNS->Get();
 601
 602     TopoDS_Shape aContextShape = aContext->shape()->impl<TopoDS_Shape>();
 603     Selector_Selector aSelector(aSelLab, baseDocumentLab());
 604     aResult = aSelector.restore(aContextShape);
 605     setInvalidIfFalse(aSelLab, aResult);
 606
 607     TopoDS_Shape aNewShape;
 608     if (aSelLab.FindAttribute(TNaming_NamedShape::GetID(), aNS))
 609       aNewShape = aNS->Get();
 610
 611     if (anOldShape.IsNull() || aNewShape.IsNull() || !anOldShape.IsEqual(aNewShape)) {
 612       // shape type should not be changed: if shape becomes compound of such shapes, then split
 613       if (myParent && !anOldShape.IsNull() && !aNewShape.IsNull() &&
 614           anOldShape.ShapeType() != aNewShape.ShapeType() &&
 615           aNewShape.ShapeType() == TopAbs_COMPOUND) {
 616         split(aContext, aNewShape, anOldShape.ShapeType());
 617       }
 618       owner()->data()->sendAttributeUpdated(this);  // send updated if shape is changed
 619     }
 620     return aResult;
 621   }
 622
 623   if (aContext->groupName() == ModelAPI_ResultConstruction::group()) {
 624     bool aResult = true;
 625     std::shared_ptr<Model_ResultConstruction> aConstructionContext =
 626       std::dynamic_pointer_cast<Model_ResultConstruction>(aContext);
 627     if (!aConstructionContext->isInfinite()) {
 628       TopoDS_Shape aContextShape = aContext->shape()->impl<TopoDS_Shape>();
 629       Selector_Selector aSelector(aSelLab, baseDocumentLab());
 630       TopoDS_Shape anOldShape = aSelector.value();
 631       aResult = aSelector.restore(aContextShape);
 632       setInvalidIfFalse(aSelLab, aResult);
 633       if (aResult && !anOldShape.IsEqual(aSelector.value()))
 634         owner()->data()->sendAttributeUpdated(this);  // send updated if shape is changed
 635     } else {
 636       owner()->data()->sendAttributeUpdated(this);  // send updated if "update" called in any way
 637     }
 638     return aResult;
 639   }
 640   return setInvalidIfFalse(aSelLab, false); // unknown case
 641 }
 642
 643 void Model_AttributeSelection::selectBody(
 644   const ResultPtr& theContext, const std::shared_ptr<GeomAPI_Shape>& theSubShape)
 645 {
 646   // perform the selection
 647   TopoDS_Shape aContext;
 648
 649   ResultPtr aBody = std::dynamic_pointer_cast<ModelAPI_Result>(theContext);
 650   if (aBody) {
 651     aContext = aBody->shape()->impl<TopoDS_Shape>();
 652   } else {
 653     ResultPtr aResult =
 654       std::dynamic_pointer_cast<ModelAPI_Result>(myRef.value());
 655     if (aResult) {
 656       aContext = aResult->shape()->impl<TopoDS_Shape>();
 657     } else {
 658       Events_InfoMessage("Model_AttributeSelection", "A result with shape is expected").send();
 659       return;
 660     }
 661   }
 662
 663   if (!aContext.IsNull()) {
 664     TDF_Label aSelLab = selectionLabel();
 665     TopoDS_Shape aNewSub = theSubShape->impl<TopoDS_Shape>();
 666
 667     bool aSelectorOk = true;
 668     Selector_Selector aSelector(aSelLab, baseDocumentLab());
 669     try {
 670       aSelectorOk = aSelector.select(aContext, aNewSub, myIsGeometricalSelection);
 671       if (aSelectorOk) {
 672         aSelectorOk = aSelector.store(aContext);
 673       }
 674     } catch(...) {
 675       aSelectorOk = false;
 676     }
 677     if (aSelectorOk) {
 678       TopoDS_Shape aShape = aSelector.value();
 679       aSelectorOk = !aShape.IsNull() && aShape.ShapeType() == aNewSub.ShapeType();
 680     }
 681     setInvalidIfFalse(aSelLab, aSelectorOk);
 682   }
 683 }
 684
 685 bool Model_AttributeSelection::selectPart(
 686   const ResultPtr& theContext, const std::shared_ptr<GeomAPI_Shape>& theSubShape,
 687   const bool theUpdate)
 688 {
 689   ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(theContext);
 690   if (!aPart.get() || !aPart->isActivated())
 691     return true; // postponed naming
 692   if (theUpdate) {
 693     Handle(TDataStd_Integer) anIndex;
 694     if (selectionLabel().FindAttribute(TDataStd_Integer::GetID(), anIndex)) {
 695       // by internal selection
 696       if (anIndex->Get() > 0) {
 697         // update the selection by index
 698         return aPart->updateInPart(anIndex->Get());
 699       } else {
 700         return true; // nothing to do, referencing just by name
 701       }
 702     }
 703     return true; // nothing to do, referencing just by name
 704   }
 705   // store the shape (in case part is not loaded it should be useful
 706   TopoDS_Shape aShape;
 707   std::string aName = theContext->data()->name();
 708   if (!theSubShape.get() || theSubShape->isNull()) {// the whole part shape is selected
 709     aShape = theContext->shape()->impl<TopoDS_Shape>();
 710   } else {
 711     aShape = theSubShape->impl<TopoDS_Shape>();
 712     int anIndex;
 713     aName += "/" + aPart->nameInPart(theSubShape, anIndex);
 714     TDataStd_Integer::Set(selectionLabel(), anIndex);
 715   }
 716   TNaming_Builder aBuilder(selectionLabel());
 717   aBuilder.Select(aShape, aShape);
 718   // identify by name in the part
 719   TDataStd_Name::Set(selectionLabel(), aName.c_str());
 720   return !aName.empty();
 721 }
 722
 723 TDF_Label Model_AttributeSelection::selectionLabel()
 724 {
 725   return myRef.myRef->Label().FindChild(1);
 726 }
 727
 728 /// prefixes of the shape names with centers defined
 729 static std::map<ModelAPI_AttributeSelection::CenterType, std::string> kCENTERS_PREFIX;
 730
 731 /// returns the map that contains all possible prefixes of the center-names
 732 static std::map<ModelAPI_AttributeSelection::CenterType, std::string>& centersMap()
 733 {
 734   if (kCENTERS_PREFIX.empty()) { // fill map by initial values
 735     kCENTERS_PREFIX[ModelAPI_AttributeSelection::CIRCLE_CENTER] = "__cc";
 736     kCENTERS_PREFIX[ModelAPI_AttributeSelection::ELLIPSE_FIRST_FOCUS] = "__eff";
 737     kCENTERS_PREFIX[ModelAPI_AttributeSelection::ELLIPSE_SECOND_FOCUS] = "__esf";
 738   }
 739   return kCENTERS_PREFIX;
 740 }
 741
 742 std::string Model_AttributeSelection::namingName(const std::string& theDefaultName)
 743 {
 744   std::string aName("");
 745   if(!this->isInitialized())
 746     return !theDefaultName.empty() ? theDefaultName : aName;
 747
 748   // not argument has not parametric name (filters)
 749   if (!this->isArgument() || (myParent && !myParent->isArgument())) {
 750     GeomShapePtr aShape = value();
 751     if (!aShape.get() && context().get())
 752       aShape = context()->shape();
 753     std::string aName;
 754     if (aShape.get()) {
 755       aName = aShape->shapeTypeStr();
 756       if (myParent) {
 757         aName += std::string("_") +
 758           TCollection_AsciiString(selectionLabel().Father().Tag()).ToCString();
 759       }
 760     }
 761     return aName;
 762   }
 763
 764   CenterType aCenterType = NOT_CENTER;
 765   std::shared_ptr<GeomAPI_Shape> aSubSh = internalValue(aCenterType);
 766   ResultPtr aCont = context();
 767
 768   if (!aCont.get() ||
 769       (aCont->groupName() == ModelAPI_ResultConstruction::group() && contextFeature().get())) {
 770     // selection of a full feature
 771     FeaturePtr aFeatureCont = contextFeature();
 772     if (aFeatureCont.get()) {
 773       return kWHOLE_FEATURE + aFeatureCont->name();
 774     }
 775     // in case of selection of removed result
 776     return "";
 777   }
 778
 779   TDF_Label aSelLab = selectionLabel();
 780   if (aSelLab.IsAttribute(kSIMPLE_REF_ID)) { // whole context, no value
 781     return contextName(aCont);
 782   }
 783
 784   // if it is in result of another part
 785   if (aCont->groupName() == ModelAPI_ResultPart::group()) {
 786     ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(aCont);
 787     int anIndex;
 788     std::string aResult = aSubSh.get() ?
 789       aPart->data()->name() + "/" + aPart->nameInPart(aSubSh, anIndex) : aPart->data()->name();
 790     if (aCenterType != NOT_CENTER)
 791       aResult += centersMap()[aCenterType];
 792     return aResult;
 793   }
 794
 795   // whole infinitive construction
 796   if (aCont->groupName() == ModelAPI_ResultConstruction::group()) {
 797     ResultConstructionPtr aConstr = std::dynamic_pointer_cast<Model_ResultConstruction>(aCont);
 798     if (aConstr->isInfinite()) {
 799       return contextName(aCont);
 800     }
 801   }
 802
 803   Selector_Selector aSelector(aSelLab, baseDocumentLab());
 804   std::string aResult;
 805   if (aCont->shape().get() && aSelector.restore(aCont->shape()->impl<TopoDS_Shape>()))
 806     aResult = aSelector.name(this);
 807   if (aCenterType != NOT_CENTER) {
 808     aResult += centersMap()[aCenterType];
 809   }
 810   return aResult;
 811 }
 812
 813 // returns the center type and modifies the shape name if this name is center-name
 814 static ModelAPI_AttributeSelection::CenterType centerTypeByName(std::string& theShapeName)
 815 {
 816   std::map<ModelAPI_AttributeSelection::CenterType, std::string>::iterator aPrefixIter =
 817     centersMap().begin();
 818   for(; aPrefixIter != centersMap().end(); aPrefixIter++) {
 819     std::size_t aFound = theShapeName.find(aPrefixIter->second);
 820     if (aFound != std::string::npos &&
 821         aFound == theShapeName.size() - aPrefixIter->second.size()) {
 822       theShapeName = theShapeName.substr(0, aFound);
 823       return aPrefixIter->first;
 824     }
 825   }
 826   return ModelAPI_AttributeSelection::NOT_CENTER;
 827 }
 828
 829 // type ::= COMP | COMS | SOLD | SHEL | FACE | WIRE | EDGE | VERT
 830 void Model_AttributeSelection::selectSubShape(
 831   const std::string& theType, const std::string& theSubShapeName)
 832 {
 833   if(theSubShapeName.empty() || theType.empty()) return;
 834
 835   std::string aSubShapeName = theSubShapeName;
 836   CenterType aCenterType = theType[0] == 'v' || theType[0] == 'V' ? // only for vertex-type
 837     centerTypeByName(aSubShapeName) : NOT_CENTER;
 838   std::string aType = aCenterType == NOT_CENTER ? theType : "EDGE"; // search for edge now
 839   static const GeomShapePtr anEmptyShape;
 840
 841   // first iteration is selection by name without center prefix, second - in case of problem,
 842   // try with initial name
 843   for(int aUseCenter = 1; aUseCenter >= 0; aUseCenter--)  {
 844     if (aUseCenter == 0 && aCenterType != NOT_CENTER) {
 845       aSubShapeName = theSubShapeName;
 846       aCenterType = NOT_CENTER;
 847       aType = theType;
 848     } else if (aUseCenter != 1) continue;
 849     TopAbs_ShapeEnum aShapeType =  TopAbs_ShapeEnum(GeomAPI_Shape::shapeTypeByStr(aType));
 850
 851     std::shared_ptr<Model_Document> aDoc =
 852       std::dynamic_pointer_cast<Model_Document>(owner()->document());
 853     // check this is Part-name: 2 delimiters in the name
 854     std::size_t aPartEnd = aSubShapeName.find('/');
 855     if (aPartEnd != std::string::npos) {
 856       std::string aPartName = aSubShapeName.substr(0, aPartEnd);
 857       DocumentPtr aRootDoc = ModelAPI_Session::get()->moduleDocument();
 858       if (aPartName == aRootDoc->kind()) {
 859         aDoc = std::dynamic_pointer_cast<Model_Document>(aRootDoc);
 860         aSubShapeName = aSubShapeName.substr(aPartEnd + 1);
 861       } else {
 862         ObjectPtr aFound =
 863           owner()->document()->objectByName(ModelAPI_ResultPart::group(), aPartName);
 864         if (aFound.get()) { // found such part, so asking it for the name
 865           ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(aFound);
 866           std::string aNameInPart = aSubShapeName.substr(aPartEnd + 1);
 867           if (aNameInPart.empty()) { // whole part
 868             setValue(aPart, anEmptyShape);
 869             return;
 870           }
 871           int anIndex;
 872           std::shared_ptr<GeomAPI_Shape> aSelected =
 873             aPart->shapeInPart(aNameInPart, aType, anIndex);
 874           if (aSelected.get()) {
 875             if (aCenterType != NOT_CENTER) {
 876               if (!aSelected->isEdge())
 877                 continue;
 878               std::shared_ptr<GeomAPI_Edge> aSelectedEdge(new GeomAPI_Edge(aSelected));
 879               setValueCenter(aPart, aSelectedEdge, aCenterType);
 880             } else
 881               setValue(aPart, aSelected);
 882             TDataStd_Integer::Set(selectionLabel(), anIndex);
 883             return;
 884           }
 885         }
 886       }
 887     }
 888
 889     // check this is a whole feature context
 890     if (aSubShapeName.size() > kWHOLE_FEATURE.size() &&
 891       aSubShapeName.substr(0, kWHOLE_FEATURE.size()) == kWHOLE_FEATURE) {
 892       std::string aFeatureName = aSubShapeName.substr(kWHOLE_FEATURE.size());
 893       ObjectPtr anObj = aDoc->objectByName(ModelAPI_Feature::group(), aFeatureName);
 894       if (anObj.get()) {
 895         setValue(anObj, anEmptyShape);
 896         return;
 897       }
 898     }
 899
 900     // the whole result selection check
 901     if (aSubShapeName.find('/') == std::string::npos) {
 902       ObjectPtr aRes = aDoc->objectByName(ModelAPI_ResultConstruction::group(), aSubShapeName);
 903       if (!aRes.get()) {
 904         aRes = aDoc->objectByName(ModelAPI_ResultBody::group(), aSubShapeName);
 905         if (!aRes.get())
 906           aRes = aDoc->objectByName(ModelAPI_ResultGroup::group(), aSubShapeName);
 907       }
 908       if (aRes.get()) {
 909         setValue(aRes, anEmptyShape);
 910         return;
 911       }
 912     }
 913
 914     Selector_Selector aSelector(selectionLabel(), baseDocumentLab());
 915     myRestoreDocument = aDoc;
 916     TDF_Label aContextLabel = aSelector.restoreByName(
 917       aSubShapeName, aShapeType, this, myIsGeometricalSelection);
 918     myRestoreDocument.reset();
 919     if (!aContextLabel.IsNull()) {
 920       ResultPtr aContext = aDoc->resultByLab(aContextLabel); // any label for document access
 921       if (aContext.get() && aContext->shape().get()) {
 922         TopoDS_Shape aContextShape = aContext->shape()->impl<TopoDS_Shape>();
 923         if (aSelector.solve(aContextShape)) {
 924           TopoDS_Shape aSelectorShape = aSelector.value();
 925           GeomShapePtr aShapeToBeSelected(new GeomAPI_Shape);
 926           aShapeToBeSelected->setImpl<TopoDS_Shape>(new TopoDS_Shape(aSelectorShape));
 927           // make the context result the latest existing
 928           aContext = newestContext(aContext, aShapeToBeSelected);
 929           if (myIsGeometricalSelection || aCenterType == NOT_CENTER) {
 930             // store the currently generated name
 931             selectionLabel().ForgetAllAttributes(true);
 932             bool aToUnblock = false;
 933             aToUnblock = !owner()->data()->blockSendAttributeUpdated(true);
 934             myRef.setValue(aContext);
 935             aSelector.store(aContextShape);
 936             owner()->data()->sendAttributeUpdated(this);
 937             if (aToUnblock)
 938               owner()->data()->blockSendAttributeUpdated(false);
 939             return;
 940           } else { // re-select center of circle/arc by context and value
 941             if (!aShapeToBeSelected->isEdge())
 942               continue;
 943             std::shared_ptr<GeomAPI_Edge> aSelectedEdge(new GeomAPI_Edge(aShapeToBeSelected));
 944             setValueCenter(aContext, aSelectedEdge, aCenterType);
 945           }
 946           return;
 947         }
 948       }
 949     }
 950   }
 951   // invalid
 952   TDF_Label aSelLab = selectionLabel();
 953   setInvalidIfFalse(aSelLab, false);
 954   reset();
 955 }
 956
 957 void Model_AttributeSelection::selectSubShape(const std::string& theType,
 958                                               const GeomPointPtr& thePoint)
 959 {
 960   if (theType.empty() || !thePoint)
 961     return;
 962
 963   // list of parent features
 964   FeaturePtr anOwner = ModelAPI_Feature::feature(owner());
 965   std::set<FeaturePtr> aParents = ModelAPI_Tools::getParents(anOwner);
 966
 967   int aSelectionIndex = 0;
 968   GeomAPI_Shape::ShapeType aType = GeomAPI_Shape::shapeTypeByStr(theType);
 969   if (aType == GeomAPI_Shape::SHAPE) {
 970     // possibly, the string consists of the type and the index,
 971     // thus, try to separate them
 972     size_t aUndersporePos = theType.find_first_of('_');
 973     if (aUndersporePos != std::string::npos)
 974       aType = GeomAPI_Shape::shapeTypeByStr(theType.substr(0, aUndersporePos));
 975
 976     if (aType != GeomAPI_Shape::SHAPE) {
 977       for (std::string::const_iterator aChar = theType.begin() + aUndersporePos + 1;
 978            aChar != theType.end(); ++aChar) {
 979         if (std::isdigit(*aChar))
 980           aSelectionIndex = aSelectionIndex * 10 + (*aChar - '0');
 981         else {
 982           aSelectionIndex = 1;
 983           break;
 984         }
 985       }
 986       aSelectionIndex -= 1;
 987     }
 988   }
 989
 990   std::list<ModelGeomAlgo_Shape::SubshapeOfResult> anAppropriate;
 991
 992   // collect features from PartSet and the current part
 993   SessionPtr aSession = ModelAPI_Session::get();
 994   std::list<FeaturePtr> aFeatures = aSession->moduleDocument()->allFeatures();
 995   if (aSession->moduleDocument() != owner()->document()) {
 996     std::list<FeaturePtr> aPartFeatures = owner()->document()->allFeatures();
 997     aFeatures.insert(aFeatures.end(), aPartFeatures.begin(), aPartFeatures.end());
 998   }
 999   // Process results of all features from the last to the first
1000   // to find appropriate sub-shape
1001   for (std::list<FeaturePtr>::const_reverse_iterator anIt = aFeatures.rbegin();
1002        anIt != aFeatures.rend(); ++anIt) {
1003     // selection cannot be linked to the parent features
1004     if (aParents.find(*anIt) != aParents.end())
1005       continue;
1006     // check the feature is a part of composite feature (like sketch elements),
1007     // then do not process it, it will be processed in scope of composite feature
1008     bool isSubOfComposite = false;
1009     const std::set<AttributePtr>& aRefs = (*anIt)->data()->refsToMe();
1010     for (std::set<AttributePtr>::const_iterator aRefIt = aRefs.begin();
1011          aRefIt != aRefs.end() && !isSubOfComposite; ++aRefIt) {
1012       FeaturePtr aFeature = ModelAPI_Feature::feature((*aRefIt)->owner());
1013       CompositeFeaturePtr aCompFeature =
1014           std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aFeature);
1015       isSubOfComposite = aCompFeature && aCompFeature->isSub(*anIt);
1016     }
1017     if (isSubOfComposite)
1018       continue;
1019
1020     // process results of the current feature to find appropriate sub-shape
1021     if (ModelGeomAlgo_Shape::findSubshapeByPoint(*anIt, thePoint, aType, anAppropriate)) {
1022       std::list<ModelGeomAlgo_Shape::SubshapeOfResult>::iterator anApIt = anAppropriate.begin();
1023       for (; aSelectionIndex > 0 && anApIt != anAppropriate.end(); --aSelectionIndex)
1024         ++anApIt; // skip this shape, because one of the previous is selected
1025
1026       if (anApIt != anAppropriate.end()) {
1027         if (anApIt->myCenterType == (int)ModelAPI_AttributeSelection::NOT_CENTER)
1028           setValue(anApIt->myResult, anApIt->mySubshape);
1029         else
1030           setValueCenter(anApIt->myResult, anApIt->mySubshape->edge(),
1031                          (ModelAPI_AttributeSelection::CenterType)anApIt->myCenterType);
1032         return;
1033       }
1034     }
1035   }
1036
1037   TDF_Label aSelLab = selectionLabel();
1038   setInvalidIfFalse(aSelLab, false);
1039   reset();
1040 }
1041
1042 void Model_AttributeSelection::selectSubShape(const std::string& theType,
1043   const std::string& theContextName, const int theIndex)
1044 {
1045   // selection of context by name
1046   selectSubShape(theType, theContextName);
1047   ResultPtr aContext = context();
1048   if (aContext.get()) {
1049     GeomShapePtr aContShape = aContext->shape();
1050     if (aContShape.get()) {
1051       GeomAlgoAPI_NExplode aNExp(aContShape, GeomAPI_Shape::shapeTypeByStr(theType));
1052       GeomShapePtr aValue = aNExp.shape(theIndex);
1053       if (aValue.get())
1054         setValue(aContext, aValue);
1055     }
1056   }
1057 }
1058
1059 void Model_AttributeSelection::setId(int theID)
1060 {
1061   std::shared_ptr<GeomAPI_Shape> aSelection;
1062
1063   ResultPtr aContextRes = context();
1064   // support for compsolids:
1065   while(ModelAPI_Tools::bodyOwner(aContextRes).get()) {
1066     aContextRes = ModelAPI_Tools::bodyOwner(aContextRes);
1067   }
1068   std::shared_ptr<GeomAPI_Shape> aContext = aContextRes->shape();
1069
1070   TopoDS_Shape aMainShape = aContext->impl<TopoDS_Shape>();
1071   // searching for the latest main shape
1072   if (theID > 0 && aContext && !aContext->isNull())
1073   {
1074     std::shared_ptr<Model_Document> aDoc =
1075       std::dynamic_pointer_cast<Model_Document>(aContextRes->document());
1076     if (aDoc.get()) {
1077       Handle(TNaming_NamedShape) aNS = TNaming_Tool::NamedShape(aMainShape, aDoc->generalLabel());
1078       if (!aNS.IsNull()) {
1079         aMainShape = TNaming_Tool::CurrentShape(aNS);
1080       }
1081     }
1082
1083     TopTools_IndexedMapOfShape aSubShapesMap;
1084     TopExp::MapShapes(aMainShape, aSubShapesMap);
1085     const TopoDS_Shape& aSelShape = aSubShapesMap.FindKey(theID);
1086
1087     std::shared_ptr<GeomAPI_Shape> aResult(new GeomAPI_Shape);
1088     aResult->setImpl(new TopoDS_Shape(aSelShape));
1089
1090     aSelection = aResult;
1091   }
1092
1093   setValue(aContextRes, aSelection);
1094 }
1095
1096 std::string Model_AttributeSelection::contextName(const ResultPtr& theContext) const
1097 {
1098   std::string aResult;
1099   if (owner()->document() != theContext->document()) {
1100     if (theContext->document() == ModelAPI_Session::get()->moduleDocument()) {
1101       aResult = theContext->document()->kind() + "/";
1102     } else {
1103       ResultPtr aDocRes = ModelAPI_Tools::findPartResult(
1104         ModelAPI_Session::get()->moduleDocument(), theContext->document());
1105       if (aDocRes.get()) {
1106         aResult = aDocRes->data()->name() + "/";
1107       }
1108     }
1109   }
1110   aResult += theContext->data()->name();
1111   return aResult;
1112 }
1113
1114 void Model_AttributeSelection::computeValues(
1115   ResultPtr theOldContext, ResultPtr theNewContext, TopoDS_Shape theValShape,
1116   TopTools_ListOfShape& theShapes)
1117 {
1118   bool aWasWholeContext = theValShape.IsNull();
1119   if (aWasWholeContext) {
1120     theValShape = theOldContext->shape()->impl<TopoDS_Shape>();
1121   }
1122   TopAbs_ShapeEnum aValType = theValShape.ShapeType();
1123   TopoDS_Shape aNewContShape = theNewContext->shape()->impl<TopoDS_Shape>();
1124   // if a new value is unchanged in the new context, do nothing: value is correct
1125   TopExp_Explorer aSubExp(aNewContShape, aValType);
1126   for(; aSubExp.More(); aSubExp.Next()) {
1127     if (aSubExp.Current().IsSame(theValShape)) {
1128       theShapes.Append(theValShape);
1129       return;
1130     }
1131   }
1132   // if new context becomes compsolid, the resulting sub may be in sub-solids
1133   std::list<ResultPtr> aNewToIterate;
1134   aNewToIterate.push_back(theNewContext);
1135   ResultBodyPtr aComp = std::dynamic_pointer_cast<ModelAPI_ResultBody>(theNewContext);
1136   if (aComp.get()) {
1137     std::list<ResultPtr> allNewContextSubs;
1138     ModelAPI_Tools::allSubs(aComp, allNewContextSubs);
1139     std::list<ResultPtr>::iterator aSub = allNewContextSubs.begin();
1140     for(; aSub != allNewContextSubs.end(); aSub++) {
1141       ResultBodyPtr aBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(*aSub);
1142       if (aBody.get() && aBody->numberOfSubs() == 0) // add only lower level subs
1143         aNewToIterate.push_back(aBody);
1144     }
1145   }
1146
1147   // first iteration: searching for the whole shape appearance (like face of the box)
1148   // second iteration: searching for sub-shapes that contain the sub (like vertex on faces)
1149   int aToFindPart = 0;
1150   TopTools_DataMapOfShapeShape aNewToOld; // map from new containers to old containers (with val)
1151   TopTools_MapOfShape anOlds; // to know how many olds produced new containers
1152   for(; aToFindPart != 2 && theShapes.IsEmpty(); aToFindPart++) {
1153     std::list<ResultPtr>::iterator aNewContIter = aNewToIterate.begin();
1154     for(; aNewContIter != aNewToIterate.end(); aNewContIter++) {
1155       std::shared_ptr<Model_Data> aNewData =
1156         std::dynamic_pointer_cast<Model_Data>((*aNewContIter)->data());
1157       TDF_Label aNewLab = aNewData->shapeLab();
1158       // searching for produced sub-shape fully on some label
1159       TDF_ChildIDIterator aNSIter(aNewLab, TNaming_NamedShape::GetID());
1160       for(; aNSIter.More(); aNSIter.Next()) {
1161         Handle(TNaming_NamedShape) aNS = Handle(TNaming_NamedShape)::DownCast(aNSIter.Value());
1162         for(TNaming_Iterator aPairIter(aNS); aPairIter.More(); aPairIter.Next()) {
1163           if (aToFindPart == 0) { // search shape is fully inside
1164             if (aPairIter.OldShape().IsSame(theValShape)) {
1165               if (aPairIter.NewShape().IsNull()) {// value was removed
1166                 theShapes.Clear();
1167                 return;
1168               }
1169               // don't add edges generated from faces
1170               if (aPairIter.NewShape().ShapeType() <= aValType)
1171                 theShapes.Append(aPairIter.NewShape());
1172             }
1173           } else if (!aPairIter.OldShape().IsNull()) { // search shape that contains this sub
1174             TopExp_Explorer anExp(aPairIter.OldShape(), aValType);
1175             for(; anExp.More(); anExp.Next()) {
1176               if (anExp.Current().IsSame(theValShape)) { // found a new container
1177                 if (aPairIter.NewShape().IsNull()) // skip removed high-level shape
1178                   continue;
1179                 aNewToOld.Bind(aPairIter.NewShape(), aPairIter.OldShape());
1180                 anOlds.Add(aPairIter.OldShape());
1181                 break;
1182               }
1183             }
1184           }
1185         }
1186       }
1187     }
1188   }
1189   if (aToFindPart == 2 && !aNewToOld.IsEmpty()) {
1190     // also iterate the whole old shape to find not-modified shapes that contain this old
1191     TopoDS_Shape anOldContShape = theOldContext->shape()->impl<TopoDS_Shape>();
1192     NCollection_Map<TopAbs_ShapeEnum> aNewTypes; // types of shapes to iterate
1193     TopTools_DataMapOfShapeShape::Iterator aNewTypeIter(aNewToOld);
1194     for(; aNewTypeIter.More(); aNewTypeIter.Next()) {
1195       if (aNewTypeIter.Key().ShapeType() != aValType)
1196         aNewTypes.Add(aNewTypeIter.Key().ShapeType());
1197     }
1198     NCollection_Map<TopAbs_ShapeEnum>::Iterator aTypeIter(aNewTypes);
1199     for(; aTypeIter.More(); aTypeIter.Next()) {
1200       for(TopExp_Explorer anExp(anOldContShape, aTypeIter.Value()); anExp.More(); anExp.Next()) {
1201         TopoDS_Shape anOld = anExp.Current();
1202         if (aNewToOld.IsBound(anOld) || anOlds.Contains(anOld)) // this was modified
1203           continue;
1204         TopExp_Explorer aValExp(anOld, aValType);
1205         for(; aValExp.More(); aValExp.Next()) {
1206           const TopoDS_Shape& anUnchanged = aValExp.Current();
1207           if (anUnchanged.IsSame(theValShape)) {
1208             aNewToOld.Bind(anOld, anOld);
1209             anOlds.Add(anOld);
1210             break;
1211           }
1212         }
1213       }
1214     }
1215
1216     // map of sub-shapes -> number of occurrences of these shapes in containers
1217     NCollection_DataMap<TopoDS_Shape, TopTools_MapOfShape, TopTools_ShapeMapHasher> aSubs;
1218     TopTools_DataMapOfShapeShape::Iterator aContIter(aNewToOld);
1219     for(; aContIter.More(); aContIter.Next()) {
1220       TopExp_Explorer aSubExp(aContIter.Key(), aValType);
1221       for(; aSubExp.More(); aSubExp.Next()) {
1222         if (!aSubs.IsBound(aSubExp.Current())) {
1223           aSubs.Bind(aSubExp.Current(), TopTools_MapOfShape());
1224         }
1225         // store old to know how many olds produced this shape
1226         aSubs.ChangeFind(aSubExp.Current()).Add(aContIter.Value());
1227       }
1228     }
1229     // if sub is appeared same times in containers as the number of old shapes that contain it
1230     int aCountInOld = anOlds.Size();
1231     NCollection_DataMap<TopoDS_Shape, TopTools_MapOfShape, TopTools_ShapeMapHasher>::Iterator
1232       aSubsIter(aSubs);
1233     for(; aSubsIter.More(); aSubsIter.Next()) {
1234       if (aSubsIter.Value().Size() == aCountInOld) {
1235         TopoDS_Shape anOld = aSubsIter.Key();
1236         // check this exists in the new shape
1237         TopExp_Explorer aNew(aNewContShape, anOld.ShapeType());
1238         for (; aNew.More(); aNew.Next()) {
1239           if (aNew.Current().IsSame(anOld))
1240             break;
1241         }
1242         if (aNew.More())
1243           theShapes.Append(anOld);
1244       }
1245     }
1246   }
1247   if (theShapes.IsEmpty()) { // nothing was changed
1248     if (aWasWholeContext)
1249       theShapes.Append(TopoDS_Shape());
1250     else { // if theValShape exists in new context, add it without changes, otherwise - nothing
1251       for (TopExp_Explorer aNew(aNewContShape, aValType); aNew.More(); aNew.Next()){
1252         if (aNew.Current().IsSame(theValShape)) {
1253           theShapes.Append(theValShape);
1254           break;
1255         }
1256       }
1257     }
1258   } else if (theShapes.Size() > 1) {
1259     // check it is possible to remove extra sub-shapes:
1260     // keep only shapes with the same number of containers if possible
1261     TopAbs_ShapeEnum anAncType = TopAbs_FACE;
1262     if (aValType == TopAbs_VERTEX)
1263       anAncType = TopAbs_EDGE;
1264     TopoDS_Shape anOldContext = theOldContext->shape()->impl<TopoDS_Shape>();
1265     TopTools_IndexedDataMapOfShapeListOfShape anOldMap;
1266     TopExp::MapShapesAndUniqueAncestors(anOldContext, aValType,  anAncType, anOldMap);
1267     if (anOldMap.Contains(theValShape)) {
1268       int aNumInOld = anOldMap.FindFromKey(theValShape).Extent();
1269       TopTools_IndexedDataMapOfShapeListOfShape aNewMap;
1270       TopExp::MapShapesAndUniqueAncestors(aNewContShape, aValType,  anAncType, aNewMap);
1271       TopTools_ListOfShape aNewResults;
1272       for(TopTools_ListOfShape::Iterator aNewSubs(theShapes); aNewSubs.More(); aNewSubs.Next()) {
1273         TopoDS_Shape aCand = aNewSubs.Value();
1274         if (aNewMap.Contains(aCand) && aNewMap.FindFromKey(aCand).Extent() == aNumInOld)
1275           aNewResults.Append(aCand);
1276       }
1277       if (!aNewResults.IsEmpty() && aNewResults.Size() < theShapes.Size())
1278         theShapes = aNewResults;
1279     }
1280   }
1281 }
1282
1283
1284 void Model_AttributeSelection::concealedFeature(
1285   const FeaturePtr theFeature, const FeaturePtr theStop, std::list<FeaturePtr>& theConcealers,
1286   const ResultPtr theResultOfFeature)
1287 {
1288   std::set<FeaturePtr> alreadyProcessed;
1289   alreadyProcessed.insert(theFeature);
1290   if (theStop.get())
1291     alreadyProcessed.insert(theStop);
1292   /// iterate all results to find the concealment-attribute
1293   std::list<ResultPtr> aRootRes;
1294   if (theResultOfFeature.get()) {
1295     ResultPtr aRoot = ModelAPI_Tools::bodyOwner(theResultOfFeature, true);
1296     aRootRes.push_back(aRoot ? aRoot : theResultOfFeature);
1297   } else { // all results of a feature
1298    aRootRes = theFeature->results();
1299   }
1300   std::list<ResultPtr>::const_iterator aRootIter = aRootRes.cbegin();
1301   for(; aRootIter != aRootRes.cend(); aRootIter++) {
1302     std::list<ResultPtr> allRes;
1303     allRes.push_back(*aRootIter);
1304     ResultBodyPtr aRootBody = ModelAPI_Tools::bodyOwner(*aRootIter, true);
1305     if (!aRootBody.get())
1306       aRootBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(*aRootIter);
1307     if (aRootBody.get()) {
1308       ModelAPI_Tools::allSubs(aRootBody, allRes);
1309     }
1310     for(std::list<ResultPtr>::iterator aRIter = allRes.begin(); aRIter != allRes.end(); aRIter++) {
1311       const std::set<AttributePtr>& aRefs = (*aRIter)->data()->refsToMe();
1312       std::set<AttributePtr>::const_iterator aRef = aRefs.cbegin();
1313       for (; aRef != aRefs.cend(); aRef++) {
1314         if (!aRef->get() || !(*aRef)->owner().get())
1315           continue;
1316         // concealed attribute only
1317         FeaturePtr aRefFeat = std::dynamic_pointer_cast<ModelAPI_Feature>((*aRef)->owner());
1318         if (alreadyProcessed.find(aRefFeat) != alreadyProcessed.end()) // optimization
1319           continue;
1320         alreadyProcessed.insert(aRefFeat);
1321         if (ModelAPI_Session::get()->validators()->isConcealed(aRefFeat->getKind(), (*aRef)->id()))
1322         {
1323           // for extrusion cut in python script the nested sketch reference may be concealed before
1324           // it is nested, so, check this composite feature is valid
1325           static ModelAPI_ValidatorsFactory* aFactory = ModelAPI_Session::get()->validators();
1326           // need to be validated to update the "Apply" state if not previewed
1327           if (aFactory->validate(aRefFeat)) {
1328             if (theStop.get()) {
1329               std::shared_ptr<Model_Document> aDoc =
1330                 std::dynamic_pointer_cast<Model_Document>(theStop->document());
1331               if (!aDoc->isLaterByDep(theStop, aRefFeat)) // skip feature later than stop
1332                 continue;
1333             }
1334             theConcealers.push_back(aRefFeat);
1335           }
1336         }
1337       }
1338     }
1339   }
1340 }
1341
1342 bool Model_AttributeSelection::searchNewContext(std::shared_ptr<Model_Document> theDoc,
1343   const TopoDS_Shape theContShape, ResultPtr theContext, TopoDS_Shape theValShape,
1344   TDF_Label theAccessLabel,
1345   std::list<ResultPtr>& theResults, TopTools_ListOfShape& theValShapes)
1346 {
1347   std::set<ResultPtr> aResults; // to avoid duplicates, new context, null if deleted
1348   // iterate context and shape, but also if it is sub-shape of main shape, check also it
1349   TopTools_ListOfShape aContextList;
1350   aContextList.Append(theContShape);
1351   if (theContext.get()) {
1352     ResultPtr aComposite = ModelAPI_Tools::bodyOwner(theContext);
1353     if (aComposite.get() && aComposite->shape().get() && !aComposite->shape()->isNull())
1354       aContextList.Append(aComposite->shape()->impl<TopoDS_Shape>());
1355   }
1356   for(TopTools_ListOfShape::Iterator aContIter(aContextList); aContIter.More(); aContIter.Next()) {
1357     TNaming_SameShapeIterator aModifIter(aContIter.ChangeValue(), theAccessLabel);
1358     for(; aModifIter.More(); aModifIter.Next()) {
1359       TDF_Label anObjLab = aModifIter.Label().Father();
1360       ResultPtr aModifierObj = std::dynamic_pointer_cast<ModelAPI_Result>
1361         (theDoc->objects()->object(anObjLab));
1362       if (!aModifierObj.get()) {
1363         // #2241: shape may be sub-element of new object, not main (shell created from faces)
1364         if (!anObjLab.IsRoot())
1365           aModifierObj = std::dynamic_pointer_cast<ModelAPI_Result>
1366           (theDoc->objects()->object(anObjLab.Father()));
1367         if (!aModifierObj.get())
1368           continue;
1369       }
1370       FeaturePtr aModifierFeat = theDoc->feature(aModifierObj);
1371       if (!aModifierFeat.get())
1372         continue;
1373       FeaturePtr aThisFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
1374       if (aModifierFeat == aThisFeature || !theDoc->isLaterByDep(aThisFeature, aModifierFeat))
1375         continue; // the modifier feature is later than this, so, should not be used
1376       FeaturePtr aCurrentModifierFeat = theDoc->feature(theContext);
1377       if (aCurrentModifierFeat == aModifierFeat ||
1378         !theDoc->isLaterByDep(aModifierFeat, aCurrentModifierFeat))
1379         continue; // the current modifier is later than the found, so, useless
1380       Handle(TNaming_NamedShape) aNewNS;
1381       aModifIter.Label().FindAttribute(TNaming_NamedShape::GetID(), aNewNS);
1382       if (aNewNS->Evolution() == TNaming_MODIFY || aNewNS->Evolution() == TNaming_GENERATED) {
1383         aResults.insert(aModifierObj);
1384       } else if (aNewNS->Evolution() == TNaming_DELETE) { // a shape was deleted => result is empty
1385         aResults.insert(ResultPtr());
1386       } else { // not-processed modification => don't support it
1387         continue;
1388       }
1389     }
1390   }
1391   // if there exist context composite and sub-result(s), leave only sub(s)
1392   for(std::set<ResultPtr>::iterator aResIter = aResults.begin(); aResIter != aResults.end();) {
1393     ResultPtr aParent = ModelAPI_Tools::bodyOwner(*aResIter);
1394     for(; aParent.get(); aParent = ModelAPI_Tools::bodyOwner(aParent))
1395       if (aResults.count(aParent))
1396         break;
1397     if (aParent.get()) { // erase from set, so, restart iteration
1398       aResults.erase(aParent);
1399       aResIter = aResults.begin();
1400     } else aResIter++;
1401   }
1402
1403   if (aResults.empty()) {
1404     // check the context become concealed by operation which is earlier than this selection
1405     FeaturePtr aThisFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
1406     FeaturePtr aContextOwner = theDoc->feature(theContext);
1407     std::list<FeaturePtr> aConcealers;
1408     concealedFeature(aContextOwner, aThisFeature, aConcealers, theContext);
1409     std::list<FeaturePtr>::iterator aConcealer = aConcealers.begin();
1410     for(; aConcealer != aConcealers.end(); aConcealer++) {
1411       std::list<ResultPtr> aRefResults;
1412       ModelAPI_Tools::allResults(*aConcealer, aRefResults);
1413       std::list<ResultPtr>::iterator aRefIter = aRefResults.begin();
1414       for(; aRefIter != aRefResults.end(); aRefIter++) {
1415         ResultBodyPtr aRefBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(*aRefIter);
1416         if (!aRefBody.get() || aRefBody->numberOfSubs() != 0) // iterate only leafs
1417           continue;
1418         GeomShapePtr aRefShape = aRefBody->shape();
1419         if (!aRefShape.get() || aRefShape->isNull())
1420           continue;
1421         if (aRefShape->impl<TopoDS_Shape>().IsSame(theContShape)) {
1422           // add the new context result with the same shape
1423           aResults.insert(aRefBody);
1424         }
1425       }
1426       if (aResults.empty())
1427         return true; // feature conceals result, return true, so the context will be removed
1428     }
1429     if (aResults.empty())
1430       return false; // no modifications found, must stay the same
1431   }
1432   // iterate all results to find further modifications
1433   std::set<ResultPtr>::iterator aResIter = aResults.begin();
1434   for(; aResIter != aResults.end(); aResIter++) {
1435     if (aResIter->get() != NULL) {
1436       // compute new values by two contexts: the old and the new
1437       TopTools_ListOfShape aValShapes;
1438       computeValues(theContext, *aResIter, theValShape, aValShapes);
1439
1440       TopTools_ListIteratorOfListOfShape aNewVal(aValShapes);
1441       for(; aNewVal.More(); aNewVal.Next()) {
1442         std::list<ResultPtr> aNewRes;
1443         TopTools_ListOfShape aNewUpdatedVal;
1444         TopoDS_Shape aNewValSh = aNewVal.Value();
1445         TopoDS_Shape aNewContShape = (*aResIter)->shape()->impl<TopoDS_Shape>();
1446         if (theValShape.IsNull() && aNewContShape.IsSame(aNewValSh))
1447           aNewValSh.Nullify();
1448         if (searchNewContext(theDoc, aNewContShape, *aResIter, aNewValSh,
1449                              theAccessLabel, aNewRes, aNewUpdatedVal))
1450         {
1451           // append new results instead of the current ones
1452           std::list<ResultPtr>::iterator aNewIter = aNewRes.begin();
1453           TopTools_ListIteratorOfListOfShape aNewUpdVal(aNewUpdatedVal);
1454           for(; aNewIter != aNewRes.end(); aNewIter++, aNewUpdVal.Next()) {
1455             theResults.push_back(*aNewIter);
1456             theValShapes.Append(aNewUpdVal.Value());
1457           }
1458         } else { // the current result is good
1459           theResults.push_back(*aResIter);
1460           theValShapes.Append(aNewValSh);
1461         }
1462       }
1463     }
1464   }
1465   return true; // theResults must be empty: everything is deleted
1466 }
1467
1468 void Model_AttributeSelection::updateInHistory(bool& theRemove)
1469 {
1470   static std::shared_ptr<GeomAPI_Shape> anEmptyShape;
1471
1472   ResultPtr aContext = std::dynamic_pointer_cast<ModelAPI_Result>(myRef.value());
1473   if (!aContext.get() || (aContext->groupName() != ModelAPI_ResultBody::group() &&
1474       aContext->groupName() != ModelAPI_ResultPart::group())) {
1475     // but check the case the whole results are allowed: whole features may be selected
1476     if (myParent && myParent->isWholeResultAllowed()) {
1477       FeaturePtr aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(myRef.value());
1478       if (aFeature.get()) {
1479         FeaturePtr aThisFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
1480         std::list<FeaturePtr> aConcealers;
1481         concealedFeature(aFeature, aThisFeature, aConcealers, ResultPtr());
1482         if (aConcealers.empty())
1483           return;
1484         bool aChanged = false;
1485         std::list<FeaturePtr>::iterator aConcealer = aConcealers.begin();
1486         for(; aConcealer != aConcealers.end(); aConcealer++)
1487           if (!myParent->isInList(*aConcealer, anEmptyShape)) {// avoid addition of duplicates
1488             setValue(*aConcealer, anEmptyShape);
1489             aChanged = true;
1490           }
1491         if (aConcealer == aConcealers.end()) {
1492           if (!aChanged) // remove this
1493             theRemove = true;
1494         } else { // append new
1495           for(aConcealer++; aConcealer != aConcealers.end(); aConcealer++)
1496             if (!myParent->isInList(*aConcealer, anEmptyShape)) // avoid addition of duplicates
1497               myParent->append(*aConcealer, anEmptyShape);
1498         }
1499         if (aChanged) // searching for the further modifications
1500           updateInHistory(theRemove);
1501       }
1502     }
1503     return;// only bodies and parts may be modified later in the history, skip otherwise
1504   }
1505
1506   std::shared_ptr<Model_Document> aDoc =
1507     std::dynamic_pointer_cast<Model_Document>(aContext->document());
1508   std::shared_ptr<Model_Data> aContData = std::dynamic_pointer_cast<Model_Data>(aContext->data());
1509   if (!aContData.get() || !aContData->isValid())
1510     return;
1511   TDF_Label aContLab = aContData->shapeLab(); // named shape where the selected context is located
1512
1513   // checking this may be just a reference to another context (same shape), so use that label
1514   Handle(TNaming_NamedShape) aContNS;
1515   Handle(TDF_Reference) aRefAttr;
1516   while(!aContLab.FindAttribute(TNaming_NamedShape::GetID(), aContNS) &&
1517         aContLab.FindAttribute(TDF_Reference::GetID(), aRefAttr))
1518     aContLab = aRefAttr->Get();
1519
1520   if (!aContLab.FindAttribute(TNaming_NamedShape::GetID(), aContNS)) {
1521     bool aFoundNewContext = true;
1522     ResultPtr aNewContext = aContext;
1523     while(aFoundNewContext) {
1524       aFoundNewContext = false;
1525       // parts have no shape in result, so, trace references using the Part info
1526       if (aNewContext->groupName() == ModelAPI_ResultPart::group()) {
1527         ResultPartPtr aPartContext = std::dynamic_pointer_cast<ModelAPI_ResultPart>(aNewContext);
1528         if (aPartContext.get()) { // searching for the up to date references to the referenced cont
1529           const std::set<AttributePtr>& aRefs = aPartContext->data()->refsToMe();
1530           std::set<AttributePtr>::const_iterator aRef = aRefs.begin();
1531           for(; aRef != aRefs.end(); aRef++) {
1532             // to avoid detection of part changes by local selection only
1533             AttributeSelectionPtr aSel =
1534               std::dynamic_pointer_cast<ModelAPI_AttributeSelection>(*aRef);
1535             if (aSel.get() && aSel->value().get() &&
1536                 !aSel->value()->isSame(aSel->context()->shape()))
1537               continue;
1538
1539             FeaturePtr aRefFeat = std::dynamic_pointer_cast<ModelAPI_Feature>((*aRef)->owner());
1540
1541             if (aRefFeat.get() && aRefFeat != owner() && aRefFeat->firstResult().get()) {
1542               // check the reference is concealed: #2900
1543               ModelAPI_ValidatorsFactory* aValidators = ModelAPI_Session::get()->validators();
1544               if (!aValidators->isConcealed(aRefFeat->getKind(), (*aRef)->id()))
1545                 continue;
1546               FeaturePtr aThisFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
1547               if (!aDoc->isLaterByDep(aRefFeat, aThisFeature)) { // found better feature
1548                 aFoundNewContext = true;
1549                 aNewContext = aRefFeat->firstResult();
1550               }
1551             }
1552           }
1553         }
1554       }
1555     }
1556     if (aNewContext != aContext) {
1557       setValue(aNewContext, value());
1558     }
1559     return;
1560   }
1561   FeaturePtr aThisFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(owner());
1562   FeaturePtr aCurrentModifierFeat = aDoc->feature(aContext);
1563   // iterate the context shape modifications in order to find a feature that is upper in history
1564   // that this one and is really modifies the referenced result to refer to it
1565   ResultPtr aModifierResFound;
1566   TNaming_Iterator aPairIter(aContNS);
1567   if (!aPairIter.More())
1568     return;
1569   TopoDS_Shape aNewCShape = aPairIter.NewShape();
1570   bool anIterate = true;
1571   // trying to update also the sub-shape selected
1572   GeomShapePtr aSubShape = value();
1573   if (aSubShape.get() && aSubShape->isEqual(aContext->shape()))
1574     aSubShape.reset();
1575   TopoDS_Shape aValShape;
1576   if (aSubShape.get()) {
1577     aValShape = aSubShape->impl<TopoDS_Shape>();
1578   }
1579
1580   std::list<ResultPtr> aNewContexts;
1581   TopTools_ListOfShape aValShapes;
1582   if (searchNewContext(aDoc, aNewCShape, aContext, aValShape, aContLab, aNewContexts, aValShapes))
1583   {
1584     std::set<ResultPtr> allContexts, aSkippedContext;
1585     std::list<ResultPtr>::iterator aNewContext = aNewContexts.begin();
1586     for(; aNewContext != aNewContexts.end(); aNewContext++)
1587       allContexts.insert(*aNewContext);
1588
1589     // if there exist context composite and sub-result(s), leave only sub(s)
1590     std::set<ResultPtr>::iterator aResIter = allContexts.begin();
1591     for(; aResIter != allContexts.end(); aResIter++) {
1592       ResultPtr aParent = ModelAPI_Tools::bodyOwner(*aResIter);
1593       for(; aParent.get(); aParent = ModelAPI_Tools::bodyOwner(aParent))
1594         if (allContexts.count(aParent))
1595           aSkippedContext.insert(aParent);
1596     }
1597
1598     GeomAPI_Shape::ShapeType aListShapeType = GeomAPI_Shape::SHAPE;
1599     if (myParent) {
1600       if (myParent->selectionType() == "VERTEX" || myParent->selectionType() == "Vertices")
1601         aListShapeType = GeomAPI_Shape::VERTEX;
1602       else if (myParent->selectionType() == "EDGE" || myParent->selectionType() == "Edges")
1603         aListShapeType = GeomAPI_Shape::EDGE;
1604       else if (myParent->selectionType() == "FACE" || myParent->selectionType() == "Faces")
1605         aListShapeType = GeomAPI_Shape::FACE;
1606     }
1607
1608     // issue #3031: skip topology if there is more convenient shape type presents in the
1609     // same context as a result of this
1610     bool isWholeResult = myParent && myParent->isWholeResultAllowed() && !aSubShape.get();
1611     GeomAPI_Shape::ShapeType allowedType = GeomAPI_Shape::SHAPE;
1612     if (isWholeResult) {
1613       std::list<ResultPtr>::iterator aNewCont = aNewContexts.begin();
1614       TopTools_ListIteratorOfListOfShape aNewValues(aValShapes);
1615       for(; aNewCont != aNewContexts.end(); aNewCont++, aNewValues.Next()) {
1616         if (aNewValues.Value().IsNull()) { // only for the whole context
1617           GeomAPI_Shape::ShapeType aShapeType = (*aNewCont)->shape()->shapeType();
1618           if (allowedType == GeomAPI_Shape::SHAPE) { // just set this one
1619             allowedType = aShapeType;
1620           } else {
1621             GeomAPI_Shape::ShapeType anAllowed = allowedType;
1622             if (anAllowed != aShapeType) { // select the best, nearest to the origin
1623               GeomAPI_Shape::ShapeType anOldShapeType = aContext->shape()->shapeType();
1624               GeomAPI_Shape::ShapeType aDeltaAllowed =
1625                 (GeomAPI_Shape::ShapeType)(anOldShapeType - anAllowed);
1626               if (aDeltaAllowed < 0)
1627                 aDeltaAllowed = (GeomAPI_Shape::ShapeType)(-aDeltaAllowed);
1628               GeomAPI_Shape::ShapeType aDeltaThis =
1629                 (GeomAPI_Shape::ShapeType)(anOldShapeType - aShapeType);
1630               if (aDeltaThis < 0)
1631                 aDeltaThis = (GeomAPI_Shape::ShapeType)(-aDeltaThis);
1632               if (aDeltaThis == aDeltaAllowed) { // equal distance to context, select complicated
1633                 if (anOldShapeType < anAllowed)
1634                   allowedType = aShapeType;
1635               } else if (aDeltaAllowed > aDeltaThis) { // this wins
1636                 allowedType = aShapeType;
1637               }
1638             }
1639           }
1640         }
1641       }
1642     }
1643
1644     std::list<ResultPtr>::iterator aNewCont = aNewContexts.begin();
1645     TopTools_ListIteratorOfListOfShape aNewValues(aValShapes);
1646     bool aFirst = true; // first is set to this, next are appended to parent
1647     for(; aNewCont != aNewContexts.end(); aNewCont++, aNewValues.Next()) {
1648       if (aSkippedContext.count(*aNewCont))
1649         continue;
1650
1651       if (isWholeResult && aNewValues.Value().IsNull())
1652         if (allowedType != GeomAPI_Shape::SHAPE &&
1653             (*aNewCont)->shape()->shapeType() != allowedType)
1654           continue; // there is better result exists with the better shape type (issue #3031)
1655
1656       GeomShapePtr aValueShape;
1657       if (!aNewValues.Value().IsNull()) {
1658         aValueShape = std::make_shared<GeomAPI_Shape>();
1659         aValueShape->setImpl<TopoDS_Shape>(new TopoDS_Shape(aNewValues.Value()));
1660       }
1661       // Check that list has the same type of shape selection before adding.
1662       GeomAPI_Shape::ShapeType aShapeShapeType = GeomAPI_Shape::SHAPE;
1663       if (aValueShape.get()) {
1664         aShapeShapeType = aValueShape->shapeType();
1665       } else {
1666         aShapeShapeType = (*aNewCont)->shape()->shapeType();
1667       }
1668       if (aListShapeType != GeomAPI_Shape::SHAPE && aListShapeType != aShapeShapeType) {
1669         // exception is for whole results selected
1670         if (!isWholeResult) {
1671           continue;
1672         }
1673       }
1674
1675       ResultPtr aNewContext = *aNewCont;
1676       if (aValueShape.get()) { // #2892 if context is higher level result, search this sub in lower
1677         ResultBodyPtr aBodyContext = std::dynamic_pointer_cast<ModelAPI_ResultBody>(aNewContext);
1678         if (aBodyContext.get() && aBodyContext->numberOfSubs() != 0) {
1679           std::list<ResultPtr> aLower;
1680           ModelAPI_Tools::allSubs(aBodyContext, aLower, true);
1681           for(std::list<ResultPtr>::iterator aL = aLower.begin(); aL != aLower.end(); aL++) {
1682             GeomShapePtr aLShape = (*aL)->shape();
1683             if (aLShape.get() && !aLShape->isNull()) {
1684               if (aLShape->isSubShape(aValueShape, false)) {
1685                 aNewContext = *aL;
1686                 break;
1687               }
1688             }
1689           }
1690         }
1691       }
1692
1693       if (aFirst) {
1694         if (!myParent || !myParent->isInList(aNewContext, aValueShape)) { // avoid duplicates
1695           setValue(aNewContext, aValueShape);
1696           aFirst = false;
1697         }
1698       } else if (myParent) {
1699         if (!myParent->isInList(aNewContext, aValueShape)) // avoid addition of duplicates
1700           myParent->append(aNewContext, aValueShape);
1701       }
1702     }
1703     if (aFirst) { // nothing was added, all results were deleted
1704       if (myParent) {
1705         theRemove = true;
1706       } else {
1707         static ResultPtr anEmptyContext;
1708         setValue(anEmptyContext, anEmptyShape); // nullify the selection
1709         return;
1710       }
1711     }
1712   }
1713 }
1714
1715 void Model_AttributeSelection::setParent(Model_AttributeSelectionList* theParent)
1716 {
1717   myParent = theParent;
1718 }
1719
1720 std::string Model_AttributeSelection::contextName(const TDF_Label theSelectionLab)
1721 {
1722   std::shared_ptr<Model_Document> aDoc = myRestoreDocument.get() ? myRestoreDocument :
1723     std::dynamic_pointer_cast<Model_Document>(owner()->document());
1724   FeaturePtr aFeatureOwner = aDoc->featureByLab(theSelectionLab);
1725   bool aBaseDocumnetUsed = false;
1726   if (!aFeatureOwner.get()) { // use module document
1727     aDoc = std::dynamic_pointer_cast<Model_Document>(ModelAPI_Session::get()->moduleDocument());
1728     aFeatureOwner = aDoc->featureByLab(theSelectionLab);
1729     aBaseDocumnetUsed = true;
1730   }
1731   if (aFeatureOwner.get()) {
1732     // if it is sub-element of the sketch, the context name is the name of the sketch
1733     // searching also for result - real context
1734     ResultPtr aResult;
1735     FeaturePtr aComposite = ModelAPI_Tools::compositeOwner(aFeatureOwner);
1736     if (aComposite.get() && aComposite->results().size() == 1 &&
1737         aComposite->firstResult()->groupName() == ModelAPI_ResultConstruction::group()) {
1738       aFeatureOwner = aComposite;
1739       aResult = aFeatureOwner->firstResult();
1740     } else {
1741       aResult = aDoc->resultByLab(theSelectionLab);
1742     }
1743     if (aResult.get()) {
1744       // this is to avoid duplicated names of results problem
1745       std::string aContextName = aResult->data()->name();
1746       // myLab corresponds to the current time
1747       TDF_Label aCurrentLab = selectionLabel();
1748       while(aCurrentLab.Depth() > 3)
1749         aCurrentLab = aCurrentLab.Father();
1750
1751       int aNumInHistoryNames =
1752         aDoc->numberOfNameInHistory(aResult, aCurrentLab);
1753       while(aNumInHistoryNames > 1) { // add "_" before name the needed number of times
1754         aContextName = "_" + aContextName;
1755         aNumInHistoryNames--;
1756       }
1757       if (aBaseDocumnetUsed)
1758         aContextName = aDoc->kind() + "/" + aContextName;
1759       return aContextName;
1760     }
1761   }
1762   return ""; // invalid case
1763 }
1764
1765 /// This method restores by the context and value name the context label and
1766 /// sub-label where the value is. Returns true if it is valid.
1767 bool Model_AttributeSelection::restoreContext(std::string theName,
1768   TDF_Label& theContext, TDF_Label& theValue)
1769 {
1770   static const GeomShapePtr anEmptyShape; // to store context only
1771   std::string aName = theName;
1772   std::shared_ptr<Model_Document> aDoc = myRestoreDocument.get() ? myRestoreDocument :
1773     std::dynamic_pointer_cast<Model_Document>(owner()->document());
1774
1775   // remove the sub-value part if exists
1776   std::string aSubShapeName = aName;
1777   std::string::size_type n = aName.find('/');
1778   if (n != std::string::npos) {
1779     aName = aName.substr(0, n);
1780   }
1781
1782   if (aName.empty()) return false;
1783   bool anUniqueContext = false;
1784   ResultPtr aCont = aDoc->findByName(aName, aSubShapeName, anUniqueContext);
1785   if (!aCont.get() || !aCont->shape().get() || aCont->shape()->isNull()) {
1786     // name in PartSet?
1787     aDoc = std::dynamic_pointer_cast<Model_Document>(
1788       ModelAPI_Session::get()->moduleDocument());
1789     if (theName.find(aDoc->kind()) == 0) { // remove the document identifier from name if exists
1790       aSubShapeName = theName.substr(aDoc->kind().size() + 1);
1791       aName = aSubShapeName;
1792       std::string::size_type n = aName.find('/');
1793       if (n != std::string::npos) {
1794         aName = aName.substr(0, n);
1795       }
1796     }
1797     aCont = aDoc->findByName(aName, aSubShapeName, anUniqueContext);
1798     if (!aCont.get() || !aCont->shape().get() || aCont->shape()->isNull())
1799       return false;
1800   }
1801
1802   // searching the sub-shape
1803   static const ResultPtr anEmpty;
1804   theValue = aDoc->findNamingName(aSubShapeName, anUniqueContext ? aCont : anEmpty);
1805
1806   // sketch sub-component shape and name is located in separated feature label, try the sub-name
1807   if (theValue.IsNull() && aCont->groupName() == ModelAPI_ResultConstruction::group()) {
1808     std::string::size_type aSlash = aSubShapeName.rfind('/');
1809     if (aSlash != std::string::npos) {
1810       std::string aCompName = aSubShapeName.substr(aSlash + 1);
1811       CompositeFeaturePtr aComposite =
1812         std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aDoc->feature(aCont));
1813       if (aComposite.get() && aComposite->numberOfSubs()) {
1814         const int aSubNum = aComposite->numberOfSubs();
1815         for (int a = 0; a < aSubNum && theValue.IsNull(); a++) {
1816           FeaturePtr aSub = aComposite->subFeature(a);
1817           const std::list<std::shared_ptr<ModelAPI_Result> >& aResults = aSub->results();
1818           std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aRes = aResults.cbegin();
1819           for (; aRes != aResults.cend() && theValue.IsNull(); aRes++) {
1820             if ((*aRes)->data()->name() == aCompName) {
1821               theValue = std::dynamic_pointer_cast<Model_Data>((*aRes)->data())->shapeLab();
1822               break;
1823             } else { // any sub-label because the sketch line may be renamed, but not sub-vertices
1824               TDF_Label aLab = std::dynamic_pointer_cast<Model_Data>((*aRes)->data())->shapeLab();
1825               TDF_ChildIDIterator aSubNames(aLab, TDataStd_Name::GetID());
1826               for(; aSubNames.More(); aSubNames.Next()) {
1827                 if (Handle(TDataStd_Name)::DownCast(aSubNames.Value())->Get().
1828                   IsEqual(aCompName.c_str())) {
1829                   theValue = aSubNames.Value()->Label();
1830                   break;
1831                 }
1832               }
1833             }
1834           }
1835         }
1836       }
1837     }
1838   }
1839
1840   if (aCont.get()) {
1841     theContext = std::dynamic_pointer_cast<Model_Data>(aCont->data())->label();
1842   }
1843   return true;
1844 }
1845
1846 bool Model_AttributeSelection::isLater(
1847   const TDF_Label theResult1, const TDF_Label theResult2) const
1848 {
1849   std::shared_ptr<Model_Document> aDoc = myRestoreDocument.get() ? myRestoreDocument :
1850     std::dynamic_pointer_cast<Model_Document>(owner()->document());
1851   FeaturePtr aFeat1 = aDoc->featureByLab(theResult1);
1852   if (!aFeat1.get())
1853     return false;
1854   FeaturePtr aFeat2 = aDoc->featureByLab(theResult2);
1855   if (!aFeat2.get())
1856     return false;
1857   return aDoc->isLaterByDep(aFeat1, aFeat2);
1858 }
1859
1860 ResultPtr Model_AttributeSelection::newestContext(
1861   const ResultPtr theCurrent, const GeomShapePtr theValue)
1862 {
1863   ResultPtr aResult = theCurrent;
1864   GeomShapePtr aSelectedShape = theValue.get() ? theValue : theCurrent->shape();
1865   std::shared_ptr<Model_Document> aDoc =
1866     std::dynamic_pointer_cast<Model_Document>(owner()->document());
1867   bool aFindNewContext = true;
1868   while (aFindNewContext && aResult.get()) {
1869     aFindNewContext = false;
1870     // try to find the last context to find the up to date shape
1871     TopoDS_Shape aConShape = aResult->shape()->impl<TopoDS_Shape>();
1872     if (TNaming_Tool::HasLabel(selectionLabel(), aConShape)) {
1873       Handle(TNaming_NamedShape) aNS = TNaming_Tool::NamedShape(aConShape, selectionLabel());
1874       if (!aNS.IsNull()) {
1875         aNS = TNaming_Tool::CurrentNamedShape(aNS);
1876         if (!aNS.IsNull()
1877             && isLater(selectionLabel(), aNS->Label())
1878             && isLater(aNS->Label(),
1879                        std::dynamic_pointer_cast<Model_Data>(aResult->data())->label()))
1880         {
1881           TDF_Label aLab = aNS->Label();
1882           ResultPtr aRes = aDoc->resultByLab(aLab);
1883           if (aRes.get()) {
1884             if (aRes->shape()->isSubShape(aSelectedShape)) {
1885               aResult = aRes;
1886               aFindNewContext = true;
1887               continue;
1888             }
1889           }
1890         }
1891       }
1892     }
1893
1894     // TestFillWireVertex.py - sketch constructions for wire may participate too
1895     //if (aResult->groupName() == ModelAPI_ResultBody::group()) {
1896       // try to search newer context by the concealment references
1897       // take references to all results: root one, any sub
1898     std::list<ResultPtr> allRes;
1899     ResultPtr aCompContext;
1900     ResultBodyPtr aCompBody = ModelAPI_Tools::bodyOwner(aResult, true);
1901     if (aCompBody.get()) {
1902       ModelAPI_Tools::allSubs(aCompBody, allRes);
1903       allRes.push_back(aCompBody);
1904       aCompContext = aCompBody;
1905     }
1906     if (allRes.empty())
1907       allRes.push_back(aResult);
1908
1909     for (std::list<ResultPtr>::iterator aSub = allRes.begin(); aSub != allRes.end(); aSub++) {
1910       ResultPtr aResCont = *aSub;
1911       ResultBodyPtr aResBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(aResCont);
1912       if (aResBody.get() && aResBody->numberOfSubs() > 0 && aResBody != aCompContext)
1913         continue; // only lower and higher level subs are counted
1914       const std::set<AttributePtr>& aRefs = aResCont->data()->refsToMe();
1915       std::set<AttributePtr>::const_iterator aRef = aRefs.begin();
1916       for (; !aFindNewContext && aRef != aRefs.end(); aRef++) {
1917         if (!aRef->get() || !(*aRef)->owner().get())
1918           continue;
1919         // concealed attribute only
1920         FeaturePtr aRefFeat = std::dynamic_pointer_cast<ModelAPI_Feature>((*aRef)->owner());
1921         if (!ModelAPI_Session::get()->validators()->isConcealed(
1922           aRefFeat->getKind(), (*aRef)->id()))
1923           continue;
1924         // search the feature result that contains sub-shape selected
1925         std::list<std::shared_ptr<ModelAPI_Result> > aResults;
1926
1927         // take all sub-results or one result
1928         std::list<ResultPtr> aRefFeatResults;
1929         ModelAPI_Tools::allResults(aRefFeat, aRefFeatResults);
1930         std::list<ResultPtr>::iterator aRefResIter = aRefFeatResults.begin();
1931         for (; aRefResIter != aRefFeatResults.end(); aRefResIter++) {
1932           ResultBodyPtr aBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(*aRefResIter);
1933           if (aBody.get() && aBody->numberOfSubs() == 0) // add only lower level subs
1934             aResults.push_back(aBody);
1935         }
1936         std::list<std::shared_ptr<ModelAPI_Result> >::iterator aResIter = aResults.begin();
1937
1938         // searching by sub-shape
1939         for (; aResIter != aResults.end(); aResIter++) {
1940           if (!aResIter->get() || !(*aResIter)->data()->isValid() || (*aResIter)->isDisabled())
1941             continue;
1942           GeomShapePtr aShape = (*aResIter)->shape();
1943           if (aShape.get() && aShape->isSubShape(aSelectedShape, false)) {
1944             aResult = *aResIter; // found new context (produced from this) with same subshape
1945             aFindNewContext = true; // continue searching further
1946             break;
1947           }
1948         }
1949       }
1950     }
1951   }
1952   // if compsolid is context, try to take sub-solid as context: like in GUI and scripts
1953   ResultBodyPtr aComp = std::dynamic_pointer_cast<ModelAPI_ResultBody>(aResult);
1954   if (aComp && aComp->numberOfSubs()) {
1955     std::list<ResultPtr> allSubs;
1956     ModelAPI_Tools::allSubs(aComp, allSubs);
1957     std::list<ResultPtr>::iterator aS = allSubs.begin();
1958     for (; aS != allSubs.end(); aS++) {
1959       ResultBodyPtr aSub = std::dynamic_pointer_cast<ModelAPI_ResultBody>(*aS);
1960       if (aSub && aSub->numberOfSubs() == 0 && aSub->shape().get() &&
1961           aSub->shape()->isSubShape(aSelectedShape)) {
1962         aResult = aSub;
1963         break;
1964       }
1965     }
1966   }
1967   // in case sketch line was selected for wire, but wire was concealed and not such line anymore,
1968   // so, actually, the sketch element was selected (which is never concealed)
1969   if (aResult != theCurrent && aResult->isConcealed())
1970     aResult = theCurrent;
1971   return aResult;
1972 }
1973
1974 void Model_AttributeSelection::combineGeometrical()
1975 {
1976   if (myTmpContext.get() || myTmpSubShape.get())
1977     return;
1978   TDF_Label aSelLab = selectionLabel();
1979   if (aSelLab.IsAttribute(kINVALID_SELECTION) || !myRef.isInitialized())
1980     return;
1981
1982   if (aSelLab.IsAttribute(kCIRCLE_CENTER) || aSelLab.IsAttribute(kELLIPSE_CENTER1) ||
1983       aSelLab.IsAttribute(kELLIPSE_CENTER2) || aSelLab.IsAttribute(kSIMPLE_REF_ID))
1984     return;
1985
1986   if (aSelLab.IsAttribute(kPART_REF_ID)) {
1987     ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(context());
1988     if (!aPart.get() || !aPart->isActivated())
1989       return; // postponed naming needed
1990     Handle(TDataStd_Integer) anIndex;
1991     if (aSelLab.FindAttribute(TDataStd_Integer::GetID(), anIndex)) {
1992       if (anIndex->Get()) { // special selection attribute was created, use it
1993         std::string aNewName;
1994         aPart->combineGeometrical(anIndex->Get(), aNewName);
1995         TDataStd_Name::Set(aSelLab, aNewName.c_str());
1996       }
1997     }
1998     return;
1999   }
2000
2001   std::shared_ptr<Model_ResultConstruction> aConstr =
2002     std::dynamic_pointer_cast<Model_ResultConstruction>(context());
2003   if (aConstr.get())
2004     return;
2005   FeaturePtr aFeature = contextFeature();
2006   if (aFeature.get())
2007     return;
2008
2009   Selector_Selector aSelector(aSelLab, baseDocumentLab());
2010   TopoDS_Shape aContextShape = context()->shape()->impl<TopoDS_Shape>();
2011   if (aSelector.restore(aContextShape)) {
2012     aSelector.combineGeometrical(aContextShape);
2013   }
2014 }
2015
2016 TDF_Label Model_AttributeSelection::baseDocumentLab()
2017 {
2018   if (ModelAPI_Session::get()->moduleDocument() != owner()->document())
2019     return std::dynamic_pointer_cast<Model_Document>
2020       (ModelAPI_Session::get()->moduleDocument())->extConstructionsLabel();
2021   static TDF_Label anEmpty;
2022   return anEmpty;
2023 }