]> SALOME platform Git repositories - modules/shaper.git/commitdiff
Salome HOME
Merge remote-tracking branch 'remotes/origin/SolveSpace'
authorazv <azv@opencascade.com>
Fri, 16 May 2014 04:19:27 +0000 (08:19 +0400)
committerazv <azv@opencascade.com>
Fri, 16 May 2014 04:19:27 +0000 (08:19 +0400)
src/PartSet/PartSet_OperationSketchLine.cpp
src/SketchSolver/SketchSolver_ConstraintManager.cpp
src/SketchSolver/SketchSolver_ConstraintManager.h
src/SketchSolver/SketchSolver_Solver.cpp
src/SketchSolver/SketchSolver_Solver.h

index 97fd78c7458142d2d08ae6f017575c402ab72b4c..ba194d6ba7aa1c03f324d28f31ba363f49e29fac 100644 (file)
@@ -229,6 +229,12 @@ void PartSet_OperationSketchLine::createConstraint(boost::shared_ptr<GeomDataAPI
   boost::shared_ptr<ModelAPI_Document> aDoc = document();
   boost::shared_ptr<ModelAPI_Feature> aFeature = aDoc->addFeature("SketchConstraintCoincidence");
 
+  if (mySketch) {
+    boost::shared_ptr<SketchPlugin_Feature> aSketch = 
+                           boost::dynamic_pointer_cast<SketchPlugin_Feature>(mySketch);
+    aSketch->addSub(aFeature);
+  }
+
   boost::shared_ptr<ModelAPI_Data> aData = aFeature->data();
 
   boost::shared_ptr<ModelAPI_AttributeRefAttr> aRef1 =
index bc33e9e5cee91d58ce223d1c04ece095a345abb2..68bc58b3c9e31a5cb3a6b4f8601410837438c354 100644 (file)
 #include <SketchPlugin_Constraint.h>
 #include <SketchPlugin_ConstraintCoincidence.h>
 #include <SketchPlugin_Line.h>
+#include <SketchPlugin_Point.h>
 #include <SketchPlugin_Sketch.h>
 
+#include <math.h>
+
+/// Tolerance for value of parameters
+const double tolerance = 1.e-10;
+
+// Initialization of constraint manager self pointer
 SketchSolver_ConstraintManager* SketchSolver_ConstraintManager::_self = 0;
 
 /// Global constaint manager object
@@ -25,17 +32,6 @@ SketchSolver_ConstraintManager* myManager = SketchSolver_ConstraintManager::Inst
 /// This value is used to give unique index to the groups
 static Slvs_hGroup myGroupIndexer = 0;
 
-/** \brief Makes transformation from ModelAPI_Attribute to the list of parameters' values
- *  \remark Convertion of normal in 3D needs two attributes (coordinate axis of transversal plane)
- *  \param[in,out] theParams        list of converted values which appended to incoming list
- *  \param[in]     theAttr          attribute to be converted
- *  \param[in]     theNormExtraAttr additional attribute for conversion of a normal
- */
-static void ConvertAttributeToParamList(
-        std::vector<double>&                  theParams, 
-        boost::shared_ptr<ModelAPI_Attribute> theAttr, 
-        boost::shared_ptr<ModelAPI_Attribute> theNormExtraAttr = boost::shared_ptr<ModelAPI_Attribute>());
-
 /** \brief Search the entity/parameter with specified ID in the list of elements
  *  \param[in] theEntityID unique ID of the element
  *  \param[in] theEntities list of elements
@@ -45,6 +41,7 @@ template <typename T>
 static int Search(const uint32_t& theEntityID, const std::vector<T>& theEntities);
 
 
+
 // ========================================================
 // ========= SketchSolver_ConstraintManager ===============
 // ========================================================
@@ -72,95 +69,77 @@ SketchSolver_ConstraintManager::~SketchSolver_ConstraintManager()
 
 void SketchSolver_ConstraintManager::processEvent(const Events_Message* theMessage)
 {
-  if (theMessage->eventID() == Events_Loop::loop()->eventByName(EVENT_FEATURE_CREATED))
-  {
-    const Model_FeatureUpdatedMessage* aCreateMsg = dynamic_cast<const Model_FeatureUpdatedMessage*>(theMessage);
-
-    // Only sketches and constraints can be added by Create event
-    boost::shared_ptr<SketchPlugin_Sketch> aSketch = 
-      boost::dynamic_pointer_cast<SketchPlugin_Sketch>(aCreateMsg->feature());
-    if (aSketch)
-    {
-      addWorkplane(aSketch);
-      return ;
-    }
-    boost::shared_ptr<SketchPlugin_Constraint> aConstraint = 
-      boost::dynamic_pointer_cast<SketchPlugin_Constraint>(aCreateMsg->feature());
-    if (aConstraint)
-    {
-      addConstraint(aConstraint);
-      return ;
-    }
-  }
-  else if (theMessage->eventID() == Events_Loop::loop()->eventByName(EVENT_FEATURE_DELETED))
-  {
-    const Model_FeatureDeletedMessage* aDeleteMsg = dynamic_cast<const Model_FeatureDeletedMessage*>(theMessage);
-    /// \todo Implement deleting objects on event
-  }
-  else if (theMessage->eventID() == Events_Loop::loop()->eventByName(EVENT_FEATURE_UPDATED))
+  if (theMessage->eventID() == Events_Loop::loop()->eventByName(EVENT_FEATURE_CREATED) ||
+      theMessage->eventID() == Events_Loop::loop()->eventByName(EVENT_FEATURE_UPDATED))
   {
     const Model_FeatureUpdatedMessage* aUpdateMsg = dynamic_cast<const Model_FeatureUpdatedMessage*>(theMessage);
 
+    // Only sketches and constraints can be added by Create event
     boost::shared_ptr<SketchPlugin_Sketch> aSketch = 
       boost::dynamic_pointer_cast<SketchPlugin_Sketch>(aUpdateMsg->feature());
     if (aSketch)
     {
-      updateWorkplane(aSketch);
+      changeWorkplane(aSketch);
       return ;
     }
-
     boost::shared_ptr<SketchPlugin_Constraint> aConstraint = 
       boost::dynamic_pointer_cast<SketchPlugin_Constraint>(aUpdateMsg->feature());
     if (aConstraint)
     {
-//      updateConstraint(aConstraint);
+      changeConstraint(aConstraint);
+
+      // Solve the set of constraints
+      ResolveConstraints();
       return ;
     }
 
+    /// \todo Implement feature update handling
     boost::shared_ptr<SketchPlugin_Feature> aFeature = 
       boost::dynamic_pointer_cast<SketchPlugin_Feature>(aUpdateMsg->feature());
-//    if (aFeature)
-//      updateEntity(aFeature);
-  }
-}
-
+    if (aFeature)
+    {
+      updateEntity(aFeature);
 
-bool SketchSolver_ConstraintManager::addWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch)
-{
-  // Check the specified workplane is already used
-  std::vector<SketchSolver_ConstraintGroup>::const_iterator aGroupIter;
-  for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
-    if (aGroupIter->isBaseWorkplane(theSketch))
-      return true;
-  // Create new group for workplane
-  SketchSolver_ConstraintGroup aNewGroup(theSketch);
-  // Verify that the group is created successfully
-  if (!aNewGroup.isBaseWorkplane(theSketch))
-    return false;
-  myGroups.push_back(aNewGroup);
-  return true;
+      // Solve the set of constraints
+      ResolveConstraints();
+    }
+  }
+  else if (theMessage->eventID() == Events_Loop::loop()->eventByName(EVENT_FEATURE_DELETED))
+  {
+    const Model_FeatureDeletedMessage* aDeleteMsg = dynamic_cast<const Model_FeatureDeletedMessage*>(theMessage);
+    /// \todo Implement deleting objects on event
+  }
 }
 
-bool SketchSolver_ConstraintManager::updateWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch)
+bool SketchSolver_ConstraintManager::changeWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch)
 {
   bool aResult = true; // changed when a workplane wrongly updated
   bool isUpdated = false;
   // Try to update specified workplane in all groups
-  std::vector<SketchSolver_ConstraintGroup>::iterator aGroupIter;
+  std::vector<SketchSolver_ConstraintGroup*>::iterator aGroupIter;
   for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
-    if (aGroupIter->isBaseWorkplane(theSketch))
+    if ((*aGroupIter)->isBaseWorkplane(theSketch))
     {
       isUpdated = true;
-      if (!aGroupIter->updateWorkplane(theSketch))
+      if (!(*aGroupIter)->updateWorkplane())
         aResult = false;
     }
-  // If the workplane is not updates, so this is a new workplane
+  // If the workplane is not updated, so this is a new workplane
   if (!isUpdated)
-    return addWorkplane(theSketch);
+  {
+    SketchSolver_ConstraintGroup* aNewGroup = new SketchSolver_ConstraintGroup(theSketch);
+    // Verify that the group is created successfully
+    if (!aNewGroup->isBaseWorkplane(theSketch))
+    {
+      delete aNewGroup;
+      return false;
+    }
+    myGroups.push_back(aNewGroup);
+  }
   return aResult;
 }
 
-bool SketchSolver_ConstraintManager::addConstraint(
+bool SketchSolver_ConstraintManager::changeConstraint(
               boost::shared_ptr<SketchPlugin_Constraint> theConstraint)
 {
   // Search the groups which this constraint touchs
@@ -172,18 +151,22 @@ bool SketchSolver_ConstraintManager::addConstraint(
   { // There are no groups applicable for this constraint => create new one
     boost::shared_ptr<SketchPlugin_Sketch> aWP = findWorkplaneForConstraint(theConstraint);
     if (!aWP) return false;
-    SketchSolver_ConstraintGroup aGroup(aWP);
-    aGroup.addConstraint(theConstraint);
+    SketchSolver_ConstraintGroup* aGroup = new SketchSolver_ConstraintGroup(aWP);
+    if (!aGroup->changeConstraint(theConstraint))
+    {
+      delete aGroup;
+      return false;
+    }
     myGroups.push_back(aGroup);
     return true;
   }
   else if (aGroups.size() == 1)
   { // Only one group => add constraint into it
     Slvs_hGroup aGroupId = *(aGroups.begin());
-    std::vector<SketchSolver_ConstraintGroup>::iterator aGroupIter;
+    std::vector<SketchSolver_ConstraintGroup*>::iterator aGroupIter;
     for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
-      if (aGroupIter->getId() == aGroupId)
-        return aGroupIter->addConstraint(theConstraint);
+      if ((*aGroupIter)->getId() == aGroupId)
+        return (*aGroupIter)->changeConstraint(theConstraint);
   }
   else if (aGroups.size() > 1)
   { // Several groups applicable for this constraint => need to merge them
@@ -194,24 +177,59 @@ bool SketchSolver_ConstraintManager::addConstraint(
   return false;
 }
 
+void SketchSolver_ConstraintManager::updateEntity(boost::shared_ptr<SketchPlugin_Feature> theFeature)
+{
+  // Create list of attributes depending on type of the feature
+  std::vector<std::string> anAttrList;
+  // Point
+  boost::shared_ptr<SketchPlugin_Point> aPoint = 
+    boost::dynamic_pointer_cast<SketchPlugin_Point>(theFeature);
+  if (aPoint)
+    anAttrList.push_back(POINT_ATTR_COORD);
+  // Line
+  boost::shared_ptr<SketchPlugin_Line> aLine = 
+    boost::dynamic_pointer_cast<SketchPlugin_Line>(theFeature);
+  if (aLine)
+  {
+    anAttrList.push_back(LINE_ATTR_START);
+    anAttrList.push_back(LINE_ATTR_END);
+  }
+  /// \todo Other types of features should be implemented
+
+  // Check changing of feature's attributes (go through the groups and search usage of the attributes)
+  std::vector<std::string>::const_iterator anAttrIter;
+  for (anAttrIter = anAttrList.begin(); anAttrIter != anAttrList.end(); anAttrIter++)
+  {
+    std::vector<SketchSolver_ConstraintGroup*>::iterator aGroupIter;
+    for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
+    {
+      boost::shared_ptr<ModelAPI_Attribute> anAttribute = 
+        boost::dynamic_pointer_cast<ModelAPI_Attribute>(theFeature->data()->attribute(*anAttrIter));
+      (*aGroupIter)->updateEntityIfPossible(anAttribute);
+    }
+  }
+}
+
 
 void SketchSolver_ConstraintManager::findGroups(
               boost::shared_ptr<SketchPlugin_Constraint> theConstraint, 
               std::vector<Slvs_hGroup>&                  theGroupIDs) const
 {
-  std::vector<SketchSolver_ConstraintGroup>::const_iterator aGroupIter;
+  boost::shared_ptr<SketchPlugin_Sketch> aWP = findWorkplaneForConstraint(theConstraint);
+
+  std::vector<SketchSolver_ConstraintGroup*>::const_iterator aGroupIter;
   for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
-    if (aGroupIter->isInteract(theConstraint))
-      theGroupIDs.push_back(aGroupIter->getId());
+    if (aWP == (*aGroupIter)->getWorkplane() && (*aGroupIter)->isInteract(theConstraint))
+      theGroupIDs.push_back((*aGroupIter)->getId());
 }
 
 boost::shared_ptr<SketchPlugin_Sketch> SketchSolver_ConstraintManager::findWorkplaneForConstraint(
               boost::shared_ptr<SketchPlugin_Constraint> theConstraint) const
 {
-  std::vector<SketchSolver_ConstraintGroup>::const_iterator aGroupIter;
+  std::vector<SketchSolver_ConstraintGroup*>::const_iterator aGroupIter;
   for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
   {
-    boost::shared_ptr<SketchPlugin_Sketch> aWP = aGroupIter->getWorkplane();
+    boost::shared_ptr<SketchPlugin_Sketch> aWP = (*aGroupIter)->getWorkplane();
     boost::shared_ptr<ModelAPI_AttributeRefList> aWPFeatures = 
       boost::dynamic_pointer_cast<ModelAPI_AttributeRefList>(aWP->data()->attribute(SKETCH_ATTR_FEATURES));
     std::list< boost::shared_ptr<ModelAPI_Feature> > aFeaturesList = aWPFeatures->list();
@@ -224,6 +242,13 @@ boost::shared_ptr<SketchPlugin_Sketch> SketchSolver_ConstraintManager::findWorkp
   return boost::shared_ptr<SketchPlugin_Sketch>();
 }
 
+void SketchSolver_ConstraintManager::ResolveConstraints()
+{
+  std::vector<SketchSolver_ConstraintGroup*>::iterator aGroupIter; 
+  for (aGroupIter = myGroups.begin(); aGroupIter != myGroups.end(); aGroupIter++)
+    (*aGroupIter)->ResolveConstraints();
+}
+
 
 
 // ========================================================
@@ -236,20 +261,16 @@ SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::
     myParamMaxID(0),
     myEntityMaxID(0),
     myConstrMaxID(0),
-    myConstraintMap()
+    myConstraintMap(),
+    myNeedToSolve(false),
+    myConstrSolver()
 {
   myParams.clear();
   myEntities.clear();
   myConstraints.clear();
 
-  // Nullify all elements of the set of equations
-  myConstrSet.param = 0;
-  myConstrSet.entity = 0;
-  myConstrSet.constraint = 0;
-  myConstrSet.failed = 0;
-
   // Initialize workplane
-  myWorkplane.h = 0;
+  myWorkplane.h = SLVS_E_UNKNOWN;
   addWorkplane(theWorkplane);
 }
 
@@ -259,15 +280,6 @@ SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::~SketchSolver_Cons
   myEntities.clear();
   myConstraints.clear();
   myConstraintMap.clear();
-
-  if (myConstrSet.param)
-    delete [] myConstrSet.param;
-  if (myConstrSet.entity)
-    delete [] myConstrSet.entity;
-  if (myConstrSet.constraint)
-    delete [] myConstrSet.constraint;
-  if (myConstrSet.failed)
-    delete [] myConstrSet.failed;
 }
 
 bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::isBaseWorkplane(
@@ -279,20 +291,35 @@ bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::isBaseWorkpla
 bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::isInteract(
                 boost::shared_ptr<SketchPlugin_Constraint> theConstraint) const
 {
+  // Check the group is empty
+  if (myWorkplane.h != SLVS_E_UNKNOWN && myConstraints.empty())
+    return true;
+
   /// \todo Should be implemented
   return false;
 }
 
-bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addConstraint(
+bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::changeConstraint(
                 boost::shared_ptr<SketchPlugin_Constraint> theConstraint)
 {
   // There is no workplane yet, something wrong
-  if (myWorkplane.h == 0)
+  if (myWorkplane.h == SLVS_E_UNKNOWN)
     return false;
 
+  // Search this constraint in the current group to update it
+  std::map<boost::shared_ptr<SketchPlugin_Constraint>, Slvs_hConstraint>::const_iterator
+    aConstrMapIter = myConstraintMap.find(theConstraint);
+  std::vector<Slvs_Constraint>::iterator aConstrIter;
+  if (aConstrMapIter != myConstraintMap.end())
+  {
+    int aConstrPos = Search(aConstrMapIter->second, myConstraints);
+    aConstrIter = myConstraints.begin() + aConstrPos;
+  }
+
   // Get constraint type and verify the constraint parameters are correct
   int aConstrType = getConstraintType(theConstraint);
-  if (aConstrType == SLVS_C_UNKNOWN)
+  if (aConstrType == SLVS_C_UNKNOWN || 
+     (aConstrMapIter != myConstraintMap.end() && aConstrIter->type != aConstrType))
     return false;
 
   // Create constraint parameters
@@ -300,31 +327,55 @@ bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addConstraint
   boost::shared_ptr<ModelAPI_AttributeDouble> aDistAttr =
     boost::dynamic_pointer_cast<ModelAPI_AttributeDouble>(theConstraint->data()->attribute(CONSTRAINT_ATTR_VALUE));
   if (aDistAttr)
+  {
     aDistance = aDistAttr->value();
+    if (aConstrMapIter != myConstraintMap.end() && aConstrIter->valA != aDistance)
+    {
+      myNeedToSolve = true;
+      aConstrIter->valA = aDistance;
+    }
+  }
 
   Slvs_hEntity aConstrEnt[CONSTRAINT_ATTR_SIZE]; // parameters of the constraint
   for (unsigned int indAttr = 0; indAttr < CONSTRAINT_ATTR_SIZE; indAttr++)
   {
+    aConstrEnt[indAttr] = SLVS_E_UNKNOWN;
     boost::shared_ptr<ModelAPI_AttributeRefAttr> aConstrAttr = 
       boost::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
         theConstraint->data()->attribute(CONSTRAINT_ATTRIBUTES[indAttr])
       );
-    aConstrEnt[indAttr] = addEntity(aConstrAttr->attr());
+    if (!aConstrAttr) continue;
+    aConstrEnt[indAttr] = changeEntity(aConstrAttr->attr());
   }
 
-  // Create SolveSpace constraint structure
-  Slvs_Constraint aConstraint = 
-    Slvs_MakeConstraint(++myConstrMaxID, myID, aConstrType, myWorkplane.h, 
-                        aDistance, aConstrEnt[0], aConstrEnt[1], aConstrEnt[2], aConstrEnt[3]);
-  myConstraints.push_back(aConstraint);
-  myConstraintMap[theConstraint] = *(myConstraints.rbegin());
-
+  if (aConstrMapIter == myConstraintMap.end())
+  {
+    // Create SolveSpace constraint structure
+    Slvs_Constraint aConstraint = 
+      Slvs_MakeConstraint(++myConstrMaxID, myID, aConstrType, myWorkplane.h, 
+                          aDistance, aConstrEnt[0], aConstrEnt[1], aConstrEnt[2], aConstrEnt[3]);
+    myConstraints.push_back(aConstraint);
+    myConstraintMap[theConstraint] = aConstraint.h;
+  }
   return true;
 }
 
-Slvs_hEntity SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addEntity(
+Slvs_hEntity SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::changeEntity(
                 boost::shared_ptr<ModelAPI_Attribute> theEntity)
 {
+  // If the entity is already in the group, try to find it
+  std::map<boost::shared_ptr<ModelAPI_Attribute>, Slvs_hEntity>::const_iterator
+    aEntIter = myEntityMap.find(theEntity);
+  std::vector<Slvs_Param>::const_iterator aParamIter; // looks at first parameter of already existent entity or at the end of vector otherwise
+  if (aEntIter == myEntityMap.end()) // no such entity => should be created
+    aParamIter = myParams.end();
+  else
+  { // the entity already exists
+    int aEntPos = Search(aEntIter->second, myEntities);
+    int aParamPos = Search(myEntities[aEntPos].param[0], myParams);
+    aParamIter = myParams.begin() + aParamPos;
+  }
+
   // Look over supported types of entities
 
   // Point in 3D
@@ -332,11 +383,17 @@ Slvs_hEntity SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addEn
     boost::dynamic_pointer_cast<GeomDataAPI_Point>(theEntity);
   if (aPoint)
   {
-    Slvs_hParam aX = addParameter(aPoint->x());
-    Slvs_hParam aY = addParameter(aPoint->y());
-    Slvs_hParam aZ = addParameter(aPoint->z());
+    Slvs_hParam aX = changeParameter(aPoint->x(), aParamIter);
+    Slvs_hParam aY = changeParameter(aPoint->y(), aParamIter);
+    Slvs_hParam aZ = changeParameter(aPoint->z(), aParamIter);
+
+    if (aEntIter != myEntityMap.end()) // the entity already exists
+      return aEntIter->second;
+
+    // New entity
     Slvs_Entity aPtEntity = Slvs_MakePoint3d(++myEntityMaxID, myID, aX, aY, aZ);
     myEntities.push_back(aPtEntity);
+    myEntityMap[theEntity] = aPtEntity.h;
     return aPtEntity.h;
   }
 
@@ -346,36 +403,74 @@ Slvs_hEntity SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addEn
   if (aPoint2D)
   {
     // The 2D points are created on workplane. So, if there is no workplane yet, then error
-    if (myWorkplane.h == 0)
-      return 0;
-    Slvs_hParam aU = addParameter(aPoint2D->x());
-    Slvs_hParam aV = addParameter(aPoint2D->y());
+    if (myWorkplane.h == SLVS_E_UNKNOWN)
+      return SLVS_E_UNKNOWN;
+    Slvs_hParam aU = changeParameter(aPoint2D->x(), aParamIter);
+    Slvs_hParam aV = changeParameter(aPoint2D->y(), aParamIter);
+
+    if (aEntIter != myEntityMap.end()) // the entity already exists
+      return aEntIter->second;
+
+    // New entity
     Slvs_Entity aPt2DEntity = Slvs_MakePoint2d(++myEntityMaxID, myID, myWorkplane.h, aU, aV);
     myEntities.push_back(aPt2DEntity);
+    myEntityMap[theEntity] = aPt2DEntity.h;
     return aPt2DEntity.h;
   }
 
   /// \todo Other types of entities
   
   // Unsupported or wrong entity type
-  return 0;
+  return SLVS_E_UNKNOWN;
 }
 
-Slvs_hEntity SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addNormal(
+Slvs_hEntity SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::changeNormal(
                 boost::shared_ptr<ModelAPI_Attribute> theDirX, 
-                boost::shared_ptr<ModelAPI_Attribute> theDirY)
+                boost::shared_ptr<ModelAPI_Attribute> theDirY, 
+                boost::shared_ptr<ModelAPI_Attribute> theNorm)
 {
-  // Convert axes to the coordinates of normal
-  std::vector<double> aNormCoord;
-  ConvertAttributeToParamList(aNormCoord, theDirX, theDirY);
-  
-  // Create a normal
+  boost::shared_ptr<GeomDataAPI_Dir> aDirX = 
+    boost::dynamic_pointer_cast<GeomDataAPI_Dir>(theDirX);
+  boost::shared_ptr<GeomDataAPI_Dir> aDirY = 
+    boost::dynamic_pointer_cast<GeomDataAPI_Dir>(theDirY);
+  if (!aDirX || !aDirY || 
+     (fabs(aDirX->x()) + fabs(aDirX->y()) + fabs(aDirX->z()) < tolerance) ||
+     (fabs(aDirY->x()) + fabs(aDirY->y()) + fabs(aDirY->z()) < tolerance))
+    return SLVS_E_UNKNOWN;
+
+  // quaternion parameters of normal vector
+  double qw, qx, qy, qz;
+  Slvs_MakeQuaternion(aDirX->x(), aDirX->y(), aDirX->z(), 
+                      aDirY->x(), aDirY->y(), aDirY->z(), 
+                      &qw, &qx, &qy, &qz);
+  double aNormCoord[4] = {qw, qx, qy, qz};
+
+  // Try to find existent normal
+  std::map<boost::shared_ptr<ModelAPI_Attribute>, Slvs_hEntity>::const_iterator
+    aEntIter = myEntityMap.find(theNorm);
+  std::vector<Slvs_Param>::const_iterator aParamIter; // looks at first parameter of already existent entity or at the end of vector otherwise
+  if (aEntIter == myEntityMap.end()) // no such entity => should be created
+    aParamIter = myParams.end();
+  else
+  { // the entity already exists, update it
+    int aEntPos = Search(aEntIter->second, myEntities);
+    int aParamPos = Search(myEntities[aEntPos].param[0], myParams);
+    aParamIter = myParams.begin() + aParamPos;
+  }
+
+  // Change parameters of the normal
   Slvs_hParam aNormParams[4];
   for (int i = 0; i < 4; i++)
-    aNormParams[i] = addParameter(aNormCoord[i]);
+    aNormParams[i] = changeParameter(aNormCoord[i], aParamIter);
+
+  if (aEntIter != myEntityMap.end()) // the entity already exists
+    return aEntIter->second;
+
+  // Create a normal
   Slvs_Entity aNormal = Slvs_MakeNormal3d(++myEntityMaxID, myID, 
                 aNormParams[0], aNormParams[1], aNormParams[2], aNormParams[3]);
   myEntities.push_back(aNormal);
+  myEntityMap[theNorm] = aNormal.h;
   return aNormal.h;
 }
 
@@ -386,74 +481,54 @@ bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addWorkplane(
   if (myWorkplane.h)
     return false; // the workplane already exists
 
-  // Get parameters of workplane
-  boost::shared_ptr<ModelAPI_Attribute> aDirX    = theSketch->data()->attribute(SKETCH_ATTR_DIRX);
-  boost::shared_ptr<ModelAPI_Attribute> aDirY    = theSketch->data()->attribute(SKETCH_ATTR_DIRY);
-  boost::shared_ptr<ModelAPI_Attribute> anOrigin = theSketch->data()->attribute(SKETCH_ATTR_ORIGIN);
-  // Transform them into SolveSpace format
-  Slvs_hEntity aNormalWP = addNormal(aDirX, aDirY);
-  if (!aNormalWP) return false;
-  Slvs_hEntity anOriginWP = addEntity(anOrigin);
-  if (!anOriginWP) return false;
-  // Create workplane
-  myWorkplane = Slvs_MakeWorkplane(++myEntityMaxID, myID, anOriginWP, aNormalWP);
   mySketch = theSketch;
-  // Workplane should be added to the list of entities
-  myEntities.push_back(myWorkplane);
-  return true;
+  return updateWorkplane();
 }
 
-bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::updateWorkplane(
-                boost::shared_ptr<SketchPlugin_Sketch> theSketch)
+bool SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::updateWorkplane()
 {
-  // Renew Sketch pointer
-  mySketch = theSketch;
-
   // Get parameters of workplane
-  boost::shared_ptr<ModelAPI_Attribute> aDirX  = theSketch->data()->attribute(SKETCH_ATTR_DIRX);
-  boost::shared_ptr<ModelAPI_Attribute> aDirY  = theSketch->data()->attribute(SKETCH_ATTR_DIRY);
-  boost::shared_ptr<ModelAPI_Attribute> anOrig = theSketch->data()->attribute(SKETCH_ATTR_ORIGIN);
-  // Transform them to lists of coordinates
-  std::vector<double> aNormal;
-  ConvertAttributeToParamList(aNormal, aDirX, aDirY);
-  std::vector<double> anOrigin;
-  ConvertAttributeToParamList(anOrigin, anOrig);
-
-  // Search the normal and the origin in the parameters list and update their values.
-  // Remark: entities are sorted in the vector, so the most probable position 
-  //         of the entity is equal to identifier the entity
-  
-  // search normal
-  int aEntPos = Search(myWorkplane.normal, myEntities);
-  if (aEntPos < 0) return false;
-  // search first parameter of normal
-  int aParamPos = Search(myEntities[aEntPos].param[0], myParams);
-  if (aParamPos < 0) return false;
-  std::vector<Slvs_Param>::iterator aParamIter = myParams.begin() + aParamPos;
-  // change normal parameters
-  std::vector<double>::iterator anIter;
-  for (anIter = aNormal.begin(); anIter != aNormal.end(); anIter++, aParamIter++)
-    aParamIter->val = *anIter;
-
-  // search origin
-  aEntPos = Search(myWorkplane.point[0], myEntities);
-  if (aEntPos < 0) return false;
-  // search first parameter of origin
-  aParamPos = Search(myEntities[aEntPos].param[0], myParams);
-  if (aParamPos < 0) return false;
-  aParamIter = myParams.begin() + aParamPos;
-  // change origin's parameters
-  for (anIter = anOrigin.begin(); anIter != anOrigin.end(); anIter++, aParamIter++)
-    aParamIter->val = *anIter;
+  boost::shared_ptr<ModelAPI_Attribute> aDirX    = mySketch->data()->attribute(SKETCH_ATTR_DIRX);
+  boost::shared_ptr<ModelAPI_Attribute> aDirY    = mySketch->data()->attribute(SKETCH_ATTR_DIRY);
+  boost::shared_ptr<ModelAPI_Attribute> aNorm    = mySketch->data()->attribute(SKETCH_ATTR_NORM);
+  boost::shared_ptr<ModelAPI_Attribute> anOrigin = mySketch->data()->attribute(SKETCH_ATTR_ORIGIN);
+  // Transform them into SolveSpace format
+  Slvs_hEntity aNormalWP = changeNormal(aDirX, aDirY, aNorm);
+  if (!aNormalWP) return false;
+  Slvs_hEntity anOriginWP = changeEntity(anOrigin);
+  if (!anOriginWP) return false;
 
+  if (!myWorkplane.h)
+  {
+    // Create workplane
+    myWorkplane = Slvs_MakeWorkplane(++myEntityMaxID, myID, anOriginWP, aNormalWP);
+    // Workplane should be added to the list of entities
+    myEntities.push_back(myWorkplane);
+  }
   return true;
 }
 
 
-Slvs_hParam SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::addParameter(double theParam)
+Slvs_hParam SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::changeParameter(
+                const double&                            theParam, 
+                std::vector<Slvs_Param>::const_iterator& thePrmIter)
 {
+  if (thePrmIter != myParams.end())
+  { // Parameter should be updated
+    if (thePrmIter->val != theParam)
+      myNeedToSolve = true; // parameter is changed, need to resolve constraints
+    int aParamPos = thePrmIter - myParams.begin();
+    myParams[aParamPos].val = theParam;
+    thePrmIter++;
+    return myParams[aParamPos].h;
+  }
+
+  // Newly created parameter
   Slvs_Param aParam = Slvs_MakeParam(++myParamMaxID, myID, theParam);
   myParams.push_back(aParam);
+  myNeedToSolve = true;
+  // The list of parameters is changed, move iterator to the end of the list to avoid problems
+  thePrmIter = myParams.end();
   return aParam.h;
 }
 
@@ -474,6 +549,7 @@ int SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::getConstraintT
         boost::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
           theConstraint->data()->attribute(CONSTRAINT_ATTRIBUTES[indAttr])
         );
+      if (!anAttr) continue;
       // Verify the attribute is a 2D point
       boost::shared_ptr<GeomDataAPI_Point2D> aPoint2D = 
         boost::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttr->attr());
@@ -506,59 +582,80 @@ int SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::getConstraintT
   return SLVS_C_UNKNOWN;
 }
 
+void SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::ResolveConstraints()
+{
+  if (!myNeedToSolve)
+    return;
+
+  myConstrSolver.setGroupID(myID);
+  myConstrSolver.setParameters(myParams);
+  myConstrSolver.setEntities(myEntities);
+  myConstrSolver.setConstraints(myConstraints);
+
+  if (myConstrSolver.solve() == SLVS_RESULT_OKAY)
+  { // solution succeeded, store results into correspondent attributes
+    // Obtain result into the same list of parameters
+    if (!myConstrSolver.getResult(myParams))
+      return;
+
+    std::map<boost::shared_ptr<ModelAPI_Attribute>, Slvs_hEntity>::iterator
+      anEntIter = myEntityMap.begin();
+    for ( ; anEntIter != myEntityMap.end(); anEntIter++)
+      updateAttribute(anEntIter->first, anEntIter->second);
+  }
+  /// \todo Implement error handling
 
+  myNeedToSolve = false;
+}
 
-// ========================================================
-// =========      Auxiliary functions       ===============
-// ========================================================
 
-void ConvertAttributeToParamList(
-                std::vector<double>&                  theParams, 
-                boost::shared_ptr<ModelAPI_Attribute> theAttr, 
-                boost::shared_ptr<ModelAPI_Attribute> theNormExtraAttr)
+void SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::updateAttribute(
+                boost::shared_ptr<ModelAPI_Attribute> theAttribute, 
+                const Slvs_hEntity&                   theEntityID)
 {
+  // Search the position of the first parameter of the entity
+  int anEntPos = Search(theEntityID, myEntities);
+  int aFirstParamPos = Search(myEntities[anEntPos].param[0], myParams);
+
+  // Look over supported types of entities
+
   // Point in 3D
   boost::shared_ptr<GeomDataAPI_Point> aPoint = 
-    boost::dynamic_pointer_cast<GeomDataAPI_Point>(theAttr);
+    boost::dynamic_pointer_cast<GeomDataAPI_Point>(theAttribute);
   if (aPoint)
   {
-    theParams.push_back(aPoint->x());
-    theParams.push_back(aPoint->y());
-    theParams.push_back(aPoint->z());
+    aPoint->setValue(myParams[aFirstParamPos].val,
+                     myParams[aFirstParamPos+1].val,
+                     myParams[aFirstParamPos+2].val);
     return ;
   }
 
   // Point in 2D
   boost::shared_ptr<GeomDataAPI_Point2D> aPoint2D = 
-    boost::dynamic_pointer_cast<GeomDataAPI_Point2D>(theAttr);
+    boost::dynamic_pointer_cast<GeomDataAPI_Point2D>(theAttribute);
   if (aPoint2D)
   {
-    theParams.push_back(aPoint2D->x());
-    theParams.push_back(aPoint2D->y());
+    aPoint2D->setValue(myParams[aFirstParamPos].val,
+                       myParams[aFirstParamPos+1].val);
     return ;
   }
 
-  // Normal in 3D
-  boost::shared_ptr<GeomDataAPI_Dir> aDirX = 
-    boost::dynamic_pointer_cast<GeomDataAPI_Dir>(theAttr);
-  boost::shared_ptr<GeomDataAPI_Dir> aDirY = 
-    boost::dynamic_pointer_cast<GeomDataAPI_Dir>(theNormExtraAttr);
-  if (aDirX && aDirY)
-  {
-    // quaternion parameters of normal vector
-    double qw, qx, qy, qz;
-    Slvs_MakeQuaternion(aDirX->x(), aDirX->y(), aDirX->z(), 
-                        aDirY->x(), aDirY->y(), aDirY->z(), 
-                        &qw, &qx, &qy, &qz);
-    theParams.push_back(qw);
-    theParams.push_back(qx);
-    theParams.push_back(qy);
-    theParams.push_back(qz);
-  }
+  /// \todo Support other types of entities
+}
 
-  /// \todo Other types of entities
+void SketchSolver_ConstraintManager::SketchSolver_ConstraintGroup::updateEntityIfPossible(
+                boost::shared_ptr<ModelAPI_Attribute> theEntity)
+{
+  if (myEntityMap.find(theEntity) != myEntityMap.end())
+    changeEntity(theEntity);
 }
 
+
+
+// ========================================================
+// =========      Auxiliary functions       ===============
+// ========================================================
+
 template <typename T>
 int Search(const uint32_t& theEntityID, const std::vector<T>& theEntities)
 {
@@ -571,4 +668,3 @@ int Search(const uint32_t& theEntityID, const std::vector<T>& theEntities)
     return -1;
   return aEntIter - theEntities.begin();
 }
-
index e44cb6444b040a1b3b6e1b789c903e3043b1bea2..95512cf3bcd2d30d4c49557589eaebb65ae97338 100644 (file)
@@ -6,14 +6,11 @@
 #define SketchSolver_ConstraintManager_Headerfile
 
 #include "SketchSolver.h"
+#include <SketchSolver_Solver.h>
 
 #include <Events_Listener.h>
 #include <SketchPlugin_Constraint.h>
 
-// Need to be defined before including SolveSpace to avoid additional dependances on Windows platform
-#if defined(WIN32) && !defined(HAVE_C99_INTEGER_TYPES)
-typedef unsigned int UINT32;
-#endif
 #include <string.h>
 #include <slvs.h>
 
@@ -23,6 +20,8 @@ typedef unsigned int UINT32;
 
 // Unknown constraint (for error reporting)
 #define SLVS_C_UNKNOWN 0
+// Unknown entity
+#define SLVS_E_UNKNOWN 0
 
 /** \class   SketchSolver_ConstraintManager
  *  \ingroup DataModel
@@ -52,11 +51,11 @@ protected:
   SketchSolver_ConstraintManager();
   ~SketchSolver_ConstraintManager();
 
-  /** \brief Adds a constraint into the manager
-   *  \param[in] theConstraint constraint to be added
-   *  \return \c true if the constraint added successfully
+  /** \brief Adds or updates a constraint in the suitable group
+   *  \param[in] theConstraint constraint to be changed
+   *  \return \c true if the constraint changed successfully
    */
-  bool addConstraint(boost::shared_ptr<SketchPlugin_Constraint> theConstraint);
+  bool changeConstraint(boost::shared_ptr<SketchPlugin_Constraint> theConstraint);
 
   /** \brief Removes a constraint from the manager
    *  \param[in] theConstraint constraint to be removed
@@ -64,17 +63,11 @@ protected:
    */
   bool removeConstraint(boost::shared_ptr<SketchPlugin_Constraint> theConstraint);
 
-  /** \brief Updates a constraint
-   *  \param[in] theConstraint constraint to be updated
-   *  \return \c true if the constraint was updated
-   */
-  bool updateConstraint(boost::shared_ptr<SketchPlugin_Constraint> theConstraint);
-
-  /** \brief Adds a workplane into the manager
-   *  \param[in] theSketch the feature to create workplane
-   *  \return \c true if the workplane added successfully
+  /** \brief Adds or updates a workplane in the manager
+   *  \param[in] theSketch the feature to create or update workplane
+   *  \return \c true if the workplane cahnged successfully
    */
-  bool addWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch);
+  bool changeWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch);
 
   /** \brief Removes a workplane from the manager. 
    *         All groups based on such workplane will be removed too.
@@ -83,17 +76,14 @@ protected:
    */
   bool removeWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch);
 
