Salome HOME
Update of icons: recover, split, construction planes subtypes
[modules/shaper.git] / src / ModelAPI / ModelAPI_Tools.cpp
old mode 100644 (file)
new mode 100755 (executable)
index 8656f2f..4e094e6
-// File:        ModelAPI_Tools.cpp
-// Created:     06 Aug 2014
-// Author:      Vitaly Smetannikov
-
-#include "ModelAPI_Tools.h"
-#include <ModelAPI_ResultBody.h>
-#include <ModelAPI_ResultConstruction.h>
-#include <ModelAPI_ResultGroup.h>
-
-namespace ModelAPI_Tools {
-
-  boost::shared_ptr<GeomAPI_Shape> shape(const ResultPtr& theResult)
-  {
-
-    ResultBodyPtr aBody = boost::dynamic_pointer_cast<ModelAPI_ResultBody>(theResult);
-    if (aBody)
-      return aBody->shape();
-
-    ResultConstructionPtr aConstruct = boost::dynamic_pointer_cast<ModelAPI_ResultConstruction>(
-      theResult);
-    if (aConstruct)
-      return aConstruct->shape();
-
-    ResultGroupPtr aGroup = boost::dynamic_pointer_cast<ModelAPI_ResultGroup>(theResult);
-    if (aGroup)
-      return aGroup->shape();
-    return boost::shared_ptr<GeomAPI_Shape>();
-  }
-}
+// Copyright (C) 2014-20xx CEA/DEN, EDF R&D\r
+\r
+// File:        ModelAPI_Tools.cpp\r
+// Created:     06 Aug 2014\r
+// Author:      Vitaly Smetannikov\r
+\r
+#include "ModelAPI_Tools.h"\r
+#include <ModelAPI_Session.h>\r
+#include <ModelAPI_Document.h>\r
+#include <ModelAPI_Object.h>\r
+#include <ModelAPI_AttributeDouble.h>\r
+#include <ModelAPI_ResultParameter.h>\r
+#include <ModelAPI_ResultPart.h>\r
+#include <ModelAPI_AttributeDocRef.h>\r
+#include <ModelAPI_Validator.h>\r
+#include <list>\r
+#include <map>\r
+#include <iostream>\r
+\r
+#include <Events_Loop.h>\r
+#include <ModelAPI_Events.h>\r
+\r
+#define RECURSE_TOP_LEVEL 50\r
+\r
+//#define DEBUG_REMOVE_FEATURES\r
+//#define DEBUG_REMOVE_FEATURES_RECURSE\r
+//#define DEBUG_CYCLING_1550\r
+\r
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE\r
+#include <sstream>\r
+std::string getFeatureInfo(FeaturePtr theFeature)\r
+{\r
+  if (!theFeature.get())\r
+    return "";\r
+  //std::ostringstream aPtrStr;\r
+  //aPtrStr << "[" << theFeature.get() << "] ";\r
+  std::string aFeatureInfo = /*aPtrStr.str() + */theFeature->name();\r
+  CompositeFeaturePtr aComposite = ModelAPI_Tools::compositeOwner(theFeature);\r
+  if (aComposite.get()) {\r
+      aFeatureInfo = aFeatureInfo + "[in " + aComposite->name() + "]";\r
+  }\r
+  return aFeatureInfo;\r
+}\r
+#endif\r
+\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+void printMapInfo(const std::map<FeaturePtr, std::set<FeaturePtr> >& theMainList,\r
+                  const std::string& thePrefix)\r
+{\r
+  std::map<FeaturePtr, std::set<FeaturePtr> >::const_iterator aMainIt = theMainList.begin(),\r
+                                                              aMainLast = theMainList.end();\r
+  std::string anInfo;\r
+  for (; aMainIt != aMainLast; aMainIt++) {\r
+    FeaturePtr aMainListFeature = aMainIt->first;\r
+    std::set<FeaturePtr> aMainRefList = aMainIt->second;\r
+    std::set<FeaturePtr>::const_iterator anIt = aMainRefList.begin(), aLast = aMainRefList.end();\r
+    std::string aRefsInfo;\r
+    for (; anIt != aLast; anIt++) {\r
+      aRefsInfo += (*anIt)->name().c_str();\r
+      if (anIt != aLast)\r
+        aRefsInfo += ", ";\r
+    }\r
+    if (!aRefsInfo.empty()) {\r
+      anInfo = anInfo + aMainListFeature->name().c_str() + ": " + aRefsInfo + "\n";\r
+    }\r
+  }\r
+  std::cout << thePrefix.c_str() << " [feature: references to]: \n" << anInfo.c_str() << std::endl;\r
+}\r
+\r
+void printListInfo(const std::set<FeaturePtr>& theMainList,\r
+                  const std::string& thePrefix)\r
+{\r
+  std::set<FeaturePtr>::const_iterator aMainIt = theMainList.begin(),\r
+                                       aMainLast = theMainList.end();\r
+  std::string anInfo;\r
+  for (; aMainIt != aMainLast; aMainIt++) {\r
+    FeaturePtr aRefFeature = *aMainIt;\r
+    anInfo += aRefFeature->name().c_str();\r
+    if (aMainIt != aMainLast)\r
+      anInfo += ", ";\r
+  }\r
+  std::cout << thePrefix.c_str() << ": " << anInfo.c_str() << std::endl;\r
+}\r
+#endif\r
+\r
+namespace ModelAPI_Tools {\r
+\r
+std::shared_ptr<GeomAPI_Shape> shape(const ResultPtr& theResult)\r
+{\r
+  return theResult->shape();\r
+}\r
+\r
+void shapesOfType(const FeaturePtr& theFeature,\r
+                  const GeomAPI_Shape::ShapeType& theType,\r
+                  std::set<GeomShapePtr>& theShapes)\r
+{\r
+  std::list<ResultPtr> aResults = theFeature->results();\r
+  std::list<ResultPtr>::const_iterator aRIter = aResults.cbegin();\r
+  for (; aRIter != aResults.cend(); aRIter++) {\r
+    ResultPtr aResult = *aRIter;\r
+    GeomShapePtr aShape = aResult->shape();\r
+    if (aShape.get() && aShape->shapeType() == theType)\r
+      theShapes.insert(aShape);\r
+  }\r
+}\r
+\r
+const char* toString(ModelAPI_ExecState theExecState) \r
+{\r
+#define TO_STRING(__NAME__) case __NAME__: return #__NAME__;\r
+  switch (theExecState) {\r
+  TO_STRING(ModelAPI_StateDone)\r
+  TO_STRING(ModelAPI_StateMustBeUpdated)\r
+  TO_STRING(ModelAPI_StateExecFailed)\r
+  TO_STRING(ModelAPI_StateInvalidArgument)\r
+  TO_STRING(ModelAPI_StateNothing)\r
+  default: return "Unknown ExecState.";\r
+  }\r
+#undef TO_STRING\r
+}\r
+\r
+std::string getFeatureError(const FeaturePtr& theFeature)\r
+{\r
+  std::string anError;\r
+  if (!theFeature.get() || !theFeature->data()->isValid() || theFeature->isAction())\r
+    return anError;\r
+\r
+  // to be removed later, this error should be got from the feature\r
+  if (theFeature->data()->execState() == ModelAPI_StateDone ||\r
+      theFeature->data()->execState() == ModelAPI_StateMustBeUpdated)\r
+    return anError;\r
+\r
+  // set error indication\r
+  anError = theFeature->error();\r
+  if (anError.empty()) {\r
+    bool isDone = ( theFeature->data()->execState() == ModelAPI_StateDone\r
+                 || theFeature->data()->execState() == ModelAPI_StateMustBeUpdated );\r
+    if (!isDone) {\r
+      anError = toString(theFeature->data()->execState());\r
+      // If the feature is Composite and error is StateInvalidArgument,\r
+      // error text should include error of first invalid sub-feature. Otherwise\r
+      // it is not clear what is the reason of the invalid argument.\r
+      if (theFeature->data()->execState() == ModelAPI_StateInvalidArgument) {\r
+        CompositeFeaturePtr aComposite =\r
+                    std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(theFeature);\r
+        if (aComposite) {\r
+          for (int i = 0, aSize = aComposite->numberOfSubs(); i < aSize; i++) {\r
+            FeaturePtr aSubFeature = aComposite->subFeature(i);\r
+            std::string aSubFeatureError = getFeatureError(aSubFeature);\r
+            if (!aSubFeatureError.empty()) {\r
+              anError = anError + " in " + aSubFeature->getKind() + ".\n" + aSubFeatureError;\r
+              break;\r
+            }\r
+          }\r
+        }\r
+      }\r
+    }\r
+  }\r
+\r
+  return anError;\r
+}\r
+\r
+ObjectPtr objectByName(const DocumentPtr& theDocument, const std::string& theGroup, const std::string& theName)\r
+{\r
+  for (int anIndex = 0; anIndex < theDocument->size(theGroup); ++anIndex) {\r
+    ObjectPtr anObject = theDocument->object(theGroup, anIndex);\r
+    if (anObject->data()->name() == theName)\r
+      return anObject;\r
+  }\r
+  // not found\r
+  return ObjectPtr();\r
+}\r
+\r
+bool findVariable(const DocumentPtr& theDocument, FeaturePtr theSearcher,\r
+                  const std::string& theName, double& outValue, ResultParameterPtr& theParam)\r
+{\r
+  ObjectPtr aParamObj = objectByName(theDocument, ModelAPI_ResultParameter::group(), theName);\r
+  theParam = std::dynamic_pointer_cast<ModelAPI_ResultParameter>(aParamObj);\r
+  if (!theParam.get())\r
+    return false;\r
+  // avoid usage of parameters created later than the initial parameter\r
+  if (theSearcher.get() && theDocument->isLater(theDocument->feature(theParam), theSearcher))\r
+    return false;\r
+  AttributeDoublePtr aValueAttribute = theParam->data()->real(ModelAPI_ResultParameter::VALUE());\r
+  outValue = aValueAttribute->value();\r
+  return true;\r
+}\r
+\r
+bool findVariable(FeaturePtr theSearcher, const std::string& theName, double& outValue, ResultParameterPtr& theParam,\r
+                  const DocumentPtr& theDocument)\r
+{\r
+  SessionPtr aSession = ModelAPI_Session::get();\r
+  std::list<DocumentPtr> aDocList;\r
+  DocumentPtr aDocument = theDocument.get() ? theDocument : aSession->activeDocument();\r
+  if (findVariable(aDocument, theSearcher, theName, outValue, theParam))\r
+    return true;\r
+  DocumentPtr aRootDocument = aSession->moduleDocument();\r
+  if (aDocument != aRootDocument) {\r
+    // any parameters in PartSet is okindependently on the Part position (issu #1504)\r
+    if (findVariable(aRootDocument, FeaturePtr(), theName, outValue, theParam))\r
+      return true;\r
+  }\r
+  return false;\r
+}\r
+\r
+static std::map<int, std::vector<int> > myColorMap;\r
+\r
+void appendValues(std::vector<int>& theRGB, const int theRed, const int theGreen, const int theBlue)\r
+{\r
+  theRGB.push_back(theRed);\r
+  theRGB.push_back(theGreen);\r
+  theRGB.push_back(theBlue);\r
+}\r
+\r
+bool containsValues(std::map<int, std::vector<int> >& theColorMap, std::vector<int>& theValues)\r
+{\r
+  std::map<int, std::vector<int> >::const_iterator anIt = theColorMap.begin(), aLast = theColorMap.end();\r
+  bool isFound = false;\r
+  for (; anIt != aLast && !isFound; anIt++) {\r
+    std::vector<int> aValues = anIt->second;\r
+    isFound = aValues[0] == theValues[0] &&\r
+              aValues[1] == theValues[1] &&\r
+              aValues[2] == theValues[2];\r
+  }\r
+  return isFound;\r
+}\r
+\r
+std::vector<int> HSVtoRGB(int theH, int theS, int theV)\r
+{\r
+  std::vector<int> aRGB;\r
+  if (theH < 0 || theH > 360 ||\r
+      theS < 0 || theS > 100 ||\r
+      theV < 0 || theV > 100)\r
+    return aRGB;\r
+\r
+  int aHi = (int)theH/60;\r
+\r
+  double aV = theV;\r
+  double aVmin = (100 - theS)*theV/100;\r
+\r
+  double anA = (theV - aVmin)* (theH % 60) / 60;\r
+\r
+  double aVinc = aVmin + anA;\r
+  double aVdec = theV - anA;\r
+\r
+  double aPercentToValue = 255./100;\r
+  int aV_int    = (int)(aV*aPercentToValue);\r
+  int aVinc_int = (int)(aVinc*aPercentToValue);\r
+  int aVmin_int = (int)(aVmin*aPercentToValue);\r
+  int aVdec_int = (int)(aVdec*aPercentToValue);\r
+\r
+  switch(aHi) {\r
+    case 0: appendValues(aRGB, aV_int,    aVinc_int, aVmin_int); break;\r
+    case 1: appendValues(aRGB, aVdec_int, aV_int,    aVmin_int); break;\r
+    case 2: appendValues(aRGB, aVmin_int, aV_int,    aVinc_int); break;\r
+    case 3: appendValues(aRGB, aVmin_int, aVdec_int, aV_int); break;\r
+    case 4: appendValues(aRGB, aVinc_int, aVmin_int, aV_int); break;\r
+    case 5: appendValues(aRGB, aV_int,    aVmin_int, aVdec_int); break;\r
+    default: break;\r
+  }\r
+  return aRGB;\r
+}\r
+\r
+\r
+void fillColorMap()\r
+{\r
+  if (!myColorMap.empty())\r
+    return;\r
+\r
+  int i = 0;\r
+  for (int s = 100; s > 0; s = s - 50)\r
+  {\r
+    for (int v = 100; v >= 40; v = v - 20)\r
+    {\r
+      for (int h = 0; h < 359 ; h = h + 60)\r
+      {\r
+        std::vector<int> aColor = HSVtoRGB(h, s, v);\r
+        if (containsValues(myColorMap, aColor))\r
+          continue;\r
+        myColorMap[i] = aColor;\r
+        i++;\r
+      }\r
+    }\r
+  }\r
+}\r
+\r
+void findRandomColor(std::vector<int>& theValues)\r
+{\r
+  theValues.clear();\r
+  if (myColorMap.empty()) {\r
+    fillColorMap();\r
+  }\r
+\r
+  size_t aSize = myColorMap.size();\r
+  int anIndex = rand() % aSize;\r
+  if (myColorMap.find(anIndex) != myColorMap.end()) {\r
+    theValues = myColorMap.at(anIndex);\r
+  }\r
+}\r
+\r
+ResultPtr findPartResult(const DocumentPtr& theMain, const DocumentPtr& theSub)\r
+{\r
+  if (theMain != theSub) { // to optimize and avoid of crash on partset document close (don't touch the sub-document structure)\r
+    for (int a = theMain->size(ModelAPI_ResultPart::group()) - 1; a >= 0; a--) {\r
+      ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(\r
+          theMain->object(ModelAPI_ResultPart::group(), a));\r
+      if (aPart && aPart->data()->document(ModelAPI_ResultPart::DOC_REF())->value() == theSub) {\r
+        return aPart;\r
+      }\r
+    }\r
+  }\r
+  return ResultPtr();\r
+}\r
+\r
+FeaturePtr findPartFeature(const DocumentPtr& theMain, const DocumentPtr& theSub)\r
+{\r
+  if (theMain != theSub) { // to optimize and avoid of crash on partset document close (don't touch the sub-document structure)\r
+    for (int a = theMain->size(ModelAPI_Feature::group()) - 1; a >= 0; a--) {\r
+      FeaturePtr aPartFeat = std::dynamic_pointer_cast<ModelAPI_Feature>(\r
+          theMain->object(ModelAPI_Feature::group(), a));\r
+      if (aPartFeat.get()) {\r
+        const std::list<std::shared_ptr<ModelAPI_Result> >& aResList = aPartFeat->results();\r
+        std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aRes = aResList.begin();\r
+        for(; aRes != aResList.end(); aRes++) {\r
+          ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(*aRes);\r
+          if (aPart.get()) {\r
+            if (aPart->isActivated() && aPart->partDoc() == theSub)\r
+              return aPartFeat;\r
+          } else break; // if the first is not Part, others are also not\r
+        }\r
+      }\r
+    }\r
+  }\r
+  return FeaturePtr();\r
+}\r
+\r
+CompositeFeaturePtr compositeOwner(const FeaturePtr& theFeature)\r
+{\r
+  if (theFeature.get() && theFeature->data()->isValid()) {\r
+    const std::set<std::shared_ptr<ModelAPI_Attribute> >& aRefs = theFeature->data()->refsToMe();\r
+    std::set<std::shared_ptr<ModelAPI_Attribute> >::const_iterator aRefIter = aRefs.begin();\r
+    for(; aRefIter != aRefs.end(); aRefIter++) {\r
+      CompositeFeaturePtr aComp = std::dynamic_pointer_cast<ModelAPI_CompositeFeature>\r
+        ((*aRefIter)->owner());\r
+      if (aComp.get() && aComp->data()->isValid() && aComp->isSub(theFeature))\r
+        return aComp;\r
+    }\r
+  }\r
+  return CompositeFeaturePtr(); // not found\r
+}\r
+\r
+ResultCompSolidPtr compSolidOwner(const ResultPtr& theSub)\r
+{\r
+  ResultBodyPtr aBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(theSub);\r
+  if (aBody.get()) {\r
+    FeaturePtr aFeatureOwner = aBody->document()->feature(aBody);\r
+    if (aFeatureOwner.get()) {\r
+      std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aResIter = \r
+        aFeatureOwner->results().cbegin();\r
+      for(; aResIter != aFeatureOwner->results().cend(); aResIter++) {\r
+        ResultCompSolidPtr aComp = std::dynamic_pointer_cast<ModelAPI_ResultCompSolid>(*aResIter);\r
+        if (aComp && aComp->isSub(aBody))\r
+          return aComp;\r
+      }\r
+    }\r
+  }\r
+  return ResultCompSolidPtr(); // not found\r
+}\r
+\r
+bool hasSubResults(const ResultPtr& theResult)\r
+{\r
+  ResultCompSolidPtr aCompSolid = std::dynamic_pointer_cast<ModelAPI_ResultCompSolid>(theResult);\r
+  return aCompSolid.get() && aCompSolid->numberOfSubs() > 0;\r
+}\r
+\r
+void allResults(const FeaturePtr& theFeature, std::list<ResultPtr>& theResults)\r
+{\r
+  if (!theFeature.get()) // safety: for empty feature no results\r
+    return;\r
+  const std::list<std::shared_ptr<ModelAPI_Result> >& aResults = theFeature->results();\r
+  std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aRIter = aResults.begin();\r
+  for (; aRIter != aResults.cend(); aRIter++) {\r
+    theResults.push_back(*aRIter);\r
+    // iterate sub-bodies of compsolid\r
+    ResultCompSolidPtr aComp = std::dynamic_pointer_cast<ModelAPI_ResultCompSolid>(*aRIter);\r
+    if (aComp.get()) {\r
+      int aNumSub = aComp->numberOfSubs();\r
+      for(int a = 0; a < aNumSub; a++) {\r
+        theResults.push_back(aComp->subResult(a));\r
+      }\r
+    }\r
+  }\r
+}\r
+\r
+//******************************************************************\r
+bool allDocumentsActivated(std::string& theNotActivatedNames)\r
+{\r
+  theNotActivatedNames = "";\r
+  bool anAllPartActivated = true;\r
+\r
+  DocumentPtr aRootDoc = ModelAPI_Session::get()->moduleDocument();\r
+  int aSize = aRootDoc->size(ModelAPI_ResultPart::group());\r
+  for (int i = 0; i < aSize; i++) {\r
+    ObjectPtr aObject = aRootDoc->object(ModelAPI_ResultPart::group(), i);\r
+    ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(aObject);\r
+    if (!aPart->isActivated()) {\r
+      anAllPartActivated = false;\r
+      if (!theNotActivatedNames.empty())\r
+        theNotActivatedNames += ", ";\r
+      theNotActivatedNames += aObject->data()->name().c_str();\r
+    }\r
+  }\r
+  return anAllPartActivated;\r
+}\r
+\r
+bool removeFeaturesAndReferences(const std::set<FeaturePtr>& theFeatures,\r
+                                 const bool theFlushRedisplay,\r
+                                 const bool theUseComposite,\r
+                                 const bool theUseRecursion)\r
+{\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+  printListInfo(theFeatures, "selection: ");\r
+#endif\r
+\r
+  std::map<FeaturePtr, std::set<FeaturePtr> > aReferences;\r
+  ModelAPI_Tools::findAllReferences(theFeatures, aReferences, theUseComposite, theUseRecursion);\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+  printMapInfo(aReferences, "allDependencies: ");\r
+#endif\r
+\r
+  std::set<FeaturePtr> aFeaturesRefsTo;\r
+  ModelAPI_Tools::findRefsToFeatures(theFeatures, aReferences, aFeaturesRefsTo);\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+  printListInfo(aFeaturesRefsTo, "references: ");\r
+#endif\r
+\r
+  std::set<FeaturePtr> aFeatures = theFeatures;\r
+  if (!aFeaturesRefsTo.empty())\r
+    aFeatures.insert(aFeaturesRefsTo.begin(), aFeaturesRefsTo.end());\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+  printListInfo(aFeatures, "removeFeatures: ");\r
+#endif\r
+\r
+  return ModelAPI_Tools::removeFeatures(aFeatures, false);\r
+}\r
+\r
+bool removeFeatures(const std::set<FeaturePtr>& theFeatures,\r
+                    const bool theFlushRedisplay)\r
+{\r
+  bool isDone = false;\r
+  std::set<FeaturePtr>::const_iterator anIt = theFeatures.begin(),\r
+                                       aLast = theFeatures.end();\r
+  for (; anIt != aLast; anIt++) {\r
+    FeaturePtr aFeature = *anIt;\r
+    if (aFeature.get()) {\r
+      DocumentPtr aDoc = aFeature->document();\r
+      // flush REDISPLAY signal after remove feature\r
+      aDoc->removeFeature(aFeature);\r
+      isDone = true;\r
+    }\r
+  }\r
+  if (isDone && theFlushRedisplay) {\r
+    // the redisplay signal should be flushed in order to erase the feature presentation in the viewer\r
+    // if should be done after removeFeature() of document\r
+    Events_Loop::loop()->flush(Events_Loop::loop()->eventByName(EVENT_OBJECT_TO_REDISPLAY));\r
+  }\r
+  return true;\r
+}\r
+\r
+// Fills the references list by all references of the feature from the references map.\r
+// This is a recusive method to find references by next found feature in the map of references.\r
+// \param theFeature a feature to find references\r
+// \param theReferencesMap a map of references\r
+// \param theReferences an out container of references\r
+void addRefsToFeature(const FeaturePtr& theFeature,\r
+                      const std::map<FeaturePtr, std::set<FeaturePtr> >& theReferencesMap,\r
+                      std::map<FeaturePtr, std::set<FeaturePtr> >& theProcessedReferences,\r
+                      int theRecLevel,\r
+                      std::set<FeaturePtr>& theReferences)\r
+{\r
+  if (theRecLevel > RECURSE_TOP_LEVEL)\r
+    return;\r
+  theRecLevel++;\r
+\r
+  // if the feature is already processed, get the ready references from the map\r
+  if (theProcessedReferences.find(theFeature) != theProcessedReferences.end()) {\r
+    std::set<FeaturePtr> aReferences = theProcessedReferences.at(theFeature);\r
+    theReferences.insert(aReferences.begin(), aReferences.end());\r
+    return;\r
+  }\r
+\r
+  if (theReferencesMap.find(theFeature) == theReferencesMap.end())\r
+    return; // this feature is not in the selection list, so exists without references to it\r
+  std::set<FeaturePtr> aMainReferences = theReferencesMap.at(theFeature);\r
+\r
+  std::set<FeaturePtr>::const_iterator anIt = aMainReferences.begin(),\r
+                                       aLast = aMainReferences.end();\r
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE\r
+  std::string aSpacing;\r
+  for (int i = 0; i < theRecLevel; i++)\r
+    aSpacing.append(" ");\r
+#endif\r
+\r
+  for (; anIt != aLast; anIt++) {\r
+    FeaturePtr aRefFeature = *anIt;\r
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE\r
+  std::cout << aSpacing << " Ref: " << getFeatureInfo(aRefFeature) << std::endl;\r
+#endif\r
+    if (theReferences.find(aRefFeature) == theReferences.end())\r
+      theReferences.insert(aRefFeature);\r
+    addRefsToFeature(aRefFeature, theReferencesMap, theProcessedReferences, theRecLevel, theReferences);\r
+  }\r
+}\r
+\r
+// For each feature from the feature list it searches references to the feature and append them\r
+// to the references map. This is a recusive method.\r
+// \param theFeature a feature to find references\r
+// \param theReferencesMap a map of references\r
+// \param theReferences an out container of references\r
+void findReferences(const std::set<FeaturePtr>& theFeatures,\r
+                    std::map<FeaturePtr, std::set<FeaturePtr> >& theReferences,\r
+                    const bool theUseComposite, const bool theUseRecursion, int theRecLevel)\r
+{\r
+  if (theRecLevel > RECURSE_TOP_LEVEL)\r
+    return;\r
+  theRecLevel++;\r
+  std::set<FeaturePtr>::const_iterator anIt = theFeatures.begin(),\r
+                                        aLast = theFeatures.end();\r
+  for (; anIt != aLast; anIt++) {\r
+    FeaturePtr aFeature = *anIt;\r
+    if (aFeature.get() && theReferences.find(aFeature) == theReferences.end()) {\r
+      DocumentPtr aSelFeatureDoc = aFeature->document();\r
+      std::set<FeaturePtr> aSelRefFeatures;\r
+      aSelFeatureDoc->refsToFeature(aFeature, aSelRefFeatures, false/*do not emit signals*/);\r
+      if (theUseComposite) { // do not filter selection\r
+        theReferences[aFeature] = aSelRefFeatures;\r
+      }\r
+      else { // filter references to skip composition features of the current feature\r
+        std::set<FeaturePtr> aFilteredFeatures;\r
+        std::set<FeaturePtr>::const_iterator anIt = aSelRefFeatures.begin(),\r
+                                             aLast = aSelRefFeatures.end();\r
+        for (; anIt != aLast; anIt++) {\r
+          FeaturePtr aCFeature = *anIt;\r
+          CompositeFeaturePtr aComposite = std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aCFeature);\r
+          if (aComposite.get() && aComposite->isSub(aFeature))\r
+            continue; /// composite of the current feature should be skipped\r
+          aFilteredFeatures.insert(aCFeature);\r
+        }\r
+        theReferences[aFeature] = aFilteredFeatures;\r
+      }\r
+      if (theUseRecursion) {\r
+#ifdef DEBUG_CYCLING_1550\r
+        findReferences(aSelRefFeatures, theReferences, theUseComposite, theUseRecursion, theRecLevel);\r
+#else\r
+        findReferences(theReferences[aFeature], theReferences, theUseComposite, theUseRecursion,\r
+                       theRecLevel);\r
+#endif\r
+      }\r
+    }\r
+  }\r
+}\r
+\r
+void findAllReferences(const std::set<FeaturePtr>& theFeatures,\r
+                       std::map<FeaturePtr, std::set<FeaturePtr> >& theReferences,\r
+                       const bool theUseComposite,\r
+                       const bool theUseRecursion)\r
+{\r
+  // For dependencies, find main_list:\r
+  // sk_1(ext_1, vertex_1)\r
+  // ext_1(bool_1, sk_3)\r
+  // vertex_1()\r
+  // sk_2(ext_2)\r
+  // ext_2(bool_2)\r
+  // sk_3()\r
+  // Information: bool_1 is not selected, ext_2(bool_2) exists\r
+  // find all referenced features\r
+  std::map<FeaturePtr, std::set<FeaturePtr> > aMainList;\r
+  int aRecLevel = 0;\r
+  findReferences(theFeatures, aMainList, theUseComposite, theUseRecursion, aRecLevel);\r
+\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+  printMapInfo(aMainList, "firstDependencies");\r
+#endif\r
+  // find all dependencies for each object:\r
+  // sk_1(ext_1, vertex_1) + (sk_3, bool_1)\r
+  // ext_1(bool_1, sk_3)\r
+  // vertex_1()\r
+  // sk_2(ext_2) + (bool_1)\r
+  // ext_2(bool_1)\r
+  // sk_3()\r
+  std::map<FeaturePtr, std::set<FeaturePtr> >::const_iterator aMainIt = aMainList.begin(),\r
+                                                              aMainLast = aMainList.end();\r
+  for (; aMainIt != aMainLast; aMainIt++) {\r
+    FeaturePtr aMainListFeature = aMainIt->first;\r
+\r
+    //std::string aName = aMainListFeature->name();\r
+    std::set<FeaturePtr> aMainRefList = aMainIt->second;\r
+\r
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE\r
+    char aBuf[50];\r
+    int n = sprintf(aBuf, "%d", aMainRefList.size());\r
+    std::string aSize(aBuf);\r
+    std::cout << "_findAllReferences for the Feature: " << getFeatureInfo(aMainListFeature)\r
+              << ", references size = " << aSize << std::endl;\r
+#endif\r
+    std::set<FeaturePtr>::const_iterator anIt = aMainRefList.begin(),\r
+                                         aLast = aMainRefList.end();\r
+    std::set<FeaturePtr> aResultRefList;\r
+    aResultRefList.insert(aMainRefList.begin(), aMainRefList.end());\r
+    for (; anIt != aLast; anIt++) {\r
+      FeaturePtr aFeature = *anIt;\r
+      int aRecLevel = 0;\r
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE\r
+      std::cout << " Ref: " << getFeatureInfo(aFeature) << std::endl;\r
+#endif\r
+      aRecLevel++;\r
+      addRefsToFeature(aFeature, aMainList, theReferences, aRecLevel, aResultRefList/*aMainRefList*/);\r
+    }\r
+    theReferences[aMainListFeature] = aResultRefList;\r
+  }\r
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE\r
+    std::cout << std::endl;\r
+#endif\r
+\r
+#ifdef DEBUG_REMOVE_FEATURES\r
+  printMapInfo(theReferences, "allDependencies");\r
+#endif\r
+}\r
+\r
+void findRefsToFeatures(const std::set<FeaturePtr>& theFeatures,\r
+                        const std::map<FeaturePtr, std::set<FeaturePtr> >& theReferences,\r
+                        std::set<FeaturePtr>& theFeaturesRefsTo)\r
+{\r
+  std::set<FeaturePtr>::const_iterator anIt = theFeatures.begin(),\r
+                                       aLast = theFeatures.end();\r
+  for (; anIt != aLast; anIt++) {\r
+    FeaturePtr aFeature = *anIt;\r
+    if (theReferences.find(aFeature) == theReferences.end())\r
+      continue;\r
+    std::set<FeaturePtr> aRefList = theReferences.at(aFeature);\r
+    std::set<FeaturePtr>::const_iterator aRefIt = aRefList.begin(), aRefLast = aRefList.end();\r
+    for (; aRefIt != aRefLast; aRefIt++) {\r
+      FeaturePtr aRefFeature = *aRefIt;\r
+      CompositeFeaturePtr aComposite = std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aRefFeature);\r
+      if (aComposite.get() && aComposite->isSub(aFeature))\r
+        continue; /// composite of the current feature should not be removed\r
+\r
+      if (theFeatures.find(aRefFeature) == theFeatures.end() && // it is not selected\r
+          theFeaturesRefsTo.find(aRefFeature) == theFeaturesRefsTo.end()) // it is not added\r
+        theFeaturesRefsTo.insert(aRefFeature);\r
+    }\r
+  }\r
+}\r
+\r
+void getConcealedResults(const FeaturePtr& theFeature,\r
+                         std::list<std::shared_ptr<ModelAPI_Result> >& theResults)\r
+{\r
+  SessionPtr aSession = ModelAPI_Session::get();\r
+\r
+  std::list<std::pair<std::string, std::list<std::shared_ptr<ModelAPI_Object> > > > aRefs;\r
+  theFeature->data()->referencesToObjects(aRefs);\r
+  std::list<std::pair<std::string, std::list<ObjectPtr> > >::const_iterator\r
+                                                  anIt = aRefs.begin(), aLast = aRefs.end();\r
+  std::set<ResultPtr> alreadyThere; // to avoid duplications\r
+  for (; anIt != aLast; anIt++) {\r
+    if (!aSession->validators()->isConcealed(theFeature->getKind(), anIt->first))\r
+      continue; // use only concealed attributes\r
+    std::list<ObjectPtr> anObjects = (*anIt).second;\r
+    std::list<ObjectPtr>::const_iterator anOIt = anObjects.begin(), anOLast = anObjects.end();\r
+    for (; anOIt != anOLast; anOIt++) {\r
+      ResultPtr aResult = std::dynamic_pointer_cast<ModelAPI_Result>(*anOIt);\r
+      if (aResult && aResult->isConcealed()) {\r
+        if (alreadyThere.find(aResult) == alreadyThere.end()) // issue 1712, avoid duplicates\r
+          alreadyThere.insert(aResult);\r
+        else continue;\r
+        theResults.push_back(aResult);\r
+      }\r
+    }\r
+  }\r
+}\r
+\r
+} // namespace ModelAPI_Tools\r