src/SketchSolver/SketchSolver_Group.cpp

   1 // Copyright (C) 2014-20xx CEA/DEN, EDF R&D
   2
   3 // File:    SketchSolver_Group.cpp
   4 // Created: 27 May 2014
   5 // Author:  Artem ZHIDKOV
   6
   7 #include "SketchSolver_Group.h"
   8
   9 #include <SketchSolver_Builder.h>
  10 #include <SketchSolver_Constraint.h>
  11 #include <SketchSolver_ConstraintCoincidence.h>
  12 #include <SketchSolver_Error.h>
  13
  14 #include <Events_Error.h>
  15 #include <Events_Loop.h>
  16 #include <GeomAPI_XY.h>
  17 #include <GeomAPI_Dir2d.h>
  18 #include <GeomAPI_Pnt2d.h>
  19 #include <GeomDataAPI_Dir.h>
  20 #include <GeomDataAPI_Point.h>
  21 #include <GeomDataAPI_Point2D.h>
  22 #include <ModelAPI_AttributeDouble.h>
  23 #include <ModelAPI_Document.h>
  24 #include <ModelAPI_Events.h>
  25 #include <ModelAPI_ResultConstruction.h>
  26
  27 #include <SketchPlugin_Constraint.h>
  28 #include <SketchPlugin_ConstraintEqual.h>
  29 #include <SketchPlugin_ConstraintFillet.h>
  30 #include <SketchPlugin_ConstraintLength.h>
  31 #include <SketchPlugin_ConstraintCoincidence.h>
  32 #include <SketchPlugin_ConstraintMirror.h>
  33 #include <SketchPlugin_ConstraintRigid.h>
  34 #include <SketchPlugin_ConstraintTangent.h>
  35 #include <SketchPlugin_Feature.h>
  36
  37 #include <SketchPlugin_Arc.h>
  38 #include <SketchPlugin_Circle.h>
  39 #include <SketchPlugin_Line.h>
  40 #include <SketchPlugin_Point.h>
  41 #include <SketchPlugin_Sketch.h>
  42
  43 #include <math.h>
  44 #include <assert.h>
  45
  46
  47 /// \brief This class is used to give unique index to the groups
  48 class GroupIndexer
  49 {
  50 public:
  51   /// \brief Return vacant index
  52   static Slvs_hGroup NEW_GROUP() { return ++myGroupIndex; }
  53   /// \brief Removes the index
  54   static void REMOVE_GROUP(const Slvs_hGroup& theIndex) {
  55     if (myGroupIndex == theIndex)
  56       myGroupIndex--;
  57   }
  58
  59 private:
  60   GroupIndexer() {};
  61
  62   static Slvs_hGroup myGroupIndex; ///< index of the group
  63 };
  64
  65 Slvs_hGroup GroupIndexer::myGroupIndex = 0;
  66
  67
  68
  69 // ========================================================
  70 // =========  SketchSolver_Group  ===============
  71 // ========================================================
  72
  73 SketchSolver_Group::SketchSolver_Group(
  74     std::shared_ptr<ModelAPI_CompositeFeature> theWorkplane)
  75     : myID(GroupIndexer::NEW_GROUP())
  76 {
  77   // Initialize workplane
  78   myWorkplaneID = SLVS_E_UNKNOWN;
  79 #ifndef NDEBUG
  80   assert(addWorkplane(theWorkplane));
  81 #else
  82   addWorkplane(theWorkplane);
  83 #endif
  84 }
  85
  86 SketchSolver_Group::~SketchSolver_Group()
  87 {
  88   myConstraints.clear();
  89   GroupIndexer::REMOVE_GROUP(myID);
  90 }
  91
  92 // ============================================================================
  93 //  Function: isBaseWorkplane
  94 //  Class:    SketchSolver_Group
  95 //  Purpose:  verify the group is based on the given workplane
  96 // ============================================================================
  97 bool SketchSolver_Group::isBaseWorkplane(CompositeFeaturePtr theWorkplane) const
  98 {
  99   return theWorkplane == mySketch;
 100 }
 101
 102 // ============================================================================
 103 //  Function: isInteract
 104 //  Class:    SketchSolver_Group
 105 //  Purpose:  verify are there any entities in the group used by given constraint
 106 // ============================================================================
 107 bool SketchSolver_Group::isInteract(
 108     std::shared_ptr<SketchPlugin_Feature> theFeature) const
 109 {
 110   // Empty group interacts with everything
 111   if (isEmpty()) return true;
 112   ConstraintPtr aConstraint = std::dynamic_pointer_cast<SketchPlugin_Constraint>(theFeature);
 113   if (aConstraint)
 114     return myFeatureStorage->isInteract(aConstraint);
 115   return myFeatureStorage->isInteract(std::dynamic_pointer_cast<ModelAPI_Feature>(theFeature));
 116 }
 117
 118 // ============================================================================
 119 //  Function: getFeatureId
 120 //  Class:    SketchSolver_Group
 121 //  Purpose:  Find the identifier of the feature, if it already exists in the group
 122 // ============================================================================
 123 Slvs_hEntity SketchSolver_Group::getFeatureId(FeaturePtr theFeature) const
 124 {
 125   Slvs_hEntity aResult = SLVS_E_UNKNOWN;
 126   if (!myFeatureStorage)
 127     return aResult;
 128   std::set<ConstraintPtr> aConstraints = myFeatureStorage->getConstraints(theFeature);
 129   if (aConstraints.empty())
 130     return aResult;
 131   std::set<ConstraintPtr>::iterator aConstrIter = aConstraints.begin();
 132   for (; aConstrIter != aConstraints.end(); aConstrIter++) {
 133     ConstraintConstraintMap::const_iterator aCIter = myConstraints.find(*aConstrIter);
 134     if (aCIter == myConstraints.end())
 135       continue;
 136     aResult = aCIter->second->getId(theFeature);
 137     if (aResult != SLVS_E_UNKNOWN)
 138       return aResult;
 139   }
 140   return SLVS_E_UNKNOWN;
 141 }
 142
 143 // ============================================================================
 144 //  Function: getAttributeId
 145 //  Class:    SketchSolver_Group
 146 //  Purpose:  Find the identifier of the attribute, if it already exists in the group
 147 // ============================================================================
 148 Slvs_hEntity SketchSolver_Group::getAttributeId(AttributePtr theAttribute) const
 149 {
 150   Slvs_hEntity aResult = SLVS_E_UNKNOWN;
 151   if (!myFeatureStorage)
 152     return aResult;
 153   std::set<ConstraintPtr> aConstraints = myFeatureStorage->getConstraints(theAttribute);
 154   if (aConstraints.empty())
 155     return aResult;
 156   std::set<ConstraintPtr>::iterator aConstrIter = aConstraints.begin();
 157   for (; aConstrIter != aConstraints.end(); aConstrIter++) {
 158     ConstraintConstraintMap::const_iterator aCIter = myConstraints.find(*aConstrIter);
 159     if (aCIter == myConstraints.end())
 160       continue;
 161     aResult = aCIter->second->getId(theAttribute);
 162     if (aResult != SLVS_E_UNKNOWN)
 163       return aResult;
 164   }
 165   return SLVS_E_UNKNOWN;
 166 }
 167
 168 // ============================================================================
 169 //  Function: changeConstraint
 170 //  Class:    SketchSolver_Group
 171 //  Purpose:  create/update the constraint in the group
 172 // ============================================================================
 173 bool SketchSolver_Group::changeConstraint(
 174     std::shared_ptr<SketchPlugin_Constraint> theConstraint)
 175 {
 176   // There is no workplane yet, something wrong
 177   if (myWorkplaneID == SLVS_E_UNKNOWN)
 178     return false;
 179
 180   if (!theConstraint)
 181     return false;
 182
 183   bool isNewConstraint = myConstraints.find(theConstraint) == myConstraints.end();
 184   if (isNewConstraint) {
 185     // Add constraint to the current group
 186     SolverConstraintPtr aConstraint =
 187         SketchSolver_Builder::getInstance()->createConstraint(theConstraint);
 188     if (!aConstraint)
 189       return false;
 190     aConstraint->setGroup(this);
 191     aConstraint->setStorage(myStorage);
 192     if (!aConstraint->error().empty()) {
 193       if (aConstraint->error() == SketchSolver_Error::NOT_INITIALIZED())
 194         return false; // some attribute are not initialized yet, don't show message
 195       Events_Error::send(aConstraint->error(), this);
 196     }
 197
 198     // Additional verification of coincidence of several points
 199     if (theConstraint->getKind() == SketchPlugin_ConstraintCoincidence::ID()) {
 200       ConstraintConstraintMap::iterator aCIter = myConstraints.begin();
 201       for (; aCIter != myConstraints.end(); aCIter++) {
 202         std::shared_ptr<SketchSolver_ConstraintCoincidence> aCoincidence =
 203           std::dynamic_pointer_cast<SketchSolver_ConstraintCoincidence>(aCIter->second);
 204         if (!aCoincidence)
 205           continue;
 206         std::shared_ptr<SketchSolver_ConstraintCoincidence> aCoinc2 =
 207           std::dynamic_pointer_cast<SketchSolver_ConstraintCoincidence>(aConstraint);
 208         if (aCoincidence != aCoinc2 && aCoincidence->isCoincide(aCoinc2)) {
 209           aCoincidence->attach(aCoinc2);
 210           aConstraint = aCoincidence;
 211         }
 212       }
 213     }
 214     myConstraints[theConstraint] = aConstraint;
 215   }
 216   else
 217     myConstraints[theConstraint]->update();
 218
 219   // Fix base features for fillet
 220   if (isNewConstraint && theConstraint->getKind() == SketchPlugin_ConstraintFillet::ID()) {
 221     std::list<AttributePtr> anAttrList =
 222         theConstraint->data()->attributes(ModelAPI_AttributeRefAttr::typeId());
 223     std::list<AttributePtr>::iterator anAttrIter = anAttrList.begin();
 224     for (; anAttrIter != anAttrList.end(); anAttrIter++) {
 225       AttributeRefAttrPtr aRefAttr =
 226           std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*anAttrIter);
 227       if (!aRefAttr || !aRefAttr->isObject())
 228         continue;
 229       FeaturePtr aFeature = ModelAPI_Feature::feature(aRefAttr->object());
 230       SolverConstraintPtr aConstraint =
 231           SketchSolver_Builder::getInstance()->createRigidConstraint(aFeature);
 232       if (!aConstraint)
 233         continue;
 234       aConstraint->setGroup(this);
 235       aConstraint->setStorage(myStorage);
 236       setTemporary(aConstraint);
 237     }
 238   }
 239   //// Fix base features for mirror
 240   //if (theConstraint->getKind() == SketchPlugin_ConstraintMirror::ID()) {
 241   //  AttributeRefListPtr aRefList = std::dynamic_pointer_cast<ModelAPI_AttributeRefList>(
 242   //      theConstraint->attribute(SketchPlugin_ConstraintMirror::ENTITY_B()));
 243   //  fixFeaturesList(aRefList);
 244   //}
 245
 246   if (!myFeatureStorage)
 247     myFeatureStorage = FeatureStoragePtr(new SketchSolver_FeatureStorage);
 248   myFeatureStorage->changeConstraint(theConstraint);
 249
 250   return true;
 251 }
 252
 253
 254 bool SketchSolver_Group::updateFeature(std::shared_ptr<SketchPlugin_Feature> theFeature)
 255 {
 256   std::set<ConstraintPtr> aConstraints =
 257       myFeatureStorage->getConstraints(std::dynamic_pointer_cast<ModelAPI_Feature>(theFeature));
 258   if (aConstraints.empty())
 259     return false;
 260   std::set<ConstraintPtr>::iterator aCIter = aConstraints.begin();
 261   for (; aCIter != aConstraints.end(); aCIter++) {
 262     ConstraintConstraintMap::iterator aSolConIter = myConstraints.find(*aCIter);
 263     if (aSolConIter == myConstraints.end() || !aSolConIter->first->data() ||
 264         !aSolConIter->first->data()->isValid())
 265       continue;
 266     myFeatureStorage->changeFeature(theFeature, aSolConIter->first);
 267     aSolConIter->second->addFeature(theFeature);
 268     aSolConIter->second->update();
 269   }
 270   return true;
 271 }
 272
 273 void SketchSolver_Group::moveFeature(std::shared_ptr<SketchPlugin_Feature> theFeature)
 274 {
 275   updateFeature(theFeature);
 276   // Temporary rigid constraint
 277   SolverConstraintPtr aConstraint =
 278       SketchSolver_Builder::getInstance()->createRigidConstraint(theFeature);
 279   if (!aConstraint)
 280     return;
 281   aConstraint->setGroup(this);
 282   aConstraint->setStorage(myStorage);
 283   if (aConstraint->error().empty())
 284     setTemporary(aConstraint);
 285 }
 286
 287 // ============================================================================
 288 //  Function: fixFeaturesList
 289 //  Class:    SketchSolver_Group
 290 //  Purpose:  Apply temporary rigid constraints for the list of features
 291 // ============================================================================
 292 void SketchSolver_Group::fixFeaturesList(AttributeRefListPtr theList)
 293 {
 294   std::list<ObjectPtr> aList = theList->list();
 295   std::list<ObjectPtr>::iterator anIt = aList.begin();
 296   std::list<FeaturePtr> aFeatures;
 297   // Sort features, at begining there are features used by Equal constraint
 298   for (; anIt != aList.end(); anIt++) {
 299     if (!(*anIt))
 300       continue;
 301     FeaturePtr aFeature = ModelAPI_Feature::feature(*anIt);
 302     std::set<ConstraintPtr> aConstraints = myFeatureStorage->getConstraints(aFeature);
 303     std::set<ConstraintPtr>::iterator aCIter = aConstraints.begin();
 304     for (; aCIter != aConstraints.end(); aCIter++)
 305       if ((*aCIter)->getKind() == SketchPlugin_ConstraintEqual::ID())
 306         break;
 307     if (aCIter != aConstraints.end())
 308       aFeatures.push_front(aFeature);
 309     else
 310       aFeatures.push_back(aFeature);
 311   }
 312
 313   std::list<FeaturePtr>::iterator aFeatIter = aFeatures.begin();
 314   for (; aFeatIter != aFeatures.end(); aFeatIter++) {
 315     SolverConstraintPtr aConstraint =
 316         SketchSolver_Builder::getInstance()->createRigidConstraint(*aFeatIter);
 317     if (!aConstraint)
 318       continue;
 319     aConstraint->setGroup(this);
 320     aConstraint->setStorage(myStorage);
 321     setTemporary(aConstraint);
 322   }
 323 }
 324
 325 // ============================================================================
 326 //  Function: addWorkplane
 327 //  Class:    SketchSolver_Group
 328 //  Purpose:  create workplane for the group
 329 // ============================================================================
 330 bool SketchSolver_Group::addWorkplane(CompositeFeaturePtr theSketch)
 331 {
 332   if (myWorkplaneID != SLVS_E_UNKNOWN || theSketch->getKind() != SketchPlugin_Sketch::ID())
 333     return false;  // the workplane already exists or the function parameter is not Sketch
 334
 335   mySketch = theSketch;
 336   updateWorkplane();
 337   return true;
 338 }
 339
 340 // ============================================================================
 341 //  Function: updateWorkplane
 342 //  Class:    SketchSolver_Group
 343 //  Purpose:  update parameters of workplane
 344 // ============================================================================
 345 bool SketchSolver_Group::updateWorkplane()
 346 {
 347   if (!myStorage) // Create storage if not exists
 348     myStorage = StoragePtr(new SketchSolver_Storage);
 349   SketchSolver_Builder* aBuilder = SketchSolver_Builder::getInstance();
 350
 351   std::vector<Slvs_Entity> anEntities;
 352   std::vector<Slvs_Param> aParams;
 353   if (!aBuilder->createWorkplane(mySketch, anEntities, aParams))
 354     return false;
 355
 356   if (myWorkplaneID == SLVS_E_UNKNOWN) {
 357     myWorkplaneID = anEntities.back().h;
 358     // Add new workplane elements
 359     std::vector<Slvs_Param>::iterator aParIter = aParams.begin();
 360     for (; aParIter != aParams.end(); aParIter++) {
 361       aParIter->h = SLVS_E_UNKNOWN; // the ID should be generated by storage
 362       aParIter->group = myID;
 363       aParIter->h = myStorage->addParameter(*aParIter);
 364     }
 365     std::vector<Slvs_Entity>::iterator anEntIter = anEntities.begin();
 366     for (; anEntIter != anEntities.end(); anEntIter++) {
 367       anEntIter->h = SLVS_E_UNKNOWN; // the ID should be generated by storage
 368       anEntIter->group = myID;
 369       anEntIter->wrkpl = myWorkplaneID;
 370       for (int i = 0; i < 4; i++)
 371         if (anEntIter->param[i] != SLVS_E_UNKNOWN)
 372           anEntIter->param[i] = aParams[anEntIter->param[i]-1].h;
 373       for (int i = 0; i < 4; i++)
 374         if (anEntIter->point[i] != SLVS_E_UNKNOWN)
 375           anEntIter->point[i] = anEntities[anEntIter->point[i]-1].h;
 376       anEntIter->h = myStorage->addEntity(*anEntIter);
 377     }
 378   } else {
 379     // Update existent workplane
 380     const Slvs_Entity& aWP = myStorage->getEntity(myWorkplaneID);
 381     const Slvs_Entity& anOrigin = myStorage->getEntity(aWP.point[0]);
 382     const Slvs_Entity& aNormal = myStorage->getEntity(aWP.normal);
 383     // Get parameters and update them
 384     Slvs_hParam aWPParams[7] = {
 385         anOrigin.param[0], anOrigin.param[1], anOrigin.param[2],
 386         aNormal.param[0], aNormal.param[1], aNormal.param[2], aNormal.param[3]
 387       };
 388     std::vector<Slvs_Param>::iterator aParIter = aParams.begin();
 389     for (int i = 0; aParIter != aParams.end(); aParIter++, i++) {
 390       Slvs_Param aParam = myStorage->getParameter(aWPParams[i]);
 391       aParam.val = aParIter->val;
 392       myStorage->updateParameter(aParam);
 393     }
 394   }
 395   return myWorkplaneID > 0;
 396 }
 397
 398 // ============================================================================
 399 //  Function: resolveConstraints
 400 //  Class:    SketchSolver_Group
 401 //  Purpose:  solve the set of constraints for the current group
 402 // ============================================================================
 403 bool SketchSolver_Group::resolveConstraints()
 404 {
 405   bool aResolved = false;
 406   if (myStorage->isNeedToResolve() && !isEmpty()) {
 407     myConstrSolver.setGroupID(myID);
 408     myStorage->initializeSolver(myConstrSolver);
 409
 410     int aResult = SLVS_RESULT_OKAY;
 411     try {
 412       if (myStorage->hasDuplicatedConstraint())
 413         aResult = SLVS_RESULT_INCONSISTENT;
 414       else {
 415         // To avoid overconstraint situation, we will remove temporary constraints one-by-one
 416         // and try to find the case without overconstraint
 417         int aNbTemp = (int)myTempConstraints.size();
 418         while (true) {
 419           aResult = myConstrSolver.solve();
 420           if (aResult == SLVS_RESULT_OKAY || aNbTemp <= 0)
 421             break;
 422           aNbTemp = myStorage->removeFirstTemporaryConstraint();
 423           myStorage->initializeSolver(myConstrSolver);
 424         }
 425       }
 426     } catch (...) {
 427       Events_Error::send(SketchSolver_Error::SOLVESPACE_CRASH(), this);
 428       return false;
 429     }
 430     if (aResult == SLVS_RESULT_OKAY) {  // solution succeeded, store results into correspondent attributes
 431       myFeatureStorage->blockEvents(true);
 432       ConstraintConstraintMap::iterator aConstrIter = myConstraints.begin();
 433       for (; aConstrIter != myConstraints.end(); aConstrIter++)
 434         aConstrIter->second->refresh();
 435       myFeatureStorage->blockEvents(false);
 436     } else if (!myConstraints.empty())
 437       Events_Error::send(SketchSolver_Error::CONSTRAINTS(), this);
 438
 439     aResolved = true;
 440   }
 441   removeTemporaryConstraints();
 442   myStorage->setNeedToResolve(false);
 443   return aResolved;
 444 }
 445
 446 // ============================================================================
 447 //  Function: mergeGroups
 448 //  Class:    SketchSolver_Group
 449 //  Purpose:  append specified group to the current group
 450 // ============================================================================
 451 void SketchSolver_Group::mergeGroups(const SketchSolver_Group& theGroup)
 452 {
 453   // If specified group is empty, no need to merge
 454   if (theGroup.isEmpty())
 455     return;
 456   if (!myFeatureStorage)
 457     myFeatureStorage = FeatureStoragePtr(new SketchSolver_FeatureStorage);
 458
 459   std::vector<ConstraintPtr> aComplexConstraints;
 460   ConstraintConstraintMap::const_iterator aConstrIter = theGroup.myConstraints.begin();
 461   // append simple constraints
 462   for (; aConstrIter != theGroup.myConstraints.end(); aConstrIter++)
 463     if (isComplexConstraint(aConstrIter->first))
 464       aComplexConstraints.push_back(aConstrIter->first);
 465     else
 466       changeConstraint(aConstrIter->first);
 467   // append complex constraints
 468   std::vector<ConstraintPtr>::iterator aComplexIter = aComplexConstraints.begin();
 469   for (; aComplexIter != aComplexConstraints.end(); aComplexIter++)
 470       changeConstraint(*aComplexIter);
 471 }
 472
 473 // ============================================================================
 474 //  Function: splitGroup
 475 //  Class:    SketchSolver_Group
 476 //  Purpose:  divide the group into several subgroups
 477 // ============================================================================
 478 void SketchSolver_Group::splitGroup(std::vector<SketchSolver_Group*>& theCuts)
 479 {
 480   // Obtain constraints, which should be separated
 481   FeatureStoragePtr aNewFeatStorage(new SketchSolver_FeatureStorage);
 482   std::vector<ConstraintPtr> anUnusedConstraints;
 483   ConstraintConstraintMap::iterator aCIter = myConstraints.begin();
 484   for ( ; aCIter != myConstraints.end(); aCIter++) {
 485     std::list<ConstraintPtr> aBaseConstraints = aCIter->second->constraints();
 486     std::list<ConstraintPtr>::iterator anIter = aBaseConstraints.begin();
 487     for (; anIter != aBaseConstraints.end(); anIter++)
 488       if (aNewFeatStorage->isInteract(*anIter)) {
 489         aNewFeatStorage->changeConstraint(*anIter);
 490       } else
 491         anUnusedConstraints.push_back(*anIter);
 492   }
 493
 494   // Check the unused constraints once again, because they may become interacted with new storage since adding constraints
 495   std::vector<ConstraintPtr>::iterator aUnuseIt = anUnusedConstraints.begin();
 496   while (aUnuseIt != anUnusedConstraints.end()) {
 497     if (aNewFeatStorage->isInteract(*aUnuseIt)) {
 498       size_t aShift = aUnuseIt - anUnusedConstraints.begin();
 499       anUnusedConstraints.erase(aUnuseIt);
 500       aUnuseIt = anUnusedConstraints.begin() + aShift;
 501       continue;
 502     }
 503     aUnuseIt++;
 504   }
 505
 506   std::vector<SketchSolver_Group*>::iterator aCutsIter;
 507   aUnuseIt = anUnusedConstraints.begin();
 508   for ( ; aUnuseIt != anUnusedConstraints.end(); aUnuseIt++) {
 509     // Remove unused constraints
 510     removeConstraint(*aUnuseIt);
 511     // Try to append constraint to already existent group
 512     for (aCutsIter = theCuts.begin(); aCutsIter != theCuts.end(); aCutsIter++)
 513       if ((*aCutsIter)->isInteract(*aUnuseIt)) {
 514         (*aCutsIter)->changeConstraint(*aUnuseIt);
 515         break;
 516       }
 517     if (aCutsIter == theCuts.end()) {
 518       // Add new group
 519       SketchSolver_Group* aGroup = new SketchSolver_Group(mySketch);
 520       aGroup->changeConstraint(*aUnuseIt);
 521       theCuts.push_back(aGroup);
 522     }
 523   }
 524 }
 525
 526 // ============================================================================
 527 //  Function: isConsistent
 528 //  Class:    SketchSolver_Group
 529 //  Purpose:  search removed entities and constraints
 530 // ============================================================================
 531 bool SketchSolver_Group::isConsistent()
 532 {
 533   if (!myFeatureStorage) // no one constraint is initialized yet
 534     return true;
 535
 536   bool aResult = myFeatureStorage->isConsistent();
 537   if (!aResult) {
 538     // remove invalid entities
 539     ConstraintConstraintMap::iterator aCIter = myConstraints.begin();
 540     while (aCIter != myConstraints.end()) {
 541       std::list<ConstraintPtr> aConstraints = aCIter->second->constraints();
 542       std::list<ConstraintPtr>::iterator anIt = aConstraints.begin();
 543       for (; anIt != aConstraints.end(); anIt++)
 544         if (!(*anIt)->data() || !(*anIt)->data()->isValid())
 545           if (aCIter->second->remove(*anIt)) {
 546             // the constraint is fully removed, detach it from the list
 547             ConstraintConstraintMap::iterator aTmpIt = aCIter++;
 548             myFeatureStorage->removeConstraint(aTmpIt->first);
 549             myConstraints.erase(aTmpIt);
 550             break;
 551           }
 552       if (anIt == aConstraints.end())
 553         aCIter++;
 554     }
 555   }
 556   return aResult;
 557 }
 558
 559 // ============================================================================
 560 //  Function: removeTemporaryConstraints
 561 //  Class:    SketchSolver_Group
 562 //  Purpose:  remove all transient SLVS_C_WHERE_DRAGGED constraints after
 563 //            resolving the set of constraints
 564 // ============================================================================
 565 void SketchSolver_Group::removeTemporaryConstraints()
 566 {
 567   myTempConstraints.clear();
 568   myStorage->removeTemporaryConstraints();
 569   // Clean lists of removed entities in the storage
 570   std::set<Slvs_hParam> aRemPar;
 571   std::set<Slvs_hEntity> aRemEnt;
 572   std::set<Slvs_hConstraint> aRemCon;
 573   myStorage->getRemoved(aRemPar, aRemEnt, aRemCon);
 574   myStorage->setNeedToResolve(false);
 575 }
 576
 577 // ============================================================================
 578 //  Function: removeConstraint
 579 //  Class:    SketchSolver_Group
 580 //  Purpose:  remove constraint and all unused entities
 581 // ============================================================================
 582 void SketchSolver_Group::removeConstraint(ConstraintPtr theConstraint)
 583 {
 584   myFeatureStorage->removeConstraint(theConstraint);
 585   ConstraintConstraintMap::iterator aCIter = myConstraints.begin();
 586   for (; aCIter != myConstraints.end(); aCIter++)
 587     if (aCIter->second->hasConstraint(theConstraint)) {
 588       if (!aCIter->second->remove(theConstraint)) // the constraint is not fully removed
 589         aCIter = myConstraints.end();
 590       break;
 591     }
 592   if (aCIter != myConstraints.end())
 593     myConstraints.erase(aCIter);
 594 }
 595
 596 // ============================================================================
 597 //  Function: isComplexConstraint
 598 //  Class:    SketchSolver_Group
 599 //  Purpose:  verifies the constraint is complex, i.e. it needs another constraints to be created before
 600 // ============================================================================
 601 bool SketchSolver_Group::isComplexConstraint(FeaturePtr theConstraint)
 602 {
 603   return theConstraint->getKind() == SketchPlugin_ConstraintFillet::ID() ||
 604          theConstraint->getKind() == SketchPlugin_ConstraintMirror::ID() ||
 605          theConstraint->getKind() == SketchPlugin_ConstraintTangent::ID();
 606 }
 607
 608 // ============================================================================
 609 //  Function: setTemporary
 610 //  Class:    SketchSolver_Group
 611 //  Purpose:  append given constraint to th group of temporary constraints
 612 // ============================================================================
 613 void SketchSolver_Group::setTemporary(SolverConstraintPtr theConstraint)
 614 {
 615   theConstraint->makeTemporary();
 616   myTempConstraints.insert(theConstraint);
 617 }
 618