-  /** \brief Updates a workplane
-   *  \param[in] theSketch workplane to be updated
-   *  \return \c true if the workplane was updated
-   */
-  bool updateWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch);
-
   /** \brief Updates entity which is neither workplane nor constraint
    *  \param[in] theFeature entity to be updated
-   *  \return \c true if the entity updated successfully
    */
-  bool updateEntity(boost::shared_ptr<SketchPlugin_Feature> theFeature);
+  void updateEntity(boost::shared_ptr<SketchPlugin_Feature> theFeature);
+
+  /** \brief Goes through the list of groups and solve the constraints
+   */
+  void ResolveConstraints();
 
 private:
   class SketchSolver_ConstraintGroup;
@@ -114,7 +104,7 @@ private:
 
 private:
   static SketchSolver_ConstraintManager*    _self;    ///< Self pointer to implement singleton functionality
-  std::vector<SketchSolver_ConstraintGroup> myGroups; ///< Groups of constraints
+  std::vector<SketchSolver_ConstraintGroup*> myGroups; ///< Groups of constraints
 };
 
 
@@ -135,11 +125,11 @@ public:
   const Slvs_hGroup& getId() const
   {return myID;}
 
-  /** \brief Adds a constraint into the group
-   *  \param[in] theConstraint constraint to be added
-   *  \return \c true if the constraint added successfully
+  /** \brief Adds or updates a constraint in the group
+   *  \param[in] theConstraint constraint to be changed
+   *  \return \c true if the constraint added or updated successfully
    */
