Salome HOME
CEA : Lot 2 - Auto Color for Groups
[modules/shaper.git] / src / ModelAPI / ModelAPI_Tools.cpp
old mode 100755 (executable)
new mode 100644 (file)
index 36fd94b..0dd1c1c
-// 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
-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
+// Copyright (C) 2014-2021  CEA/DEN, EDF R&D
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
+
+#include <ModelAPI_AttributeBoolean.h>
+#include <ModelAPI_AttributeDocRef.h>
+#include <ModelAPI_AttributeDouble.h>
+#include <ModelAPI_AttributeIntArray.h>
+#include <ModelAPI_AttributeSelectionList.h>
+#include <ModelAPI_CompositeFeature.h>
+#include <ModelAPI_Document.h>
+#include <ModelAPI_Events.h>
+#include <ModelAPI_Object.h>
+#include <ModelAPI_ResultBody.h>
+#include <ModelAPI_ResultConstruction.h>
+#include <ModelAPI_ResultGroup.h>
+#include <ModelAPI_ResultParameter.h>
+#include <ModelAPI_ResultPart.h>
+#include <ModelAPI_Session.h>
+#include "ModelAPI_Tools.h"
+#include <ModelAPI_Validator.h>
+
+#include <Config_Translator.h>
+#include <Events_Loop.h>
+#include <Locale_Convert.h>
+
+#include <GeomAPI_ShapeHierarchy.h>
+#include <GeomAPI_ShapeIterator.h>
+
+#include <algorithm>
+#include <iostream>
+#include <list>
+#include <map>
+#include <sstream>
+
+#define RECURSE_TOP_LEVEL 50
+
+//#define DEBUG_REMOVE_FEATURES
+//#define DEBUG_REMOVE_FEATURES_RECURSE
+//#define DEBUG_CYCLING_1550
+
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE
+#include <sstream>
+std::string getFeatureInfo(FeaturePtr theFeature)
+{
+  if (!theFeature.get())
+    return "";
+  //std::ostringstream aPtrStr;
+  //aPtrStr << "[" << theFeature.get() << "] ";
+  std::string aFeatureInfo = /*aPtrStr.str() + */theFeature->name();
+  CompositeFeaturePtr aComposite = ModelAPI_Tools::compositeOwner(theFeature);
+  if (aComposite.get()) {
+      aFeatureInfo = aFeatureInfo + "[in " + aComposite->name() + "]";
+  }
+  return aFeatureInfo;
+}
+#endif
+
+#ifdef DEBUG_REMOVE_FEATURES
+void printMapInfo(const std::map<FeaturePtr, std::set<FeaturePtr> >& theMainList,
+                  const std::string& thePrefix)
+{
+  std::map<FeaturePtr, std::set<FeaturePtr> >::const_iterator aMainIt = theMainList.begin(),
+                                                              aMainLast = theMainList.end();
+  std::string anInfo;
+  for (; aMainIt != aMainLast; aMainIt++) {
+    FeaturePtr aMainListFeature = aMainIt->first;
+    std::set<FeaturePtr> aMainRefList = aMainIt->second;
+    std::set<FeaturePtr>::const_iterator anIt = aMainRefList.begin(), aLast = aMainRefList.end();
+    std::string aRefsInfo;
+    for (; anIt != aLast; anIt++) {
+      aRefsInfo += (*anIt)->name().c_str();
+      if (anIt != aLast)
+        aRefsInfo += ", ";
+    }
+    if (!aRefsInfo.empty()) {
+      anInfo = anInfo + aMainListFeature->name().c_str() + ": " + aRefsInfo + "\n";
+    }
+  }
+  std::cout << thePrefix.c_str() << " [feature: references to]: \n" << anInfo.c_str() << std::endl;
+}
+
+void printListInfo(const std::set<FeaturePtr>& theMainList,
+                  const std::string& thePrefix)
+{
+  std::set<FeaturePtr>::const_iterator aMainIt = theMainList.begin(),
+                                       aMainLast = theMainList.end();
+  std::string anInfo;
+  for (; aMainIt != aMainLast; aMainIt++) {
+    FeaturePtr aRefFeature = *aMainIt;
+    anInfo += aRefFeature->name().c_str();
+    if (aMainIt != aMainLast)
+      anInfo += ", ";
+  }
+  std::cout << thePrefix.c_str() << ": " << anInfo.c_str() << std::endl;
+}
+#endif
+
+namespace ModelAPI_Tools {
+
+std::shared_ptr<GeomAPI_Shape> shape(const ResultPtr& theResult)
+{
+  return theResult->shape();
+}
+
+// LCOV_EXCL_START
+const char* toString(ModelAPI_ExecState theExecState)
+{
+  switch (theExecState) {
+  case ModelAPI_StateDone: return "Done";
+  case ModelAPI_StateMustBeUpdated: return "Must be updated";
+  case ModelAPI_StateExecFailed: return "Execution failed";
+  case ModelAPI_StateInvalidArgument: return "Invalid argument";
+  case ModelAPI_StateNothing: return "Empty state";
+  default: return "Unknown ExecState.";
+  }
+}
+
+std::string getFeatureError(const FeaturePtr& theFeature)
+{
+  std::string anError;
+  if (!theFeature.get() || !theFeature->data()->isValid() || theFeature->isAction())
+    return anError;
+
+  // to be removed later, this error should be got from the feature
+  if (theFeature->data()->execState() == ModelAPI_StateDone ||
+      theFeature->data()->execState() == ModelAPI_StateMustBeUpdated)
+    return anError;
+
+  // set error indication
+  anError = theFeature->error();
+  if (anError.empty()) {
+    bool isDone = ( theFeature->data()->execState() == ModelAPI_StateDone
+                 || theFeature->data()->execState() == ModelAPI_StateMustBeUpdated );
+    if (!isDone) {
+      anError = toString(theFeature->data()->execState());
+      // If the feature is Composite and error is StateInvalidArgument,
+      // error text should include error of first invalid sub-feature. Otherwise
+      // it is not clear what is the reason of the invalid argument.
+      if (theFeature->data()->execState() == ModelAPI_StateInvalidArgument) {
+        CompositeFeaturePtr aComposite =
+                    std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(theFeature);
+        if (aComposite) {
+          bool aHasSubError = false;
+          for (int i = 0, aSize = aComposite->numberOfSubs(); i < aSize; i++) {
+            FeaturePtr aSubFeature = aComposite->subFeature(i);
+            std::string aSubFeatureError = getFeatureError(aSubFeature);
+            if (!aSubFeatureError.empty()) {
+              anError = anError + " in " + aSubFeature->getKind() + ".\n" + aSubFeatureError;
+              aHasSubError = true;
+              break;
+            }
+          }
+          if (!aHasSubError) { // #24260: error not in the sub-features, but in the argument
+            if (aComposite->getKind() == "Sketch" &&
+                aComposite->selection("External")->isInvalid()) {
+              std::string aMsg = "The sketch base plane is invalid, ";
+              aMsg += "please push 'Change sketch plane' button to reselect it.";
+              anError = Config_Translator::translate(aComposite->getKind(), aMsg);
+            }
+          }
+        }
+      }
+    }
+  }
+
+  return anError;
+}
+// LCOV_EXCL_STOP
+
+ObjectPtr objectByName(const DocumentPtr& theDocument, const std::string& theGroup,
+                       const std::wstring& theName)
+{
+  for (int anIndex = 0; anIndex < theDocument->size(theGroup); ++anIndex) {
+    ObjectPtr anObject = theDocument->object(theGroup, anIndex);
+    if (anObject->data()->name() == theName)
+      return anObject;
+  }
+  // not found
+  return ObjectPtr();
+}
+
+bool findVariable(const DocumentPtr& theDocument, FeaturePtr theSearcher,
+                  const std::wstring& theName, double& outValue, ResultParameterPtr& theParam)
+{
+  ObjectPtr aParamObj = objectByName(theDocument, ModelAPI_ResultParameter::group(), theName);
+  theParam = std::dynamic_pointer_cast<ModelAPI_ResultParameter>(aParamObj);
+  if (!theParam.get())
+    return false;
+  // avoid usage of parameters created later than the initial parameter
+
+  if (theSearcher.get()) {
+    FeaturePtr aParamFeat = theDocument->feature(theParam);
+    if (aParamFeat == theSearcher || theDocument->isLater(aParamFeat, theSearcher))
+      return false;
+  }
+  AttributeDoublePtr aValueAttribute = theParam->data()->real(ModelAPI_ResultParameter::VALUE());
+  outValue = aValueAttribute->value();
+  return true;
+}
+
+bool findVariable(FeaturePtr theSearcher, const std::wstring& theName, double& outValue,
+                  ResultParameterPtr& theParam, const DocumentPtr& theDocument)
+{
+  SessionPtr aSession = ModelAPI_Session::get();
+  std::list<DocumentPtr> aDocList;
+  DocumentPtr aDocument = theDocument.get() ? theDocument : aSession->activeDocument();
+  if (findVariable(aDocument, theSearcher, theName, outValue, theParam))
+    return true;
+  DocumentPtr aRootDocument = aSession->moduleDocument();
+  if (aDocument != aRootDocument) {
+    // any parameters in PartSet is okindependently on the Part position (issu #1504)
+    if (findVariable(aRootDocument, FeaturePtr(), theName, outValue, theParam))
+      return true;
+  }
+  return false;
+}
+
+ResultPtr findPartResult(const DocumentPtr& theMain, const DocumentPtr& theSub)
+{
+  // to optimize and avoid of crash on partset document close
+  // (don't touch the sub-document structure)
+  if (theMain != theSub) {
+    for (int a = theMain->size(ModelAPI_ResultPart::group()) - 1; a >= 0; a--) {
+      ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(
+          theMain->object(ModelAPI_ResultPart::group(), a));
+      if (aPart && aPart->data()->document(ModelAPI_ResultPart::DOC_REF())->value() == theSub) {
+        return aPart;
+      }
+    }
+  }
+  return ResultPtr();
+}
+
+FeaturePtr findPartFeature(const DocumentPtr& theMain, const DocumentPtr& theSub)
+{
+  // to optimize and avoid of crash on partset document close
+  // (don't touch the sub-document structure)
+  if (theMain != theSub) {
+    // iteration from top to bottom to avoid finding the movement documents before the original
+    int aSize = theMain->size(ModelAPI_Feature::group());
+    for (int a = 0; a < aSize; a++) {
+      FeaturePtr aPartFeat = std::dynamic_pointer_cast<ModelAPI_Feature>(
+          theMain->object(ModelAPI_Feature::group(), a));
+      if (aPartFeat.get()) {
+        const std::list<std::shared_ptr<ModelAPI_Result> >& aResList = aPartFeat->results();
+        std::list<std::shared_ptr<ModelAPI_Result> >::const_iterator aRes = aResList.begin();
+        for(; aRes != aResList.end(); aRes++) {
+          ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(*aRes);
+          if (aPart.get()) {
+            if (aPart->isActivated() && aPart->partDoc() == theSub)
+              return aPartFeat;
+          } else break; // if the first is not Part, others are also not
+        }
+      }
+    }
+  }
+  return FeaturePtr();
+}
+
+CompositeFeaturePtr compositeOwner(const FeaturePtr& theFeature)
+{
+  if (theFeature.get() && theFeature->data() && theFeature->data()->isValid()) {
+    const std::set<std::shared_ptr<ModelAPI_Attribute> >& aRefs = theFeature->data()->refsToMe();
+    std::set<std::shared_ptr<ModelAPI_Attribute> >::const_iterator aRefIter = aRefs.begin();
+    for(; aRefIter != aRefs.end(); aRefIter++) {
+      CompositeFeaturePtr aComp = std::dynamic_pointer_cast<ModelAPI_CompositeFeature>
+        ((*aRefIter)->owner());
+      if (aComp.get() && aComp->data()->isValid() && aComp->isSub(theFeature))
+        return aComp;
+    }
+  }
+  return CompositeFeaturePtr(); // not found
+}
+
+ResultBodyPtr bodyOwner(const ResultPtr& theSub, const bool theRoot)
+{
+  if (theSub.get()) {
+    ObjectPtr aParent = theSub->document()->parent(theSub);
+    if (aParent.get()) {
+      if (theRoot) { // try to find parent of parent
+        ResultPtr aResultParent = std::dynamic_pointer_cast<ModelAPI_Result>(aParent);
+        ResultBodyPtr aGrandParent = bodyOwner(aResultParent, true);
+        if (aGrandParent.get())
+          aParent = aGrandParent;
+      }
+      return std::dynamic_pointer_cast<ModelAPI_ResultBody>(aParent);
+    }
+  }
+  return ResultBodyPtr(); // not found
+}
+
+int bodyIndex(const ResultPtr& theSub)
+{
+  int anIndex = -1;
+  ResultBodyPtr aParent = bodyOwner(theSub);
+  if (aParent.get()) {
+    ResultBodyPtr aBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(theSub);
+    if (aBody.get() && aParent->isSub(aBody, anIndex))
+      return anIndex;
+  }
+  return anIndex; // not found
+}
+
+bool hasSubResults(const ResultPtr& theResult)
+{
+  ResultBodyPtr aCompSolid = std::dynamic_pointer_cast<ModelAPI_ResultBody>(theResult);
+  return aCompSolid.get() && aCompSolid->numberOfSubs() > 0;
+}
+
+void allSubs(const ResultBodyPtr& theResult, std::list<ResultPtr>& theResults,
+             const bool theLowerOnly) {
+  // iterate sub-bodies of compsolid
+  ResultBodyPtr aComp = std::dynamic_pointer_cast<ModelAPI_ResultBody>(theResult);
+  if (aComp.get()) {
+    int aNumSub = aComp->numberOfSubs();
+    for (int a = 0; a < aNumSub; a++) {
+      ResultBodyPtr aSub = aComp->subResult(a);
+      if (!theLowerOnly || aSub->numberOfSubs() == 0)
+        theResults.push_back(aSub);
+      allSubs(aSub, theResults);
+    }
+  }
+}
+
+void allResults(const FeaturePtr& theFeature, std::list<ResultPtr>& theResults)
+{
+  if (!theFeature.get()) // safety: for empty feature no results
+    return;
+  const std::list<std::shared_ptr<ModelAPI_Result> >& aResults = theFeature->results();
+  std::list<ResultPtr>::const_iterator aRIter = aResults.begin();
+  for (; aRIter != aResults.cend(); aRIter++) {
+    theResults.push_back(*aRIter);
+    ResultBodyPtr aResult = std::dynamic_pointer_cast<ModelAPI_ResultBody>(*aRIter);
+    allSubs(aResult, theResults);
+  }
+}
+
+//******************************************************************
+bool allDocumentsActivated(std::wstring& theNotActivatedNames)
+{
+  theNotActivatedNames = L"";
+  bool anAllPartActivated = true;
+
+  DocumentPtr aRootDoc = ModelAPI_Session::get()->moduleDocument();
+  int aSize = aRootDoc->size(ModelAPI_ResultPart::group());
+  for (int i = 0; i < aSize; i++) {
+    ObjectPtr aObject = aRootDoc->object(ModelAPI_ResultPart::group(), i);
+    ResultPartPtr aPart = std::dynamic_pointer_cast<ModelAPI_ResultPart>(aObject);
+    if (!aPart->isActivated()) {
+      anAllPartActivated = false;
+      if (!theNotActivatedNames.empty())
+        theNotActivatedNames += L", ";
+      theNotActivatedNames += aObject->data()->name().c_str();
+    }
+  }
+  return anAllPartActivated;
+}
+
+bool removeFeaturesAndReferences(const std::set<FeaturePtr>& theFeatures,
+                                 const bool /*theFlushRedisplay*/,
+                                 const bool theUseComposite,
+                                 const bool theUseRecursion)
+{
+#ifdef DEBUG_REMOVE_FEATURES
+  printListInfo(theFeatures, "selection: ");
+#endif
+
+  std::map<FeaturePtr, std::set<FeaturePtr> > aReferences;
+  ModelAPI_Tools::findAllReferences(theFeatures, aReferences, theUseComposite, theUseRecursion);
+#ifdef DEBUG_REMOVE_FEATURES
+  printMapInfo(aReferences, "allDependencies: ");
+#endif
+
+  std::set<FeaturePtr> aFeaturesRefsTo;
+  ModelAPI_Tools::findRefsToFeatures(theFeatures, aReferences, aFeaturesRefsTo);
+#ifdef DEBUG_REMOVE_FEATURES
+  printListInfo(aFeaturesRefsTo, "references: ");
+#endif
+
+  std::set<FeaturePtr> aFeatures = theFeatures;
+  if (!aFeaturesRefsTo.empty())
+    aFeatures.insert(aFeaturesRefsTo.begin(), aFeaturesRefsTo.end());
+#ifdef DEBUG_REMOVE_FEATURES
+  printListInfo(aFeatures, "removeFeatures: ");
+#endif
+
+  return ModelAPI_Tools::removeFeatures(aFeatures, false);
+}
+
+//***********************************************************************
+bool removeFeatures(const std::set<FeaturePtr>& theFeatures,
+                    const bool theFlushRedisplay)
+{
+  bool isDone = false;
+  std::set<FeaturePtr>::const_iterator anIt = theFeatures.begin(),
+                                       aLast = theFeatures.end();
+  for (; anIt != aLast; anIt++) {
+    FeaturePtr aFeature = *anIt;
+    if (aFeature.get()) {
+      DocumentPtr aDoc = aFeature->document();
+      // flush REDISPLAY signal after remove feature
+      aDoc->removeFeature(aFeature);
+      isDone = true;
+    }
+  }
+  if (isDone && theFlushRedisplay) {
+    // the redisplay signal should be flushed in order to erase
+    // the feature presentation in the viewer
+    // if should be done after removeFeature() of document
+    Events_Loop::loop()->flush(Events_Loop::loop()->eventByName(EVENT_OBJECT_TO_REDISPLAY));
+  }
+  return true;
+}
+
+//***********************************************************************
+// Fills the references list by all references of the feature from the references map.
+// This is a recusive method to find references by next found feature in the map of references.
+// \param theFeature a feature to find references
+// \param theReferencesMap a map of references
+// \param theReferences an out container of references
+void addRefsToFeature(const FeaturePtr& theFeature,
+                      const std::map<FeaturePtr, std::set<FeaturePtr> >& theReferencesMap,
+                      int theRecLevel,
+                      std::set<FeaturePtr>& theReferences)
+{
+  if (theRecLevel > RECURSE_TOP_LEVEL)
+    return;
+  theRecLevel++;
+
+  if (theReferencesMap.find(theFeature) == theReferencesMap.end())
+    return; // this feature is not in the selection list, so exists without references to it
+  std::set<FeaturePtr> aMainReferences = theReferencesMap.at(theFeature);
+
+  std::set<FeaturePtr>::const_iterator anIt = aMainReferences.begin(),
+                                       aLast = aMainReferences.end();
+  for (; anIt != aLast; anIt++) {
+    FeaturePtr aRefFeature = *anIt;
+    if (theReferences.find(aRefFeature) == theReferences.end()) {
+      addRefsToFeature(aRefFeature, theReferencesMap, theRecLevel, theReferences);
+      theReferences.insert(aRefFeature);
+    }
+  }
+}
+
+// For each feature from the feature list it searches references to the feature and append them
+// to the references map. This is a recusive method.
+// \param theFeature a feature to find references
+// \param theReferencesMap a map of references
+// \param theReferences an out container of references
+void findReferences(const std::set<FeaturePtr>& theFeatures,
+                    std::map<FeaturePtr, std::set<FeaturePtr> >& theReferences,
+                    const bool theUseComposite, const bool theUseRecursion, int theRecLevel)
+{
+  if (theRecLevel > RECURSE_TOP_LEVEL)
+    return;
+  theRecLevel++;
+  std::set<FeaturePtr>::const_iterator anIt = theFeatures.begin(),
+                                        aLast = theFeatures.end();
+  for (; anIt != aLast; anIt++) {
+    FeaturePtr aFeature = *anIt;
+    if (aFeature.get() && theReferences.find(aFeature) == theReferences.end()) {
+      DocumentPtr aSelFeatureDoc = aFeature->document();
+      std::set<FeaturePtr> aSelRefFeatures;
+      aSelFeatureDoc->refsToFeature(aFeature, aSelRefFeatures, false/*do not emit signals*/);
+      if (theUseComposite) { // do not filter selection
+        theReferences[aFeature] = aSelRefFeatures;
+      }
+      else { // filter references to skip composition features of the current feature
+        std::set<FeaturePtr> aFilteredFeatures;
+        std::set<FeaturePtr>::const_iterator aRefIt = aSelRefFeatures.begin(),
+                                             aRefLast = aSelRefFeatures.end();
+        for (; aRefIt != aRefLast; aRefIt++) {
+          FeaturePtr aCFeature = *aRefIt;
+          CompositeFeaturePtr aComposite =
+            std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aCFeature);
+          if (aComposite.get() && aComposite->isSub(aFeature))
+            continue; /// composite of the current feature should be skipped
+          aFilteredFeatures.insert(aCFeature);
+        }
+        theReferences[aFeature] = aFilteredFeatures;
+      }
+      if (theUseRecursion) {
+#ifdef DEBUG_CYCLING_1550
+        findReferences(aSelRefFeatures, theReferences, theUseComposite,
+                       theUseRecursion, theRecLevel);
+#else
+        findReferences(theReferences[aFeature], theReferences, theUseComposite, theUseRecursion,
+                       theRecLevel);
+#endif
+      }
+    }
+  }
+}
+
+void findAllReferences(const std::set<FeaturePtr>& theFeatures,
+                       std::map<FeaturePtr, std::set<FeaturePtr> >& theReferences,
+                       const bool theUseComposite,
+                       const bool theUseRecursion)
+{
+  // For dependencies, find main_list:
+  // sk_1(ext_1, vertex_1)
+  // ext_1(bool_1, sk_3)
+  // vertex_1()
+  // sk_2(ext_2)
+  // ext_2(bool_2)
+  // sk_3()
+  // Information: bool_1 is not selected, ext_2(bool_2) exists
+  // find all referenced features
+  std::map<FeaturePtr, std::set<FeaturePtr> > aMainList;
+  int aRecLevel = 0;
+  findReferences(theFeatures, aMainList, theUseComposite, theUseRecursion, aRecLevel);
+
+#ifdef DEBUG_REMOVE_FEATURES
+  printMapInfo(aMainList, "firstDependencies");
+#endif
+  // find all dependencies for each object:
+  // sk_1(ext_1, vertex_1) + (sk_3, bool_1)
+  // ext_1(bool_1, sk_3)
+  // vertex_1()
+  // sk_2(ext_2) + (bool_1)
+  // ext_2(bool_1)
+  // sk_3()
+  std::map<FeaturePtr, std::set<FeaturePtr> >::const_iterator aMainIt = aMainList.begin(),
+                                                              aMainLast = aMainList.end();
+  for (; aMainIt != aMainLast; aMainIt++) {
+    FeaturePtr aMainListFeature = aMainIt->first;
+
+    //std::string aName = aMainListFeature->name();
+    std::set<FeaturePtr> aMainRefList = aMainIt->second;
+
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE
+    char aBuf[50];
+    int n = sprintf(aBuf, "%d", aMainRefList.size());
+    std::string aSize(aBuf);
+    std::cout << "_findAllReferences for the Feature: " << getFeatureInfo(aMainListFeature)
+              << ", references size = " << aSize << std::endl;
+#endif
+    std::set<FeaturePtr>::const_iterator anIt = aMainRefList.begin(),
+                                         aLast = aMainRefList.end();
+    std::set<FeaturePtr> aResultRefList;
+    aResultRefList.insert(aMainRefList.begin(), aMainRefList.end());
+    for (; anIt != aLast; anIt++) {
+      FeaturePtr aFeature = *anIt;
+      aRecLevel = 0;
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE
+      std::cout << " Ref: " << getFeatureInfo(aFeature) << std::endl;
+#endif
+      aRecLevel++;
+      addRefsToFeature(aFeature, aMainList,
+                       aRecLevel, aResultRefList/*aMainRefList*/);
+    }
+    theReferences[aMainListFeature] = aResultRefList;
+  }
+#ifdef DEBUG_REMOVE_FEATURES_RECURSE
+    std::cout << std::endl;
+#endif
+
+#ifdef DEBUG_REMOVE_FEATURES
+  printMapInfo(theReferences, "allDependencies");
+#endif
+}
+
+void findRefsToFeatures(const std::set<FeaturePtr>& theFeatures,
+                        const std::map<FeaturePtr, std::set<FeaturePtr> >& theReferences,
+                        std::set<FeaturePtr>& theFeaturesRefsTo)
+{
+  std::set<FeaturePtr>::const_iterator anIt = theFeatures.begin(),
+                                       aLast = theFeatures.end();
+  for (; anIt != aLast; anIt++) {
+    FeaturePtr aFeature = *anIt;
+    if (theReferences.find(aFeature) == theReferences.end())
+      continue;
+    std::set<FeaturePtr> aRefList = theReferences.at(aFeature);
+    std::set<FeaturePtr>::const_iterator aRefIt = aRefList.begin(), aRefLast = aRefList.end();
+    for (; aRefIt != aRefLast; aRefIt++) {
+      FeaturePtr aRefFeature = *aRefIt;
+      CompositeFeaturePtr aComposite =
+        std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aRefFeature);
+      if (aComposite.get() && aComposite->isSub(aFeature))
+        continue; /// composite of the current feature should not be removed
+
+      if (theFeatures.find(aRefFeature) == theFeatures.end() && // it is not selected
+          theFeaturesRefsTo.find(aRefFeature) == theFeaturesRefsTo.end()) // it is not added
+        theFeaturesRefsTo.insert(aRefFeature);
+    }
+  }
+}
+
+void getConcealedResults(const FeaturePtr& theFeature,
+                         std::list<std::shared_ptr<ModelAPI_Result> >& theResults)
+{
+  SessionPtr aSession = ModelAPI_Session::get();
+
+  std::list<std::pair<std::string, std::list<std::shared_ptr<ModelAPI_Object> > > > aRefs;
+  theFeature->data()->referencesToObjects(aRefs);
+  std::list<std::pair<std::string, std::list<ObjectPtr> > >::const_iterator
+                                                  anIt = aRefs.begin(), aLast = aRefs.end();
+  std::set<ResultPtr> alreadyThere; // to avoid duplications
+  for (; anIt != aLast; anIt++) {
+    if (!aSession->validators()->isConcealed(theFeature->getKind(), anIt->first))
+      continue; // use only concealed attributes
+    std::list<ObjectPtr> anObjects = (*anIt).second;
+    std::list<ObjectPtr>::const_iterator anOIt = anObjects.begin(), anOLast = anObjects.end();
+    for (; anOIt != anOLast; anOIt++) {
+      ResultPtr aResult = std::dynamic_pointer_cast<ModelAPI_Result>(*anOIt);
+      if (aResult && aResult->isConcealed()) {
+        if (alreadyThere.find(aResult) == alreadyThere.end()) // issue 1712, avoid duplicates
+          alreadyThere.insert(aResult);
+        else continue;
+        theResults.push_back(aResult);
+      }
+    }
+  }
+}
+
+std::pair<std::wstring, bool> getDefaultName(const std::shared_ptr<ModelAPI_Result>& theResult,
+                                             const bool theInherited,
+                                             const bool theRecursive)
+{
+  typedef std::list< std::pair < std::string, std::list<ObjectPtr> > > ListOfReferences;
+
+  SessionPtr aSession = ModelAPI_Session::get();
+
+  ResultBodyPtr anOwnerRes = bodyOwner(theResult);
+  if (anOwnerRes) {
+    // names of sub-solids in CompSolid should be default (for example,
+    // result of boolean operation 'Boolean_1_1' is a CompSolid which is renamed to 'MyBOOL',
+    // however, sub-elements of 'MyBOOL' should be named 'Boolean_1_1_1', 'Boolean_1_1_2' etc.)
+    std::wostringstream aDefaultName;
+    aDefaultName << getDefaultName(anOwnerRes).first;
+    aDefaultName << "_" << (bodyIndex(theResult) + 1);
+    return std::pair<std::wstring, bool>(aDefaultName.str(), false);
+  }
+
+  FeaturePtr anOwner = ModelAPI_Feature::feature(theResult->data()->owner());
+  DataPtr aData = anOwner->data();
+
+  ListOfReferences aReferences;
+  // find first result with user-defined name
+  ListOfReferences::const_iterator aFoundRef = aReferences.end();
+  if (theInherited) {
+    aData->referencesToObjects(aReferences);
+
+    for (ListOfReferences::const_iterator aRefIt = aReferences.begin();
+         aRefIt != aReferences.end(); ++aRefIt) {
+      bool isConcealed = aSession->validators()->isConcealed(anOwner->getKind(), aRefIt->first);
+      bool isMainArg = isConcealed &&
+                       aSession->validators()->isMainArgument(anOwner->getKind(), aRefIt->first);
+      if (isConcealed) {
+        // check the referred object is a Body
+        // (for example, ExtrusionCut has a sketch as a first attribute which is concealing)
+        bool isBody = aRefIt->second.size() > 1 || (aRefIt->second.size() == 1 &&
+                      aRefIt->second.front().get() &&
+                      aRefIt->second.front()->groupName() == ModelAPI_ResultBody::group());
+        if (isBody && (isMainArg || aFoundRef == aReferences.end() ||
+            aData->isPrecedingAttribute(aRefIt->first, aFoundRef->first)))
+          aFoundRef = aRefIt;
+
+        if (isMainArg)
+          break;
+      }
+    }
+  }
+  // get the result number in the feature
+  int anIndexInOwner = 0;
+  const std::list<ResultPtr>& anOwnerResults = anOwner->results();
+  std::list<ResultPtr>::const_iterator aResIt = anOwnerResults.cbegin();
+  for(; aResIt != anOwnerResults.cend(); aResIt++) {
+    if(*aResIt == theResult)
+      break;
+    anIndexInOwner++;
+  }
+
+  // find an object which is concealed by theResult
+  if (aFoundRef != aReferences.end() && !aFoundRef->second.empty()) {
+    // store number of references for each object
+    std::map<ResultPtr, int> aNbRefToObject;
+    // search the object by result index
+    std::list<ObjectPtr>::const_iterator anObjIt = aFoundRef->second.begin();
+    int aResultIndex = anIndexInOwner;
+    while (--aResultIndex >= 0) {
+      ResultPtr aCurRes = std::dynamic_pointer_cast<ModelAPI_Result>(*anObjIt);
+      ResultBodyPtr aParentBody = ModelAPI_Tools::bodyOwner(aCurRes);
+      if (aParentBody)
+        aCurRes = aParentBody;
+      if (aNbRefToObject.find(aCurRes) == aNbRefToObject.end())
+        aNbRefToObject[aCurRes] = 1;
+      else
+        aNbRefToObject[aCurRes] += 1;
+
+      ++anObjIt;
+      if (anObjIt == aFoundRef->second.end()) {
+        anObjIt = aFoundRef->second.begin();
+        break;
+      }
+    }
+    // check the result is a Body
+    if (anObjIt->get() && (*anObjIt)->groupName() == ModelAPI_ResultBody::group()) {
+      // check the result is part of CompSolid
+      ResultPtr anObjRes = std::dynamic_pointer_cast<ModelAPI_Result>(*anObjIt);
+      ResultBodyPtr aParentBody = ModelAPI_Tools::bodyOwner(anObjRes);
+      if (aParentBody)
+        anObjRes = aParentBody;
+
+      // return name of reference result only if it has been renamed by the user,
+      // in other case compose a default name
+      if (anObjRes->data()->hasUserDefinedName() ||
+          (theRecursive && anObjRes->data()->name() != getDefaultName(anObjRes).first)) {
+        std::wstringstream aName;
+        aName << anObjRes->data()->name();
+        std::map<ResultPtr, int>::iterator aFound = aNbRefToObject.find(anObjRes);
+        if (aFound != aNbRefToObject.end()) {
+          // to generate unique name, add suffix if there are several results
+          // referring to the same shape
+          aName << "_" << aFound->second + 1;
+        }
+        return std::pair<std::wstring, bool>(aName.str(), true);
+      }
+    }
+  }
+
+  // compose default name by the name of the feature and the index of result
+  std::wstringstream aDefaultName;
+  aDefaultName << anOwner->name();
+  // if there are several results (issue #899: any number of result),
+  // add unique prefix starting from second
+  if (anIndexInOwner > 0 || theResult->groupName() == ModelAPI_ResultBody::group())
+    aDefaultName << "_" << anIndexInOwner + 1;
+  return std::pair<std::wstring, bool>(aDefaultName.str(), false);
+}
+
+std::set<FeaturePtr> getParents(const FeaturePtr& theFeature)
+{
+  std::set<FeaturePtr> aParents;
+  for (FeaturePtr aCurFeat = theFeature; aCurFeat; ) {
+    CompositeFeaturePtr aFoundComposite;
+    const std::set<AttributePtr>& aRefs = aCurFeat->data()->refsToMe();
+    for (std::set<AttributePtr>::const_iterator anIt = aRefs.begin();
+      anIt != aRefs.end(); ++anIt) {
+      FeaturePtr aF = ModelAPI_Feature::feature((*anIt)->owner());
+      aFoundComposite = std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(aF);
+      if (aFoundComposite && aFoundComposite->isSub(aCurFeat))
+        break;
+      else
+        aFoundComposite = CompositeFeaturePtr();
+    }
+
+    if (aFoundComposite) {
+      aParents.insert(aFoundComposite);
+      aCurFeat = aFoundComposite;
+    }
+    else {
+      // add the part containing high-level feature
+      SessionPtr aSession = ModelAPI_Session::get();
+      DocumentPtr aPartSetDoc = aSession->moduleDocument();
+      std::list<FeaturePtr> aPartSetFeatures = aPartSetDoc->allFeatures();
+      for (std::list<FeaturePtr>::const_iterator anIt = aPartSetFeatures.begin();
+        anIt != aPartSetFeatures.end(); ++anIt) {
+        aFoundComposite = std::dynamic_pointer_cast<ModelAPI_CompositeFeature>(*anIt);
+        if (aFoundComposite && aFoundComposite->isSub(aCurFeat)) {
+          aParents.insert(aFoundComposite);
+          break;
+        }
+      }
+
+      aCurFeat = FeaturePtr();
+    }
+  }
+  return aParents;
+}
+
+void fillShapeHierarchy(const GeomShapePtr& theShape,
+                        const ResultPtr& theContext,
+                        GeomAPI_ShapeHierarchy& theHierarchy)
+{
+  ResultBodyPtr aResCompSolidPtr = ModelAPI_Tools::bodyOwner(theContext);
+  if (aResCompSolidPtr.get()) {
+    std::shared_ptr<GeomAPI_Shape> aContextShape = aResCompSolidPtr->shape();
+    if (aContextShape->shapeType() <= GeomAPI_Shape::COMPSOLID) {
+      theHierarchy.addParent(theShape, aContextShape);
+      fillShapeHierarchy(aContextShape, aResCompSolidPtr, theHierarchy);
+    }
+  }
+}
+
+
+void removeResults(const std::list<ResultPtr>& theResults)
+{
+  // collect all documents where the results must be removed
+  std::map<DocumentPtr, std::list<ResultPtr> > aDocs;
+
+  std::list<ResultPtr>::const_iterator aResIter = theResults.cbegin();
+  for(; aResIter != theResults.cend(); aResIter++) {
+    DocumentPtr aDoc = (*aResIter)->document();
+    if (!aDocs.count(aDoc))
+      aDocs[aDoc] = std::list<ResultPtr>();
+    aDocs[aDoc].push_back(*aResIter);
+  }
+  // create a "remove" feature in each doc
+  std::map<DocumentPtr, std::list<ResultPtr> >::iterator aDoc = aDocs.begin();
+  for(; aDoc != aDocs.end(); aDoc++) {
+    FeaturePtr aRemove = aDoc->first->addFeature("RemoveResults");
+    if (aRemove) {
+      for(aResIter = aDoc->second.cbegin(); aResIter != aDoc->second.cend(); aResIter++)
+        aRemove->selectionList("results")->append(*aResIter, GeomShapePtr());
+    }
+  }
+}
+
+// used by GUI only
+// LCOV_EXCL_START
+
+//**************************************************************
+void setDeflection(ResultPtr theResult, const double theDeflection)
+{
+  if (!theResult.get())
+    return;
+
+  AttributeDoublePtr aDeflectionAttr = theResult->data()->real(ModelAPI_Result::DEFLECTION_ID());
+  if (aDeflectionAttr.get() != NULL) {
+    aDeflectionAttr->setValue(theDeflection);
+  }
+}
+
+double getDeflection(const std::shared_ptr<ModelAPI_Result>& theResult)
+{
+  double aDeflection = -1;
+  // get deflection from the attribute of the result
+  if (theResult.get() != NULL &&
+    theResult->data()->attribute(ModelAPI_Result::DEFLECTION_ID()).get() != NULL) {
+    AttributeDoublePtr aDoubleAttr = theResult->data()->real(ModelAPI_Result::DEFLECTION_ID());
+    if (aDoubleAttr.get() && aDoubleAttr->isInitialized()) {
+      double aValue = aDoubleAttr->value();
+      if (aValue > 0) /// zero value should not be used as a deflection(previous studies)
+        aDeflection = aDoubleAttr->value();
+    }
+  }
+  return aDeflection;
+}
+
+//******************************************************
+void setColor(ResultPtr theResult, const std::vector<int>& theColor)
+{
+  if (!theResult.get())
+    return;
+
+  AttributeIntArrayPtr aColorAttr = theResult->data()->intArray(ModelAPI_Result::COLOR_ID());
+  if (aColorAttr.get() != NULL) {
+    if (!aColorAttr->size()) {
+      aColorAttr->setSize(3);
+    }
+    aColorAttr->setValue(0, theColor[0]);
+    aColorAttr->setValue(1, theColor[1]);
+    aColorAttr->setValue(2, theColor[2]);
+  }
+}
+
+void getColor(const std::shared_ptr<ModelAPI_Result>& theResult, std::vector<int>& theColor)
+{
+  theColor.clear();
+  // get color from the attribute of the result
+  if (theResult.get() != NULL &&
+    theResult->data()->attribute(ModelAPI_Result::COLOR_ID()).get() != NULL) {
+    AttributeIntArrayPtr aColorAttr = theResult->data()->intArray(ModelAPI_Result::COLOR_ID());
+    if (aColorAttr.get() && aColorAttr->size()) {
+      theColor.push_back(aColorAttr->value(0));
+      theColor.push_back(aColorAttr->value(1));
+      theColor.push_back(aColorAttr->value(2));
+    }
+  }
+}
+
+//******************************************************
+void getIsoLines(const std::shared_ptr<ModelAPI_Result>& theResult,
+  bool& isVisible, std::vector<int>& theNbLines)
+{
+  theNbLines.clear();
+  isVisible = false;
+  if (!theResult.get())
+    return;
+  if (theResult->groupName() == ModelAPI_ResultConstruction::group()) {
+    theNbLines.push_back(0);
+    theNbLines.push_back(0);
+  }
+  else {
+    // get color from the attribute of the result
+    AttributeIntArrayPtr aAttr = theResult->data()->intArray(ModelAPI_Result::ISO_LINES_ID());
+    if (aAttr.get()) {
+      if (aAttr->size()) {
+        theNbLines.push_back(aAttr->value(0));
+        theNbLines.push_back(aAttr->value(1));
+      }
+    }
+    AttributeBooleanPtr aBoolAttr =
+      theResult->data()->boolean(ModelAPI_Result::SHOW_ISO_LINES_ID());
+    if (aBoolAttr.get()) {
+      isVisible = aBoolAttr->value();
+    }
+  }
+}
+
+//******************************************************
+void setIsoLines(ResultPtr theResult, const std::vector<int>& theIso)
+{
+  if (!theResult.get())
+    return;
+
+  AttributeIntArrayPtr aAttr = theResult->data()->intArray(ModelAPI_Result::ISO_LINES_ID());
+  if (aAttr.get() != NULL) {
+    if (!aAttr->size()) {
+      aAttr->setSize(2);
+    }
+    aAttr->setValue(0, theIso[0]);
+    aAttr->setValue(1, theIso[1]);
+  }
+}
+
+//******************************************************
+void showIsoLines(std::shared_ptr<ModelAPI_Result> theResult, bool theShow)
+{
+  if (!theResult.get())
+    return;
+
+  AttributeBooleanPtr aAttr = theResult->data()->boolean(ModelAPI_Result::SHOW_ISO_LINES_ID());
+  if (aAttr.get() != NULL) {
+    aAttr->setValue(theShow);
+  }
+}
+
+//******************************************************
+bool isShownIsoLines(std::shared_ptr<ModelAPI_Result> theResult)
+{
+  if (!theResult.get())
+    return false;
+
+  AttributeBooleanPtr aAttr = theResult->data()->boolean(ModelAPI_Result::SHOW_ISO_LINES_ID());
+  if (aAttr.get() != NULL) {
+    return aAttr->value();
+  }
+  return false;
+}
+
+//**************************************************************
+void setTransparency(ResultPtr theResult, double theTransparency)
+{
+  if (!theResult.get())
+    return;
+
+  AttributeDoublePtr anAttribute = theResult->data()->real(ModelAPI_Result::TRANSPARENCY_ID());
+  if (anAttribute.get() != NULL) {
+    anAttribute->setValue(theTransparency);
+  }
+}
+
+double getTransparency(const std::shared_ptr<ModelAPI_Result>& theResult)
+{
+  double aTransparency = -1;
+  // get transparency from the attribute of the result
+  if (theResult.get() != NULL &&
+    theResult->data()->attribute(ModelAPI_Result::TRANSPARENCY_ID()).get() != NULL) {
+    AttributeDoublePtr aDoubleAttr = theResult->data()->real(ModelAPI_Result::TRANSPARENCY_ID());
+    if (aDoubleAttr.get() && aDoubleAttr->isInitialized()) {
+      aTransparency = aDoubleAttr->value();
+    }
+  }
+  return aTransparency;
+}
+
+void copyVisualizationAttrs(
+  std::shared_ptr<ModelAPI_Result> theSource, std::shared_ptr<ModelAPI_Result> theDest)
+{
+  // color
+  AttributeIntArrayPtr aSourceColor = theSource->data()->intArray(ModelAPI_Result::COLOR_ID());
+  if (aSourceColor.get() && aSourceColor->isInitialized() && aSourceColor->size()) {
+    AttributeIntArrayPtr aDestColor = theDest->data()->intArray(ModelAPI_Result::COLOR_ID());
+    if (aDestColor.get()) {
+      aDestColor->setSize(aSourceColor->size());
+      for(int a = 0; a < aSourceColor->size(); a++)
+        aDestColor->setValue(a, aSourceColor->value(a));
+    }
+  }
+  // Iso-lines
+  AttributeIntArrayPtr aSource = theSource->data()->intArray(ModelAPI_Result::ISO_LINES_ID());
+  if (aSource.get() && aSource->isInitialized() && aSource->size()) {
+    AttributeIntArrayPtr aDest = theDest->data()->intArray(ModelAPI_Result::ISO_LINES_ID());
+    if (aDest.get()) {
+      aDest->setSize(aSource->size());
+      for(int a = 0; a < aSource->size(); a++)
+        aDest->setValue(a, aSource->value(a));
+    }
+  }
+  // deflection
+  AttributeDoublePtr aSourceDefl = theSource->data()->real(ModelAPI_Result::DEFLECTION_ID());
+  if (aSourceDefl.get() && aSourceDefl->isInitialized()) {
+    AttributeDoublePtr aDestDefl = theDest->data()->real(ModelAPI_Result::DEFLECTION_ID());
+    if (aDestDefl.get()) {
+      aDestDefl->setValue(aSourceDefl->value());
+    }
+  }
+  // transparency
+  AttributeDoublePtr aSourceTransp = theSource->data()->real(ModelAPI_Result::TRANSPARENCY_ID());
+  if (aSourceTransp.get() && aSourceTransp->isInitialized()) {
+    AttributeDoublePtr aDestTransp = theDest->data()->real(ModelAPI_Result::TRANSPARENCY_ID());
+    if (aDestTransp.get()) {
+      aDestTransp->setValue(aSourceTransp->value());
+    }
+  }
+}
+
+std::list<FeaturePtr> referencedFeatures(
+  std::shared_ptr<ModelAPI_Result> theTarget, const std::string& theFeatureKind,
+  const bool theSortResults)
+{
+  std::set<FeaturePtr> aResSet; // collect in the set initially to avoid duplicates
+  std::list<ResultPtr> allSubRes;
+  allSubRes.push_back(theTarget);
+  ResultBodyPtr aBody = std::dynamic_pointer_cast<ModelAPI_ResultBody>(theTarget);
+  if (aBody.get())
+    allSubs(aBody, allSubRes);
+  std::list<ResultPtr>::iterator aSub = allSubRes.begin();
+  for(; aSub != allSubRes.end(); aSub++) {
+    const std::set<AttributePtr>& aRefs = (*aSub)->data()->refsToMe();
+    std::set<AttributePtr>::const_iterator aRef = aRefs.cbegin();
+    for(; aRef != aRefs.cend(); aRef++) {
+      FeaturePtr aFeat = std::dynamic_pointer_cast<ModelAPI_Feature>((*aRef)->owner());
+      if (aFeat.get() && (theFeatureKind.empty() || aFeat->getKind() == theFeatureKind))
+        aResSet.insert(aFeat);
+    }
+  }
+  // add also feature of the target that may be referenced as a whole
+  FeaturePtr aTargetFeature = theTarget->document()->feature(theTarget);
+  const std::set<AttributePtr>& aRefs = aTargetFeature->data()->refsToMe();
+  std::set<AttributePtr>::const_iterator aRef = aRefs.cbegin();
+  for(; aRef != aRefs.cend(); aRef++) {
+    FeaturePtr aFeat = std::dynamic_pointer_cast<ModelAPI_Feature>((*aRef)->owner());
+    if (aFeat.get() && (theFeatureKind.empty() || aFeat->getKind() == theFeatureKind))
+      aResSet.insert(aFeat);
+  }
+  // check also Group-operations that may refer to groups - add them for theFeatureKind "Group"
+  if (theFeatureKind == "Group") {
+    std::set<FeaturePtr> aGroupOperations;
+    for(bool aNeedIterate = true; aNeedIterate; ) {
+      std::set<FeaturePtr>::iterator aResIter = aResSet.begin();
+      for(; aResIter != aResSet.end(); aResIter++) {
+        std::list<ResultPtr>::const_iterator aGroupRes = (*aResIter)->results().cbegin();
+        for(; aGroupRes != (*aResIter)->results().cend(); aGroupRes++) {
+          const std::set<AttributePtr>& aGroupRefs = (*aGroupRes)->data()->refsToMe();
+          std::set<AttributePtr>::const_iterator aRefIt = aGroupRefs.cbegin();
+          for(; aRefIt != aGroupRefs.cend(); aRefIt++) {
+            FeaturePtr aFeat = std::dynamic_pointer_cast<ModelAPI_Feature>((*aRefIt)->owner());
+            if (aFeat.get() && !aGroupOperations.count(aFeat) && !aFeat->results().empty() &&
+                aFeat->firstResult()->groupName() == ModelAPI_ResultGroup::group()) {
+              // iterate results of this group operation because it may be without theTarget shape
+              GeomShapePtr aTargetShape = theTarget->shape();
+              bool anIsIn = false;
+              std::list<ResultPtr>::const_iterator anOpRes = aFeat->results().cbegin();
+              for(; anOpRes != aFeat->results().cend() && !anIsIn; anOpRes++) {
+                GeomShapePtr anOpShape = (*anOpRes)->shape();
+                if (!anOpShape.get() || anOpShape->isNull())
+                  continue;
+                for(GeomAPI_ShapeIterator aSubIt(anOpShape); aSubIt.more(); aSubIt.next()) {
+                  if (aTargetShape->isSubShape(aSubIt.current(), false)) {
+                    anIsIn = true;
+                    break;
+                  }
+                }
+              }
+              if (anIsIn)
+                aGroupOperations.insert(aFeat);
+            }
+          }
+        }
+      }
+      // insert all new group operations into result and if they are, check for next dependencies
+      aNeedIterate = false;
+      std::set<FeaturePtr>::iterator aGroupOpIter = aGroupOperations.begin();
+      for(; aGroupOpIter != aGroupOperations.end(); aGroupOpIter++) {
+        if (aResSet.find(*aGroupOpIter) == aResSet.end()) {
+          aResSet.insert(*aGroupOpIter);
+          aNeedIterate = true;
+        }
+      }
+    }
+  }
+
+  std::list<FeaturePtr> aResList;
+  std::set<FeaturePtr>::iterator aResIter = aResSet.begin();
+  for(; aResIter != aResSet.end(); aResIter++) {
+    if (theSortResults) { // sort results by creation-order
+      std::list<FeaturePtr>::iterator aListIter = aResList.begin();
+      for(; aListIter != aResList.end(); aListIter++) {
+        if ((*aResIter)->document()->isLater(*aListIter, *aResIter))
+          break;
+      }
+      if (aListIter == aResList.end()) // goes to the end
+        aResList.push_back(*aResIter);
+      else
+        aResList.insert(aListIter, *aResIter);
+    } else //just push to the end unsorted
+      aResList.push_back(*aResIter);
+  }
+  return aResList;
+}
+
+void setValues(std::vector<int>& theRGB, const int theRed, const int theGreen, const int theBlue)
+{
+  theRGB.push_back(theRed);
+  theRGB.push_back(theGreen);
+  theRGB.push_back(theBlue);
+}
+
+std::vector<int> HSVtoRGB(int theH, int theS, int theV)
+{
+  std::vector<int> aRGB;
+  if (theH < 0 || theH > 360 ||
+      theS < 0 || theS > 100 ||
+      theV < 0 || theV > 100)
+    return aRGB;
+
+  int aHi = (int)theH/60;
+  double aV = theV;
+  double aVmin = (100 - theS)*theV/100;
+  double anA = (theV - aVmin)* (theH % 60) / 60;
+  double aVinc = aVmin + anA;
+  double aVdec = theV - anA;
+  double aPercentToValue = 255./100;
+  int aV_int    = (int)(aV*aPercentToValue);
+  int aVinc_int = (int)(aVinc*aPercentToValue);
+  int aVmin_int = (int)(aVmin*aPercentToValue);
+  int aVdec_int = (int)(aVdec*aPercentToValue);
+
+  switch(aHi) {
+    case 0: setValues(aRGB, aV_int,    aVinc_int, aVmin_int); break;
+    case 1: setValues(aRGB, aVdec_int, aV_int,    aVmin_int); break;
+    case 2: setValues(aRGB, aVmin_int, aV_int,    aVinc_int); break;
+    case 3: setValues(aRGB, aVmin_int, aVdec_int, aV_int); break;
+    case 4: setValues(aRGB, aVinc_int, aVmin_int, aV_int); break;
+    case 5: setValues(aRGB, aV_int,    aVmin_int, aVdec_int); break;
+    default: break;
+  }
+  return aRGB;
+}
+
+std::array<std::vector<int>, 10> myColorTab = {
+  std::vector<int> {255, 0, 0},
+  std::vector<int> {0, 255, 0},
+  std::vector<int> {0, 0, 255},
+  std::vector<int> {255, 255, 0},
+  std::vector<int> {0, 255, 255},
+  std::vector<int> {255, 0, 255},
+  std::vector<int> {255, 94, 0},
+  std::vector<int> {132, 255, 0},
+  std::vector<int> {132, 0, 255},
+  std::vector<int> {0, 0, 0},
+};
+
+void findRandomColor(std::vector<int>& theValues, bool theReset)
+{
+  static int i = 0;
+  static std::vector<std::vector<int>> usedGeneratedColor;
+
+  // True when disabling auto-color
+  if ( theReset ) {
+    i = 0;
+    return;
+  }
+
+  theValues.clear();
+  if (i < myColorTab.size()) {
+    theValues = myColorTab[i++];
+  } else {
+      int timeout = 0;
+      std::vector<int> aHSVColor;
+      std::vector<int> aRGBColor;
+      do {
+        aHSVColor = {rand() % 360 , rand() % (100 - 50 + 1) + 50, rand() % (100 - 50 + 1) + 50};
+        aRGBColor = HSVtoRGB(aHSVColor[0], aHSVColor[1], aHSVColor[2]);
+        timeout++;
+      } while (
+        timeout < 20 &&
+        std::find(usedGeneratedColor.begin(), usedGeneratedColor.end(), aHSVColor)
+        != usedGeneratedColor.end() &&
+        std::find(myColorTab.begin(), myColorTab.end(), aRGBColor) != myColorTab.end());
+      usedGeneratedColor.push_back(aHSVColor);
+      theValues = aRGBColor;
+  }
+}
+
+// LCOV_EXCL_STOP
+
+} // namespace ModelAPI_Tools