src/SketchPlugin/SketchPlugin_Split.cpp

   1 // Copyright (C) 2014-2017  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or
  18 // email : webmaster.salome@opencascade.com<mailto:webmaster.salome@opencascade.com>
  19 //
  20
  21 #include "SketchPlugin_Split.h"
  22
  23 #include <Events_Message.h>
  24
  25 #include <GeomAPI_Dir2d.h>
  26 #include <GeomAPI_Edge.h>
  27 #include <GeomAPI_Pnt2d.h>
  28 #include <GeomAPI_XY.h>
  29 #include <GeomDataAPI_Point2D.h>
  30 #include <GeomAlgoAPI_ShapeTools.h>
  31 #include <GeomAlgoAPI_CompoundBuilder.h>
  32
  33 #include <ModelAPI_AttributeBoolean.h>
  34 #include <ModelAPI_AttributeDouble.h>
  35 #include <ModelAPI_AttributeRefAttr.h>
  36 #include <ModelAPI_AttributeReference.h>
  37 #include <ModelAPI_AttributeString.h>
  38 #include <ModelAPI_Events.h>
  39 #include <ModelAPI_Validator.h>
  40 #include <ModelAPI_Session.h>
  41 #include <ModelAPI_Tools.h>
  42
  43 #include <ModelGeomAlgo_Shape.h>
  44
  45 #include <SketchPlugin_Arc.h>
  46 #include <SketchPlugin_Circle.h>
  47 #include <SketchPlugin_ConstraintCoincidence.h>
  48 #include <SketchPlugin_ConstraintEqual.h>
  49 #include <SketchPlugin_ConstraintLength.h>
  50 #include <SketchPlugin_ConstraintMiddle.h>
  51 #include <SketchPlugin_ConstraintMirror.h>
  52 #include <SketchPlugin_ConstraintParallel.h>
  53 #include <SketchPlugin_ConstraintTangent.h>
  54 #include <SketchPlugin_Line.h>
  55 #include <SketchPlugin_MultiRotation.h>
  56 #include <SketchPlugin_MultiTranslation.h>
  57 #include <SketchPlugin_Point.h>
  58 #include <SketchPlugin_Tools.h>
  59
  60 #include <ModelGeomAlgo_Point2D.h>
  61 #include <ModelAPI_EventReentrantMessage.h>
  62 #include <Events_Loop.h>
  63
  64 #include <cmath>
  65
  66 //#define CREATE_CONSTRAINTS
  67
  68 #ifdef DEBUG_SPLIT
  69 #include <iostream>
  70 #endif
  71
  72 static const double PI = 3.141592653589793238463;
  73
  74 SketchPlugin_Split::SketchPlugin_Split()
  75 {
  76 }
  77
  78 void SketchPlugin_Split::initAttributes()
  79 {
  80   data()->addAttribute(SELECTED_OBJECT(), ModelAPI_AttributeReference::typeId());
  81   data()->addAttribute(SELECTED_POINT(), GeomDataAPI_Point2D::typeId());
  82
  83   data()->addAttribute(PREVIEW_POINT(), GeomDataAPI_Point2D::typeId());
  84   data()->addAttribute(PREVIEW_OBJECT(), ModelAPI_AttributeReference::typeId());
  85
  86   data()->attribute(PREVIEW_POINT())->setIsArgument(false);
  87   data()->attribute(SELECTED_POINT())->setIsArgument(false);
  88   data()->attribute(PREVIEW_OBJECT())->setIsArgument(false);
  89
  90   ModelAPI_Session::get()->validators()->registerNotObligatory(getKind(), PREVIEW_POINT());
  91   ModelAPI_Session::get()->validators()->registerNotObligatory(getKind(), PREVIEW_OBJECT());
  92
  93   // TODO: remove
  94   //data()->addAttribute(SketchPlugin_Constraint::VALUE(), ModelAPI_AttributeReference::typeId());
  95   //data()->addAttribute(SketchPlugin_Constraint::ENTITY_A(), ModelAPI_AttributeRefAttr::typeId());
  96   //data()->addAttribute(SketchPlugin_Constraint::ENTITY_B(), ModelAPI_AttributeRefAttr::typeId());
  97 }
  98
  99 void SketchPlugin_Split::execute()
 100 {
 101   std::shared_ptr<ModelAPI_Data> aData = data();
 102
 103   // Check the base objects are initialized.
 104   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 105                                                            data()->attribute(SELECTED_OBJECT()));
 106   //ObjectPtr aBaseObject = anObjectAttr->value();
 107   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 108   //                                          aData->attribute(SketchPlugin_Constraint::VALUE()));
 109   if(!aBaseObjectAttr->isInitialized()) {
 110     setError("Error: Base object is not initialized.");
 111     return;
 112   }
 113   ObjectPtr aBaseObject = aBaseObjectAttr->value();
 114   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
 115   //  getPointOfRefAttr(aData->attribute(SketchPlugin_Constraint::ENTITY_A()));
 116   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
 117   //  getPointOfRefAttr(aData->attribute(SketchPlugin_Constraint::ENTITY_B()));
 118   if (!aFirstPointAttrOfSplit.get() || !aFirstPointAttrOfSplit->isInitialized() ||
 119       !aSecondPointAttrOfSplit.get() || !aSecondPointAttrOfSplit->isInitialized()) {
 120     setError("Error: Sub-shape is not initialized.");
 121     return;
 122   }
 123
 124   /// Remove reference of this feature to feature used in preview, it is not necessary anymore
 125   /// as trim will be removed after execute
 126   AttributeReferencePtr aPreviewObjectAttr =
 127                      std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 128                      data()->attribute(PREVIEW_OBJECT()));
 129
 130   ObjectPtr aPreviewObject = aPreviewObjectAttr->value();
 131   AttributePoint2DPtr aPoint = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 132                                            data()->attribute(PREVIEW_POINT()));
 133   std::shared_ptr<GeomAPI_Pnt2d> aPreviewPnt2d = aPoint->pnt();
 134   // nullify pointer of preview attribute
 135   aPreviewObjectAttr->setValue(ResultPtr());
 136   bool anIsEqualPreviewAndSelected = aPreviewObject == aBaseObject;
 137
 138   // Wait all constraints being created, then send update events
 139   static Events_ID anUpdateEvent = Events_Loop::eventByName(EVENT_OBJECT_UPDATED);
 140   bool isUpdateFlushed = Events_Loop::loop()->isFlushed(anUpdateEvent);
 141   if (isUpdateFlushed)
 142     Events_Loop::loop()->setFlushed(anUpdateEvent, false);
 143
 144   // Find feature constraints
 145   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObject);
 146   ResultPtr aBaseFeatureResult = getFeatureResult(aBaseFeature);
 147   std::set<FeaturePtr> aFeaturesToDelete, aFeaturesToUpdate;
 148
 149   //std::map<FeaturePtr, IdToPointPair> aTangentFeatures;
 150   std::map<FeaturePtr, IdToPointPair> aCoincidenceToFeature;
 151   getConstraints(aFeaturesToDelete, aFeaturesToUpdate, /*aTangentFeatures, */
 152                  aCoincidenceToFeature);
 153
 154   std::map<AttributePtr, std::list<AttributePtr> > aBaseRefAttributes;
 155   std::list<AttributePtr> aRefsToFeature;
 156   getRefAttributes(aBaseFeature, aBaseRefAttributes, aRefsToFeature);
 157
 158   std::map<AttributePtr, AttributePtr> aBasePointModifiedAttributes;
 159
 160 #ifdef DEBUG_SPLIT
 161   std::cout << std::endl;
 162   std::cout << "SketchPlugin_Split::execute()" << std::endl;
 163   std::cout << std::endl;
 164
 165   SketchPlugin_Sketch* aSketch = sketch();
 166   std::cout << "SKETCH FEATURES (before split) [" << aSketch->numberOfSubs() << "]:" << std::endl;
 167   for (int i = 0, aNbSubs = aSketch->numberOfSubs(); i < aNbSubs; i++) {
 168     std::cout << getFeatureInfo(aSketch->subFeature(i), false) << std::endl;
 169   }
 170
 171   std::cout << std::endl;
 172   std::cout << "---- IN PARAMETERS ----" << std::endl;
 173   std::cout << "Base feature:" << getFeatureInfo(aBaseFeature) << std::endl;
 174   std::cout << std::endl;
 175
 176   if (!aCoincidenceToFeature.empty()) {
 177     std::cout << "Coincidences to base feature[" <<
 178       aCoincidenceToFeature.size() << "]: " << std::endl;
 179     std::map<FeaturePtr, IdToPointPair>::const_iterator anIt = aCoincidenceToFeature.begin(),
 180                                                         aLast = aCoincidenceToFeature.end();
 181     for (int i = 1; anIt != aLast; anIt++, i++) {
 182       FeaturePtr aFeature = (*anIt).first;
 183       std::string anAttributeId = (*anIt).second.first;
 184       std::shared_ptr<GeomDataAPI_Point2D> aPointAttr = (*anIt).second.second;
 185
 186       std::cout << i << "-" << getFeatureInfo(aFeature) << std::endl;
 187       std::cout <<     " -Attribute to correct:" << anAttributeId << std::endl;
 188       std::cout <<     " -Point attribute:" <<
 189         ModelGeomAlgo_Point2D::getPointAttributeInfo(aPointAttr) << std::endl;
 190     }
 191   }
 192
 193   std::map<AttributePtr, std::list<AttributePtr> >::const_iterator
 194     aRefIt = aBaseRefAttributes.begin(), aRefLast = aBaseRefAttributes.end();
 195   std::cout << std::endl << "References to attributes of base feature [" <<
 196     aBaseRefAttributes.size() << "]" << std::endl;
 197   for (; aRefIt != aRefLast; aRefIt++) {
 198     AttributePtr aBaseAttr = aRefIt->first;
 199     std::list<AttributePtr> aRefAttributes = aRefIt->second;
 200     std::string aRefsInfo;
 201     std::list<AttributePtr>::const_iterator aRefAttrIt = aRefAttributes.begin(),
 202                                             aRefAttrLast = aRefAttributes.end();
 203     for (; aRefAttrIt != aRefAttrLast; aRefAttrIt++) {
 204       if (!aRefsInfo.empty())
 205         aRefsInfo.append(",");
 206       AttributePtr aRAttr = *aRefAttrIt;
 207       aRefsInfo.append(aRAttr->id());
 208       FeaturePtr aRFeature = ModelAPI_Feature::feature(aRAttr->owner());
 209       aRefsInfo.append("(" + aRFeature->name() + ") ");
 210     }
 211     std::shared_ptr<GeomDataAPI_Point2D> aPointAttr =
 212       std::dynamic_pointer_cast<GeomDataAPI_Point2D>(aBaseAttr);
 213     std::cout << aPointAttr->id().c_str() <<
 214       ": " << "[" << aRefAttributes.size() << "] " << aRefsInfo << std::endl;
 215   }
 216   std::cout << std::endl;
 217   std::cout << std::endl << "References to base feature [" <<
 218     aRefsToFeature.size() << "]" << std::endl;
 219   std::list<AttributePtr>::const_iterator aRefAttrIt = aRefsToFeature.begin(),
 220                                           aRefAttrLast = aRefsToFeature.end();
 221   std::string aRefsInfo;
 222   for (; aRefAttrIt != aRefAttrLast; aRefAttrIt++) {
 223     if (!aRefsInfo.empty())
 224       aRefsInfo.append(",");
 225     AttributePtr aRAttr = *aRefAttrIt;
 226     aRefsInfo.append(aRAttr->id());
 227     FeaturePtr aRFeature = ModelAPI_Feature::feature(aRAttr->owner());
 228     aRefsInfo.append("(" + aRFeature->name() + ") ");
 229   }
 230   std::cout << "[" << aRefsToFeature.size() << "] " << aRefsInfo << std::endl;
 231
 232
 233   std::cout << std::endl;
 234   std::cout << "---- SPLIT ----" << std::endl;
 235   std::cout << std::endl;
 236 #endif
 237
 238   std::string aFeatureKind = aBaseFeature->getKind();
 239   FeaturePtr aSplitFeature, anAfterFeature;
 240   std::set<AttributePoint2DPtr> aFurtherCoincidences;
 241   std::set<FeaturePtr> aCreatedFeatures;
 242   std::set<std::pair<AttributePtr, AttributePtr>> aModifiedAttributes;
 243   FeaturePtr aReplacingFeature, aNewFeature;
 244   if (aFeatureKind == SketchPlugin_Line::ID())
 245     aNewFeature = splitLine(aSplitFeature, aBaseFeature, anAfterFeature,
 246                             aFurtherCoincidences, aCreatedFeatures, aModifiedAttributes);
 247   else if (aFeatureKind == SketchPlugin_Arc::ID())
 248     aNewFeature = splitArc(aSplitFeature, aBaseFeature, anAfterFeature, aFurtherCoincidences,
 249                            aCreatedFeatures, aModifiedAttributes);
 250   if (aFeatureKind == SketchPlugin_Circle::ID()) {
 251     FeaturePtr aCircleFeature = aBaseFeature;
 252     aReplacingFeature = splitCircle(aSplitFeature, aBaseFeature, anAfterFeature,
 253                                     aFurtherCoincidences, aCreatedFeatures, aModifiedAttributes);
 254
 255     updateRefFeatureConstraints(getFeatureResult(aBaseFeature), aRefsToFeature);
 256
 257     AttributePtr aCenterAttr = aCircleFeature->attribute(SketchPlugin_Circle::CENTER_ID());
 258     aFeaturesToDelete.insert(aCircleFeature);
 259     // as circle is removed, temporary fill this attribute*/
 260     aBaseObjectAttr->setObject(ResultPtr());
 261   }
 262
 263 #ifdef DEBUG_SPLIT
 264   std::cout << "---- OUT PARAMETERS ----" << std::endl;
 265   std::cout << "Base modified feature:" << getFeatureInfo(aBaseFeature) << std::endl;
 266   std::cout << "Split feature:" << getFeatureInfo(aSplitFeature) << std::endl;
 267   std::cout << "After feature:" << getFeatureInfo(anAfterFeature) << std::endl;
 268   std::cout << std::endl;
 269
 270   std::cout << "Created features by split:[" << aCreatedFeatures.size() << "]" << std::endl;
 271   std::set<FeaturePtr>::const_iterator aFIt = aCreatedFeatures.begin(),
 272                                        aFLast = aCreatedFeatures.end();
 273   for (; aFIt != aFLast; aFIt++) {
 274     std::cout << getFeatureInfo(*aFIt) << std::endl;
 275   }
 276   std::cout << std::endl;
 277
 278   std::cout << "Attributes for further Coincidences:" << std::endl;
 279   std::set<AttributePoint2DPtr>::const_iterator anIt = aFurtherCoincidences.begin(),
 280                                                 aLast = aFurtherCoincidences.end();
 281   for (; anIt != aLast; anIt++) {
 282     AttributePtr anAttribute = *anIt;
 283     FeaturePtr aFeature = ModelAPI_Feature::feature(anAttribute->owner());
 284     std::cout << ModelGeomAlgo_Point2D::getPointAttributeInfo(anAttribute)
 285               << " [" << getFeatureInfo(aFeature, false) << "]" << std::endl;
 286   }
 287
 288   std::cout << "Modifed attributes (constraints to attributes are moved here):" << std::endl;
 289   std::set<std::pair<AttributePtr, AttributePtr> >::const_iterator
 290     aPIt = aModifiedAttributes.begin(), aPLast = aModifiedAttributes.end();
 291   std::string aResInfo;
 292   for (; aPIt != aPLast; aPIt++) {
 293     if (!aResInfo.empty())
 294       aResInfo += "\n";
 295
 296     std::pair<AttributePtr, AttributePtr> aPair = *aPIt;
 297
 298     AttributePtr anAttr = aPair.first;
 299     aResInfo.append(anAttr->id());
 300     FeaturePtr aFeature = ModelAPI_Feature::feature(anAttr->owner());
 301     aResInfo.append("(" + aFeature->name() + ") ");
 302
 303     aResInfo.append("  - is modified to -  ");
 304
 305     anAttr = aPair.second;
 306     aResInfo.append(anAttr->id());
 307     aFeature = ModelAPI_Feature::feature(anAttr->owner());
 308     aResInfo.append("(" + aFeature->name() + ") ");
 309   }
 310   std::cout << aResInfo << std::endl;
 311 #endif
 312
 313   std::set<ResultPtr> aFeatureResults;
 314   aFeatureResults.insert(getFeatureResult(aBaseFeature));
 315   if (anAfterFeature.get() && anAfterFeature != aBaseFeature)
 316     aFeatureResults.insert(getFeatureResult(anAfterFeature));
 317
 318   // coincidence to feature
 319   updateCoincidenceConstraintsToFeature(aCoincidenceToFeature, aFurtherCoincidences,
 320                                         aFeatureResults, aSplitFeature, aFeaturesToDelete);
 321
 322   updateRefAttConstraints(aBaseRefAttributes, aModifiedAttributes);
 323
 324   // delete constraints
 325 #ifdef DEBUG_SPLIT
 326   std::cout << "remove features and references:" << std::endl;
 327   std::set<FeaturePtr>::const_iterator aDIt = aFeaturesToDelete.begin(),
 328                                        aDLast = aFeaturesToDelete.end();
 329   for (; aDIt != aDLast; aDIt++) {
 330     std::cout << getFeatureInfo(*aDIt, false) << std::endl;
 331     std::cout << std::endl;
 332   }
 333 #endif
 334   ModelAPI_Tools::removeFeaturesAndReferences(aFeaturesToDelete);
 335   Events_Loop::loop()->flush(Events_Loop::eventByName(EVENT_OBJECT_DELETED));
 336
 337 #ifdef DEBUG_SPLIT
 338   std::cout << "update features after split:" << std::endl;
 339   std::set<FeaturePtr>::const_iterator anUIt = aFeaturesToUpdate.begin(),
 340                                        anULast = aFeaturesToUpdate.end();
 341   for (; anUIt != anULast; anUIt++) {
 342     std::cout << getFeatureInfo(*anUIt, false) << std::endl;
 343     std::cout << std::endl;
 344   }
 345 #endif
 346   updateFeaturesAfterSplit(aFeaturesToUpdate);
 347
 348   // Send events to update the sub-features by the solver.
 349   if(isUpdateFlushed) {
 350     Events_Loop::loop()->setFlushed(anUpdateEvent, true);
 351   }
 352
 353     if (anIsEqualPreviewAndSelected) {
 354     // equal preview and selected objects
 355     // nothing to do if the preview and selected objects are different
 356     ResultPtr aReplacingResult;
 357     if (aReplacingFeature.get()) {
 358       aReplacingFeature->execute(); // need it to obtain result
 359       aReplacingResult = getFeatureResult(aReplacingFeature);
 360     }
 361     if (aReplacingResult.get()) { // base object was removed
 362       aPreviewObject = aReplacingResult;
 363       //aMessage->setSelectedObject(aReplacingResult);
 364
 365       //GeomShapePtr aSelectedShape = aReplacingResult->shape();
 366       //std::shared_ptr<GeomAPI_Pnt> aPreviewPnt = sketch()->to3D(aPreviewPnt2d->x(),
 367       //                                                          aPreviewPnt2d->y());
 368       //std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 369       //if (ModelGeomAlgo_Point2D::isPointOnEdge(aSelectedShape, aPreviewPnt, aProjectedPoint)) {
 370         //bool aValue = true;
 371       //}
 372       //aBaseShape = aShape;
 373
 374 #ifdef DEBUG_SPLIT
 375       if (!aSelectedShape.get())
 376         std::cout << "Set empty selected object" << std::endl;
 377       else
 378         std::cout << "Set shape with ShapeType: " << aSelectedShape->shapeTypeStr() << std::endl;
 379 #endif
 380       //bool aValue = true;
 381     }
 382     else {
 383       aPreviewObject = ObjectPtr();
 384
 385       aBaseFeature->execute(); // should recompute shapes of result to do not check obsolete one
 386       aBaseObject = getFeatureResult(aBaseFeature);
 387       std::shared_ptr<GeomAPI_Pnt> aPreviewPnt = sketch()->to3D(aPreviewPnt2d->x(),
 388                                                                 aPreviewPnt2d->y());
 389       ResultPtr aBaseResult = std::dynamic_pointer_cast<ModelAPI_Result>(aBaseObject);
 390       if (aBaseResult) {
 391         GeomShapePtr aShape = aBaseResult->shape();
 392         std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 393         if (ModelGeomAlgo_Point2D::isPointOnEdge(aShape, aPreviewPnt, aProjectedPoint))
 394           aPreviewObject = aBaseResult;
 395       }
 396       if (!aPreviewObject.get() && aNewFeature.get()) {
 397         ResultPtr aNewFeatureResult = getFeatureResult(aNewFeature);
 398         if (aNewFeatureResult.get()) {
 399           GeomShapePtr aShape = aNewFeatureResult->shape();
 400           std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 401           if (ModelGeomAlgo_Point2D::isPointOnEdge(aShape, aPreviewPnt, aProjectedPoint))
 402             aPreviewObject = aNewFeatureResult;
 403         }
 404       }
 405     }
 406   }
 407   if (aPreviewObject.get()) {
 408     std::shared_ptr<ModelAPI_EventReentrantMessage> aMessage = std::shared_ptr
 409       <ModelAPI_EventReentrantMessage>(new ModelAPI_EventReentrantMessage(
 410                                            ModelAPI_EventReentrantMessage::eventId(), this));
 411     aMessage->setSelectedObject(aPreviewObject);
 412     Events_Loop::loop()->send(aMessage);
 413   }
 414
 415
 416 #ifdef DEBUG_SPLIT
 417   std::cout << "SKETCH FEATURES (after split) [" << aSketch->numberOfSubs() << "]:" << std::endl;
 418   for (int i = 0, aNbSubs = aSketch->numberOfSubs(); i < aNbSubs; i++) {
 419     std::cout << getFeatureInfo(aSketch->subFeature(i), false) << std::endl;
 420   }
 421 #endif
 422 }
 423
 424 std::string SketchPlugin_Split::processEvent(const std::shared_ptr<Events_Message>& theMessage)
 425 {
 426 #ifdef DEBUG_SPLIT
 427   std::cout << "SketchPlugin_Split::processEvent:" << data()->name() << std::endl;
 428 #endif
 429   std::string aFilledAttributeName;
 430
 431   std::shared_ptr<ModelAPI_EventReentrantMessage> aMessage =
 432         std::dynamic_pointer_cast<ModelAPI_EventReentrantMessage>(theMessage);
 433   if (aMessage.get()) {
 434     ObjectPtr anObject = aMessage->selectedObject();
 435     std::shared_ptr<GeomAPI_Pnt2d> aPoint = aMessage->clickedPoint();
 436
 437     if (anObject.get() && aPoint.get()) {
 438       //if (myCashedShapes.find(anObject) == myCashedShapes.end())
 439       //  fillObjectShapes(anObject, sketch()->data()->owner(), myCashedShapes, myObjectToPoints);
 440       if (myCashedShapes.find(anObject) == myCashedShapes.end())
 441         fillObjectShapes(anObject, sketch()->data()->owner());
 442       const std::set<GeomShapePtr>& aShapes = myCashedShapes[anObject];
 443       if (aShapes.size() > 1) {
 444         std::shared_ptr<ModelAPI_AttributeReference> aRefSelectedAttr =
 445                               std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 446                               data()->attribute(SELECTED_OBJECT()));
 447         std::shared_ptr<ModelAPI_AttributeReference> aRefPreviewAttr =
 448                               std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 449                               data()->attribute(PREVIEW_OBJECT()));
 450         aRefSelectedAttr->setValue(anObject);
 451         aRefPreviewAttr->setValue(anObject);
 452
 453         std::shared_ptr<GeomDataAPI_Point2D> aPointSelectedAttr =
 454                               std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 455                               data()->attribute(SELECTED_POINT()));
 456         std::shared_ptr<GeomDataAPI_Point2D> aPointPreviewAttr =
 457                               std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 458                               data()->attribute(PREVIEW_POINT()));
 459         aPointSelectedAttr->setValue(aPoint);
 460         aPointPreviewAttr->setValue(aPoint);
 461
 462         Events_Loop::loop()->flush(Events_Loop::eventByName(EVENT_OBJECT_UPDATED));
 463
 464         GeomShapePtr aSelectedShape = getSubShape(SELECTED_OBJECT(), SELECTED_POINT());
 465         if (aSelectedShape.get()) {
 466           aFilledAttributeName = SELECTED_OBJECT();
 467         }
 468         else {
 469           // #2480 - sub shape is not initialized when split sketch
 470           // If restarted operation use some selection on the shape that is split and
 471           // result selectiona can not participate in new split(checked shape above is null),
 472           // reset filled values of selection set in this method above
 473           aRefSelectedAttr->setValue(ResultPtr());
 474           aRefPreviewAttr->setValue(ResultPtr());
 475         }
 476   #ifdef DEBUG_SPLIT
 477         if (!aSelectedShape.get())
 478           std::cout << "Set empty selected object" << std::endl;
 479         else
 480           std::cout << "Set shape with ShapeType: " << aSelectedShape->shapeTypeStr() << std::endl;
 481   #endif
 482       }
 483     }
 484   }
 485   return aFilledAttributeName;
 486 }
 487
 488 AISObjectPtr SketchPlugin_Split::getAISObject(AISObjectPtr thePrevious)
 489 {
 490   /*AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 491                                            data()->attribute(SketchPlugin_Constraint::VALUE()));
 492   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
 493
 494   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointOfRefAttr(
 495                                         data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
 496   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointOfRefAttr(
 497                                         data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
 498
 499   if (aBaseObjectAttr->isInitialized() && aBaseFeature.get() &&
 500       aFirstPointAttrOfSplit->isInitialized() &&
 501       aSecondPointAttrOfSplit->isInitialized()) {
 502
 503     ResultPtr aResult = getFeatureResult(aBaseFeature);
 504     GeomShapePtr aBaseShape = aResult->shape();
 505     std::list<std::shared_ptr<GeomAPI_Pnt> > aPoints;
 506
 507     std::shared_ptr<GeomAPI_Pnt2d> aStartPnt2d = aFirstPointAttrOfSplit->pnt();
 508     std::shared_ptr<GeomAPI_Pnt> aStartPoint = sketch()->to3D(aStartPnt2d->x(), aStartPnt2d->y());
 509     aPoints.push_back(aStartPoint);
 510
 511     std::shared_ptr<GeomAPI_Pnt2d> aSecondPnt2d = aSecondPointAttrOfSplit->pnt();
 512     std::shared_ptr<GeomAPI_Pnt> aSecondPoint =
 513       sketch()->to3D(aSecondPnt2d->x(), aSecondPnt2d->y());
 514     aPoints.push_back(aSecondPoint);
 515
 516     std::set<std::shared_ptr<GeomAPI_Shape> > aSplitShapes;
 517
 518     GeomAlgoAPI_ShapeTools::splitShape_p(aBaseShape, aPoints, aSplitShapes);
 519     std::shared_ptr<GeomAPI_Shape> aShape =
 520       GeomAlgoAPI_ShapeTools::findShape(aPoints, aSplitShapes);
 521
 522     AISObjectPtr anAIS = thePrevious;
 523     if (aShape) {
 524       if (!anAIS)
 525         anAIS = AISObjectPtr(new GeomAPI_AISObject);
 526       anAIS->createShape(aShape);
 527       std::shared_ptr<ModelAPI_AttributeBoolean> anAuxiliaryAttr =
 528              aBaseFeature->data()->boolean(SketchPlugin_SketchEntity::AUXILIARY_ID());
 529
 530       bool isConstruction = anAuxiliaryAttr.get() != NULL && anAuxiliaryAttr->value();
 531
 532       std::vector<int> aColor;
 533       double aWidth = SketchPlugin_SketchEntity::SKETCH_LINE_WIDTH();
 534       int aLineStyle = SketchPlugin_SketchEntity::SKETCH_LINE_STYLE();
 535       if (isConstruction) {
 536         aColor = Config_PropManager::color("Visualization", "sketch_auxiliary_color");
 537         aWidth = SketchPlugin_SketchEntity::SKETCH_LINE_WIDTH_AUXILIARY();
 538         aLineStyle = SketchPlugin_SketchEntity::SKETCH_LINE_STYLE_AUXILIARY();
 539       }
 540       else {
 541         aColor = Config_PropManager::color("Visualization", "sketch_entity_color");
 542       }
 543       anAIS->setColor(aColor[0], aColor[1], aColor[2]);
 544       anAIS->setWidth(aWidth + 1);
 545       anAIS->setLineStyle(aLineStyle);
 546     }
 547     return anAIS;
 548   }
 549   return AISObjectPtr();*/
 550 #ifdef DEBUG_SPLIT
 551   std::cout << "SketchPlugin_Split::getAISObject: " << data()->name() << std::endl;
 552 #endif
 553
 554   AISObjectPtr anAIS = thePrevious;
 555
 556   std::list<std::shared_ptr<GeomAPI_Shape> > aShapes;
 557   GeomShapePtr aPreviewShape = getSubShape(PREVIEW_OBJECT(), PREVIEW_POINT());
 558   if (aPreviewShape.get())
 559     aShapes.push_back(aPreviewShape);
 560   GeomShapePtr aSelectedShape = getSubShape(SELECTED_OBJECT(), SELECTED_POINT());
 561   if (aSelectedShape.get())
 562     aShapes.push_back(aSelectedShape);
 563
 564   if (aShapes.empty())
 565     return AISObjectPtr();
 566
 567   GeomShapePtr aBaseShape = GeomAlgoAPI_CompoundBuilder::compound(aShapes);
 568   if (!aBaseShape.get())
 569     return AISObjectPtr();
 570
 571   if (aBaseShape.get()) {
 572     if (!anAIS)
 573       anAIS = AISObjectPtr(new GeomAPI_AISObject);
 574     anAIS->createShape(aBaseShape);
 575
 576     std::vector<int> aColor;
 577     aColor = Config_PropManager::color("Visualization", "operation_remove_feature_color");
 578     double aWidth = SketchPlugin_SketchEntity::SKETCH_LINE_WIDTH();
 579     int aLineStyle = SketchPlugin_SketchEntity::SKETCH_LINE_STYLE();
 580     anAIS->setColor(aColor[0], aColor[1], aColor[2]);
 581     // width when there is not base object should be extened in several points
 582     // in order to see this preview over highlight
 583     anAIS->setWidth(aWidth+4);
 584     anAIS->setLineStyle(aLineStyle);
 585   }
 586   else
 587     anAIS = AISObjectPtr();
 588   return anAIS;
 589 }
 590
 591 //********************************************************************
 592 void SketchPlugin_Split::fillObjectShapes(const ObjectPtr& theObject,
 593                                           const ObjectPtr& theSketch)
 594 {
 595   std::set<std::shared_ptr<GeomAPI_Shape> > aShapes;
 596   std::map<std::shared_ptr<GeomDataAPI_Point2D>, std::shared_ptr<GeomAPI_Pnt> > aPointToAttributes;
 597   std::set<std::shared_ptr<GeomDataAPI_Point2D> > aRefAttributes;
 598   // current feature
 599   FeaturePtr aFeature = ModelAPI_Feature::feature(theObject);
 600   // edges on feature
 601   std::set<ResultPtr> anEdgeResults;
 602   ModelGeomAlgo_Shape::shapesOfType(aFeature, GeomAPI_Shape::EDGE, anEdgeResults);
 603   if (!anEdgeResults.empty()) {
 604     GeomShapePtr aFeatureShape = (*anEdgeResults.begin())->shape();
 605
 606     // coincidences to the feature
 607     ModelGeomAlgo_Point2D::getPointsOfReference(aFeature, SketchPlugin_ConstraintCoincidence::ID(),
 608                          aRefAttributes, SketchPlugin_Point::ID(), SketchPlugin_Point::COORD_ID());
 609     // layed on feature coincidences to divide it on several shapes
 610     std::shared_ptr<ModelAPI_Data> aData = theSketch->data();
 611     std::shared_ptr<GeomDataAPI_Point> aC = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 612         aData->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 613     std::shared_ptr<GeomDataAPI_Dir> aX = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 614         aData->attribute(SketchPlugin_Sketch::DIRX_ID()));
 615     std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 616         aData->attribute(SketchPlugin_Sketch::NORM_ID()));
 617     std::shared_ptr<GeomAPI_Dir> aY(new GeomAPI_Dir(aNorm->dir()->cross(aX->dir())));
 618     std::list<std::shared_ptr<GeomAPI_Pnt> > aPoints;
 619     ModelGeomAlgo_Point2D::getPointsInsideShape_p(aFeatureShape, aRefAttributes, aC->pnt(),
 620                                                 aX->dir(), aY, aPoints, aPointToAttributes);
 621
 622     GeomAlgoAPI_ShapeTools::splitShape_p(aFeatureShape, aPoints, aShapes);
 623   }
 624   myCashedShapes[theObject] = aShapes;
 625   myCashedReferences[theObject] = aPointToAttributes;
 626 }
 627
 628 GeomShapePtr SketchPlugin_Split::getSubShape(const std::string& theObjectAttributeId,
 629                                              const std::string& thePointAttributeId)
 630 {
 631   GeomShapePtr aBaseShape;
 632
 633   AttributeReferencePtr anObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 634                                                        data()->attribute(theObjectAttributeId));
 635   ObjectPtr aBaseObject = anObjectAttr->value();
 636   if (!aBaseObject.get())
 637     return aBaseShape;
 638
 639   // point on feature
 640   AttributePoint2DPtr aPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 641                                            data()->attribute(thePointAttributeId));
 642   std::shared_ptr<GeomAPI_Pnt2d> anAttributePnt2d = aPointAttr->pnt();
 643   std::shared_ptr<GeomAPI_Pnt> anAttributePnt = sketch()->to3D(anAttributePnt2d->x(),
 644                                                                anAttributePnt2d->y());
 645
 646 #ifdef TRIM_SHAPE
 647   if (myCashedShapes.find(aBaseObject) == myCashedShapes.end())
 648     fillObjectShapes(aBaseObject, sketch()->data()->owner(), myCashedShapes, myObjectToPoints);
 649
 650   const std::set<GeomShapePtr>& aShapes = myCashedShapes[aBaseObject];
 651   if (!aShapes.empty()) {
 652     std::set<GeomShapePtr>::const_iterator anIt = aShapes.begin(), aLast = aShapes.end();
 653     for (; anIt != aLast; anIt++) {
 654       GeomShapePtr aShape = *anIt;
 655       std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 656       if (ModelGeomAlgo_Point2D::isPointOnEdge(aShape, anAttributePnt, aProjectedPoint))
 657         aBaseShape = aShape;
 658     }
 659   }
 660 #else
 661   if (myCashedShapes.find(aBaseObject) == myCashedShapes.end())
 662     fillObjectShapes(aBaseObject, sketch()->data()->owner());
 663
 664   std::shared_ptr<GeomAPI_Pnt> aStartPoint;
 665   std::shared_ptr<GeomAPI_Pnt> aSecondPoint;
 666   const std::set<GeomShapePtr>& aShapes = myCashedShapes[aBaseObject];
 667   std::set<GeomShapePtr>::const_iterator anIt = aShapes.begin(), aLast = aShapes.end();
 668   for (; anIt != aLast; anIt++) {
 669     GeomShapePtr aCurrentShape = *anIt;
 670     std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 671     if (ModelGeomAlgo_Point2D::isPointOnEdge(aCurrentShape, anAttributePnt, aProjectedPoint)) {
 672       if (aCurrentShape->shapeType() == GeomAPI_Shape::EDGE) {
 673         std::shared_ptr<GeomAPI_Edge> anEdge(new GeomAPI_Edge(aCurrentShape));
 674         aStartPoint = anEdge->firstPoint();
 675         aSecondPoint = anEdge->lastPoint();
 676       }
 677       break;
 678     }
 679   }
 680
 681   if (!aStartPoint.get() || !aSecondPoint.get())
 682     return aBaseShape;
 683
 684   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 685   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
 686   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObject/*aBaseObjectAttr->value()*/);
 687
 688   //AttributePoint2DPtr aFirstPointAttrOfSplit = getPointOfRefAttr(
 689   //                                      data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
 690   //AttributePoint2DPtr aSecondPointAttrOfSplit = getPointOfRefAttr(
 691   //                                      data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
 692   if (anObjectAttr->isInitialized() && aBaseFeature.get() && aPointAttr->isInitialized()) {
 693       //aFirstPointAttrOfSplit->isInitialized() &&
 694       //aSecondPointAttrOfSplit->isInitialized()) {
 695     ResultPtr aResult = getFeatureResult(aBaseFeature);
 696     GeomShapePtr aResultShape = aResult->shape();
 697     std::list<std::shared_ptr<GeomAPI_Pnt> > aPoints;
 698
 699     //std::shared_ptr<GeomAPI_Pnt2d> aStartPnt2d = aFirstPointAttrOfSplit->pnt();
 700     //std::shared_ptr<GeomAPI_Pnt> aStartPoint = sketch()->to3D(aStartPnt2d->x(), aStartPnt2d->y());
 701     aPoints.push_back(aStartPoint);
 702
 703     //std::shared_ptr<GeomAPI_Pnt2d> aSecondPnt2d = aSecondPointAttrOfSplit->pnt();
 704     //std::shared_ptr<GeomAPI_Pnt> aSecondPoint =
 705     //  sketch()->to3D(aSecondPnt2d->x(), aSecondPnt2d->y());
 706     aPoints.push_back(aSecondPoint);
 707
 708     std::set<std::shared_ptr<GeomAPI_Shape> > aSplitShapes;
 709     GeomAlgoAPI_ShapeTools::splitShape_p(aResultShape, aPoints, aSplitShapes);
 710     aBaseShape = GeomAlgoAPI_ShapeTools::findShape(aPoints, aSplitShapes);
 711 #endif
 712   }
 713   return aBaseShape;
 714 }
 715
 716 void SketchPlugin_Split::getFeaturePoints(const FeaturePtr& theFeature,
 717                                                     AttributePoint2DPtr& theStartPointAttr,
 718                                                     AttributePoint2DPtr& theEndPointAttr)
 719 {
 720   std::string aFeatureKind = theFeature->getKind();
 721   std::string aStartAttributeName, anEndAttributeName;
 722   if (aFeatureKind == SketchPlugin_Line::ID()) {
 723     aStartAttributeName = SketchPlugin_Line::START_ID();
 724     anEndAttributeName = SketchPlugin_Line::END_ID();
 725   }
 726   else if (aFeatureKind == SketchPlugin_Arc::ID()) {
 727     aStartAttributeName = SketchPlugin_Arc::START_ID();
 728     anEndAttributeName = SketchPlugin_Arc::END_ID();
 729   }
 730   if (!aStartAttributeName.empty() && !anEndAttributeName.empty()) {
 731     theStartPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 732                                          theFeature->attribute(aStartAttributeName));
 733     theEndPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 734                                          theFeature->attribute(anEndAttributeName));
 735   }
 736 }
 737
 738 void SketchPlugin_Split::getConstraints(std::set<FeaturePtr>& theFeaturesToDelete,
 739                                     std::set<FeaturePtr>& theFeaturesToUpdate,
 740                                     std::map<FeaturePtr, IdToPointPair>& theCoincidenceToFeature)
 741 {
 742   std::shared_ptr<ModelAPI_Data> aData = data();
 743
 744   // Check the base objects are initialized.
 745   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 746                                                            data()->attribute(SELECTED_OBJECT()));
 747   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 748   //                                          aData->attribute(SketchPlugin_Constraint::VALUE()));
 749   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
 750   ResultPtr aBaseFeatureResult = getFeatureResult(aBaseFeature);
 751
 752   std::set<AttributePtr> aRefsList = aBaseFeatureResult->data()->refsToMe();
 753   std::set<AttributePtr> aFRefsList = aBaseFeature->data()->refsToMe();
 754   aRefsList.insert(aFRefsList.begin(), aFRefsList.end());
 755
 756   std::set<AttributePtr>::const_iterator aIt;
 757   for (aIt = aRefsList.cbegin(); aIt != aRefsList.cend(); ++aIt) {
 758     std::shared_ptr<ModelAPI_Attribute> aAttr = (*aIt);
 759     FeaturePtr aRefFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aAttr->owner());
 760     std::string aRefFeatureKind = aRefFeature->getKind();
 761     if (aRefFeatureKind == SketchPlugin_ConstraintMirror::ID() ||
 762         aRefFeatureKind == SketchPlugin_MultiRotation::ID() ||
 763         aRefFeatureKind == SketchPlugin_MultiTranslation::ID() ||
 764         aRefFeatureKind == SketchPlugin_ConstraintMiddle::ID())
 765       theFeaturesToDelete.insert(aRefFeature);
 766     else if (aRefFeatureKind == SketchPlugin_ConstraintLength::ID())
 767       theFeaturesToUpdate.insert(aRefFeature);
 768     else if (aRefFeatureKind == SketchPlugin_ConstraintCoincidence::ID()) {
 769       std::string anAttributeToBeModified;
 770       AttributePoint2DPtr aCoincidentPoint;
 771       AttributeRefAttrPtr anAttrA = aRefFeature->refattr(SketchPlugin_Constraint::ENTITY_A());
 772       AttributeRefAttrPtr anAttrB = aRefFeature->refattr(SketchPlugin_Constraint::ENTITY_B());
 773       bool isToFeature = false;
 774       if (anAttrA->isObject() || anAttrB->isObject()) { /// coincidence to base feature
 775         FeaturePtr aFeature = anAttrA->isObject() ? ModelAPI_Feature::feature(anAttrA->object())
 776                                                   : FeaturePtr();
 777         isToFeature = aFeature.get() && aFeature == aBaseFeature;
 778         anAttributeToBeModified = anAttrA->id();
 779         if (!isToFeature) {
 780           aFeature = anAttrB->isObject() ? ModelAPI_Feature::feature(anAttrB->object())
 781                                          : FeaturePtr();
 782           isToFeature = aFeature.get() && aFeature == aBaseFeature;
 783           anAttributeToBeModified = anAttrB->id();
 784         }
 785         if (isToFeature)
 786           aCoincidentPoint = SketchPlugin_ConstraintCoincidence::getPoint(aRefFeature);
 787       }
 788       if (!isToFeature) { /// coincidence to point on base feature
 789         AttributePtr anAttribute;
 790
 791         if (!anAttrA->isObject()) {
 792           AttributePtr aCurAttribute = anAttrA->attr();
 793           if (aCurAttribute.get()) {
 794             FeaturePtr aCurFeature = ModelAPI_Feature::feature(aCurAttribute->owner());
 795             if (aCurFeature.get() && aCurFeature == aBaseFeature) {
 796               anAttribute = anAttrB->attr();
 797               anAttributeToBeModified = anAttrA->id();
 798             }
 799           }
 800         }
 801         if (!anAttribute.get() && !anAttrB->isObject()) {
 802           AttributePtr aCurAttribute = anAttrB->attr();
 803           if (aCurAttribute.get()) {
 804             FeaturePtr aCurFeature = ModelAPI_Feature::feature(aCurAttribute->owner());
 805             if (aCurFeature.get() && aCurFeature == aBaseFeature) {
 806               anAttribute = anAttrA->attr();
 807               anAttributeToBeModified = anAttrB->id();
 808             }
 809           }
 810         }
 811         if (anAttribute.get())
 812           aCoincidentPoint = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttribute);
 813       }
 814       if (aCoincidentPoint.get() && isToFeature)
 815         theCoincidenceToFeature[aRefFeature] = std::make_pair(anAttributeToBeModified,
 816                                                               aCoincidentPoint);
 817     }
 818   }
 819 }
 820
 821 void SketchPlugin_Split::getRefAttributes(const FeaturePtr& theFeature,
 822                                     std::map<AttributePtr, std::list<AttributePtr> >& theRefs,
 823                                     std::list<AttributePtr>& theRefsToFeature)
 824 {
 825   theRefs.clear();
 826
 827   std::list<AttributePtr> aPointAttributes =
 828     theFeature->data()->attributes(GeomDataAPI_Point2D::typeId());
 829   std::set<AttributePtr> aPointAttributesSet;
 830
 831   std::list<AttributePtr>::const_iterator aPIt =
 832     aPointAttributes.begin(), aPLast = aPointAttributes.end();
 833   for (; aPIt != aPLast; aPIt++)
 834     aPointAttributesSet.insert(*aPIt);
 835
 836   std::set<AttributePtr> aRefsAttributes = getFeatureResult(theFeature)->data()->refsToMe();
 837   std::set<AttributePtr> aFRefsList = theFeature->data()->refsToMe();
 838   aRefsAttributes.insert(aFRefsList.begin(), aFRefsList.end());
 839
 840   std::set<AttributePtr>::const_iterator aIt;
 841   for (aIt = aRefsAttributes.cbegin(); aIt != aRefsAttributes.cend(); ++aIt) {
 842     AttributePtr anAttr = (*aIt);
 843     FeaturePtr anAttrFeature = ModelAPI_Feature::feature(anAttr->owner());
 844     if (anAttrFeature.get() != this &&
 845         anAttr.get() && anAttr->attributeType() == ModelAPI_AttributeRefAttr::typeId()) {
 846       AttributeRefAttrPtr aRefAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(anAttr);
 847       if (!aRefAttr->isObject()) { /// find attributes referenced to feature point attributes
 848         AttributePtr anAttrInRef = aRefAttr->attr();
 849         if (anAttrInRef.get() &&
 850             aPointAttributesSet.find(anAttrInRef) != aPointAttributesSet.end()) {
 851           if (theRefs.find(anAttrInRef) != theRefs.end())
 852             theRefs[anAttrInRef].push_back(aRefAttr);
 853           else {
 854             std::list<AttributePtr> anAttrList;
 855             anAttrList.push_back(aRefAttr);
 856             theRefs[anAttrInRef] = anAttrList;
 857           }
 858         }
 859       }
 860       else { /// find attributes referenced to feature itself
 861         theRefsToFeature.push_back(anAttr);
 862       }
 863     }
 864   }
 865 }
 866
 867 void SketchPlugin_Split::updateCoincidenceConstraintsToFeature(
 868       const std::map<std::shared_ptr<ModelAPI_Feature>, IdToPointPair>& theCoincidenceToFeature,
 869       const std::set<std::shared_ptr<GeomDataAPI_Point2D> >& theFurtherCoincidences,
 870       const std::set<ResultPtr>& theFeatureResults,
 871       const FeaturePtr& theSplitFeature,
 872       std::set<FeaturePtr>& theFeaturesToDelete)
 873 {
 874   if (theCoincidenceToFeature.empty())
 875     return;
 876
 877   // we should build coincidence constraints to end of the split feature
 878   std::set<std::shared_ptr<GeomDataAPI_Point2D> > aNewCoincidencesToSplitFeature;
 879   AttributePoint2DPtr aStartPointAttr, anEndPointAttr;
 880   getFeaturePoints(theSplitFeature, aStartPointAttr, anEndPointAttr);
 881   if (theFurtherCoincidences.find(aStartPointAttr) == theFurtherCoincidences.end())
 882     aNewCoincidencesToSplitFeature.insert(aStartPointAttr);
 883   if (theFurtherCoincidences.find(anEndPointAttr) == theFurtherCoincidences.end())
 884     aNewCoincidencesToSplitFeature.insert(anEndPointAttr);
 885
 886   std::map<FeaturePtr, IdToPointPair>::const_iterator aCIt = theCoincidenceToFeature.begin(),
 887                                                             aCLast = theCoincidenceToFeature.end();
 888 #ifdef DEBUG_SPLIT
 889   std::cout << std::endl;
 890   std::cout << "Coincidences to feature(modified):"<< std::endl;
 891 #endif
 892   for (; aCIt != aCLast; aCIt++) {
 893     FeaturePtr aCoincFeature = aCIt->first;
 894     std::string anAttributeId = aCIt->second.first;
 895     std::string aSecondAttribute = anAttributeId == SketchPlugin_Constraint::ENTITY_A() ?
 896         SketchPlugin_Constraint::ENTITY_B() : SketchPlugin_Constraint::ENTITY_A();
 897
 898     AttributePoint2DPtr aCoincPoint = aCIt->second.second;
 899     std::set<AttributePoint2DPtr>::const_iterator aFCIt = theFurtherCoincidences.begin(),
 900                                                   aFCLast = theFurtherCoincidences.end();
 901     std::shared_ptr<GeomAPI_Pnt2d> aCoincPnt = aCoincPoint->pnt();
 902     AttributePoint2DPtr aFeaturePointAttribute;
 903     for (; aFCIt != aFCLast && !aFeaturePointAttribute.get(); aFCIt++) {
 904       AttributePoint2DPtr aFCAttribute = *aFCIt;
 905       if (aCoincPnt->isEqual(aFCAttribute->pnt()))
 906         aFeaturePointAttribute = aFCAttribute;
 907     }
 908     if (aFeaturePointAttribute.get()) {
 909       // create new constraint and remove the current
 910       aCoincFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
 911           SketchPlugin_ConstraintCoincidence::ID(),
 912           aFeaturePointAttribute, aCoincFeature->refattr(aSecondAttribute)->attr());
 913       theFeaturesToDelete.insert(aCIt->first);
 914       // create new coincidences to split feature points
 915       std::set<AttributePoint2DPtr>::const_iterator aSFIt = aNewCoincidencesToSplitFeature.begin(),
 916                                                     aSFLast = aNewCoincidencesToSplitFeature.end();
 917       for (; aSFIt != aSFLast; aSFIt++) {
 918         AttributePoint2DPtr aSFAttribute = *aSFIt;
 919         if (aCoincPnt->isEqual(aSFAttribute->pnt())) {
 920           SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
 921                            SketchPlugin_ConstraintCoincidence::ID(),
 922                            aSFAttribute, aCoincFeature->refattr(aSecondAttribute)->attr());
 923         }
 924       }
 925     }
 926     else {
 927       /// find feature by shape intersected the point
 928       ResultPtr aResultForCoincidence = *(theFeatureResults.begin());
 929
 930       if (theFeatureResults.size() > 1) { // try to find point on additional feature
 931         ResultPtr anAddtionalResult = *(theFeatureResults.begin()++);
 932         GeomShapePtr aShape = anAddtionalResult->shape();
 933
 934         std::shared_ptr<GeomAPI_Pnt2d> aPnt2d = aCoincPoint->pnt();
 935         std::shared_ptr<GeomAPI_Pnt> aPoint = sketch()->to3D(aPnt2d->x(), aPnt2d->y());
 936
 937         std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 938         if (ModelGeomAlgo_Point2D::isPointOnEdge(aShape, aPoint, aProjectedPoint))
 939           aResultForCoincidence = anAddtionalResult;
 940       }
 941       aCoincFeature->refattr(anAttributeId)->setObject(aResultForCoincidence);
 942     }
 943 #ifdef DEBUG_SPLIT
 944   std::cout << " -" << getFeatureInfo(aCoincFeature) << std::endl;
 945 #endif
 946   }
 947 }
 948
 949 void SketchPlugin_Split::updateRefFeatureConstraints(
 950                                                   const ResultPtr& theFeatureBaseResult,
 951                                                   const std::list<AttributePtr>& theRefsToFeature)
 952 {
 953   std::list<AttributePtr>::const_iterator anIt = theRefsToFeature.begin(),
 954                                           aLast = theRefsToFeature.end();
 955   for (; anIt != aLast; anIt++) {
 956     AttributeRefAttrPtr aRefAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*anIt);
 957     if (aRefAttr.get())
 958       aRefAttr->setObject(theFeatureBaseResult);
 959   }
 960 }
 961
 962 void SketchPlugin_Split::updateRefAttConstraints(
 963                     const std::map<AttributePtr, std::list<AttributePtr> >& theBaseRefAttributes,
 964                     const std::set<std::pair<AttributePtr, AttributePtr> >& theModifiedAttributes)
 965 {
 966 #ifdef DEBUG_SPLIT
 967   std::cout << "SketchPlugin_Split::updateRefAttConstraints" << std::endl;
 968 #endif
 969
 970   std::set<std::pair<AttributePtr, AttributePtr> >::const_iterator
 971     anIt = theModifiedAttributes.begin(),  aLast = theModifiedAttributes.end();
 972   for (; anIt != aLast; anIt++) {
 973     AttributePtr anAttribute = anIt->first;
 974
 975     /// not found in references
 976     if (theBaseRefAttributes.find(anAttribute) == theBaseRefAttributes.end())
 977       continue;
 978     std::list<AttributePtr> aRefAttributes = theBaseRefAttributes.at(anAttribute);
 979     std::list<AttributePtr>::const_iterator aRefIt = aRefAttributes.begin(),
 980                                             aRLast = aRefAttributes.end();
 981
 982     AttributePtr aNewAttribute = anIt->second;
 983     for (; aRefIt != aRLast; aRefIt++) {
 984       AttributeRefAttrPtr aRefAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*aRefIt);
 985       if (aRefAttr.get()) {
 986         aRefAttr->setAttr(aNewAttribute);
 987 #ifdef DEBUG_SPLIT
 988         FeaturePtr aFeature = ModelAPI_Feature::feature(aRefAttr->owner());
 989         std::cout << " -" << getFeatureInfo(aFeature) << std::endl;
 990 #endif
 991       }
 992     }
 993   }
 994 }
 995
 996 FeaturePtr SketchPlugin_Split::splitLine(FeaturePtr& theSplitFeature,
 997                                              FeaturePtr& theBaseFeatureModified,
 998                                              FeaturePtr& theAfterFeature,
 999                                              std::set<AttributePoint2DPtr>& thePoints,
1000                                              std::set<FeaturePtr>& theCreatedFeatures,
1001                  std::set<std::pair<AttributePtr, AttributePtr>>& theModifiedAttributes)
1002 {
1003   FeaturePtr anNewFeature;
1004
1005   std::set<FeaturePtr> aCreatedFeatures;
1006   FeaturePtr aConstraintFeature;
1007   theBaseFeatureModified = FeaturePtr(); // it will contain modified base feature
1008
1009   SketchPlugin_Sketch* aSketch = sketch();
1010   if (!aSketch)
1011     return anNewFeature;
1012
1013   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1014                                                            data()->attribute(SELECTED_OBJECT()));
1015   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1016   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
1017   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
1018   std::string aFeatureKind = aBaseFeature->getKind();
1019   if (aFeatureKind != SketchPlugin_Line::ID())
1020     return anNewFeature;
1021
1022   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
1023     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
1024   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
1025     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
1026   AttributePoint2DPtr aStartPointAttrOfBase, anEndPointAttrOfBase;
1027
1028   getFeaturePoints(aBaseFeature, aStartPointAttrOfBase, anEndPointAttrOfBase);
1029   if (!aStartPointAttrOfBase.get() && !anEndPointAttrOfBase.get()) {
1030     setError("Error: Feature has no start and end points.");
1031     return anNewFeature;
1032   }
1033
1034   arrangePointsOnLine(aStartPointAttrOfBase, anEndPointAttrOfBase,
1035                       aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1036
1037 #ifdef DEBUG_SPLIT
1038   std::cout << "Arranged points (to build split between 1st and 2nd points:" << std::endl;
1039   std::cout << "Start point: " <<
1040     ModelGeomAlgo_Point2D::getPointAttributeInfo(aStartPointAttrOfBase) << std::endl;
1041   std::cout << "1st point:   " <<
1042     ModelGeomAlgo_Point2D::getPointAttributeInfo(aFirstPointAttrOfSplit) << std::endl;
1043   std::cout << "2nd point:   " <<
1044     ModelGeomAlgo_Point2D::getPointAttributeInfo(aSecondPointAttrOfSplit) << std::endl;
1045   std::cout << "End point:   " <<
1046     ModelGeomAlgo_Point2D::getPointAttributeInfo(anEndPointAttrOfBase) << std::endl;
1047 #endif
1048
1049   /// create a split feature
1050   theSplitFeature =
1051     createLineFeature(aBaseFeature, aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1052   theCreatedFeatures.insert(theSplitFeature);
1053
1054   // before split feature
1055   if (aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1056     theModifiedAttributes.insert(std::make_pair(aStartPointAttrOfBase,
1057                                         theSplitFeature->attribute(SketchPlugin_Line::START_ID())));
1058   }
1059   else {
1060     theBaseFeatureModified = aBaseFeature; ///< use base feature to store all constraints here
1061     /// move end arc point to start of split
1062   }
1063
1064   // after split feature
1065   if (!aSecondPointAttrOfSplit->pnt()->isEqual(anEndPointAttrOfBase->pnt())) {
1066     FeaturePtr aFeature;
1067     if (!theBaseFeatureModified.get()) {
1068       aFeature = aBaseFeature; ///< use base feature to store all constraints here
1069       fillAttribute(aFeature->attribute(SketchPlugin_Line::START_ID()), aSecondPointAttrOfSplit);
1070       aFeature->execute(); // to update result
1071     }
1072     else {
1073       aFeature = createLineFeature(aBaseFeature, aSecondPointAttrOfSplit, anEndPointAttrOfBase);
1074       theCreatedFeatures.insert(aFeature);
1075       theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1076                                              aFeature->attribute(SketchPlugin_Line::END_ID())));
1077       anNewFeature = aFeature;
1078     }
1079     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1080                      SketchPlugin_ConstraintCoincidence::ID(),
1081                      theSplitFeature->attribute(SketchPlugin_Line::END_ID()),
1082                      aFeature->attribute(SketchPlugin_Line::START_ID()));
1083     theCreatedFeatures.insert(aConstraintFeature);
1084
1085     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1086                                 (aFeature->attribute(SketchPlugin_Line::START_ID())));
1087     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1088                                 (aFeature->attribute(SketchPlugin_Line::END_ID())));
1089
1090     if (!theBaseFeatureModified.get())
1091       theBaseFeatureModified = aFeature;
1092     else
1093       theAfterFeature = aFeature;
1094   }
1095   else {
1096     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1097                                   (theSplitFeature->attribute(SketchPlugin_Line::END_ID())));
1098     theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1099                                    theSplitFeature->attribute(SketchPlugin_Line::END_ID())));
1100   }
1101   // base split, that is defined before split feature should be changed at end
1102   // (after the after feature creation). Otherwise modified value will be used in after feature
1103   // before split feature
1104   if (!aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1105     /// move end arc point to start of split
1106     fillAttribute(theBaseFeatureModified->attribute(SketchPlugin_Line::END_ID()),
1107                                                     aFirstPointAttrOfSplit);
1108     theBaseFeatureModified->execute(); // to update result
1109     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1110                      SketchPlugin_ConstraintCoincidence::ID(),
1111                      theBaseFeatureModified->attribute(SketchPlugin_Line::END_ID()),
1112                      theSplitFeature->attribute(SketchPlugin_Line::START_ID()));
1113     theCreatedFeatures.insert(aConstraintFeature);
1114
1115     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1116                              (theBaseFeatureModified->attribute(SketchPlugin_Line::START_ID())));
1117     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1118                                (theBaseFeatureModified->attribute(SketchPlugin_Line::END_ID())));
1119   }
1120   else
1121     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1122                                        (theSplitFeature->attribute(SketchPlugin_Line::START_ID())));
1123
1124 #ifdef CREATE_CONSTRAINTS
1125   // additional constraints between split and base features
1126   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1127           SketchPlugin_ConstraintParallel::ID(),
1128           getFeatureResult(aBaseFeature),
1129           getFeatureResult(theSplitFeature));
1130   theCreatedFeatures.insert(aConstraintFeature);
1131   if (theAfterFeature.get()) {
1132     aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1133           SketchPlugin_ConstraintParallel::ID(),
1134           getFeatureResult(aBaseFeature),
1135           getFeatureResult(theAfterFeature));
1136     theCreatedFeatures.insert(aConstraintFeature);
1137   }
1138 #endif
1139   return anNewFeature;
1140 }
1141
1142 FeaturePtr SketchPlugin_Split::splitArc(FeaturePtr& theSplitFeature,
1143                                             FeaturePtr& theBaseFeatureModified,
1144                                             FeaturePtr& theAfterFeature,
1145                                             std::set<AttributePoint2DPtr>& thePoints,
1146                                             std::set<FeaturePtr>& theCreatedFeatures,
1147                  std::set<std::pair<AttributePtr, AttributePtr>>& theModifiedAttributes)
1148 {
1149   FeaturePtr anNewFeature;
1150
1151   std::set<FeaturePtr> aCreatedFeatures;
1152   FeaturePtr aConstraintFeature;
1153   theBaseFeatureModified = FeaturePtr(); // it will contain modified base feature
1154
1155   SketchPlugin_Sketch* aSketch = sketch();
1156   if (!aSketch)
1157     return anNewFeature;
1158
1159   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1160                                                            data()->attribute(SELECTED_OBJECT()));
1161   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1162   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
1163   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
1164   std::string aFeatureKind = aBaseFeature->getKind();
1165   if (aFeatureKind != SketchPlugin_Arc::ID())
1166     return anNewFeature;
1167
1168   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
1169     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
1170   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
1171     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
1172   AttributePoint2DPtr aStartPointAttrOfBase, anEndPointAttrOfBase;
1173   getFeaturePoints(aBaseFeature, aStartPointAttrOfBase, anEndPointAttrOfBase);
1174   if (!aStartPointAttrOfBase.get() && !anEndPointAttrOfBase.get()) {
1175     setError("Error: Feature has no start and end points.");
1176     return anNewFeature;
1177   }
1178
1179   arrangePointsOnArc(aBaseFeature, aStartPointAttrOfBase, anEndPointAttrOfBase,
1180                      aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1181 #ifdef DEBUG_SPLIT
1182   std::cout << "Arranged points (to build split between 1st and 2nd points:" << std::endl;
1183   std::cout << "Start point: " <<
1184     ModelGeomAlgo_Point2D::getPointAttributeInfo(aStartPointAttrOfBase) << std::endl;
1185   std::cout << "1st point:   " <<
1186     ModelGeomAlgo_Point2D::getPointAttributeInfo(aFirstPointAttrOfSplit) << std::endl;
1187   std::cout << "2nd point:   " <<
1188     ModelGeomAlgo_Point2D::getPointAttributeInfo(aSecondPointAttrOfSplit) << std::endl;
1189   std::cout << "End point:   " <<
1190     ModelGeomAlgo_Point2D::getPointAttributeInfo(anEndPointAttrOfBase) << std::endl;
1191 #endif
1192
1193   /// split feature
1194   theSplitFeature = createArcFeature(aBaseFeature, aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1195   theCreatedFeatures.insert(theSplitFeature);
1196
1197   // before split feature
1198   if (aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1199     theModifiedAttributes.insert(std::make_pair(aStartPointAttrOfBase,
1200                                   theSplitFeature->attribute(SketchPlugin_Arc::START_ID())));
1201   }
1202   else {
1203     theBaseFeatureModified = aBaseFeature; ///< use base feature to store all constraints here
1204     /// move end arc point to start of split
1205   }
1206
1207   // after split feature
1208   if (!aSecondPointAttrOfSplit->pnt()->isEqual(anEndPointAttrOfBase->pnt())) {
1209     FeaturePtr aFeature;
1210     if (!theBaseFeatureModified.get()) {
1211       aFeature = aBaseFeature; ///< use base feature to store all constraints here
1212       fillAttribute(aFeature->attribute(SketchPlugin_Arc::START_ID()), aSecondPointAttrOfSplit);
1213       aFeature->execute(); // to update result
1214     }
1215     else {
1216       aFeature = createArcFeature(aBaseFeature, aSecondPointAttrOfSplit, anEndPointAttrOfBase);
1217       theCreatedFeatures.insert(aFeature);
1218       theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1219                                                   aFeature->attribute(SketchPlugin_Arc::END_ID())));
1220       anNewFeature = aFeature;
1221     }
1222     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1223                      SketchPlugin_ConstraintCoincidence::ID(),
1224                      theSplitFeature->attribute(SketchPlugin_Arc::END_ID()),
1225                      aFeature->attribute(SketchPlugin_Arc::START_ID()));
1226     theCreatedFeatures.insert(aConstraintFeature);
1227
1228     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1229                                 (aFeature->attribute(SketchPlugin_Arc::START_ID())));
1230     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1231                                 (aFeature->attribute(SketchPlugin_Arc::END_ID())));
1232
1233     if (!theBaseFeatureModified.get())
1234       theBaseFeatureModified = aFeature;
1235     else
1236       theAfterFeature = aFeature;
1237   }
1238   else {
1239     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1240                                   (theSplitFeature->attribute(SketchPlugin_Arc::END_ID())));
1241     theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1242                                    theSplitFeature->attribute(SketchPlugin_Arc::END_ID())));
1243   }
1244   // base split, that is defined before split feature should be changed at end
1245   // (after the after feature creation). Otherwise modified value will be used in after feature
1246   // before split feature
1247   if (!aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1248     /// move end arc point to start of split
1249     fillAttribute(theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID()),
1250                                                     aFirstPointAttrOfSplit);
1251     theBaseFeatureModified->execute(); // to update result
1252     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1253                      SketchPlugin_ConstraintCoincidence::ID(),
1254                      theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID()),
1255                      theSplitFeature->attribute(SketchPlugin_Arc::START_ID()));
1256     theCreatedFeatures.insert(aConstraintFeature);
1257
1258     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1259                              (theBaseFeatureModified->attribute(SketchPlugin_Arc::START_ID())));
1260     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1261                                (theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID())));
1262   }
1263   else
1264     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1265                                        (theSplitFeature->attribute(SketchPlugin_Arc::START_ID())));
1266
1267   // additional constraints between split and base features
1268 #ifdef CREATE_CONSTRAINTS
1269   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1270           SketchPlugin_ConstraintEqual::ID(),
1271           getFeatureResult(aBaseFeature),
1272           getFeatureResult(theSplitFeature));
1273   theCreatedFeatures.insert(aConstraintFeature);
1274   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1275           SketchPlugin_ConstraintTangent::ID(),
1276           getFeatureResult(theSplitFeature),
1277           getFeatureResult(aBaseFeature));
1278   theCreatedFeatures.insert(aConstraintFeature);
1279   if (theAfterFeature.get()) {
1280     aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1281                 SketchPlugin_ConstraintEqual::ID(),
1282                 getFeatureResult(aBaseFeature),
1283                 getFeatureResult(theAfterFeature));
1284     theCreatedFeatures.insert(aConstraintFeature);
1285     aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1286                 SketchPlugin_ConstraintTangent::ID(),
1287                 getFeatureResult(theSplitFeature),
1288                 getFeatureResult(theAfterFeature));
1289     theCreatedFeatures.insert(aConstraintFeature);
1290   }
1291 #endif
1292   return anNewFeature;
1293 }
1294
1295 FeaturePtr SketchPlugin_Split::splitCircle(FeaturePtr& theSplitFeature,
1296                                                FeaturePtr& theBaseFeatureModified,
1297                                                FeaturePtr& theAfterFeature,
1298                                                std::set<AttributePoint2DPtr>& thePoints,
1299                                                std::set<FeaturePtr>& theCreatedFeatures,
1300                  std::set<std::pair<AttributePtr, AttributePtr>>& theModifiedAttributes)
1301 {
1302   FeaturePtr anNewFeature;
1303
1304   std::set<FeaturePtr> aCreatedFeatures;
1305   FeaturePtr aConstraintFeature;
1306   theBaseFeatureModified = FeaturePtr(); // it will contain modified base feature
1307
1308   SketchPlugin_Sketch* aSketch = sketch();
1309   if (!aSketch)
1310     return anNewFeature;
1311
1312   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1313                                                            data()->attribute(SELECTED_OBJECT()));
1314   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1315   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
1316   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
1317   std::string aFeatureKind = aBaseFeature->getKind();
1318   if (aFeatureKind != SketchPlugin_Circle::ID())
1319     return anNewFeature;
1320
1321   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
1322     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
1323   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
1324     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
1325
1326   /// split feature
1327   theSplitFeature =
1328     createArcFeature(aBaseFeature, aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1329   bool aSplitReversed = std::dynamic_pointer_cast<SketchPlugin_Arc>(theSplitFeature)->isReversed();
1330   theCreatedFeatures.insert(theSplitFeature);
1331
1332   /// base feature is a left part of the circle
1333   theBaseFeatureModified = createArcFeature(aBaseFeature,
1334     aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1335   anNewFeature = theBaseFeatureModified;
1336   std::dynamic_pointer_cast<SketchPlugin_Arc>(
1337     theBaseFeatureModified)->setReversed(!aSplitReversed);
1338   theBaseFeatureModified->execute();
1339
1340   theModifiedAttributes.insert(
1341     std::make_pair(aBaseFeature->attribute(SketchPlugin_Circle::CENTER_ID()),
1342                   theBaseFeatureModified->attribute(SketchPlugin_Arc::CENTER_ID())));
1343
1344   theCreatedFeatures.insert(theBaseFeatureModified);
1345
1346   thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1347                              (theBaseFeatureModified->attribute(SketchPlugin_Arc::START_ID())));
1348   thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1349                              (theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID())));
1350
1351   // additional constraints between split and base features
1352   aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1353                      SketchPlugin_ConstraintCoincidence::ID(),
1354                      theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID()),
1355                      theSplitFeature->attribute(SketchPlugin_Arc::END_ID()));
1356   theCreatedFeatures.insert(aConstraintFeature);
1357   aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1358                      SketchPlugin_ConstraintCoincidence::ID(),
1359                      theBaseFeatureModified->attribute(SketchPlugin_Arc::START_ID()),
1360                      theSplitFeature->attribute(SketchPlugin_Arc::START_ID()));
1361   theCreatedFeatures.insert(aConstraintFeature);
1362
1363 #ifdef CREATE_CONSTRAINTS
1364   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1365                      SketchPlugin_ConstraintTangent::ID(),
1366                      getFeatureResult(theSplitFeature),
1367                      getFeatureResult(theBaseFeatureModified));
1368   theCreatedFeatures.insert(aConstraintFeature);
1369 #endif
1370   return anNewFeature;
1371 }
1372
1373 void SketchPlugin_Split::arrangePointsOnLine(
1374     const AttributePoint2DPtr& theStartPointAttr,
1375     const AttributePoint2DPtr& theEndPointAttr,
1376     AttributePoint2DPtr& theFirstPointAttr,
1377     AttributePoint2DPtr& theLastPointAttr) const
1378 {
1379   // if first point is closer to last point, swap first and last values
1380   if (theStartPointAttr->pnt()->distance(theFirstPointAttr->pnt()) >
1381       theStartPointAttr->pnt()->distance(theLastPointAttr->pnt())) {
1382     AttributePoint2DPtr aTmpPoint = theFirstPointAttr;
1383     theFirstPointAttr = theLastPointAttr;
1384     theLastPointAttr = aTmpPoint;
1385   }
1386 }
1387
1388 void SketchPlugin_Split::arrangePointsOnArc(
1389     const FeaturePtr& theArc,
1390     const std::shared_ptr<GeomDataAPI_Point2D>& theStartPointAttr,
1391     const std::shared_ptr<GeomDataAPI_Point2D>& theEndPointAttr,
1392     std::shared_ptr<GeomDataAPI_Point2D>& theFirstPointAttr,
1393     std::shared_ptr<GeomDataAPI_Point2D>& theSecondPointAttr) const
1394 {
1395   static const double anAngleTol = 1.e-12;
1396
1397   std::shared_ptr<GeomAPI_Pnt2d> aCenter = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
1398       theArc->attribute(SketchPlugin_Arc::CENTER_ID()))->pnt();
1399   bool isReversed = theArc->boolean(SketchPlugin_Arc::REVERSED_ID())->value();
1400
1401   // collect directions to each point
1402   std::shared_ptr<GeomAPI_Dir2d> aStartDir(
1403       new GeomAPI_Dir2d(theStartPointAttr->pnt()->xy()->decreased(aCenter->xy())));
1404   std::shared_ptr<GeomAPI_Dir2d> aFirstPtDir(
1405       new GeomAPI_Dir2d(theFirstPointAttr->pnt()->xy()->decreased(aCenter->xy())));
1406   std::shared_ptr<GeomAPI_Dir2d> aSecondPtDir(
1407       new GeomAPI_Dir2d(theSecondPointAttr->pnt()->xy()->decreased(aCenter->xy())));
1408
1409   // sort points by their angular values
1410   double aFirstPtAngle = aStartDir->angle(aFirstPtDir);
1411   double aSecondPtAngle = aStartDir->angle(aSecondPtDir);
1412   double aPeriod = isReversed ? -2.0 * PI : 2.0 * PI;
1413   if (fabs(aFirstPtAngle) > anAngleTol && isReversed == (aFirstPtAngle > 0.))
1414     aFirstPtAngle += aPeriod;
1415   if (fabs(aSecondPtAngle) > anAngleTol && isReversed == (aSecondPtAngle > 0.))
1416     aSecondPtAngle += aPeriod;
1417
1418   if (fabs(aFirstPtAngle) > fabs(aSecondPtAngle)) {
1419     AttributePoint2DPtr aTmpPoint = theFirstPointAttr;
1420     theFirstPointAttr = theSecondPointAttr;
1421     theSecondPointAttr = aTmpPoint;
1422   }
1423 }
1424
1425 void SketchPlugin_Split::fillAttribute(const AttributePtr& theModifiedAttribute,
1426                                                  const AttributePtr& theSourceAttribute)
1427 {
1428   std::string anAttributeType = theModifiedAttribute->attributeType();
1429   if (anAttributeType == GeomDataAPI_Point2D::typeId()) {
1430     AttributePoint2DPtr aModifiedAttribute = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
1431                                               theModifiedAttribute);
1432     AttributePoint2DPtr aSourceAttribute = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
1433                                               theSourceAttribute);
1434
1435     if (aModifiedAttribute.get() && aSourceAttribute.get())
1436       aModifiedAttribute->setValue(aSourceAttribute->pnt());
1437   }
1438   else if (anAttributeType == ModelAPI_AttributeBoolean::typeId()) {
1439     AttributeBooleanPtr aModifiedAttribute = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
1440                                               theModifiedAttribute);
1441     AttributeBooleanPtr aSourceAttribute = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
1442                                               theSourceAttribute);
1443
1444     if (aModifiedAttribute.get() && aSourceAttribute.get())
1445       aModifiedAttribute->setValue(aSourceAttribute->value());
1446   }
1447 }
1448
1449 FeaturePtr SketchPlugin_Split::createLineFeature(const FeaturePtr& theBaseFeature,
1450                                                            const AttributePtr& theFirstPointAttr,
1451                                                            const AttributePtr& theSecondPointAttr)
1452 {
1453   FeaturePtr aFeature;
1454   SketchPlugin_Sketch* aSketch = sketch();
1455   if (!aSketch || !theBaseFeature.get())
1456     return aFeature;
1457
1458   aFeature = aSketch->addFeature(SketchPlugin_Line::ID());
1459
1460   fillAttribute(aFeature->attribute(SketchPlugin_Line::START_ID()), theFirstPointAttr);
1461   fillAttribute(aFeature->attribute(SketchPlugin_Line::END_ID()), theSecondPointAttr);
1462
1463   fillAttribute(aFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()),
1464                 theBaseFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()));
1465
1466   aFeature->execute(); // to obtain result
1467
1468   return aFeature;
1469 }
1470
1471 FeaturePtr SketchPlugin_Split::createArcFeature(const FeaturePtr& theBaseFeature,
1472                                                           const AttributePtr& theFirstPointAttr,
1473                                                           const AttributePtr& theSecondPointAttr)
1474 {
1475   FeaturePtr aFeature;
1476   SketchPlugin_Sketch* aSketch = sketch();
1477   if (!aSketch || !theBaseFeature.get())
1478     return aFeature;
1479
1480   std::string aCenterAttributeId;
1481   if (theBaseFeature->getKind() == SketchPlugin_Arc::ID())
1482     aCenterAttributeId = SketchPlugin_Arc::CENTER_ID();
1483   else if (theBaseFeature->getKind() == SketchPlugin_Circle::ID())
1484     aCenterAttributeId = SketchPlugin_Circle::CENTER_ID();
1485
1486   if (aCenterAttributeId.empty())
1487     return aFeature;
1488
1489   aFeature = aSketch->addFeature(SketchPlugin_Arc::ID());
1490   // update fillet arc: make the arc correct for sure, so, it is not needed to process
1491   // the "attribute updated"
1492   // by arc; moreover, it may cause cyclicity in hte mechanism of updater
1493   bool aWasBlocked = aFeature->data()->blockSendAttributeUpdated(true);
1494
1495   fillAttribute(aFeature->attribute(SketchPlugin_Arc::CENTER_ID()),
1496                 theBaseFeature->attribute(aCenterAttributeId));
1497   fillAttribute(aFeature->attribute(SketchPlugin_Arc::START_ID()), theFirstPointAttr);
1498   fillAttribute(aFeature->attribute(SketchPlugin_Arc::END_ID()), theSecondPointAttr);
1499
1500   fillAttribute(aFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()),
1501                 theBaseFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()));
1502
1503   /// fill referersed state of created arc as it is on the base arc
1504   if (theBaseFeature->getKind() == SketchPlugin_Arc::ID()) {
1505     bool aReversed = theBaseFeature->boolean(SketchPlugin_Arc::REVERSED_ID())->value();
1506     aFeature->boolean(SketchPlugin_Arc::REVERSED_ID())->setValue(aReversed);
1507   }
1508   aFeature->data()->blockSendAttributeUpdated(aWasBlocked);
1509   aFeature->execute(); // to obtain result
1510
1511   return aFeature;
1512 }
1513
1514 void SketchPlugin_Split::updateFeaturesAfterSplit(
1515                                                    const std::set<FeaturePtr>& theFeaturesToUpdate)
1516 {
1517   std::set<FeaturePtr>::const_iterator anIt = theFeaturesToUpdate.begin(),
1518                                        aLast = theFeaturesToUpdate.end();
1519   for (; anIt != aLast; anIt++) {
1520     FeaturePtr aRefFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(*anIt);
1521     std::string aRefFeatureKind = aRefFeature->getKind();
1522     if (aRefFeatureKind == SketchPlugin_ConstraintLength::ID()) {
1523       std::shared_ptr<SketchPlugin_ConstraintLength> aLenghtFeature =
1524                               std::dynamic_pointer_cast<SketchPlugin_ConstraintLength>(*anIt);
1525       if (aLenghtFeature.get()) {
1526         std::shared_ptr<ModelAPI_AttributeDouble> aValueAttr = std::dynamic_pointer_cast<
1527             ModelAPI_AttributeDouble>(aLenghtFeature->attribute(SketchPlugin_Constraint::VALUE()));
1528         double aValue;
1529         if (aLenghtFeature->computeLenghtValue(aValue) && aValueAttr.get())
1530           aValueAttr->setValue(aValue);
1531       }
1532     }
1533   }
1534 }
1535
1536 std::shared_ptr<ModelAPI_Result> SketchPlugin_Split::getFeatureResult(
1537                                     const std::shared_ptr<ModelAPI_Feature>& theFeature)
1538 {
1539   std::shared_ptr<ModelAPI_Result> aResult;
1540
1541   std::string aFeatureKind = theFeature->getKind();
1542   if (aFeatureKind == SketchPlugin_Line::ID())
1543     aResult = theFeature->firstResult();
1544   else if (aFeatureKind == SketchPlugin_Arc::ID())
1545     aResult = theFeature->lastResult();
1546   else if (aFeatureKind == SketchPlugin_Circle::ID())
1547     aResult = theFeature->lastResult();
1548
1549   return aResult;
1550 }
1551
1552 #ifdef _DEBUG
1553 std::set<std::shared_ptr<ModelAPI_Attribute> > SketchPlugin_Split::getEdgeAttributes(
1554                                            const std::shared_ptr<ModelAPI_Feature>& theFeature)
1555 {
1556   std::set<std::shared_ptr<ModelAPI_Attribute> > anAttributes;
1557
1558   std::string aFeatureKind = theFeature->getKind();
1559   if (aFeatureKind == SketchPlugin_Line::ID()) {
1560     anAttributes.insert(theFeature->attribute(SketchPlugin_Line::START_ID()));
1561     anAttributes.insert(theFeature->attribute(SketchPlugin_Line::END_ID()));
1562   }
1563   else if (aFeatureKind == SketchPlugin_Arc::ID()) {
1564     anAttributes.insert(theFeature->attribute(SketchPlugin_Arc::START_ID()));
1565     anAttributes.insert(theFeature->attribute(SketchPlugin_Arc::END_ID()));
1566   }
1567   else if (aFeatureKind == SketchPlugin_Circle::ID()) {
1568   }
1569
1570   return anAttributes;
1571 }
1572 #endif
1573
1574 std::shared_ptr<GeomDataAPI_Point2D> SketchPlugin_Split::getPointAttribute
1575                                                               (const bool isFirstAttribute)
1576 {
1577   std::shared_ptr<GeomDataAPI_Point2D> anAttribute;
1578
1579   GeomShapePtr aSelectedShape = getSubShape(SELECTED_OBJECT(), SELECTED_POINT());
1580   if (!aSelectedShape.get())
1581     return anAttribute;
1582
1583   if (aSelectedShape->shapeType() != GeomAPI_Shape::EDGE)
1584     return anAttribute;
1585
1586   AttributeReferencePtr anObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1587                                                        data()->attribute(SELECTED_OBJECT()));
1588   ObjectPtr aBaseObject = anObjectAttr->value();
1589   if (!aBaseObject.get())
1590     return anAttribute;
1591
1592   std::shared_ptr<GeomAPI_Edge> anEdge(new GeomAPI_Edge(aSelectedShape));
1593
1594   std::shared_ptr<GeomAPI_Pnt> aFirstPnt = anEdge->firstPoint();
1595   std::shared_ptr<GeomAPI_Pnt> aLastPnt = anEdge->lastPoint();
1596
1597   std::shared_ptr<GeomDataAPI_Point2D> aFirstPointAttr, aLastPointAttr;
1598   /// find the points in feature attributes
1599   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObject);
1600   std::list<AttributePtr> a2DPointAttributes = aBaseFeature->data()->attributes(
1601                                                     GeomDataAPI_Point2D::typeId());
1602   std::list<AttributePtr>::const_iterator anIt = a2DPointAttributes.begin(),
1603                                           aLast = a2DPointAttributes.end();
1604   for (; anIt != aLast; anIt++) {
1605     std::shared_ptr<GeomDataAPI_Point2D> anAttributePoint =
1606                                   std::dynamic_pointer_cast<GeomDataAPI_Point2D>(*anIt);
1607     std::shared_ptr<GeomAPI_Pnt2d> aPoint2D = anAttributePoint->pnt();
1608     std::shared_ptr<GeomAPI_Pnt> aPoint3D = sketch()->to3D(aPoint2D->x(), aPoint2D->y());
1609     if (aFirstPnt->isEqual(aPoint3D))
1610       aFirstPointAttr = anAttributePoint;
1611     else if (aLastPnt->isEqual(aPoint3D))
1612       aLastPointAttr = anAttributePoint;
1613   }
1614
1615   /// find the points in coincident features
1616   PntToAttributesMap aRefAttributes = myCashedReferences[aBaseObject];
1617   PntToAttributesMap::const_iterator
1618     aRIt = aRefAttributes.begin(), aRLast = aRefAttributes.end();
1619   for (; aRIt != aRLast; aRIt++) {
1620     std::shared_ptr<GeomDataAPI_Point2D> anAttribute = aRIt->first;
1621     std::shared_ptr<GeomAPI_Pnt> aPoint = aRIt->second;
1622     if (!aFirstPointAttr.get() && aFirstPnt->isEqual(aPoint))
1623       aFirstPointAttr = anAttribute;
1624     if (!aLastPointAttr.get() && aLastPnt->isEqual(aPoint))
1625       aLastPointAttr = anAttribute;
1626     if (aFirstPointAttr.get() && aLastPointAttr.get())
1627       break;
1628   }
1629   if (!aFirstPointAttr.get() || !aLastPointAttr)
1630     return anAttribute;
1631
1632   return isFirstAttribute ? aFirstPointAttr : aLastPointAttr;
1633 }
1634
1635 #ifdef _DEBUG
1636 std::string SketchPlugin_Split::getFeatureInfo(const std::shared_ptr<ModelAPI_Feature>& theFeature,
1637                                                const bool isUseAttributesInfo)
1638 {
1639   std::string anInfo;
1640   if (!theFeature.get()) {
1641     return "none";
1642   }
1643
1644   if (theFeature->data()->isValid())
1645     anInfo.append(theFeature->data()->name().c_str());
1646
1647   if (isUseAttributesInfo) {
1648     std::string aPointsInfo = ModelGeomAlgo_Point2D::getPontAttributesInfo(theFeature,
1649                                                              getEdgeAttributes(theFeature));
1650     /// processing of feature with point 2d attributes, like line, arc, circle
1651     if (!aPointsInfo.empty()) {
1652       anInfo += ": ";
1653       anInfo += "\n";
1654       anInfo += aPointsInfo;
1655     }
1656     else { /// process constraint coincidence, find points in ref attr attributes
1657       std::list<AttributePtr> anAttrs = theFeature->data()->attributes(
1658                                                        ModelAPI_AttributeRefAttr::typeId());
1659       std::list<AttributePtr>::const_iterator anIt = anAttrs.begin(), aLast = anAttrs.end();
1660       std::string anAttributesInfo;
1661       for(; anIt != aLast; anIt++) {
1662         if (!anAttributesInfo.empty()) {
1663           anAttributesInfo.append(", ");
1664           anAttributesInfo += "\n";
1665         }
1666         AttributePtr anAttr = *anIt;
1667         std::string aValue = "not defined";
1668         std::string aType = anAttr->attributeType();
1669         if (aType == ModelAPI_AttributeRefAttr::typeId()) {
1670           std::shared_ptr<ModelAPI_AttributeRefAttr> aRefAttr =
1671                              std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(anAttr);
1672           if (aRefAttr.get()) {
1673             if (aRefAttr->isObject()) {
1674               FeaturePtr aFeature = ModelAPI_Feature::feature(aRefAttr->object());
1675               aValue = "<object:>" + getFeatureInfo(aFeature, false);
1676             }
1677             else {
1678               AttributePtr anAttribute = aRefAttr->attr();
1679               if (anAttribute.get()) {
1680                 FeaturePtr aFeature = ModelAPI_Feature::feature(anAttribute->owner());
1681                 aValue = "<attr:>" + ModelGeomAlgo_Point2D::getPointAttributeInfo(anAttribute) +
1682                          " [" + getFeatureInfo(aFeature, false) + "]";
1683               }
1684             }
1685           }
1686         }
1687         anAttributesInfo.append("    " + anAttr->id() + ": " + aValue);
1688       }
1689       if (!anAttributesInfo.empty())
1690         anInfo = anInfo + "\n" + anAttributesInfo;
1691     }
1692   }
1693   return anInfo;
1694 }
1695 #endif