-  bool addConstraint(boost::shared_ptr<SketchPlugin_Constraint> theConstraint);
+  bool changeConstraint(boost::shared_ptr<SketchPlugin_Constraint> theConstraint);
 
   /** \brief Removes a constraint into the group
    *  \param[in] theConstraint constraint to be removed
@@ -162,24 +152,32 @@ public:
   boost::shared_ptr<SketchPlugin_Sketch> getWorkplane() const
   { return mySketch; }
 
-  /** \brief Update parameters of workplane. Should be called when Update event is coming
-   *  \param[in] theWorkplane workplane to be updated
-   *  \return \c true if workplane updated successfully
+  /** \brief Update parameters of workplane. Should be called when Update event is coming.
+   *  \return \c true if workplane updated successfully, \c false if workplane parameters are not consistent
    */
-  bool updateWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theWorkplane);
+  bool updateWorkplane();
+
+  /** \brief If the entity is in this group it will updated
+   *  \param[in] theEntity attribute, which values should update SolveSpace entity
+   */
+  void updateEntityIfPossible(boost::shared_ptr<ModelAPI_Attribute> theEntity);
+
+  /** \brief Start solution procedure if necessary and update attributes of features
+   */
+  void ResolveConstraints();
 
 protected:
-  /** \brief Adds an entity into the group
+  /** \brief Adds or updates an entity in the group
    *
    *  The parameters of entity will be parsed and added to the list of SolveSpace parameters.
    *  Parameters of certain entity will be placed sequentially in the list.
    *
    *  \param[in] theEntity the object of constraint
-   *  \return identifier of created entity or 0 if entity was not added
+   *  \return identifier of changed entity or 0 if entity could not be changed
    */
