Salome HOME
SketchSolver library refactoring
[modules/shaper.git] / src / SketchSolver / SketchSolver_Constraint.cpp
index caccc5223b1303460ed3fd6dd265e4cd728c4c9c..2cdd24eb796e3840682658844ca98449624b604f 100644 (file)
-// File:    SketchSolver_Constraint.cpp
-// Created: 27 May 2014
-// Author:  Artem ZHIDKOV
-
-#include "SketchSolver_Constraint.h"
-#include <SketchSolver_Solver.h>
+#include <SketchSolver_Constraint.h>
+#include <SketchSolver_Group.h>
+#include <SketchSolver_Error.h>
 
+#include <SketchPlugin_Arc.h>
+#include <SketchPlugin_Circle.h>
 #include <SketchPlugin_Line.h>
 #include <SketchPlugin_Point.h>
-#include <SketchPlugin_Circle.h>
-#include <SketchPlugin_Arc.h>
-#include <SketchPlugin_ConstraintCoincidence.h>
-#include <SketchPlugin_ConstraintDistance.h>
-#include <SketchPlugin_ConstraintLength.h>
-#include <SketchPlugin_ConstraintParallel.h>
-#include <SketchPlugin_ConstraintPerpendicular.h>
-#include <SketchPlugin_ConstraintRadius.h>
-#include <SketchPlugin_ConstraintRigid.h>
-
-#include <ModelAPI_AttributeRefAttr.h>
-#include <ModelAPI_Data.h>
-#include <ModelAPI_Document.h>
-#include <ModelAPI_Object.h>
-#include <ModelAPI_ResultConstruction.h>
 
 #include <GeomDataAPI_Point.h>
 #include <GeomDataAPI_Point2D.h>
-#include <GeomAPI_Edge.h>
-#include <GeomAPI_Shape.h>
+#include <ModelAPI_AttributeDouble.h>
+#include <ModelAPI_ResultConstruction.h>
+
 
-/// Possible types of attributes (used to determine constraint type)
-enum AttrType
+SketchSolver_Constraint::SketchSolver_Constraint(
+    ConstraintPtr  theConstraint)
+  : myBaseConstraint(theConstraint),
+    myGroup(0)
+{
+}
+
+SketchSolver_Constraint::~SketchSolver_Constraint()
 {
-  UNKNOWN,  // Something wrong during type determination
-  POINT2D,
-  POINT3D,
-  LINE,
-  CIRCLE,
-  ARC
-};
+  std::map<AttributePtr, Slvs_hParam>::const_iterator anIt1 = myValueMap.begin();
+  for (; anIt1 != myValueMap.end(); anIt1++)
+    myStorage->removeParameter(anIt1->second);
+  myValueMap.clear();
 
-/// Calculate type of the attribute
-static AttrType typeOfAttribute(std::shared_ptr<ModelAPI_Attribute> theAttribute);
+  std::map<AttributePtr, Slvs_hEntity>::const_iterator anIt2 = myAttributeMap.begin();
+  for (; anIt2 != myAttributeMap.end(); anIt2++)
+    myStorage->removeEntity(anIt2->second);
+  myAttributeMap.clear();
 
-SketchSolver_Constraint::SketchSolver_Constraint()
-    : myConstraint(std::shared_ptr<SketchPlugin_Constraint>()),
-      myType(SLVS_C_UNKNOWN),
-      myAttributesList()
+  std::map<FeaturePtr, Slvs_hEntity>::const_iterator anIt3 =  myFeatureMap.begin();
+  for (; anIt3 != myFeatureMap.end(); anIt3++)
+    myStorage->removeEntity(anIt3->second);
+  myFeatureMap.clear();
+
+  std::vector<Slvs_hConstraint>::const_iterator anIt4 = mySlvsConstraints.begin();
+  for (; anIt4 != mySlvsConstraints.end(); anIt4++)
+    myStorage->removeConstraint(*anIt4);
+  mySlvsConstraints.clear();
+}
+
+void SketchSolver_Constraint::setStorage(StoragePtr theStorage)
 {
+  myStorage = theStorage;
+  process();
 }
 
