1 #include <SketchSolver_Constraint.h>
2 #include <SketchSolver_Group.h>
3 #include <SketchSolver_Error.h>
5 #include <SketchPlugin_Arc.h>
6 #include <SketchPlugin_Circle.h>
7 #include <SketchPlugin_Line.h>
8 #include <SketchPlugin_Point.h>
10 #include <GeomAPI_Dir2d.h>
11 #include <GeomDataAPI_Point.h>
12 #include <GeomDataAPI_Point2D.h>
13 #include <ModelAPI_AttributeDouble.h>
14 #include <ModelAPI_ResultConstruction.h>
18 SketchSolver_Constraint::SketchSolver_Constraint(
19 ConstraintPtr theConstraint)
20 : myBaseConstraint(theConstraint),
25 SketchSolver_Constraint::~SketchSolver_Constraint()
27 std::map<AttributePtr, Slvs_hParam>::const_iterator anIt1 = myValueMap.begin();
28 for (; anIt1 != myValueMap.end(); anIt1++)
29 myStorage->removeParameter(anIt1->second);
32 std::map<AttributePtr, Slvs_hEntity>::const_iterator anIt2 = myAttributeMap.begin();
33 for (; anIt2 != myAttributeMap.end(); anIt2++)
34 myStorage->removeEntity(anIt2->second);
35 myAttributeMap.clear();
37 std::map<FeaturePtr, Slvs_hEntity>::iterator anIt3 = myFeatureMap.begin();
38 while (!myFeatureMap.empty()) {
39 std::shared_ptr<SketchPlugin_Feature> aFeature =
40 std::dynamic_pointer_cast<SketchPlugin_Feature>(anIt3->first);
41 Slvs_hEntity anEnt = anIt3->second;
42 std::map<FeaturePtr, Slvs_hEntity>::iterator aRemIt = anIt3++;
43 myFeatureMap.erase(aRemIt);
44 if (!myGroup->isInteract(aFeature))
45 myStorage->removeEntity(anEnt);
48 std::vector<Slvs_hConstraint>::const_iterator anIt4 = mySlvsConstraints.begin();
49 for (; anIt4 != mySlvsConstraints.end(); anIt4++)
50 myStorage->removeConstraint(*anIt4);
51 mySlvsConstraints.clear();
54 void SketchSolver_Constraint::setStorage(StoragePtr theStorage)
56 myStorage = theStorage;
60 void SketchSolver_Constraint::setGroup(SketchSolver_Group* theGroup)
66 void SketchSolver_Constraint::addFeature(FeaturePtr theFeature)
69 changeEntity(theFeature, aType);
73 void SketchSolver_Constraint::process()
76 if (!myBaseConstraint || !myStorage || myGroup == 0) {
77 /// TODO: Put error message here
80 if (!mySlvsConstraints.empty()) // some data is changed, update constraint
81 update(myBaseConstraint);
83 int aConstrType = getType();
85 std::vector<Slvs_hEntity> anAttributes;
86 getAttributes(aValue, anAttributes);
87 if (!myErrorMsg.empty())
89 if (aConstrType == SLVS_C_UNKNOWN)
90 aConstrType = getType();
92 Slvs_hGroup aGroupID = myGroup->getId();
93 Slvs_hEntity aWorkplaneID = myGroup->getWorkplaneId();
94 Slvs_Constraint aConstraint;
95 if (mySlvsConstraints.empty())
96 aConstraint = Slvs_MakeConstraint(SLVS_C_UNKNOWN, aGroupID, aConstrType, aWorkplaneID,
97 aValue, anAttributes[0], anAttributes[1], anAttributes[2], anAttributes[3]);
99 aConstraint = myStorage->getConstraint(mySlvsConstraints[0]);
100 aConstraint.valA = aValue;
101 static const int aNbAttrs = 6;
102 Slvs_hEntity* aConstrAttrs[aNbAttrs] = {
103 &aConstraint.ptA, &aConstraint.ptB,
104 &aConstraint.entityA, &aConstraint.entityB,
105 &aConstraint.entityC, &aConstraint.entityD};
106 std::vector<Slvs_hEntity>::const_iterator anIter = anAttributes.begin();
107 for (int i = 0; i < aNbAttrs && anIter != anAttributes.end(); i++, anIter++)
108 *(aConstrAttrs[i]) = *anIter;
111 Slvs_hConstraint anID = myStorage->addConstraint(aConstraint);
112 if (mySlvsConstraints.empty())
113 mySlvsConstraints.push_back(anID);
115 mySlvsConstraints[0] = anID;
119 bool SketchSolver_Constraint::checkAttributesChanged(ConstraintPtr theConstraint)
121 // Check the attrbutes of constraint are changed
122 ConstraintPtr aConstraint = theConstraint ? theConstraint : myBaseConstraint;
123 std::list<AttributePtr> anAttrList = aConstraint->data()->attributes(std::string());
124 std::list<AttributePtr>::iterator anAttrIter = anAttrList.begin();
125 for (; anAttrIter != anAttrList.end(); anAttrIter++) {
126 AttributeRefAttrPtr aRefAttr =
127 std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*anAttrIter);
129 if (aRefAttr->isObject()) {
130 FeaturePtr aFeature = ModelAPI_Feature::feature(aRefAttr->object());
131 if (aFeature && myFeatureMap.find(aFeature) == myFeatureMap.end())
133 } else if (aRefAttr->attr() &&
134 myAttributeMap.find(aRefAttr->attr()) == myAttributeMap.end())
137 AttributeRefListPtr aRefList =
138 std::dynamic_pointer_cast<ModelAPI_AttributeRefList>(*anAttrIter);
140 std::list<ObjectPtr> anItems = aRefList->list();
141 std::list<ObjectPtr>::iterator anIt = anItems.begin();
142 for (; anIt != anItems.end(); anIt++) {
143 FeaturePtr aFeature = ModelAPI_Feature::feature(*anIt);
144 if (aFeature && myFeatureMap.find(aFeature) == myFeatureMap.end())
152 void SketchSolver_Constraint::update(ConstraintPtr theConstraint)
155 bool needToRebuild = (theConstraint && theConstraint != myBaseConstraint);
157 needToRebuild = checkAttributesChanged(theConstraint);
159 if (theConstraint && theConstraint->getKind() != myBaseConstraint->getKind())
161 remove(myBaseConstraint);
163 myBaseConstraint = theConstraint;
168 // Update all attributes
170 std::map<Slvs_hEntity, Slvs_hEntity> aRelocationMap;
171 std::map<FeaturePtr, Slvs_hEntity>::iterator aFeatIter = myFeatureMap.begin();
172 for (; aFeatIter != myFeatureMap.end(); aFeatIter++) {
173 Slvs_hEntity aPrevID = aFeatIter->second;
174 aFeatIter->second = changeEntity(aFeatIter->first, aType);
175 if (aFeatIter->second != aPrevID)
176 aRelocationMap[aPrevID] = aFeatIter->second;
178 std::map<AttributePtr, Slvs_hEntity>::iterator anAttrIter = myAttributeMap.begin();
179 for (; anAttrIter != myAttributeMap.end(); anAttrIter++) {
180 Slvs_hEntity aPrevID = anAttrIter->second;
181 anAttrIter->second = changeEntity(anAttrIter->first, aType);
182 if (anAttrIter->second != aPrevID)
183 aRelocationMap[aPrevID] = anAttrIter->second;
187 DataPtr aData = myBaseConstraint->data();
189 AttributeDoublePtr aValueAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
190 myBaseConstraint->attribute(SketchPlugin_Constraint::VALUE()));
191 double aValue = aValueAttr ? aValueAttr->value() : 0.0;
194 std::vector<Slvs_hConstraint>::iterator aCIter = mySlvsConstraints.begin();
195 for (; aCIter != mySlvsConstraints.end(); aCIter++) {
196 Slvs_Constraint aConstraint = myStorage->getConstraint(*aCIter);
198 aConstraint.valA = aValue;
199 Slvs_hEntity* aCoeffs[6] = {
200 &aConstraint.ptA, &aConstraint.ptB,
201 &aConstraint.entityA, &aConstraint.entityB,
202 &aConstraint.entityC, &aConstraint.entityD};
203 for (int i = 0; i < 6; i++) {
204 if (*(aCoeffs[i]) == SLVS_E_UNKNOWN)
206 std::map<Slvs_hEntity, Slvs_hEntity>::iterator aFound = aRelocationMap.find(*(aCoeffs[i]));
207 if (aFound != aRelocationMap.end())
208 *(aCoeffs[i]) = aFound->second;
210 *aCIter = myStorage->addConstraint(aConstraint);
215 bool SketchSolver_Constraint::remove(ConstraintPtr theConstraint)
218 if (theConstraint && theConstraint != myBaseConstraint)
220 if (mySlvsConstraints.empty())
222 bool isFullyRemoved = myStorage->removeConstraint(mySlvsConstraints.front());
223 if (isFullyRemoved) {
224 myFeatureMap.clear();
225 myAttributeMap.clear();
228 cleanRemovedEntities();
229 mySlvsConstraints.clear();
233 void SketchSolver_Constraint::cleanRemovedEntities()
235 std::set<Slvs_hParam> aRemovedParams;
236 std::set<Slvs_hEntity> aRemovedEntities;
237 std::set<Slvs_hConstraint> aRemovedConstraints;
238 myStorage->getRemoved(aRemovedParams, aRemovedEntities, aRemovedConstraints);
239 std::map<FeaturePtr, Slvs_hEntity>::iterator aFeatIt = myFeatureMap.begin();
240 while (aFeatIt != myFeatureMap.end()) {
241 if (aRemovedEntities.find(aFeatIt->second) == aRemovedEntities.end()) {
245 std::map<FeaturePtr, Slvs_hEntity>::iterator aTmpIter = aFeatIt++;
246 myFeatureMap.erase(aTmpIter);
248 std::map<AttributePtr, Slvs_hEntity>::iterator anAttrIt = myAttributeMap.begin();
249 while (anAttrIt != myAttributeMap.end()) {
250 if (aRemovedEntities.find(anAttrIt->second) == aRemovedEntities.end()) {
254 std::map<AttributePtr, Slvs_hEntity>::iterator aTmpIter = anAttrIt++;
255 myAttributeMap.erase(aTmpIter);
257 std::map<AttributePtr, Slvs_hParam>::iterator aValIt = myValueMap.begin();
258 while (aValIt != myValueMap.end()) {
259 if (aRemovedParams.find(aValIt->second) == aRemovedParams.end()) {
263 std::map<AttributePtr, Slvs_hParam>::iterator aTmpIter = aValIt++;
264 myValueMap.erase(aTmpIter);
268 void SketchSolver_Constraint::getAttributes(
270 std::vector<Slvs_hEntity>& theAttributes)
272 static const int anInitNbOfAttr = 4;
273 theAttributes.assign(anInitNbOfAttr, SLVS_E_UNKNOWN);
275 DataPtr aData = myBaseConstraint->data();
277 AttributeDoublePtr aValueAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
278 aData->attribute(SketchPlugin_Constraint::VALUE()));
279 theValue = aValueAttr ? aValueAttr->value() : 0.0;
281 int aPtInd = 0; // index of first point in the list of attributes
282 int aEntInd = 2; // index of first antity in the list of attributes
283 std::list<AttributePtr> aConstrAttrs = aData->attributes(ModelAPI_AttributeRefAttr::typeId());
284 std::list<AttributePtr>::iterator anIter = aConstrAttrs.begin();
285 for (; anIter != aConstrAttrs.end(); anIter++) {
286 AttributeRefAttrPtr aRefAttr =
287 std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*anIter);
288 if (!aRefAttr || !aRefAttr->isInitialized()) {
289 myErrorMsg = SketchSolver_Error::NOT_INITIALIZED();
293 int aType = SLVS_E_UNKNOWN; // type of created entity
294 Slvs_hEntity anEntity = myGroup->getAttributeId(aRefAttr);
295 if (anEntity == SLVS_E_UNKNOWN)
296 anEntity = changeEntity(aRefAttr, aType);
298 Slvs_Entity anEnt = myStorage->getEntity(anEntity);
302 if (aType == SLVS_E_UNKNOWN)
304 else if (aType == SLVS_E_POINT_IN_2D || aType == SLVS_E_POINT_IN_3D)
305 theAttributes[aPtInd++] = anEntity; // the point is created
306 else { // another entity (not a point) is created
307 if (aEntInd < anInitNbOfAttr)
308 theAttributes[aEntInd] = anEntity;
310 theAttributes.push_back(anEntity);
316 Slvs_hEntity SketchSolver_Constraint::changeEntity(AttributeRefAttrPtr theAttribute, int& theType)
318 // Convert the object of the attribute to the feature
320 if (theAttribute->isObject() && theAttribute->object()) {
321 ResultConstructionPtr aRC = std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(
322 theAttribute->object());
324 myErrorMsg = SketchSolver_Error::NOT_INITIALIZED();
325 return SLVS_E_UNKNOWN;
327 std::shared_ptr<ModelAPI_Document> aDoc = aRC->document();
328 aFeature = aDoc->feature(aRC);
330 return changeEntity(aFeature, theType);
333 return changeEntity(theAttribute->attr(), theType);
336 Slvs_hEntity SketchSolver_Constraint::changeEntity(AttributePtr theEntity, int& theType)
338 Slvs_hEntity aResult = SLVS_E_UNKNOWN;
339 if (!theEntity || !isInitialized(theEntity))
340 return SLVS_E_UNKNOWN;
342 // If the entity is already in the group, try to find it
343 std::map<std::shared_ptr<ModelAPI_Attribute>, Slvs_hEntity>::const_iterator anEntIter =
344 myAttributeMap.find(theEntity);
345 Slvs_Entity aCurrentEntity;
346 aCurrentEntity.h = SLVS_E_UNKNOWN;
347 if (anEntIter != myAttributeMap.end())
348 aCurrentEntity = myStorage->getEntity(anEntIter->second);
350 aResult = myGroup->getAttributeId(theEntity);
351 if (aResult != SLVS_E_UNKNOWN) {
352 Slvs_Entity anEnt = myStorage->getEntity(aResult);
353 theType = anEnt.type;
354 myAttributeMap[theEntity] = aResult;
359 Slvs_hGroup aGroupID = myGroup->getId();
361 std::shared_ptr<GeomDataAPI_Point> aPoint =
362 std::dynamic_pointer_cast<GeomDataAPI_Point>(theEntity);
364 double aXYZ[3] = {aPoint->x(), aPoint->y(), aPoint->z()};
365 Slvs_hParam aParams[3];
366 for (int i = 0; i < 3; i++) {
367 Slvs_Param aPar = aCurrentEntity.h != SLVS_E_UNKNOWN ?
368 myStorage->getParameter(aCurrentEntity.param[i]) :
369 Slvs_MakeParam(SLVS_E_UNKNOWN, aGroupID, 0.0);
371 aParams[i] = myStorage->addParameter(aPar);
374 if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
375 aCurrentEntity = Slvs_MakePoint3d(SLVS_E_UNKNOWN, aGroupID, aParams[0], aParams[1], aParams[2]);
376 else { // update entity data
377 for (int i = 0; i < 3; i++)
378 aCurrentEntity.param[i] = aParams[i];
380 aResult = myStorage->addEntity(aCurrentEntity);
382 // All entities except 3D points are created on workplane. So, if there is no workplane yet, then error
383 Slvs_hEntity aWorkplaneID = myGroup->getWorkplaneId();
384 if (aWorkplaneID == SLVS_E_UNKNOWN)
385 return SLVS_E_UNKNOWN;
388 std::shared_ptr<GeomDataAPI_Point2D> aPoint2D =
389 std::dynamic_pointer_cast<GeomDataAPI_Point2D>(theEntity);
391 double aXY[2] = {aPoint2D->x(), aPoint2D->y()};
392 Slvs_hParam aParams[2];
393 for (int i = 0; i < 2; i++) {
394 Slvs_Param aPar = aCurrentEntity.h != SLVS_E_UNKNOWN ?
395 myStorage->getParameter(aCurrentEntity.param[i]) :
396 Slvs_MakeParam(SLVS_E_UNKNOWN, aGroupID, 0.0);
398 aParams[i] = myStorage->addParameter(aPar);
401 if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
402 aCurrentEntity = Slvs_MakePoint2d(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID, aParams[0], aParams[1]);
403 else { // update entity data
404 for (int i = 0; i < 2; i++)
405 aCurrentEntity.param[i] = aParams[i];
407 aResult = myStorage->addEntity(aCurrentEntity);
409 // Scalar value (used for the distance entities)
410 AttributeDoublePtr aScalar = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(theEntity);
412 Slvs_Param aParam = aCurrentEntity.h != SLVS_E_UNKNOWN ?
413 myStorage->getParameter(aCurrentEntity.param[0]) :
414 Slvs_MakeParam(SLVS_E_UNKNOWN, aGroupID, 0.0);
415 aParam.val = aScalar->value();
416 Slvs_hParam aValue = myStorage->addParameter(aParam);
418 if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
419 aCurrentEntity = Slvs_MakeDistance(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID, aValue);
421 aCurrentEntity.param[0] = aValue;
422 aResult = myStorage->addEntity(aCurrentEntity);
427 myAttributeMap[theEntity] = aResult;
428 theType = aCurrentEntity.type;
432 Slvs_hEntity SketchSolver_Constraint::changeEntity(FeaturePtr theEntity, int& theType)
434 Slvs_hEntity aResult = SLVS_E_UNKNOWN;
435 if (!theEntity || !theEntity->data() || !theEntity->data()->isValid())
436 return SLVS_E_UNKNOWN;
437 // If the entity is already in the group, try to find it
438 std::map<FeaturePtr, Slvs_hEntity>::const_iterator anEntIter = myFeatureMap.find(theEntity);
439 Slvs_Entity aCurrentEntity;
440 aCurrentEntity.h = SLVS_E_UNKNOWN;
441 if (anEntIter != myFeatureMap.end())
442 aCurrentEntity = myStorage->getEntity(anEntIter->second);
444 aResult = myGroup->getFeatureId(theEntity);
445 if (aResult != SLVS_E_UNKNOWN) {
446 Slvs_Entity anEnt = myStorage->getEntity(aResult);
447 theType = anEnt.type;
448 myFeatureMap[theEntity] = aResult;
453 Slvs_hGroup aGroupID = myGroup->getId();
454 Slvs_hEntity aWorkplaneID = myGroup->getWorkplaneId();
455 DataPtr aData = theEntity->data();
457 // SketchPlugin features
458 const std::string& aFeatureKind = theEntity->getKind();
459 AttributePtr anAttribute;
462 if (aFeatureKind == SketchPlugin_Line::ID()) {
463 anAttribute = aData->attribute(SketchPlugin_Line::START_ID());
464 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
465 Slvs_hEntity aStart = changeEntity(anAttribute, anAttrType);
467 anAttribute = aData->attribute(SketchPlugin_Line::END_ID());
468 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
469 Slvs_hEntity aEnd = changeEntity(anAttribute, anAttrType);
471 if (aCurrentEntity.h == SLVS_E_UNKNOWN) // New entity
472 aCurrentEntity = Slvs_MakeLineSegment(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID, aStart, aEnd);
474 aCurrentEntity.point[0] = aStart;
475 aCurrentEntity.point[1] = aEnd;
477 aResult = myStorage->addEntity(aCurrentEntity);
480 else if (aFeatureKind == SketchPlugin_Circle::ID()) {
481 anAttribute = aData->attribute(SketchPlugin_Circle::CENTER_ID());
482 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
483 Slvs_hEntity aCenter = changeEntity(anAttribute, anAttrType);
485 anAttribute = aData->attribute(SketchPlugin_Circle::RADIUS_ID());
486 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
487 Slvs_hEntity aRadius = changeEntity(anAttribute, anAttrType);
489 if (aCurrentEntity.h == SLVS_E_UNKNOWN) { // New entity
490 Slvs_Entity aWorkplane = myStorage->getEntity(aWorkplaneID);
491 aCurrentEntity = Slvs_MakeCircle(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID,
492 aCenter, aWorkplane.normal, aRadius);
494 aCurrentEntity.point[0] = aCenter;
495 aCurrentEntity.distance = aRadius;
497 aResult = myStorage->addEntity(aCurrentEntity);
500 else if (aFeatureKind == SketchPlugin_Arc::ID()) {
501 anAttribute = aData->attribute(SketchPlugin_Arc::CENTER_ID());
502 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
503 Slvs_hEntity aCenter = changeEntity(anAttribute, anAttrType);
505 anAttribute = aData->attribute(SketchPlugin_Arc::START_ID());
506 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
507 Slvs_hEntity aStart = changeEntity(anAttribute, anAttrType);
509 anAttribute = aData->attribute(SketchPlugin_Arc::END_ID());
510 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
511 Slvs_hEntity aEnd = changeEntity(anAttribute, anAttrType);
513 if (aCurrentEntity.h == SLVS_E_UNKNOWN) { // New entity
514 Slvs_Entity aWorkplane = myStorage->getEntity(aWorkplaneID);
515 aCurrentEntity = Slvs_MakeArcOfCircle(SLVS_E_UNKNOWN, aGroupID, aWorkplaneID,
516 aWorkplane.normal, aCenter, aStart, aEnd);
518 aCurrentEntity.point[0] = aCenter;
519 aCurrentEntity.point[1] = aStart;
520 aCurrentEntity.point[2] = aEnd;
522 aResult = myStorage->addEntity(aCurrentEntity);
524 // Point (it has low probability to be an attribute of constraint, so it is checked at the end)
525 else if (aFeatureKind == SketchPlugin_Point::ID()) {
526 anAttribute = aData->attribute(SketchPlugin_Point::COORD_ID());
527 if (!isInitialized(anAttribute)) return SLVS_E_UNKNOWN;
528 // Both the sketch point and its attribute (coordinates) link to the same SolveSpace point identifier
529 aResult = changeEntity(anAttribute, anAttrType);
530 aCurrentEntity.type = SLVS_E_POINT_IN_3D;
533 if (aResult != SLVS_E_UNKNOWN) {
534 myFeatureMap[theEntity] = aResult;
535 theType = aCurrentEntity.type;
540 std::list<ConstraintPtr> SketchSolver_Constraint::constraints() const
542 std::list<ConstraintPtr> aConstraints;
543 aConstraints.push_back(myBaseConstraint);
547 void SketchSolver_Constraint::refresh()
550 std::map<AttributePtr, Slvs_hEntity>::iterator anAttrIter = myAttributeMap.begin();
551 for (; anAttrIter != myAttributeMap.end(); anAttrIter++) {
552 std::shared_ptr<GeomDataAPI_Point> aPoint =
553 std::dynamic_pointer_cast<GeomDataAPI_Point>(anAttrIter->first);
555 Slvs_Entity anEntity = myStorage->getEntity(anAttrIter->second);
557 for (int i = 0; i < 3; i++) {
558 Slvs_Param aPar = myStorage->getParameter(anEntity.param[i]);
561 if (fabs(aPoint->x() - aXYZ[0]) > tolerance ||
562 fabs(aPoint->y() - aXYZ[1]) > tolerance ||
563 fabs(aPoint->z() - aXYZ[2]) > tolerance)
564 aPoint->setValue(aXYZ[0], aXYZ[1], aXYZ[2]);
567 std::shared_ptr<GeomDataAPI_Point2D> aPoint2D =
568 std::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttrIter->first);
570 Slvs_Entity anEntity = myStorage->getEntity(anAttrIter->second);
572 for (int i = 0; i < 2; i++) {
573 Slvs_Param aPar = myStorage->getParameter(anEntity.param[i]);
576 if (fabs(aPoint2D->x() - aXY[0]) > tolerance ||
577 fabs(aPoint2D->y() - aXY[1]) > tolerance)
578 aPoint2D->setValue(aXY[0], aXY[1]);
580 // Scalar value (used for the distance entities)
581 AttributeDoublePtr aScalar =
582 std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(anAttrIter->first);
584 Slvs_Entity anEntity = myStorage->getEntity(anAttrIter->second);
585 Slvs_Param aPar = myStorage->getParameter(anEntity.param[0]);
586 if (fabs(aScalar->value() - aPar.val) > tolerance)
587 aScalar->setValue(aPar.val);
593 std::map<AttributePtr, Slvs_hParam>::iterator aValIter = myValueMap.begin();
594 for (; aValIter != myValueMap.end(); aValIter++) {
595 AttributeDoublePtr aScalar =
596 std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(anAttrIter->first);
598 Slvs_Param aPar = myStorage->getParameter(anAttrIter->second);
599 aScalar->setValue(aPar.val);
604 Slvs_hEntity SketchSolver_Constraint::getId(FeaturePtr theFeature) const
606 std::map<FeaturePtr, Slvs_hEntity>::const_iterator aFIter = myFeatureMap.find(theFeature);
607 if (aFIter == myFeatureMap.end())
608 return SLVS_E_UNKNOWN;
609 // check the Feature is really in the storage
610 Slvs_Entity anEntity = myStorage->getEntity(aFIter->second);
611 if (anEntity.h == SLVS_E_UNKNOWN) {
614 anEntity.h = const_cast<SketchSolver_Constraint*>(this)->changeEntity(aFIter->first, aType);
615 const_cast<SketchSolver_Constraint*>(this)->myFeatureMap[theFeature] = anEntity.h;
620 Slvs_hEntity SketchSolver_Constraint::getId(AttributePtr theAttribute) const
622 std::map<AttributePtr, Slvs_hEntity>::const_iterator anAttrIter = myAttributeMap.find(theAttribute);
623 if (anAttrIter == myAttributeMap.end())
624 return SLVS_E_UNKNOWN;
625 return anAttrIter->second;
628 bool SketchSolver_Constraint::isInitialized(AttributePtr theAttribute)
630 if (theAttribute->isInitialized())
632 myErrorMsg = SketchSolver_Error::NOT_INITIALIZED();
637 void SketchSolver_Constraint::calculateMiddlePoint(
638 const Slvs_Entity& theEntity, double theCoeff, double& theX, double& theY) const
640 if (theEntity.type == SLVS_E_LINE_SEGMENT) {
641 double aStartEndXY[2][2];
643 for (int i = 0; i < 2; i++) {
644 aPoint = myStorage->getEntity(theEntity.point[i]);
645 for (int j = 0; j < 2; j++)
646 aStartEndXY[i][j] = myStorage->getParameter(aPoint.param[j]).val;
648 theX = (1.0 - theCoeff) * aStartEndXY[0][0] + theCoeff * aStartEndXY[1][0];
649 theY = (1.0 - theCoeff) * aStartEndXY[0][1] + theCoeff * aStartEndXY[1][1];
650 } else if (theEntity.type == SLVS_E_ARC_OF_CIRCLE) {
651 double anArcPoint[3][2];
653 for (int i = 0; i < 3; i++) {
654 aPoint = myStorage->getEntity(theEntity.point[i]);
655 for (int j = 0; j < 2; j++)
656 anArcPoint[i][j] = myStorage->getParameter(aPoint.param[j]).val;
658 // project last point of arc on the arc
659 double x = anArcPoint[1][0] - anArcPoint[0][0];
660 double y = anArcPoint[1][1] - anArcPoint[0][1];
661 double aRad = sqrt(x*x + y*y);
662 x = anArcPoint[2][0] - anArcPoint[0][0];
663 y = anArcPoint[2][1] - anArcPoint[0][1];
664 double aNorm = sqrt(x*x + y*y);
665 if (aNorm >= tolerance) {
666 anArcPoint[2][0] = x * aRad / aNorm;
667 anArcPoint[2][1] = y * aRad / aNorm;
669 anArcPoint[1][0] -= anArcPoint[0][0];
670 anArcPoint[1][1] -= anArcPoint[0][1];
671 if (theCoeff < tolerance) {
672 theX = anArcPoint[0][0] + anArcPoint[1][0];
673 theY = anArcPoint[0][1] + anArcPoint[1][1];
675 } else if (1 - theCoeff < tolerance) {
676 theX = anArcPoint[0][0] + anArcPoint[2][0];
677 theY = anArcPoint[0][1] + anArcPoint[2][1];
681 std::shared_ptr<GeomAPI_Dir2d> aStartDir(new GeomAPI_Dir2d(anArcPoint[1][0], anArcPoint[1][1]));
682 std::shared_ptr<GeomAPI_Dir2d> aEndDir(new GeomAPI_Dir2d(anArcPoint[2][0], anArcPoint[2][1]));
683 double anAngle = aStartDir->angle(aEndDir);
687 double aCos = cos(anAngle);
688 double aSin = sin(anAngle);
689 theX = anArcPoint[0][0] + anArcPoint[1][0] * aCos - anArcPoint[1][1] * aSin;
690 theY = anArcPoint[0][1] + anArcPoint[1][0] * aSin + anArcPoint[1][1] * aCos;
694 void SketchSolver_Constraint::makeTemporary() const
696 std::vector<Slvs_hConstraint>::const_iterator anIt = mySlvsConstraints.begin();
697 for (; anIt != mySlvsConstraints.end(); anIt++)
698 myStorage->addTemporaryConstraint(*anIt);