-  Slvs_hEntity addEntity(boost::shared_ptr<ModelAPI_Attribute> theEntity);
+  Slvs_hEntity changeEntity(boost::shared_ptr<ModelAPI_Attribute> theEntity);
 
-  /** \brief Adds a normal into the group
+  /** \brief Adds or updates a normal in the group
    *
    *  Normal is a special entity in SolveSpace, which defines a direction in 3D and 
    *  a rotation about this direction. So, SolveSpace represents normals as unit quaternions.
@@ -189,16 +187,21 @@ protected:
    *
    *  \param[in] theDirX first coordinate axis of the plane
    *  \param[in] theDirY second coordinate axis of the plane
-   *  \return identifier of created normal
+   *  \param[in] theNorm attribute for the normal (used to identify newly created entity)
+   *  \return identifier of created or updated normal
    */
-  Slvs_hEntity addNormal(boost::shared_ptr<ModelAPI_Attribute> theDirX, 
-                         boost::shared_ptr<ModelAPI_Attribute> theDirY);
+  Slvs_hEntity changeNormal(boost::shared_ptr<ModelAPI_Attribute> theDirX, 
+                            boost::shared_ptr<ModelAPI_Attribute> theDirY, 
+                            boost::shared_ptr<ModelAPI_Attribute> theNorm);
 