-SketchSolver_Constraint::SketchSolver_Constraint(
-    std::shared_ptr<SketchPlugin_Constraint> theConstraint)
-    : myConstraint(theConstraint),
-      myAttributesList()
-{
-  myType = getType(myConstraint);
-}
-
-const int& SketchSolver_Constraint::getType(
-    std::shared_ptr<SketchPlugin_Constraint> theConstraint)
-{
-  myType = SLVS_C_UNKNOWN;
-  if (!theConstraint)
-    return getType();
-
-  DataPtr aConstrData = theConstraint->data();
-  if (!aConstrData || !aConstrData->isValid())
-    return getType();
-
-  // Assign empty names of attributes
-  myAttributesList.clear();
-  for (int i = 0; i < CONSTRAINT_ATTR_SIZE; i++)
-    myAttributesList.push_back(std::string());
-
-  const std::string& aConstraintKind = theConstraint->getKind();
-  // Constraint for coincidence of two points
-  if (aConstraintKind.compare(SketchPlugin_ConstraintCoincidence::ID()) == 0) {
-    int anAttrPos = 0;
-    // Verify the constraint has only two attributes and they are points
-    int aPt2d = 0;  // bit-mapped field, each bit indicates whether the attribute is 2D point
-    int aPt3d = 0;  // bit-mapped field, the same information for 3D points
-    for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++) {
-      std::shared_ptr<ModelAPI_Attribute> anAttr = 
-          aConstrData->attribute(SketchPlugin_Constraint::ATTRIBUTE(indAttr));
-      if (!anAttr)
+void SketchSolver_Constraint::setGroup(SketchSolver_Group* theGroup)
+{
+  myGroup = theGroup;
+  process();
+}
+
+
+void SketchSolver_Constraint::process()
+{
+  cleanErrorMsg();
+  if (!myBaseConstraint || !myStorage || myGroup == 0) {
+    /// TODO: Put error message here
+    return;
+  }
+  if (!mySlvsConstraints.empty()) // some data is changed, update constraint
+    update(myBaseConstraint);
+
+  int aConstrType = getType();
+  double aValue = 0.0;
+  std::vector<Slvs_hEntity> anAttributes;
+  getAttributes(aValue, anAttributes);
+  if (!myErrorMsg.empty())
+    return;
+
+  Slvs_hGroup aGroupID = myGroup->getId();
+  Slvs_hEntity aWorkplaneID = myGroup->getWorkplaneId();
+  Slvs_Constraint aConstraint;
+  if (mySlvsConstraints.empty())
+    aConstraint = Slvs_MakeConstraint(SLVS_C_UNKNOWN, aGroupID, aConstrType, aWorkplaneID,
+        aValue, anAttributes[0], anAttributes[1], anAttributes[2], anAttributes[3]);
+  else {
+    aConstraint = myStorage->getConstraint(mySlvsConstraints[0]);
+    aConstraint.valA = aValue;
+    static const int aNbAttrs = 6;
+    Slvs_hEntity* aConstrAttrs[aNbAttrs] = {
+        &aConstraint.ptA, &aConstraint.ptB,
+        &aConstraint.entityA, &aConstraint.entityB,
+        &aConstraint.entityC, &aConstraint.entityD};
+    std::vector<Slvs_hEntity>::const_iterator anIter = anAttributes.begin();
+    for (int i = 0; i < aNbAttrs && anIter != anAttributes.end(); i++, anIter++)
+      *(aConstrAttrs[i]) = *anIter;
+  }
+
+  Slvs_hConstraint anID = myStorage->addConstraint(aConstraint);
+  if (mySlvsConstraints.empty())
+    mySlvsConstraints.push_back(anID);
+  else
+    mySlvsConstraints[0] = anID;
+  adjustConstraint();
+}
+
+void SketchSolver_Constraint::update(ConstraintPtr theConstraint)
+{
+  cleanErrorMsg();
+  if (theConstraint && theConstraint != myBaseConstraint) {
+    if (theConstraint->getKind() != myBaseConstraint->getKind())
+      return;
+    remove(myBaseConstraint);
+    myBaseConstraint = theConstraint;
+    process();
+  }
+
+  // Update all attributes
+  int aType;
+  std::map<Slvs_hEntity, Slvs_hEntity> aRelocationMap;
+  std::map<FeaturePtr, Slvs_hEntity>::iterator aFeatIter = myFeatureMap.begin();
+  for (; aFeatIter != myFeatureMap.end(); aFeatIter++) {
+    Slvs_hEntity aPrevID = aFeatIter->second;
+    aFeatIter->second = changeEntity(aFeatIter->first, aType);
+    if (aFeatIter->second != aPrevID)
+      aRelocationMap[aPrevID] = aFeatIter->second;
+  }
+  std::map<AttributePtr, Slvs_hEntity>::iterator anAttrIter = myAttributeMap.begin();
+  for (; anAttrIter != myAttributeMap.end(); anAttrIter++) {
+    Slvs_hEntity aPrevID = anAttrIter->second;
+    anAttrIter->second = changeEntity(anAttrIter->first, aType);
+    if (anAttrIter->second != aPrevID)
+      aRelocationMap[aPrevID] = anAttrIter->second;
+  }
+
+  // Value if exists
+  AttributeDoublePtr aValueAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
+    myBaseConstraint->data()->attribute(SketchPlugin_Constraint::VALUE()));
+  double aValue = aValueAttr ? aValueAttr->value() : 0.0;
+
+  // Update constraint
+  std::vector<Slvs_hConstraint>::iterator aCIter = mySlvsConstraints.begin();
+  for (; aCIter != mySlvsConstraints.end(); aCIter++) {
+    Slvs_Constraint aConstraint = myStorage->getConstraint(*aCIter);
+    aConstraint.valA = aValue;
+    Slvs_hEntity* aCoeffs[6] = {
+        &aConstraint.ptA, &aConstraint.ptB,
+        &aConstraint.entityA, &aConstraint.entityB,
+        &aConstraint.entityC, &aConstraint.entityD};
+    for (int i = 0; i < 6; i++) {
+      if (*(aCoeffs[i]) == SLVS_E_UNKNOWN)
         continue;
-      switch (typeOfAttribute(anAttr)) {
-        case POINT2D:  // the attribute is a 2D point
-          aPt2d |= (1 << indAttr);
-          myAttributesList[anAttrPos++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-          break;
-        case POINT3D:  // the attribute is a 3D point
-          aPt3d |= (1 << indAttr);
-          myAttributesList[anAttrPos++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-          break;
-        default:
-          // Attribute neither 2D nor 3D point is not supported by this type of constraint
-          return getType();
-      }
+      std::map<Slvs_hEntity, Slvs_hEntity>::iterator aFound = aRelocationMap.find(*(aCoeffs[i]));
+      if (aFound != aRelocationMap.end())
+        *(aCoeffs[i]) = aFound->second;
     }
-    // The constrained points should be in first and second positions,
-    // so the expected value of aPt2d or aPt3d is 3
-    if ((aPt2d == 3 && aPt3d == 0) || (aPt2d == 0 && aPt3d == 3))
-      myType = SLVS_C_POINTS_COINCIDENT;
-    // Constraint parameters are wrong
-    return getType();
-  }
-
-  // Constraint for distance between point and another entity
-  if (aConstraintKind.compare(SketchPlugin_ConstraintDistance::ID()) == 0) {
-    int aNbPoints = 0;
-    int aNbEntities = 0;
-    for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++) {
-      std::shared_ptr<ModelAPI_Attribute> anAttr = 
-          aConstrData->attribute(SketchPlugin_Constraint::ATTRIBUTE(indAttr));
-      switch (typeOfAttribute(anAttr)) {
-        case POINT2D:
-        case POINT3D:
-          myAttributesList[aNbPoints++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-          break;
-        case LINE:
-          // entities are placed starting from SketchPlugin_Constraint::ENTITY_C() attribute
-          myAttributesList[2 + aNbEntities++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-          myType = SLVS_C_PT_LINE_DISTANCE;
-          break;
-      }
+    *aCIter = myStorage->addConstraint(aConstraint);
+  }
+  adjustConstraint();
+}
+
+bool SketchSolver_Constraint::remove(ConstraintPtr theConstraint)
+{
+  cleanErrorMsg();
+  if (theConstraint && theConstraint != myBaseConstraint)
+    return false;
+  bool isFullyRemoved = myStorage->removeConstraint(mySlvsConstraints.front());
+  if (isFullyRemoved) {
+    myFeatureMap.clear();
+    myAttributeMap.clear();
+    myValueMap.clear();
+  } else
+    cleanRemovedEntities();
+  return true;
+}
+
+void SketchSolver_Constraint::cleanRemovedEntities()
+{
+  std::set<Slvs_hParam> aRemovedParams;
+  std::set<Slvs_hEntity> aRemovedEntities;
+  std::set<Slvs_hConstraint> aRemovedConstraints;
+  myStorage->getRemoved(aRemovedParams, aRemovedEntities, aRemovedConstraints);
+  std::map<FeaturePtr, Slvs_hEntity>::iterator aFeatIt = myFeatureMap.begin();
+  while (aFeatIt != myFeatureMap.end()) {
+    if (aRemovedEntities.find(aFeatIt->second) == aRemovedEntities.end()) {
+      aFeatIt++;
+      continue;
+    }
+    std::map<FeaturePtr, Slvs_hEntity>::iterator aTmpIter = aFeatIt++;
+    myFeatureMap.erase(aTmpIter);
+  }
+  std::map<AttributePtr, Slvs_hEntity>::iterator anAttrIt = myAttributeMap.begin();
+  while (anAttrIt != myAttributeMap.end()) {
+    if (aRemovedEntities.find(anAttrIt->second) == aRemovedEntities.end()) {
+      anAttrIt++;
+      continue;
+    }
+    std::map<AttributePtr, Slvs_hEntity>::iterator aTmpIter = anAttrIt++;
+    myAttributeMap.erase(aTmpIter);
+  }
+  std::map<AttributePtr, Slvs_hParam>::iterator aValIt = myValueMap.begin();
+  while (aValIt != myValueMap.end()) {
+    if (aRemovedParams.find(aValIt->second) == aRemovedParams.end()) {
+      aValIt++;
+      continue;
     }
-    // Verify the correctness of constraint arguments
-    if (aNbPoints == 2 && aNbEntities == 0)
-      myType = SLVS_C_PT_PT_DISTANCE;
-    else if (aNbPoints != 1 || aNbEntities != 1)
-      myType = SLVS_C_UNKNOWN;
-    return getType();
-  }
-
-  // Constraint for the given length of a line
-  if (aConstraintKind.compare(SketchPlugin_ConstraintLength::ID()) == 0) {
-    int aNbLines = 0;
-    for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++) {
-      std::shared_ptr<ModelAPI_Attribute> anAttr = 
-          aConstrData->attribute(SketchPlugin_Constraint::ATTRIBUTE(indAttr));
-      if (typeOfAttribute(anAttr) == LINE)
-        myAttributesList[aNbLines++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-    }
-    if (aNbLines == 1)
-      myType = SLVS_C_PT_PT_DISTANCE;
-    return getType();
-  }
-
-  // Constraint for two parallel/perpendicular lines
-  bool isParallel = (aConstraintKind.compare(SketchPlugin_ConstraintParallel::ID()) == 0);
-  bool isPerpendicular = (aConstraintKind.compare(SketchPlugin_ConstraintPerpendicular::ID()) == 0);
-  if (isParallel || isPerpendicular) {
-    int aNbEntities = 2;  // lines in SolveSpace constraints should start from SketchPlugin_Constraint::ENTITY_C() attribute
-    for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++) {
-      std::shared_ptr<ModelAPI_Attribute> anAttr = 
-          aConstrData->attribute(SketchPlugin_Constraint::ATTRIBUTE(indAttr));
-      if (typeOfAttribute(anAttr) == LINE)
-        myAttributesList[aNbEntities++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-    }
-    if (aNbEntities == 4)
-      myType = isParallel ? SLVS_C_PARALLEL : SLVS_C_PERPENDICULAR;
-    return getType();
-  }
-
-  // Constraint for radius of a circle or an arc of circle
-  if (aConstraintKind.compare(SketchPlugin_ConstraintRadius::ID()) == 0) {
-    int aNbEntities = 2;  // lines in SolveSpace constraints should started from SketchPlugin_Constraint::ENTITY_C() attribute
-    for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++) {
-      std::shared_ptr<ModelAPI_Attribute> anAttr = 
-          aConstrData->attribute(SketchPlugin_Constraint::ATTRIBUTE(indAttr));
-      AttrType aType = typeOfAttribute(anAttr);
-      if (aType == CIRCLE || aType == ARC)
-        myAttributesList[aNbEntities++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-    }
-    if (aNbEntities == 3)
-      myType = SLVS_C_DIAMETER;
-    return getType();
-  }
-
-  // Constraint for fixed entity
-  if (aConstraintKind.compare(SketchPlugin_ConstraintRigid::ID()) == 0) {
-    // Verify that only one entity is filled
-    int aNbAttrs = 0;
-    for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++) {
-      std::shared_ptr<ModelAPI_Attribute> anAttr = 
-          aConstrData->attribute(SketchPlugin_Constraint::ATTRIBUTE(indAttr));
-      AttrType aType = typeOfAttribute(anAttr);
-      if (aType != UNKNOWN)
-        myAttributesList[aNbAttrs++] = SketchPlugin_Constraint::ATTRIBUTE(indAttr);
-    }
-    if (aNbAttrs == 1)
-      myType = SLVS_C_WHERE_DRAGGED;
-    return getType();
-  }
-
-  /// \todo Implement other kind of constraints
-
-  return getType();
-}
-
-// ================= Auxiliary functions ==============================
-AttrType typeOfAttribute(std::shared_ptr<ModelAPI_Attribute> theAttribute)
-{
-  std::shared_ptr<ModelAPI_AttributeRefAttr> anAttrRef = std::dynamic_pointer_cast<
-      ModelAPI_AttributeRefAttr>(theAttribute);
-  if (!anAttrRef)
-    return UNKNOWN;
-
-  if (anAttrRef->isObject()) {
+    std::map<AttributePtr, Slvs_hParam>::iterator aTmpIter = aValIt++;
+    myValueMap.erase(aTmpIter);
+  }
+}
+
+void SketchSolver_Constraint::getAttributes(
+    double& theValue,
+    std::vector<Slvs_hEntity>& theAttributes)
+{
+  static const int anInitNbOfAttr = 4;
+  theAttributes.assign(anInitNbOfAttr, SLVS_E_UNKNOWN);
+
+  DataPtr aData = myBaseConstraint->data();
+
+  AttributeDoublePtr aValueAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
+    aData->attribute(SketchPlugin_Constraint::VALUE()));
+  theValue = aValueAttr ? aValueAttr->value() : 0.0;
+
+  int aPtInd = 0; // index of first point in the list of attributes
+  int aEntInd = 2; // index of first antity in the list of attributes
+  std::list<AttributePtr> aConstrAttrs = aData->attributes(ModelAPI_AttributeRefAttr::typeId());
+  std::list<AttributePtr>::iterator anIter = aConstrAttrs.begin();
+  for (; anIter != aConstrAttrs.end(); anIter++) {
+    AttributeRefAttrPtr aRefAttr =
+        std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*anIter);
+    if (!aRefAttr || !aRefAttr->isInitialized()) {
+      myErrorMsg = SketchSolver_Error::NOT_INITIALIZED();
+      return;
+    }
+
+    int aType = SLVS_E_UNKNOWN; // type of created entity
+    Slvs_hEntity anEntity = myGroup->getAttributeId(aRefAttr);
+    if (anEntity == SLVS_E_UNKNOWN)
+      anEntity = changeEntity(aRefAttr, aType);
+    else {
+      Slvs_Entity anEnt = myStorage->getEntity(anEntity);
+      aType = anEnt.type;
+    }
+
+    if (aType == SLVS_E_UNKNOWN)
+      continue;
+    else if (aType == SLVS_E_POINT_IN_2D || aType == SLVS_E_POINT_IN_3D)
+      theAttributes[aPtInd++] = anEntity; // the point is created
+    else { // another entity (not a point) is created
+      if (aEntInd < anInitNbOfAttr)
+        theAttributes[aEntInd] = anEntity;
+      else
+        theAttributes.push_back(anEntity);
+      aEntInd++;
+    }
+  }
+}
+
+Slvs_hEntity SketchSolver_Constraint::changeEntity(AttributeRefAttrPtr theAttribute, int& theType)
+{
+  // Convert the object of the attribute to the feature
+  FeaturePtr aFeature;
+  if (theAttribute->isObject() && theAttribute->object()) {
     ResultConstructionPtr aRC = std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(
-        anAttrRef->object());
-    if (!aRC || !aRC->shape())
-      return UNKNOWN;
-
-    if (aRC->shape()->isVertex())
-      return POINT3D;
-    else if (aRC->shape()->isEdge()) {
-      std::shared_ptr<GeomAPI_Edge> anEdge = std::dynamic_pointer_cast<GeomAPI_Edge>(
-          aRC->shape());
-      if (anEdge->isLine())
-        return LINE;
-      else if (anEdge->isCircle())
-        return CIRCLE;
-      else if (anEdge->isArc())
-        return ARC;
+        theAttribute->object());
+    if (!aRC) {
+      myErrorMsg = SketchSolver_Error::NOT_INITIALIZED();
+      return SLVS_E_UNKNOWN;
     }
+    std::shared_ptr<ModelAPI_Document> aDoc = aRC->document();
+    aFeature = aDoc->feature(aRC);
+
+    return changeEntity(aFeature, theType);
+  }
+
+  return changeEntity(theAttribute->attr(), theType);
+}
+
+Slvs_hEntity SketchSolver_Constraint::changeEntity(AttributePtr theEntity, int& theType)
+{
+  Slvs_hEntity aResult = SLVS_E_UNKNOWN;
+  if (!theEntity->isInitialized())
+    return SLVS_E_UNKNOWN;
+
+  // If the entity is already in the group, try to find it
+  std::map<std::shared_ptr<ModelAPI_Attribute>, Slvs_hEntity>::const_iterator anEntIter =
+      myAttributeMap.find(theEntity);
+  Slvs_Entity aCurrentEntity;
+  aCurrentEntity.h = SLVS_E_UNKNOWN;
+  if (anEntIter != myAttributeMap.end())
+    aCurrentEntity = myStorage->getEntity(anEntIter->second);
+  else {
+    aResult = myGroup->getAttributeId(theEntity);
+    if (aResult != SLVS_E_UNKNOWN) {
+      Slvs_Entity anEnt = myStorage->getEntity(aResult);
+      theType = anEnt.type;
+      return aResult;
+    }
+  }
+
+  Slvs_hGroup aGroupID = myGroup->getId();
+  // Point in 3D
+  std::shared_ptr<GeomDataAPI_Point> aPoint =
+      std::dynamic_pointer_cast<GeomDataAPI_Point>(theEntity);
+  if (aPoint) {
+    double aXYZ[3] = {aPoint->x(), aPoint->y(), aPoint->z()};
+    Slvs_hParam aParams[3];
+    for (int i = 0; i < 3; i++) {
+      Slvs_Param aPar = aCurrentEntity.h != SLVS_E_UNKNOWN ?
+          myStorage->getParameter(aCurrentEntity.param[i]) :
+          Slvs_MakeParam(SLVS_E_UNKNOWN, aGroupID, 0.0);
+      aPar.val = aXYZ[i];
+      aParams[i] = myStorage->addParameter(aPar);
+    }
+
+    if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
+      aCurrentEntity = Slvs_MakePoint3d(SLVS_E_UNKNOWN, aGroupID, aParams[0], aParams[1], aParams[2]);
+    else { // update entity data
+      for (int i = 0; i < 3; i++)
+        aCurrentEntity.param[i] = aParams[i];
+    }
+    aResult = myStorage->addEntity(aCurrentEntity);
   } else {
-    const std::string aType = anAttrRef->attr()->attributeType();
-    if (aType == GeomDataAPI_Point2D::type())
-      return POINT2D;
-    if (aType == GeomDataAPI_Point2D::type())
-      return POINT2D;
+    // All entities except 3D points are created on workplane. So, if there is no workplane yet, then error
+    Slvs_hEntity aWorkplaneID = myGroup->getWorkplaneId();
+    if (aWorkplaneID == SLVS_E_UNKNOWN)
+      return SLVS_E_UNKNOWN;
+
+    // Point in 2D
+    std::shared_ptr<GeomDataAPI_Point2D> aPoint2D =
+        std::dynamic_pointer_cast<GeomDataAPI_Point2D>(theEntity);
+    if (aPoint2D) {
+      double aXY[2] = {aPoint2D->x(), aPoint2D->y()};
+      Slvs_hParam aParams[2];
+      for (int i = 0; i < 2; i++) {
+        Slvs_Param aPar = aCurrentEntity.h != SLVS_E_UNKNOWN ?
+            myStorage->getParameter(aCurrentEntity.param[i]) :
+            Slvs_MakeParam(SLVS_E_UNKNOWN, aGroupID, 0.0);
+        aPar.val = aXY[i];
+        aParams[i] = myStorage->addParameter(aPar);
+      }
+
+      if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
+        aCurrentEntity = Slvs_MakePoint2d(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID, aParams[0], aParams[1]);
+      else { // update entity data
+        for (int i = 0; i < 2; i++)
+          aCurrentEntity.param[i] = aParams[i];
+      }
+      aResult = myStorage->addEntity(aCurrentEntity);
+    } else {
+      // Scalar value (used for the distance entities)
+      AttributeDoublePtr aScalar = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(theEntity);
+      if (aScalar) {
+        Slvs_Param aParam = aCurrentEntity.h != SLVS_E_UNKNOWN ?
+            myStorage->getParameter(aCurrentEntity.param[0]) :
+            Slvs_MakeParam(SLVS_E_UNKNOWN, aGroupID, 0.0);
+        aParam.val = aScalar->value();
+        Slvs_hParam aValue = myStorage->addParameter(aParam);
+
+        if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
+          aCurrentEntity = Slvs_MakeDistance(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID, aValue);
+        else
+          aCurrentEntity.param[0] = aValue;
+        aResult = myStorage->addEntity(aCurrentEntity);
+      }
+    }
+  }
+
+  myAttributeMap[theEntity] = aResult;
+  theType = aCurrentEntity.type;
+  return aResult;
+}
+
+Slvs_hEntity SketchSolver_Constraint::changeEntity(FeaturePtr theEntity, int& theType)
+{
+  Slvs_hEntity aResult = SLVS_E_UNKNOWN;
+  if (!theEntity->data()->isValid())
+    return SLVS_E_UNKNOWN;
+  // If the entity is already in the group, try to find it
+  std::map<FeaturePtr, Slvs_hEntity>::const_iterator anEntIter = myFeatureMap.find(theEntity);
+  Slvs_Entity aCurrentEntity;
+  aCurrentEntity.h = SLVS_E_UNKNOWN;
+  if (anEntIter != myFeatureMap.end())
+    aCurrentEntity = myStorage->getEntity(anEntIter->second);
+  else {
+    aResult = myGroup->getFeatureId(theEntity);
+    if (aResult != SLVS_E_UNKNOWN) {
+      Slvs_Entity anEnt = myStorage->getEntity(aResult);
+      theType = anEnt.type;
+      return aResult;
+    }
   }
 
-  return UNKNOWN;
+  Slvs_hGroup aGroupID = myGroup->getId();
+  Slvs_hEntity aWorkplaneID = myGroup->getWorkplaneId();
+  // SketchPlugin features
+  std::shared_ptr<SketchPlugin_Feature> aFeature = std::dynamic_pointer_cast<
+      SketchPlugin_Feature>(theEntity);
+  if (aFeature) {  // Verify the feature by its kind
+    const std::string& aFeatureKind = aFeature->getKind();
+    AttributePtr anAttribute;
+    int anAttrType;
+
+    // Line
+    if (aFeatureKind == SketchPlugin_Line::ID()) {
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Line::START_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aStart = changeEntity(anAttribute, anAttrType);
+
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Line::END_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aEnd = changeEntity(anAttribute, anAttrType);
+
+      if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
+        aCurrentEntity = Slvs_MakeLineSegment(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID, aStart, aEnd);
+      else {
+        aCurrentEntity.point[0] = aStart;
+        aCurrentEntity.point[1] = aEnd;
+      }
+      aResult = myStorage->addEntity(aCurrentEntity);
+    }
+    // Circle
+    else if (aFeatureKind == SketchPlugin_Circle::ID()) {
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Circle::CENTER_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aCenter = changeEntity(anAttribute, anAttrType);
+
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Circle::RADIUS_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aRadius = changeEntity(anAttribute, anAttrType);
+
+      if (aCurrentEntity.h == SLVS_E_UNKNOWN) { // New entity
+        Slvs_Entity aWorkplane = myStorage->getEntity(aWorkplaneID);
+        aCurrentEntity = Slvs_MakeCircle(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID,
+                                         aCenter, aWorkplane.normal, aRadius);
+      } else {
+        aCurrentEntity.point[0] = aCenter;
+        aCurrentEntity.distance = aRadius;
+      }
+      aResult = myStorage->addEntity(aCurrentEntity);
+    }
+    // Arc
+    else if (aFeatureKind == SketchPlugin_Arc::ID()) {
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Arc::CENTER_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aCenter = changeEntity(anAttribute, anAttrType);
+
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Arc::START_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aStart = changeEntity(anAttribute, anAttrType);
+
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Arc::END_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      Slvs_hEntity aEnd = changeEntity(anAttribute, anAttrType);
+
+      if (aCurrentEntity.h == SLVS_E_UNKNOWN) { // New entity
+        Slvs_Entity aWorkplane = myStorage->getEntity(aWorkplaneID);
+        aCurrentEntity = Slvs_MakeArcOfCircle(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID,
+                                              aWorkplane.normal, aCenter, aStart, aEnd);
+      } else {
+        aCurrentEntity.point[0] = aCenter;
+        aCurrentEntity.point[1] = aStart;
+        aCurrentEntity.point[2] = aEnd;
+      }
+      aResult = myStorage->addEntity(aCurrentEntity);
+    }
+    // Point (it has low probability to be an attribute of constraint, so it is checked at the end)
+    else if (aFeatureKind == SketchPlugin_Point::ID()) {
+      anAttribute = aFeature->data()->attribute(SketchPlugin_Point::COORD_ID());
+      if (!anAttribute->isInitialized()) return SLVS_E_UNKNOWN;
+      // Both the sketch point and its attribute (coordinates) link to the same SolveSpace point identifier
+      aResult = changeEntity(anAttribute, anAttrType);
+      aCurrentEntity.type = SLVS_E_POINT_IN_3D;
+    }
+  }
+
+  if (aResult != SLVS_E_UNKNOWN) {
+    myFeatureMap[theEntity] = aResult;
+    theType = aCurrentEntity.type;
+  }
+  return aResult;
+}
+
+std::list<ConstraintPtr> SketchSolver_Constraint::constraints() const
+{
+  std::list<ConstraintPtr> aConstraints;
+  aConstraints.push_back(myBaseConstraint);
+  return aConstraints;
+}
+
+void SketchSolver_Constraint::refresh()
+{
+  cleanErrorMsg();
+  std::map<AttributePtr, Slvs_hEntity>::iterator anAttrIter = myAttributeMap.begin();
+  for (; anAttrIter != myAttributeMap.end(); anAttrIter++) {
+    std::shared_ptr<GeomDataAPI_Point> aPoint =
+        std::dynamic_pointer_cast<GeomDataAPI_Point>(anAttrIter->first);
+    if (aPoint) {
+      Slvs_Entity anEntity = myStorage->getEntity(anAttrIter->second);
+      double aXYZ[3];
+      for (int i = 0; i < 3; i++) {
+        Slvs_Param aPar = myStorage->getParameter(anEntity.param[i]);
+        aXYZ[i] = aPar.val;
+      }
+      aPoint->setValue(aXYZ[0], aXYZ[1], aXYZ[2]);
+    } else {
+      // Point in 2D
+      std::shared_ptr<GeomDataAPI_Point2D> aPoint2D =
+          std::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttrIter->first);
+      if (aPoint2D) {
+        Slvs_Entity anEntity = myStorage->getEntity(anAttrIter->second);
+        double aXYZ[2];
+        for (int i = 0; i < 2; i++) {
+          Slvs_Param aPar = myStorage->getParameter(anEntity.param[i]);
+          aXYZ[i] = aPar.val;
+        }
+        aPoint2D->setValue(aXYZ[0], aXYZ[1]);
+      } else {
+        // Scalar value (used for the distance entities)
+        AttributeDoublePtr aScalar =
+            std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(anAttrIter->first);
+        if (aScalar) {
+          Slvs_Entity anEntity = myStorage->getEntity(anAttrIter->second);
+          Slvs_Param aPar = myStorage->getParameter(anEntity.param[0]);
+          aScalar->setValue(aPar.val);
+        }
+      }
+    }
+  }
+
+  std::map<AttributePtr, Slvs_hParam>::iterator aValIter = myValueMap.begin();
+  for (; aValIter != myValueMap.end(); aValIter++) {
+    AttributeDoublePtr aScalar =
+        std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(anAttrIter->first);
+    if (aScalar) {
+      Slvs_Param aPar = myStorage->getParameter(anAttrIter->second);
+      aScalar->setValue(aPar.val);
+    }
+  }
+}
+
+Slvs_hEntity SketchSolver_Constraint::getId(FeaturePtr theFeature) const
+{
+  std::map<FeaturePtr, Slvs_hEntity>::const_iterator aFIter = myFeatureMap.find(theFeature);
+  if (aFIter == myFeatureMap.end())
+    return SLVS_E_UNKNOWN;
+  return aFIter->second;
+}
+
+Slvs_hEntity SketchSolver_Constraint::getId(AttributePtr theAttribute) const
+{
+  std::map<AttributePtr, Slvs_hEntity>::const_iterator anAttrIter = myAttributeMap.find(theAttribute);
+  if (anAttrIter == myAttributeMap.end())
+    return SLVS_E_UNKNOWN;
+  return anAttrIter->second;
 }