-  /** \brief Adds a parameter into the group
-   *  \param[in] theParam parameter to be added
-   *  \return identifier of created parameter or 0 if it was not added
+  /** \brief Adds or updates a parameter in the group
+   *  \param[in] theParam   the value of parameter
+   *  \param[in] thePrmIter the cell in the list of parameters which should be changed 
+   *                        (the iterator will be increased if it does not reach the end of the list)
+   *  \return identifier of changed parameter; when the parameter cannot be created, returned ID is 0
    */
-  Slvs_hParam addParameter(double theParam);
+  Slvs_hParam changeParameter(const double& theParam, 
+                              std::vector<Slvs_Param>::const_iterator& thePrmIter);
 
   /** \brief Compute constraint type according to SolveSpace identifiers
    *         and verify that constraint parameters are correct
@@ -207,10 +210,16 @@ protected:
    */
   int getConstraintType(const boost::shared_ptr<SketchPlugin_Constraint>& theConstraint) const;
 
+  /** \brief Change values of attribute by parameters received from SolveSpace solver
+   *  \param[in,out] theAttribute pointer to the attribute to be changed
+   *  \param[in]     theEntityID  identifier of SolveSpace entity, which contains updated data
+   */
+  void updateAttribute(boost::shared_ptr<ModelAPI_Attribute> theAttribute, const Slvs_hEntity& theEntityID);
+
 private:
   /** \brief Creates a workplane from the sketch parameters
    *  \param[in] theSketch parameters of workplane are the attributes of this sketch
-   *  \return \c true if success
+   *  \return \c true if success, \c false if workplane parameters are not consistent
    */
   bool addWorkplane(boost::shared_ptr<SketchPlugin_Sketch> theSketch);
 
@@ -224,12 +233,16 @@ private:
   Slvs_hEntity                 myEntityMaxID;   ///< Actual maximal ID of entities (not equal to myEntities size)
   std::vector<Slvs_Constraint> myConstraints;   ///< List of constraints in SolveSpace format
   Slvs_hConstraint             myConstrMaxID;   ///< Actual maximal ID of constraints (not equal to myConstraints size)
-  Slvs_System                  myConstrSet;     ///< SolveSpace's set of equations obtained by constraints
+  bool                         myNeedToSolve;   ///< Indicator that something changed in the group and constraint system need to be rebuilt
+
+  SketchSolver_Solver          myConstrSolver;  ///< Solver for set of equations obtained by constraints
 
   // SketchPlugin entities
   boost::shared_ptr<SketchPlugin_Sketch> mySketch; ///< Equivalent to workplane
-  std::map<boost::shared_ptr<SketchPlugin_Constraint>, Slvs_Constraint> 
-                               myConstraintMap; ///< The map between SketchPlugin and SolveSpace constraints
+  std::map<boost::shared_ptr<SketchPlugin_Constraint>, Slvs_hConstraint> 
+                               myConstraintMap;    ///< The map between SketchPlugin and SolveSpace constraints
+  std::map<boost::shared_ptr<ModelAPI_Attribute>, Slvs_hEntity>
+                               myEntityMap;        ///< The map between parameters of constraints and their equivalent SolveSpace entities
 };
 
 #endif
index bf86e4d44f4b9ee6704d05266d556acc92b11312..f9ed8c37204c788a6de8c685faac995fe3b30af9 100644 (file)
@@ -3,3 +3,109 @@
 // Author:  Artem ZHIDKOV
 
 #include "SketchSolver_Solver.h"
+
+SketchSolver_Solver::SketchSolver_Solver()
+{
+  // Nullify all elements of the set of equations
+  myEquationsSystem.param = 0;
+  myEquationsSystem.params = 0;
+  myEquationsSystem.entity = 0;
+  myEquationsSystem.entities = 0;
+  myEquationsSystem.constraint = 0;
+  myEquationsSystem.constraints = 0;
+  myEquationsSystem.failed = 0;
+  myEquationsSystem.faileds = 0;
+
+  // If the set of constraints is inconsistent,
+  // the failed field will contain wrong constraints
+  myEquationsSystem.calculateFaileds = 1;
+}
+
+SketchSolver_Solver::~SketchSolver_Solver()
+{
+  if (myEquationsSystem.param)
+    delete [] myEquationsSystem.param;
+  if (myEquationsSystem.entity)
+    delete [] myEquationsSystem.entity;
+  if (myEquationsSystem.constraint)
+    delete [] myEquationsSystem.constraint;
+  if (myEquationsSystem.failed)
+    delete [] myEquationsSystem.failed;
+}
+
+void SketchSolver_Solver::setParameters(const std::vector<Slvs_Param>& theParameters)
+{
+  if (theParameters.size() != myEquationsSystem.params) // number of parameters was changed => reallocate the memory
+  {
+    if (myEquationsSystem.param)
+      delete [] myEquationsSystem.param;
+    myEquationsSystem.params = theParameters.size();
+    myEquationsSystem.param = new Slvs_Param[theParameters.size()];
+  }
+
+  // Copy data
+  std::vector<Slvs_Param>::const_iterator aParamIter = theParameters.begin();
+  for (int i = 0; i < myEquationsSystem.params; i++, aParamIter++)
+    myEquationsSystem.param[i] = *aParamIter;
+}
+
+void SketchSolver_Solver::setEntities(const std::vector<Slvs_Entity>& theEntities)
+{
+  if (theEntities.size() != myEquationsSystem.entities) // number of entities was changed => reallocate the memory
+  {
+    if (myEquationsSystem.entity)
+      delete [] myEquationsSystem.entity;
+    myEquationsSystem.entities = theEntities.size();
+    myEquationsSystem.entity = new Slvs_Entity[theEntities.size()];
+  }
+
+  // Copy data
+  std::vector<Slvs_Entity>::const_iterator aEntIter = theEntities.begin();
+  for (int i = 0; i < myEquationsSystem.entities; i++, aEntIter++)
+    myEquationsSystem.entity[i] = *aEntIter;
+}
+
+void SketchSolver_Solver::setConstraints(const std::vector<Slvs_Constraint>& theConstraints)
+{
+  if (theConstraints.size() != myEquationsSystem.constraints) // number of constraints was changed => reallocate the memory
+  {
+    if (myEquationsSystem.constraint)
+      delete [] myEquationsSystem.constraint;
+    myEquationsSystem.constraints = theConstraints.size();
+    myEquationsSystem.constraint = new Slvs_Constraint[theConstraints.size()];
+  }
+
+  // Copy data
+  std::vector<Slvs_Constraint>::const_iterator aConstrIter = theConstraints.begin();
+  for (int i = 0; i < myEquationsSystem.constraints; i++, aConstrIter++)
+    myEquationsSystem.constraint[i] = *aConstrIter;
+}
+
+int SketchSolver_Solver::solve()
+{
+  if (myEquationsSystem.constraints <= 0)
+    return SLVS_RESULT_EMPTY_SET;
+
+  Slvs_Solve(&myEquationsSystem, myGroupID);
+
+  return myEquationsSystem.result;
+}
+
+bool SketchSolver_Solver::getResult(std::vector<Slvs_Param>& theParameters)
+{
+  if (myEquationsSystem.result != SLVS_RESULT_OKAY)
+    return false;
+
+  if (theParameters.size() != myEquationsSystem.params)
+    return false; // number of parameters is not the same
+
+  std::vector<Slvs_Param>::iterator aParamIter = theParameters.begin();
+  for (int i = 0; i < myEquationsSystem.params; i++, aParamIter++)
+  {
+    if (myEquationsSystem.param[i].h != aParamIter->h)
+      return false; // sequence of parameters was changed
+    aParamIter->val = myEquationsSystem.param[i].val;
+  }
+
+  return true;
+}
index beb2f249f7d24b321c87bb31435745e3adc74e0e..aeb29ef5e44d5753a5febc9cad8baf999f8e97ea 100644 (file)
@@ -7,4 +7,59 @@
 
 #include "SketchSolver.h"
 
+// Need to be defined before including SolveSpace to avoid additional dependances on Windows platform
+#if defined(WIN32) && !defined(HAVE_C99_INTEGER_TYPES)
+typedef unsigned int UINT32;
+#endif
+#include <string.h>
+#include <slvs.h>
+
+#include <vector>
+
+
+#define SLVS_RESULT_EMPTY_SET -1
+
+
+class SketchSolver_Solver
+{
+public:
+  SketchSolver_Solver();
+  ~SketchSolver_Solver();
+
+  /** \brief Initialize the ID of the group
+   */
+  inline void setGroupID(Slvs_hGroup theGroupID)
+  { myGroupID = theGroupID; }
+
+  /** \brief Change array of parameters
+   *  \param[in] theParameters vector of parameters
+   */
+  void setParameters(const std::vector<Slvs_Param>& theParameters);
+
+  /** \brief Change array of entities
+   *  \param[in] theEntities vector of entities
+   */
+  void setEntities(const std::vector<Slvs_Entity>& theEntities);
+
+  /** \brief Change array of constraints
+   *  \param[in] theConstraints vector of constraints
+   */
+  void setConstraints(const std::vector<Slvs_Constraint>& theConstraints);
+
+  /** \brief Solve the set of equations
+   *  \return identifier whether solution succeeded
+   */
+  int solve();
+
+  /** \brief Updates the list of parameters by calculated values
+   *  \param[in,out] theParameters parameters to be updated
+   *  \return \c true if parameters are updated correctly
+   */
+  bool getResult(std::vector<Slvs_Param>& theParameters);
+
+private:
+  Slvs_hGroup myGroupID;         ///< identifier of the group to be solved
+  Slvs_System myEquationsSystem; ///< set of equations for solving in SolveSpace
+};
+
 #endif