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_TRIM_METHODS
 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_TRIM_METHODS
 427   std::cout << "SketchPlugin_Trim::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   #ifdef DEBUG_TRIM_METHODS
 466         if (!aSelectedShape.get())
 467           std::cout << "Set empty selected object" << std::endl;
 468         else
 469           std::cout << "Set shape with ShapeType: " << aSelectedShape->shapeTypeStr() << std::endl;
 470   #endif
 471         aFilledAttributeName = SELECTED_OBJECT();
 472       }
 473     }
 474   }
 475   return aFilledAttributeName;
 476 }
 477
 478 AISObjectPtr SketchPlugin_Split::getAISObject(AISObjectPtr thePrevious)
 479 {
 480   /*AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 481                                            data()->attribute(SketchPlugin_Constraint::VALUE()));
 482   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
 483
 484   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointOfRefAttr(
 485                                         data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
 486   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointOfRefAttr(
 487                                         data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
 488
 489   if (aBaseObjectAttr->isInitialized() && aBaseFeature.get() &&
 490       aFirstPointAttrOfSplit->isInitialized() &&
 491       aSecondPointAttrOfSplit->isInitialized()) {
 492
 493     ResultPtr aResult = getFeatureResult(aBaseFeature);
 494     GeomShapePtr aBaseShape = aResult->shape();
 495     std::list<std::shared_ptr<GeomAPI_Pnt> > aPoints;
 496
 497     std::shared_ptr<GeomAPI_Pnt2d> aStartPnt2d = aFirstPointAttrOfSplit->pnt();
 498     std::shared_ptr<GeomAPI_Pnt> aStartPoint = sketch()->to3D(aStartPnt2d->x(), aStartPnt2d->y());
 499     aPoints.push_back(aStartPoint);
 500
 501     std::shared_ptr<GeomAPI_Pnt2d> aSecondPnt2d = aSecondPointAttrOfSplit->pnt();
 502     std::shared_ptr<GeomAPI_Pnt> aSecondPoint =
 503       sketch()->to3D(aSecondPnt2d->x(), aSecondPnt2d->y());
 504     aPoints.push_back(aSecondPoint);
 505
 506     std::set<std::shared_ptr<GeomAPI_Shape> > aSplitShapes;
 507
 508     GeomAlgoAPI_ShapeTools::splitShape_p(aBaseShape, aPoints, aSplitShapes);
 509     std::shared_ptr<GeomAPI_Shape> aShape =
 510       GeomAlgoAPI_ShapeTools::findShape(aPoints, aSplitShapes);
 511
 512     AISObjectPtr anAIS = thePrevious;
 513     if (aShape) {
 514       if (!anAIS)
 515         anAIS = AISObjectPtr(new GeomAPI_AISObject);
 516       anAIS->createShape(aShape);
 517       std::shared_ptr<ModelAPI_AttributeBoolean> anAuxiliaryAttr =
 518              aBaseFeature->data()->boolean(SketchPlugin_SketchEntity::AUXILIARY_ID());
 519
 520       bool isConstruction = anAuxiliaryAttr.get() != NULL && anAuxiliaryAttr->value();
 521
 522       std::vector<int> aColor;
 523       double aWidth = SketchPlugin_SketchEntity::SKETCH_LINE_WIDTH();
 524       int aLineStyle = SketchPlugin_SketchEntity::SKETCH_LINE_STYLE();
 525       if (isConstruction) {
 526         aColor = Config_PropManager::color("Visualization", "sketch_auxiliary_color");
 527         aWidth = SketchPlugin_SketchEntity::SKETCH_LINE_WIDTH_AUXILIARY();
 528         aLineStyle = SketchPlugin_SketchEntity::SKETCH_LINE_STYLE_AUXILIARY();
 529       }
 530       else {
 531         aColor = Config_PropManager::color("Visualization", "sketch_entity_color");
 532       }
 533       anAIS->setColor(aColor[0], aColor[1], aColor[2]);
 534       anAIS->setWidth(aWidth + 1);
 535       anAIS->setLineStyle(aLineStyle);
 536     }
 537     return anAIS;
 538   }
 539   return AISObjectPtr();*/
 540 #ifdef DEBUG_TRIM_METHODS
 541   std::cout << "SketchPlugin_Trim::getAISObject: " << data()->name() << std::endl;
 542 #endif
 543
 544   AISObjectPtr anAIS = thePrevious;
 545
 546   std::list<std::shared_ptr<GeomAPI_Shape> > aShapes;
 547   GeomShapePtr aPreviewShape = getSubShape(PREVIEW_OBJECT(), PREVIEW_POINT());
 548   if (aPreviewShape.get())
 549     aShapes.push_back(aPreviewShape);
 550   GeomShapePtr aSelectedShape = getSubShape(SELECTED_OBJECT(), SELECTED_POINT());
 551   if (aSelectedShape.get())
 552     aShapes.push_back(aSelectedShape);
 553
 554   if (aShapes.empty())
 555     return AISObjectPtr();
 556
 557   GeomShapePtr aBaseShape = GeomAlgoAPI_CompoundBuilder::compound(aShapes);
 558   if (!aBaseShape.get())
 559     return AISObjectPtr();
 560
 561   if (aBaseShape.get()) {
 562     if (!anAIS)
 563       anAIS = AISObjectPtr(new GeomAPI_AISObject);
 564     anAIS->createShape(aBaseShape);
 565
 566     std::vector<int> aColor;
 567     aColor = Config_PropManager::color("Visualization", "operation_remove_feature_color");
 568     double aWidth = SketchPlugin_SketchEntity::SKETCH_LINE_WIDTH();
 569     int aLineStyle = SketchPlugin_SketchEntity::SKETCH_LINE_STYLE();
 570     anAIS->setColor(aColor[0], aColor[1], aColor[2]);
 571     // width when there is not base object should be extened in several points
 572     // in order to see this preview over highlight
 573     anAIS->setWidth(aWidth+4);
 574     anAIS->setLineStyle(aLineStyle);
 575   }
 576   else
 577     anAIS = AISObjectPtr();
 578   return anAIS;
 579 }
 580
 581 //********************************************************************
 582 void SketchPlugin_Split::fillObjectShapes(const ObjectPtr& theObject,
 583                                           const ObjectPtr& theSketch)
 584 {
 585   std::set<std::shared_ptr<GeomAPI_Shape> > aShapes;
 586   std::map<std::shared_ptr<GeomDataAPI_Point2D>, std::shared_ptr<GeomAPI_Pnt> > aPointToAttributes;
 587   std::set<std::shared_ptr<GeomDataAPI_Point2D> > aRefAttributes;
 588   // current feature
 589   FeaturePtr aFeature = ModelAPI_Feature::feature(theObject);
 590   // edges on feature
 591   std::set<ResultPtr> anEdgeResults;
 592   ModelGeomAlgo_Shape::shapesOfType(aFeature, GeomAPI_Shape::EDGE, anEdgeResults);
 593   if (!anEdgeResults.empty()) {
 594     GeomShapePtr aFeatureShape = (*anEdgeResults.begin())->shape();
 595
 596     // coincidences to the feature
 597     ModelGeomAlgo_Point2D::getPointsOfReference(aFeature, SketchPlugin_ConstraintCoincidence::ID(),
 598                          aRefAttributes, SketchPlugin_Point::ID(), SketchPlugin_Point::COORD_ID());
 599     // layed on feature coincidences to divide it on several shapes
 600     std::shared_ptr<ModelAPI_Data> aData = theSketch->data();
 601     std::shared_ptr<GeomDataAPI_Point> aC = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 602         aData->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 603     std::shared_ptr<GeomDataAPI_Dir> aX = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 604         aData->attribute(SketchPlugin_Sketch::DIRX_ID()));
 605     std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 606         aData->attribute(SketchPlugin_Sketch::NORM_ID()));
 607     std::shared_ptr<GeomAPI_Dir> aY(new GeomAPI_Dir(aNorm->dir()->cross(aX->dir())));
 608     std::list<std::shared_ptr<GeomAPI_Pnt> > aPoints;
 609     ModelGeomAlgo_Point2D::getPointsInsideShape_p(aFeatureShape, aRefAttributes, aC->pnt(),
 610                                                 aX->dir(), aY, aPoints, aPointToAttributes);
 611
 612     GeomAlgoAPI_ShapeTools::splitShape_p(aFeatureShape, aPoints, aShapes);
 613   }
 614   myCashedShapes[theObject] = aShapes;
 615   myCashedReferences[theObject] = aPointToAttributes;
 616 }
 617
 618 GeomShapePtr SketchPlugin_Split::getSubShape(const std::string& theObjectAttributeId,
 619                                              const std::string& thePointAttributeId)
 620 {
 621   GeomShapePtr aBaseShape;
 622
 623   AttributeReferencePtr anObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 624                                                        data()->attribute(theObjectAttributeId));
 625   ObjectPtr aBaseObject = anObjectAttr->value();
 626   if (!aBaseObject.get())
 627     return aBaseShape;
 628
 629   // point on feature
 630   AttributePoint2DPtr aPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 631                                            data()->attribute(thePointAttributeId));
 632   std::shared_ptr<GeomAPI_Pnt2d> anAttributePnt2d = aPointAttr->pnt();
 633   std::shared_ptr<GeomAPI_Pnt> anAttributePnt = sketch()->to3D(anAttributePnt2d->x(),
 634                                                                anAttributePnt2d->y());
 635
 636 #ifdef TRIM_SHAPE
 637   if (myCashedShapes.find(aBaseObject) == myCashedShapes.end())
 638     fillObjectShapes(aBaseObject, sketch()->data()->owner(), myCashedShapes, myObjectToPoints);
 639
 640   const std::set<GeomShapePtr>& aShapes = myCashedShapes[aBaseObject];
 641   if (!aShapes.empty()) {
 642     std::set<GeomShapePtr>::const_iterator anIt = aShapes.begin(), aLast = aShapes.end();
 643     for (; anIt != aLast; anIt++) {
 644       GeomShapePtr aShape = *anIt;
 645       std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 646       if (ModelGeomAlgo_Point2D::isPointOnEdge(aShape, anAttributePnt, aProjectedPoint))
 647         aBaseShape = aShape;
 648     }
 649   }
 650 #else
 651   if (myCashedShapes.find(aBaseObject) == myCashedShapes.end())
 652     fillObjectShapes(aBaseObject, sketch()->data()->owner());
 653
 654   std::shared_ptr<GeomAPI_Pnt> aStartPoint;
 655   std::shared_ptr<GeomAPI_Pnt> aSecondPoint;
 656   const std::set<GeomShapePtr>& aShapes = myCashedShapes[aBaseObject];
 657   std::set<GeomShapePtr>::const_iterator anIt = aShapes.begin(), aLast = aShapes.end();
 658   for (; anIt != aLast; anIt++) {
 659     GeomShapePtr aCurrentShape = *anIt;
 660     std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 661     if (ModelGeomAlgo_Point2D::isPointOnEdge(aCurrentShape, anAttributePnt, aProjectedPoint)) {
 662       if (aCurrentShape->shapeType() == GeomAPI_Shape::EDGE) {
 663         std::shared_ptr<GeomAPI_Edge> anEdge(new GeomAPI_Edge(aCurrentShape));
 664         aStartPoint = anEdge->firstPoint();
 665         aSecondPoint = anEdge->lastPoint();
 666       }
 667       break;
 668     }
 669   }
 670
 671   if (!aStartPoint.get() || !aSecondPoint.get())
 672     return aBaseShape;
 673
 674   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 675   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
 676   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObject/*aBaseObjectAttr->value()*/);
 677
 678   //AttributePoint2DPtr aFirstPointAttrOfSplit = getPointOfRefAttr(
 679   //                                      data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
 680   //AttributePoint2DPtr aSecondPointAttrOfSplit = getPointOfRefAttr(
 681   //                                      data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
 682   if (anObjectAttr->isInitialized() && aBaseFeature.get() && aPointAttr->isInitialized()) {
 683       //aFirstPointAttrOfSplit->isInitialized() &&
 684       //aSecondPointAttrOfSplit->isInitialized()) {
 685     ResultPtr aResult = getFeatureResult(aBaseFeature);
 686     GeomShapePtr aResultShape = aResult->shape();
 687     std::list<std::shared_ptr<GeomAPI_Pnt> > aPoints;
 688
 689     //std::shared_ptr<GeomAPI_Pnt2d> aStartPnt2d = aFirstPointAttrOfSplit->pnt();
 690     //std::shared_ptr<GeomAPI_Pnt> aStartPoint = sketch()->to3D(aStartPnt2d->x(), aStartPnt2d->y());
 691     aPoints.push_back(aStartPoint);
 692
 693     //std::shared_ptr<GeomAPI_Pnt2d> aSecondPnt2d = aSecondPointAttrOfSplit->pnt();
 694     //std::shared_ptr<GeomAPI_Pnt> aSecondPoint =
 695     //  sketch()->to3D(aSecondPnt2d->x(), aSecondPnt2d->y());
 696     aPoints.push_back(aSecondPoint);
 697
 698     std::set<std::shared_ptr<GeomAPI_Shape> > aSplitShapes;
 699     GeomAlgoAPI_ShapeTools::splitShape_p(aResultShape, aPoints, aSplitShapes);
 700     aBaseShape = GeomAlgoAPI_ShapeTools::findShape(aPoints, aSplitShapes);
 701 #endif
 702   }
 703   return aBaseShape;
 704 }
 705
 706 std::shared_ptr<GeomDataAPI_Point2D> SketchPlugin_Split::getPointOfRefAttr(
 707                                                       const AttributePtr& theAttribute)
 708 {
 709   AttributePoint2DPtr aPointAttribute;
 710
 711   if (theAttribute->attributeType() == ModelAPI_AttributeRefAttr::typeId()) {
 712     AttributeRefAttrPtr aRefAttr =
 713       std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(theAttribute);
 714     if (aRefAttr.get() && aRefAttr->isInitialized()) {
 715       AttributePtr anAttribute = aRefAttr->attr();
 716       if (anAttribute.get() && anAttribute->attributeType() == GeomDataAPI_Point2D::typeId())
 717         aPointAttribute = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttribute);
 718     }
 719   }
 720   return aPointAttribute;
 721 }
 722
 723 void SketchPlugin_Split::getFeaturePoints(const FeaturePtr& theFeature,
 724                                                     AttributePoint2DPtr& theStartPointAttr,
 725                                                     AttributePoint2DPtr& theEndPointAttr)
 726 {
 727   std::string aFeatureKind = theFeature->getKind();
 728   std::string aStartAttributeName, anEndAttributeName;
 729   if (aFeatureKind == SketchPlugin_Line::ID()) {
 730     aStartAttributeName = SketchPlugin_Line::START_ID();
 731     anEndAttributeName = SketchPlugin_Line::END_ID();
 732   }
 733   else if (aFeatureKind == SketchPlugin_Arc::ID()) {
 734     aStartAttributeName = SketchPlugin_Arc::START_ID();
 735     anEndAttributeName = SketchPlugin_Arc::END_ID();
 736   }
 737   if (!aStartAttributeName.empty() && !anEndAttributeName.empty()) {
 738     theStartPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 739                                          theFeature->attribute(aStartAttributeName));
 740     theEndPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 741                                          theFeature->attribute(anEndAttributeName));
 742   }
 743 }
 744
 745 void SketchPlugin_Split::getConstraints(std::set<FeaturePtr>& theFeaturesToDelete,
 746                                     std::set<FeaturePtr>& theFeaturesToUpdate,
 747                                     std::map<FeaturePtr, IdToPointPair>& theCoincidenceToFeature)
 748 {
 749   std::shared_ptr<ModelAPI_Data> aData = data();
 750
 751   // Check the base objects are initialized.
 752   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 753                                                            data()->attribute(SELECTED_OBJECT()));
 754   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
 755   //                                          aData->attribute(SketchPlugin_Constraint::VALUE()));
 756   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
 757   ResultPtr aBaseFeatureResult = getFeatureResult(aBaseFeature);
 758
 759   std::set<AttributePtr> aRefsList = aBaseFeatureResult->data()->refsToMe();
 760   std::set<AttributePtr> aFRefsList = aBaseFeature->data()->refsToMe();
 761   aRefsList.insert(aFRefsList.begin(), aFRefsList.end());
 762
 763   std::set<AttributePtr>::const_iterator aIt;
 764   for (aIt = aRefsList.cbegin(); aIt != aRefsList.cend(); ++aIt) {
 765     std::shared_ptr<ModelAPI_Attribute> aAttr = (*aIt);
 766     FeaturePtr aRefFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aAttr->owner());
 767     std::string aRefFeatureKind = aRefFeature->getKind();
 768     if (aRefFeatureKind == SketchPlugin_ConstraintMirror::ID() ||
 769         aRefFeatureKind == SketchPlugin_MultiRotation::ID() ||
 770         aRefFeatureKind == SketchPlugin_MultiTranslation::ID() ||
 771         aRefFeatureKind == SketchPlugin_ConstraintMiddle::ID())
 772       theFeaturesToDelete.insert(aRefFeature);
 773     else if (aRefFeatureKind == SketchPlugin_ConstraintLength::ID())
 774       theFeaturesToUpdate.insert(aRefFeature);
 775     else if (aRefFeatureKind == SketchPlugin_ConstraintCoincidence::ID()) {
 776       std::string anAttributeToBeModified;
 777       AttributePoint2DPtr aCoincidentPoint;
 778       AttributeRefAttrPtr anAttrA = aRefFeature->refattr(SketchPlugin_Constraint::ENTITY_A());
 779       AttributeRefAttrPtr anAttrB = aRefFeature->refattr(SketchPlugin_Constraint::ENTITY_B());
 780       bool isToFeature = false;
 781       if (anAttrA->isObject() || anAttrB->isObject()) { /// coincidence to base feature
 782         FeaturePtr aFeature = anAttrA->isObject() ? ModelAPI_Feature::feature(anAttrA->object())
 783                                                   : FeaturePtr();
 784         isToFeature = aFeature.get() && aFeature == aBaseFeature;
 785         anAttributeToBeModified = anAttrA->id();
 786         if (!isToFeature) {
 787           aFeature = anAttrB->isObject() ? ModelAPI_Feature::feature(anAttrB->object())
 788                                          : FeaturePtr();
 789           isToFeature = aFeature.get() && aFeature == aBaseFeature;
 790           anAttributeToBeModified = anAttrB->id();
 791         }
 792         if (isToFeature)
 793           aCoincidentPoint = SketchPlugin_ConstraintCoincidence::getPoint(aRefFeature);
 794       }
 795       if (!isToFeature) { /// coincidence to point on base feature
 796         AttributePtr anAttribute;
 797
 798         if (!anAttrA->isObject()) {
 799           AttributePtr aCurAttribute = anAttrA->attr();
 800           if (aCurAttribute.get()) {
 801             FeaturePtr aCurFeature = ModelAPI_Feature::feature(aCurAttribute->owner());
 802             if (aCurFeature.get() && aCurFeature == aBaseFeature) {
 803               anAttribute = anAttrB->attr();
 804               anAttributeToBeModified = anAttrA->id();
 805             }
 806           }
 807         }
 808         if (!anAttribute.get() && !anAttrB->isObject()) {
 809           AttributePtr aCurAttribute = anAttrB->attr();
 810           if (aCurAttribute.get()) {
 811             FeaturePtr aCurFeature = ModelAPI_Feature::feature(aCurAttribute->owner());
 812             if (aCurFeature.get() && aCurFeature == aBaseFeature) {
 813               anAttribute = anAttrA->attr();
 814               anAttributeToBeModified = anAttrB->id();
 815             }
 816           }
 817         }
 818         if (anAttribute.get())
 819           aCoincidentPoint = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttribute);
 820       }
 821       if (aCoincidentPoint.get() && isToFeature)
 822         theCoincidenceToFeature[aRefFeature] = std::make_pair(anAttributeToBeModified,
 823                                                               aCoincidentPoint);
 824     }
 825   }
 826 }
 827
 828 void SketchPlugin_Split::getRefAttributes(const FeaturePtr& theFeature,
 829                                     std::map<AttributePtr, std::list<AttributePtr> >& theRefs,
 830                                     std::list<AttributePtr>& theRefsToFeature)
 831 {
 832   theRefs.clear();
 833
 834   std::list<AttributePtr> aPointAttributes =
 835     theFeature->data()->attributes(GeomDataAPI_Point2D::typeId());
 836   std::set<AttributePtr> aPointAttributesSet;
 837
 838   std::list<AttributePtr>::const_iterator aPIt =
 839     aPointAttributes.begin(), aPLast = aPointAttributes.end();
 840   for (; aPIt != aPLast; aPIt++)
 841     aPointAttributesSet.insert(*aPIt);
 842
 843   std::set<AttributePtr> aRefsAttributes = getFeatureResult(theFeature)->data()->refsToMe();
 844   std::set<AttributePtr> aFRefsList = theFeature->data()->refsToMe();
 845   aRefsAttributes.insert(aFRefsList.begin(), aFRefsList.end());
 846
 847   std::set<AttributePtr>::const_iterator aIt;
 848   for (aIt = aRefsAttributes.cbegin(); aIt != aRefsAttributes.cend(); ++aIt) {
 849     AttributePtr anAttr = (*aIt);
 850     FeaturePtr anAttrFeature = ModelAPI_Feature::feature(anAttr->owner());
 851     if (anAttrFeature.get() != this &&
 852         anAttr.get() && anAttr->attributeType() == ModelAPI_AttributeRefAttr::typeId()) {
 853       AttributeRefAttrPtr aRefAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(anAttr);
 854       if (!aRefAttr->isObject()) { /// find attributes referenced to feature point attributes
 855         AttributePtr anAttrInRef = aRefAttr->attr();
 856         if (anAttrInRef.get() &&
 857             aPointAttributesSet.find(anAttrInRef) != aPointAttributesSet.end()) {
 858           if (theRefs.find(anAttrInRef) != theRefs.end())
 859             theRefs[anAttrInRef].push_back(aRefAttr);
 860           else {
 861             std::list<AttributePtr> anAttrList;
 862             anAttrList.push_back(aRefAttr);
 863             theRefs[anAttrInRef] = anAttrList;
 864           }
 865         }
 866       }
 867       else { /// find attributes referenced to feature itself
 868         theRefsToFeature.push_back(anAttr);
 869       }
 870     }
 871   }
 872 }
 873
 874 void SketchPlugin_Split::updateCoincidenceConstraintsToFeature(
 875       const std::map<std::shared_ptr<ModelAPI_Feature>, IdToPointPair>& theCoincidenceToFeature,
 876       const std::set<std::shared_ptr<GeomDataAPI_Point2D> >& theFurtherCoincidences,
 877       const std::set<ResultPtr>& theFeatureResults,
 878       const FeaturePtr& theSplitFeature,
 879       std::set<FeaturePtr>& theFeaturesToDelete)
 880 {
 881   if (theCoincidenceToFeature.empty())
 882     return;
 883
 884   // we should build coincidence constraints to end of the split feature
 885   std::set<std::shared_ptr<GeomDataAPI_Point2D> > aNewCoincidencesToSplitFeature;
 886   AttributePoint2DPtr aStartPointAttr, anEndPointAttr;
 887   getFeaturePoints(theSplitFeature, aStartPointAttr, anEndPointAttr);
 888   if (theFurtherCoincidences.find(aStartPointAttr) == theFurtherCoincidences.end())
 889     aNewCoincidencesToSplitFeature.insert(aStartPointAttr);
 890   if (theFurtherCoincidences.find(anEndPointAttr) == theFurtherCoincidences.end())
 891     aNewCoincidencesToSplitFeature.insert(anEndPointAttr);
 892
 893   std::map<FeaturePtr, IdToPointPair>::const_iterator aCIt = theCoincidenceToFeature.begin(),
 894                                                             aCLast = theCoincidenceToFeature.end();
 895 #ifdef DEBUG_SPLIT
 896   std::cout << std::endl;
 897   std::cout << "Coincidences to feature(modified):"<< std::endl;
 898 #endif
 899   for (; aCIt != aCLast; aCIt++) {
 900     FeaturePtr aCoincFeature = aCIt->first;
 901     std::string anAttributeId = aCIt->second.first;
 902     std::string aSecondAttribute = anAttributeId == SketchPlugin_Constraint::ENTITY_A() ?
 903         SketchPlugin_Constraint::ENTITY_B() : SketchPlugin_Constraint::ENTITY_A();
 904
 905     AttributePoint2DPtr aCoincPoint = aCIt->second.second;
 906     std::set<AttributePoint2DPtr>::const_iterator aFCIt = theFurtherCoincidences.begin(),
 907                                                   aFCLast = theFurtherCoincidences.end();
 908     std::shared_ptr<GeomAPI_Pnt2d> aCoincPnt = aCoincPoint->pnt();
 909     AttributePoint2DPtr aFeaturePointAttribute;
 910     for (; aFCIt != aFCLast && !aFeaturePointAttribute.get(); aFCIt++) {
 911       AttributePoint2DPtr aFCAttribute = *aFCIt;
 912       if (aCoincPnt->isEqual(aFCAttribute->pnt()))
 913         aFeaturePointAttribute = aFCAttribute;
 914     }
 915     if (aFeaturePointAttribute.get()) {
 916       // create new constraint and remove the current
 917       aCoincFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
 918           SketchPlugin_ConstraintCoincidence::ID(),
 919           aFeaturePointAttribute, aCoincFeature->refattr(aSecondAttribute)->attr());
 920       theFeaturesToDelete.insert(aCIt->first);
 921       // create new coincidences to split feature points
 922       std::set<AttributePoint2DPtr>::const_iterator aSFIt = aNewCoincidencesToSplitFeature.begin(),
 923                                                     aSFLast = aNewCoincidencesToSplitFeature.end();
 924       for (; aSFIt != aSFLast; aSFIt++) {
 925         AttributePoint2DPtr aSFAttribute = *aSFIt;
 926         if (aCoincPnt->isEqual(aSFAttribute->pnt())) {
 927           SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
 928                            SketchPlugin_ConstraintCoincidence::ID(),
 929                            aSFAttribute, aCoincFeature->refattr(aSecondAttribute)->attr());
 930         }
 931       }
 932     }
 933     else {
 934       /// find feature by shape intersected the point
 935       ResultPtr aResultForCoincidence = *(theFeatureResults.begin());
 936
 937       if (theFeatureResults.size() > 1) { // try to find point on additional feature
 938         ResultPtr anAddtionalResult = *(theFeatureResults.begin()++);
 939         GeomShapePtr aShape = anAddtionalResult->shape();
 940
 941         std::shared_ptr<GeomAPI_Pnt2d> aPnt2d = aCoincPoint->pnt();
 942         std::shared_ptr<GeomAPI_Pnt> aPoint = sketch()->to3D(aPnt2d->x(), aPnt2d->y());
 943
 944         std::shared_ptr<GeomAPI_Pnt> aProjectedPoint;
 945         if (ModelGeomAlgo_Point2D::isPointOnEdge(aShape, aPoint, aProjectedPoint))
 946           aResultForCoincidence = anAddtionalResult;
 947       }
 948       aCoincFeature->refattr(anAttributeId)->setObject(aResultForCoincidence);
 949     }
 950 #ifdef DEBUG_SPLIT
 951   std::cout << " -" << getFeatureInfo(aCoincFeature) << std::endl;
 952 #endif
 953   }
 954 }
 955
 956 void SketchPlugin_Split::updateRefFeatureConstraints(
 957                                                   const ResultPtr& theFeatureBaseResult,
 958                                                   const std::list<AttributePtr>& theRefsToFeature)
 959 {
 960   std::list<AttributePtr>::const_iterator anIt = theRefsToFeature.begin(),
 961                                           aLast = theRefsToFeature.end();
 962   for (; anIt != aLast; anIt++) {
 963     AttributeRefAttrPtr aRefAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*anIt);
 964     if (aRefAttr.get())
 965       aRefAttr->setObject(theFeatureBaseResult);
 966   }
 967 }
 968
 969 void SketchPlugin_Split::updateRefAttConstraints(
 970                     const std::map<AttributePtr, std::list<AttributePtr> >& theBaseRefAttributes,
 971                     const std::set<std::pair<AttributePtr, AttributePtr> >& theModifiedAttributes)
 972 {
 973 #ifdef DEBUG_SPLIT
 974   std::cout << "SketchPlugin_Split::updateRefAttConstraints" << std::endl;
 975 #endif
 976
 977   std::set<std::pair<AttributePtr, AttributePtr> >::const_iterator
 978     anIt = theModifiedAttributes.begin(),  aLast = theModifiedAttributes.end();
 979   for (; anIt != aLast; anIt++) {
 980     AttributePtr anAttribute = anIt->first;
 981
 982     /// not found in references
 983     if (theBaseRefAttributes.find(anAttribute) == theBaseRefAttributes.end())
 984       continue;
 985     std::list<AttributePtr> aRefAttributes = theBaseRefAttributes.at(anAttribute);
 986     std::list<AttributePtr>::const_iterator aRefIt = aRefAttributes.begin(),
 987                                             aRLast = aRefAttributes.end();
 988
 989     AttributePtr aNewAttribute = anIt->second;
 990     for (; aRefIt != aRLast; aRefIt++) {
 991       AttributeRefAttrPtr aRefAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(*aRefIt);
 992       if (aRefAttr.get()) {
 993         aRefAttr->setAttr(aNewAttribute);
 994 #ifdef DEBUG_SPLIT
 995         FeaturePtr aFeature = ModelAPI_Feature::feature(aRefAttr->owner());
 996         std::cout << " -" << getFeatureInfo(aFeature) << std::endl;
 997 #endif
 998       }
 999     }
1000   }
1001 }
1002
1003 FeaturePtr SketchPlugin_Split::splitLine(FeaturePtr& theSplitFeature,
1004                                              FeaturePtr& theBaseFeatureModified,
1005                                              FeaturePtr& theAfterFeature,
1006                                              std::set<AttributePoint2DPtr>& thePoints,
1007                                              std::set<FeaturePtr>& theCreatedFeatures,
1008                  std::set<std::pair<AttributePtr, AttributePtr>>& theModifiedAttributes)
1009 {
1010   FeaturePtr anNewFeature;
1011
1012   std::set<FeaturePtr> aCreatedFeatures;
1013   FeaturePtr aConstraintFeature;
1014   theBaseFeatureModified = FeaturePtr(); // it will contain modified base feature
1015
1016   SketchPlugin_Sketch* aSketch = sketch();
1017   if (!aSketch)
1018     return anNewFeature;
1019
1020   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1021                                                            data()->attribute(SELECTED_OBJECT()));
1022   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1023   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
1024   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
1025   std::string aFeatureKind = aBaseFeature->getKind();
1026   if (aFeatureKind != SketchPlugin_Line::ID())
1027     return anNewFeature;
1028
1029   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
1030     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
1031   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
1032     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
1033   AttributePoint2DPtr aStartPointAttrOfBase, anEndPointAttrOfBase;
1034
1035   getFeaturePoints(aBaseFeature, aStartPointAttrOfBase, anEndPointAttrOfBase);
1036   if (!aStartPointAttrOfBase.get() && !anEndPointAttrOfBase.get()) {
1037     setError("Error: Feature has no start and end points.");
1038     return anNewFeature;
1039   }
1040
1041   arrangePointsOnLine(aStartPointAttrOfBase, anEndPointAttrOfBase,
1042                       aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1043
1044 #ifdef DEBUG_SPLIT
1045   std::cout << "Arranged points (to build split between 1st and 2nd points:" << std::endl;
1046   std::cout << "Start point: " <<
1047     ModelGeomAlgo_Point2D::getPointAttributeInfo(aStartPointAttrOfBase) << std::endl;
1048   std::cout << "1st point:   " <<
1049     ModelGeomAlgo_Point2D::getPointAttributeInfo(aFirstPointAttrOfSplit) << std::endl;
1050   std::cout << "2nd point:   " <<
1051     ModelGeomAlgo_Point2D::getPointAttributeInfo(aSecondPointAttrOfSplit) << std::endl;
1052   std::cout << "End point:   " <<
1053     ModelGeomAlgo_Point2D::getPointAttributeInfo(anEndPointAttrOfBase) << std::endl;
1054 #endif
1055
1056   /// create a split feature
1057   theSplitFeature =
1058     createLineFeature(aBaseFeature, aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1059   theCreatedFeatures.insert(theSplitFeature);
1060
1061   // before split feature
1062   if (aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1063     theModifiedAttributes.insert(std::make_pair(aStartPointAttrOfBase,
1064                                         theSplitFeature->attribute(SketchPlugin_Line::START_ID())));
1065   }
1066   else {
1067     theBaseFeatureModified = aBaseFeature; ///< use base feature to store all constraints here
1068     /// move end arc point to start of split
1069   }
1070
1071   // after split feature
1072   if (!aSecondPointAttrOfSplit->pnt()->isEqual(anEndPointAttrOfBase->pnt())) {
1073     FeaturePtr aFeature;
1074     if (!theBaseFeatureModified.get()) {
1075       aFeature = aBaseFeature; ///< use base feature to store all constraints here
1076       fillAttribute(aFeature->attribute(SketchPlugin_Line::START_ID()), aSecondPointAttrOfSplit);
1077       aFeature->execute(); // to update result
1078     }
1079     else {
1080       aFeature = createLineFeature(aBaseFeature, aSecondPointAttrOfSplit, anEndPointAttrOfBase);
1081       theCreatedFeatures.insert(aFeature);
1082       theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1083                                              aFeature->attribute(SketchPlugin_Line::END_ID())));
1084       anNewFeature = aFeature;
1085     }
1086     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1087                      SketchPlugin_ConstraintCoincidence::ID(),
1088                      theSplitFeature->attribute(SketchPlugin_Line::END_ID()),
1089                      aFeature->attribute(SketchPlugin_Line::START_ID()));
1090     theCreatedFeatures.insert(aConstraintFeature);
1091
1092     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1093                                 (aFeature->attribute(SketchPlugin_Line::START_ID())));
1094     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1095                                 (aFeature->attribute(SketchPlugin_Line::END_ID())));
1096
1097     if (!theBaseFeatureModified.get())
1098       theBaseFeatureModified = aFeature;
1099     else
1100       theAfterFeature = aFeature;
1101   }
1102   else {
1103     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1104                                   (theSplitFeature->attribute(SketchPlugin_Line::END_ID())));
1105     theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1106                                    theSplitFeature->attribute(SketchPlugin_Line::END_ID())));
1107   }
1108   // base split, that is defined before split feature should be changed at end
1109   // (after the after feature creation). Otherwise modified value will be used in after feature
1110   // before split feature
1111   if (!aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1112     /// move end arc point to start of split
1113     fillAttribute(theBaseFeatureModified->attribute(SketchPlugin_Line::END_ID()),
1114                                                     aFirstPointAttrOfSplit);
1115     theBaseFeatureModified->execute(); // to update result
1116     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1117                      SketchPlugin_ConstraintCoincidence::ID(),
1118                      theBaseFeatureModified->attribute(SketchPlugin_Line::END_ID()),
1119                      theSplitFeature->attribute(SketchPlugin_Line::START_ID()));
1120     theCreatedFeatures.insert(aConstraintFeature);
1121
1122     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1123                              (theBaseFeatureModified->attribute(SketchPlugin_Line::START_ID())));
1124     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1125                                (theBaseFeatureModified->attribute(SketchPlugin_Line::END_ID())));
1126   }
1127   else
1128     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1129                                        (theSplitFeature->attribute(SketchPlugin_Line::START_ID())));
1130
1131 #ifdef CREATE_CONSTRAINTS
1132   // additional constraints between split and base features
1133   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1134           SketchPlugin_ConstraintParallel::ID(),
1135           getFeatureResult(aBaseFeature),
1136           getFeatureResult(theSplitFeature));
1137   theCreatedFeatures.insert(aConstraintFeature);
1138   if (theAfterFeature.get()) {
1139     aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1140           SketchPlugin_ConstraintParallel::ID(),
1141           getFeatureResult(aBaseFeature),
1142           getFeatureResult(theAfterFeature));
1143     theCreatedFeatures.insert(aConstraintFeature);
1144   }
1145 #endif
1146   return anNewFeature;
1147 }
1148
1149 FeaturePtr SketchPlugin_Split::splitArc(FeaturePtr& theSplitFeature,
1150                                             FeaturePtr& theBaseFeatureModified,
1151                                             FeaturePtr& theAfterFeature,
1152                                             std::set<AttributePoint2DPtr>& thePoints,
1153                                             std::set<FeaturePtr>& theCreatedFeatures,
1154                  std::set<std::pair<AttributePtr, AttributePtr>>& theModifiedAttributes)
1155 {
1156   FeaturePtr anNewFeature;
1157
1158   std::set<FeaturePtr> aCreatedFeatures;
1159   FeaturePtr aConstraintFeature;
1160   theBaseFeatureModified = FeaturePtr(); // it will contain modified base feature
1161
1162   SketchPlugin_Sketch* aSketch = sketch();
1163   if (!aSketch)
1164     return anNewFeature;
1165
1166   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1167                                                            data()->attribute(SELECTED_OBJECT()));
1168   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1169   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
1170   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
1171   std::string aFeatureKind = aBaseFeature->getKind();
1172   if (aFeatureKind != SketchPlugin_Arc::ID())
1173     return anNewFeature;
1174
1175   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
1176     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
1177   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
1178     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
1179   AttributePoint2DPtr aStartPointAttrOfBase, anEndPointAttrOfBase;
1180   getFeaturePoints(aBaseFeature, aStartPointAttrOfBase, anEndPointAttrOfBase);
1181   if (!aStartPointAttrOfBase.get() && !anEndPointAttrOfBase.get()) {
1182     setError("Error: Feature has no start and end points.");
1183     return anNewFeature;
1184   }
1185
1186   arrangePointsOnArc(aBaseFeature, aStartPointAttrOfBase, anEndPointAttrOfBase,
1187                      aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1188 #ifdef DEBUG_SPLIT
1189   std::cout << "Arranged points (to build split between 1st and 2nd points:" << std::endl;
1190   std::cout << "Start point: " <<
1191     ModelGeomAlgo_Point2D::getPointAttributeInfo(aStartPointAttrOfBase) << std::endl;
1192   std::cout << "1st point:   " <<
1193     ModelGeomAlgo_Point2D::getPointAttributeInfo(aFirstPointAttrOfSplit) << std::endl;
1194   std::cout << "2nd point:   " <<
1195     ModelGeomAlgo_Point2D::getPointAttributeInfo(aSecondPointAttrOfSplit) << std::endl;
1196   std::cout << "End point:   " <<
1197     ModelGeomAlgo_Point2D::getPointAttributeInfo(anEndPointAttrOfBase) << std::endl;
1198 #endif
1199
1200   /// split feature
1201   theSplitFeature = createArcFeature(aBaseFeature, aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1202   theCreatedFeatures.insert(theSplitFeature);
1203
1204   // before split feature
1205   if (aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1206     theModifiedAttributes.insert(std::make_pair(aStartPointAttrOfBase,
1207                                   theSplitFeature->attribute(SketchPlugin_Arc::START_ID())));
1208   }
1209   else {
1210     theBaseFeatureModified = aBaseFeature; ///< use base feature to store all constraints here
1211     /// move end arc point to start of split
1212   }
1213
1214   // after split feature
1215   if (!aSecondPointAttrOfSplit->pnt()->isEqual(anEndPointAttrOfBase->pnt())) {
1216     FeaturePtr aFeature;
1217     if (!theBaseFeatureModified.get()) {
1218       aFeature = aBaseFeature; ///< use base feature to store all constraints here
1219       fillAttribute(aFeature->attribute(SketchPlugin_Arc::START_ID()), aSecondPointAttrOfSplit);
1220       aFeature->execute(); // to update result
1221     }
1222     else {
1223       aFeature = createArcFeature(aBaseFeature, aSecondPointAttrOfSplit, anEndPointAttrOfBase);
1224       theCreatedFeatures.insert(aFeature);
1225       theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1226                                                   aFeature->attribute(SketchPlugin_Arc::END_ID())));
1227       anNewFeature = aFeature;
1228     }
1229     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1230                      SketchPlugin_ConstraintCoincidence::ID(),
1231                      theSplitFeature->attribute(SketchPlugin_Arc::END_ID()),
1232                      aFeature->attribute(SketchPlugin_Arc::START_ID()));
1233     theCreatedFeatures.insert(aConstraintFeature);
1234
1235     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1236                                 (aFeature->attribute(SketchPlugin_Arc::START_ID())));
1237     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1238                                 (aFeature->attribute(SketchPlugin_Arc::END_ID())));
1239
1240     if (!theBaseFeatureModified.get())
1241       theBaseFeatureModified = aFeature;
1242     else
1243       theAfterFeature = aFeature;
1244   }
1245   else {
1246     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1247                                   (theSplitFeature->attribute(SketchPlugin_Arc::END_ID())));
1248     theModifiedAttributes.insert(std::make_pair(anEndPointAttrOfBase,
1249                                    theSplitFeature->attribute(SketchPlugin_Arc::END_ID())));
1250   }
1251   // base split, that is defined before split feature should be changed at end
1252   // (after the after feature creation). Otherwise modified value will be used in after feature
1253   // before split feature
1254   if (!aStartPointAttrOfBase->pnt()->isEqual(aFirstPointAttrOfSplit->pnt())) {
1255     /// move end arc point to start of split
1256     fillAttribute(theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID()),
1257                                                     aFirstPointAttrOfSplit);
1258     theBaseFeatureModified->execute(); // to update result
1259     aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1260                      SketchPlugin_ConstraintCoincidence::ID(),
1261                      theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID()),
1262                      theSplitFeature->attribute(SketchPlugin_Arc::START_ID()));
1263     theCreatedFeatures.insert(aConstraintFeature);
1264
1265     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1266                              (theBaseFeatureModified->attribute(SketchPlugin_Arc::START_ID())));
1267     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1268                                (theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID())));
1269   }
1270   else
1271     thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1272                                        (theSplitFeature->attribute(SketchPlugin_Arc::START_ID())));
1273
1274   // additional constraints between split and base features
1275 #ifdef CREATE_CONSTRAINTS
1276   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1277           SketchPlugin_ConstraintEqual::ID(),
1278           getFeatureResult(aBaseFeature),
1279           getFeatureResult(theSplitFeature));
1280   theCreatedFeatures.insert(aConstraintFeature);
1281   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1282           SketchPlugin_ConstraintTangent::ID(),
1283           getFeatureResult(theSplitFeature),
1284           getFeatureResult(aBaseFeature));
1285   theCreatedFeatures.insert(aConstraintFeature);
1286   if (theAfterFeature.get()) {
1287     aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1288                 SketchPlugin_ConstraintEqual::ID(),
1289                 getFeatureResult(aBaseFeature),
1290                 getFeatureResult(theAfterFeature));
1291     theCreatedFeatures.insert(aConstraintFeature);
1292     aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1293                 SketchPlugin_ConstraintTangent::ID(),
1294                 getFeatureResult(theSplitFeature),
1295                 getFeatureResult(theAfterFeature));
1296     theCreatedFeatures.insert(aConstraintFeature);
1297   }
1298 #endif
1299   return anNewFeature;
1300 }
1301
1302 FeaturePtr SketchPlugin_Split::splitCircle(FeaturePtr& theSplitFeature,
1303                                                FeaturePtr& theBaseFeatureModified,
1304                                                FeaturePtr& theAfterFeature,
1305                                                std::set<AttributePoint2DPtr>& thePoints,
1306                                                std::set<FeaturePtr>& theCreatedFeatures,
1307                  std::set<std::pair<AttributePtr, AttributePtr>>& theModifiedAttributes)
1308 {
1309   FeaturePtr anNewFeature;
1310
1311   std::set<FeaturePtr> aCreatedFeatures;
1312   FeaturePtr aConstraintFeature;
1313   theBaseFeatureModified = FeaturePtr(); // it will contain modified base feature
1314
1315   SketchPlugin_Sketch* aSketch = sketch();
1316   if (!aSketch)
1317     return anNewFeature;
1318
1319   AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1320                                                            data()->attribute(SELECTED_OBJECT()));
1321   //AttributeReferencePtr aBaseObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1322   //                                         data()->attribute(SketchPlugin_Constraint::VALUE()));
1323   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObjectAttr->value());
1324   std::string aFeatureKind = aBaseFeature->getKind();
1325   if (aFeatureKind != SketchPlugin_Circle::ID())
1326     return anNewFeature;
1327
1328   AttributePoint2DPtr aFirstPointAttrOfSplit = getPointAttribute(true);
1329     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_A()));
1330   AttributePoint2DPtr aSecondPointAttrOfSplit = getPointAttribute(false);
1331     //getPointOfRefAttr(data()->attribute(SketchPlugin_Constraint::ENTITY_B()));
1332
1333   /// split feature
1334   theSplitFeature =
1335     createArcFeature(aBaseFeature, aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1336   bool aSplitReversed = std::dynamic_pointer_cast<SketchPlugin_Arc>(theSplitFeature)->isReversed();
1337   theCreatedFeatures.insert(theSplitFeature);
1338
1339   /// base feature is a left part of the circle
1340   theBaseFeatureModified = createArcFeature(aBaseFeature,
1341     aFirstPointAttrOfSplit, aSecondPointAttrOfSplit);
1342   anNewFeature = theBaseFeatureModified;
1343   std::dynamic_pointer_cast<SketchPlugin_Arc>(
1344     theBaseFeatureModified)->setReversed(!aSplitReversed);
1345   theBaseFeatureModified->execute();
1346
1347   theModifiedAttributes.insert(
1348     std::make_pair(aBaseFeature->attribute(SketchPlugin_Circle::CENTER_ID()),
1349                   theBaseFeatureModified->attribute(SketchPlugin_Arc::CENTER_ID())));
1350
1351   theCreatedFeatures.insert(theBaseFeatureModified);
1352
1353   thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1354                              (theBaseFeatureModified->attribute(SketchPlugin_Arc::START_ID())));
1355   thePoints.insert(std::dynamic_pointer_cast<GeomDataAPI_Point2D>
1356                              (theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID())));
1357
1358   // additional constraints between split and base features
1359   aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1360                      SketchPlugin_ConstraintCoincidence::ID(),
1361                      theBaseFeatureModified->attribute(SketchPlugin_Arc::END_ID()),
1362                      theSplitFeature->attribute(SketchPlugin_Arc::END_ID()));
1363   theCreatedFeatures.insert(aConstraintFeature);
1364   aConstraintFeature = SketchPlugin_Tools::createConstraintAttrAttr(sketch(),
1365                      SketchPlugin_ConstraintCoincidence::ID(),
1366                      theBaseFeatureModified->attribute(SketchPlugin_Arc::START_ID()),
1367                      theSplitFeature->attribute(SketchPlugin_Arc::START_ID()));
1368   theCreatedFeatures.insert(aConstraintFeature);
1369
1370 #ifdef CREATE_CONSTRAINTS
1371   aConstraintFeature = SketchPlugin_Tools::createConstraint(sketch(),
1372                      SketchPlugin_ConstraintTangent::ID(),
1373                      getFeatureResult(theSplitFeature),
1374                      getFeatureResult(theBaseFeatureModified));
1375   theCreatedFeatures.insert(aConstraintFeature);
1376 #endif
1377   return anNewFeature;
1378 }
1379
1380 void SketchPlugin_Split::arrangePointsOnLine(
1381     const AttributePoint2DPtr& theStartPointAttr,
1382     const AttributePoint2DPtr& theEndPointAttr,
1383     AttributePoint2DPtr& theFirstPointAttr,
1384     AttributePoint2DPtr& theLastPointAttr) const
1385 {
1386   // if first point is closer to last point, swap first and last values
1387   if (theStartPointAttr->pnt()->distance(theFirstPointAttr->pnt()) >
1388       theStartPointAttr->pnt()->distance(theLastPointAttr->pnt())) {
1389     AttributePoint2DPtr aTmpPoint = theFirstPointAttr;
1390     theFirstPointAttr = theLastPointAttr;
1391     theLastPointAttr = aTmpPoint;
1392   }
1393 }
1394
1395 void SketchPlugin_Split::arrangePointsOnArc(
1396     const FeaturePtr& theArc,
1397     const std::shared_ptr<GeomDataAPI_Point2D>& theStartPointAttr,
1398     const std::shared_ptr<GeomDataAPI_Point2D>& theEndPointAttr,
1399     std::shared_ptr<GeomDataAPI_Point2D>& theFirstPointAttr,
1400     std::shared_ptr<GeomDataAPI_Point2D>& theSecondPointAttr) const
1401 {
1402   static const double anAngleTol = 1.e-12;
1403
1404   std::shared_ptr<GeomAPI_Pnt2d> aCenter = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
1405       theArc->attribute(SketchPlugin_Arc::CENTER_ID()))->pnt();
1406   bool isReversed = theArc->boolean(SketchPlugin_Arc::REVERSED_ID())->value();
1407
1408   // collect directions to each point
1409   std::shared_ptr<GeomAPI_Dir2d> aStartDir(
1410       new GeomAPI_Dir2d(theStartPointAttr->pnt()->xy()->decreased(aCenter->xy())));
1411   std::shared_ptr<GeomAPI_Dir2d> aFirstPtDir(
1412       new GeomAPI_Dir2d(theFirstPointAttr->pnt()->xy()->decreased(aCenter->xy())));
1413   std::shared_ptr<GeomAPI_Dir2d> aSecondPtDir(
1414       new GeomAPI_Dir2d(theSecondPointAttr->pnt()->xy()->decreased(aCenter->xy())));
1415
1416   // sort points by their angular values
1417   double aFirstPtAngle = aStartDir->angle(aFirstPtDir);
1418   double aSecondPtAngle = aStartDir->angle(aSecondPtDir);
1419   double aPeriod = isReversed ? -2.0 * PI : 2.0 * PI;
1420   if (fabs(aFirstPtAngle) > anAngleTol && isReversed == (aFirstPtAngle > 0.))
1421     aFirstPtAngle += aPeriod;
1422   if (fabs(aSecondPtAngle) > anAngleTol && isReversed == (aSecondPtAngle > 0.))
1423     aSecondPtAngle += aPeriod;
1424
1425   if (fabs(aFirstPtAngle) > fabs(aSecondPtAngle)) {
1426     AttributePoint2DPtr aTmpPoint = theFirstPointAttr;
1427     theFirstPointAttr = theSecondPointAttr;
1428     theSecondPointAttr = aTmpPoint;
1429   }
1430 }
1431
1432 void SketchPlugin_Split::fillAttribute(const AttributePtr& theModifiedAttribute,
1433                                                  const AttributePtr& theSourceAttribute)
1434 {
1435   std::string anAttributeType = theModifiedAttribute->attributeType();
1436   if (anAttributeType == GeomDataAPI_Point2D::typeId()) {
1437     AttributePoint2DPtr aModifiedAttribute = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
1438                                               theModifiedAttribute);
1439     AttributePoint2DPtr aSourceAttribute = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
1440                                               theSourceAttribute);
1441
1442     if (aModifiedAttribute.get() && aSourceAttribute.get())
1443       aModifiedAttribute->setValue(aSourceAttribute->pnt());
1444   }
1445   else if (anAttributeType == ModelAPI_AttributeBoolean::typeId()) {
1446     AttributeBooleanPtr aModifiedAttribute = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
1447                                               theModifiedAttribute);
1448     AttributeBooleanPtr aSourceAttribute = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
1449                                               theSourceAttribute);
1450
1451     if (aModifiedAttribute.get() && aSourceAttribute.get())
1452       aModifiedAttribute->setValue(aSourceAttribute->value());
1453   }
1454 }
1455
1456 FeaturePtr SketchPlugin_Split::createLineFeature(const FeaturePtr& theBaseFeature,
1457                                                            const AttributePtr& theFirstPointAttr,
1458                                                            const AttributePtr& theSecondPointAttr)
1459 {
1460   FeaturePtr aFeature;
1461   SketchPlugin_Sketch* aSketch = sketch();
1462   if (!aSketch || !theBaseFeature.get())
1463     return aFeature;
1464
1465   aFeature = aSketch->addFeature(SketchPlugin_Line::ID());
1466
1467   fillAttribute(aFeature->attribute(SketchPlugin_Line::START_ID()), theFirstPointAttr);
1468   fillAttribute(aFeature->attribute(SketchPlugin_Line::END_ID()), theSecondPointAttr);
1469
1470   fillAttribute(aFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()),
1471                 theBaseFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()));
1472
1473   aFeature->execute(); // to obtain result
1474
1475   return aFeature;
1476 }
1477
1478 FeaturePtr SketchPlugin_Split::createArcFeature(const FeaturePtr& theBaseFeature,
1479                                                           const AttributePtr& theFirstPointAttr,
1480                                                           const AttributePtr& theSecondPointAttr)
1481 {
1482   FeaturePtr aFeature;
1483   SketchPlugin_Sketch* aSketch = sketch();
1484   if (!aSketch || !theBaseFeature.get())
1485     return aFeature;
1486
1487   std::string aCenterAttributeId;
1488   if (theBaseFeature->getKind() == SketchPlugin_Arc::ID())
1489     aCenterAttributeId = SketchPlugin_Arc::CENTER_ID();
1490   else if (theBaseFeature->getKind() == SketchPlugin_Circle::ID())
1491     aCenterAttributeId = SketchPlugin_Circle::CENTER_ID();
1492
1493   if (aCenterAttributeId.empty())
1494     return aFeature;
1495
1496   aFeature = aSketch->addFeature(SketchPlugin_Arc::ID());
1497   // update fillet arc: make the arc correct for sure, so, it is not needed to process
1498   // the "attribute updated"
1499   // by arc; moreover, it may cause cyclicity in hte mechanism of updater
1500   bool aWasBlocked = aFeature->data()->blockSendAttributeUpdated(true);
1501
1502   fillAttribute(aFeature->attribute(SketchPlugin_Arc::CENTER_ID()),
1503                 theBaseFeature->attribute(aCenterAttributeId));
1504   fillAttribute(aFeature->attribute(SketchPlugin_Arc::START_ID()), theFirstPointAttr);
1505   fillAttribute(aFeature->attribute(SketchPlugin_Arc::END_ID()), theSecondPointAttr);
1506
1507   fillAttribute(aFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()),
1508                 theBaseFeature->attribute(SketchPlugin_SketchEntity::AUXILIARY_ID()));
1509
1510   /// fill referersed state of created arc as it is on the base arc
1511   if (theBaseFeature->getKind() == SketchPlugin_Arc::ID()) {
1512     bool aReversed = theBaseFeature->boolean(SketchPlugin_Arc::REVERSED_ID())->value();
1513     aFeature->boolean(SketchPlugin_Arc::REVERSED_ID())->setValue(aReversed);
1514   }
1515   aFeature->data()->blockSendAttributeUpdated(aWasBlocked);
1516   aFeature->execute(); // to obtain result
1517
1518   return aFeature;
1519 }
1520
1521 void SketchPlugin_Split::updateFeaturesAfterSplit(
1522                                                    const std::set<FeaturePtr>& theFeaturesToUpdate)
1523 {
1524   std::set<FeaturePtr>::const_iterator anIt = theFeaturesToUpdate.begin(),
1525                                        aLast = theFeaturesToUpdate.end();
1526   for (; anIt != aLast; anIt++) {
1527     FeaturePtr aRefFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(*anIt);
1528     std::string aRefFeatureKind = aRefFeature->getKind();
1529     if (aRefFeatureKind == SketchPlugin_ConstraintLength::ID()) {
1530       std::shared_ptr<SketchPlugin_ConstraintLength> aLenghtFeature =
1531                               std::dynamic_pointer_cast<SketchPlugin_ConstraintLength>(*anIt);
1532       if (aLenghtFeature.get()) {
1533         std::shared_ptr<ModelAPI_AttributeDouble> aValueAttr = std::dynamic_pointer_cast<
1534             ModelAPI_AttributeDouble>(aLenghtFeature->attribute(SketchPlugin_Constraint::VALUE()));
1535         double aValue;
1536         if (aLenghtFeature->computeLenghtValue(aValue) && aValueAttr.get())
1537           aValueAttr->setValue(aValue);
1538       }
1539     }
1540   }
1541 }
1542
1543 std::shared_ptr<ModelAPI_Result> SketchPlugin_Split::getFeatureResult(
1544                                     const std::shared_ptr<ModelAPI_Feature>& theFeature)
1545 {
1546   std::shared_ptr<ModelAPI_Result> aResult;
1547
1548   std::string aFeatureKind = theFeature->getKind();
1549   if (aFeatureKind == SketchPlugin_Line::ID())
1550     aResult = theFeature->firstResult();
1551   else if (aFeatureKind == SketchPlugin_Arc::ID())
1552     aResult = theFeature->lastResult();
1553   else if (aFeatureKind == SketchPlugin_Circle::ID())
1554     aResult = theFeature->lastResult();
1555
1556   return aResult;
1557 }
1558
1559 std::set<std::shared_ptr<ModelAPI_Attribute> > SketchPlugin_Split::getEdgeAttributes(
1560                                            const std::shared_ptr<ModelAPI_Feature>& theFeature)
1561 {
1562   std::set<std::shared_ptr<ModelAPI_Attribute> > anAttributes;
1563
1564   std::string aFeatureKind = theFeature->getKind();
1565   if (aFeatureKind == SketchPlugin_Line::ID()) {
1566     anAttributes.insert(theFeature->attribute(SketchPlugin_Line::START_ID()));
1567     anAttributes.insert(theFeature->attribute(SketchPlugin_Line::END_ID()));
1568   }
1569   else if (aFeatureKind == SketchPlugin_Arc::ID()) {
1570     anAttributes.insert(theFeature->attribute(SketchPlugin_Arc::START_ID()));
1571     anAttributes.insert(theFeature->attribute(SketchPlugin_Arc::END_ID()));
1572   }
1573   else if (aFeatureKind == SketchPlugin_Circle::ID()) {
1574   }
1575
1576   return anAttributes;
1577 }
1578
1579 std::shared_ptr<GeomDataAPI_Point2D> SketchPlugin_Split::getPointAttribute
1580                                                               (const bool isFirstAttribute)
1581 {
1582   std::shared_ptr<GeomDataAPI_Point2D> anAttribute;
1583
1584   GeomShapePtr aSelectedShape = getSubShape(SELECTED_OBJECT(), SELECTED_POINT());
1585   if (!aSelectedShape.get())
1586     return anAttribute;
1587
1588   if (aSelectedShape->shapeType() != GeomAPI_Shape::EDGE)
1589     return anAttribute;
1590
1591   AttributeReferencePtr anObjectAttr = std::dynamic_pointer_cast<ModelAPI_AttributeReference>(
1592                                                        data()->attribute(SELECTED_OBJECT()));
1593   ObjectPtr aBaseObject = anObjectAttr->value();
1594   if (!aBaseObject.get())
1595     return anAttribute;
1596
1597   std::shared_ptr<GeomAPI_Edge> anEdge(new GeomAPI_Edge(aSelectedShape));
1598
1599   std::shared_ptr<GeomAPI_Pnt> aFirstPnt = anEdge->firstPoint();
1600   std::shared_ptr<GeomAPI_Pnt> aLastPnt = anEdge->lastPoint();
1601
1602   std::shared_ptr<GeomDataAPI_Point2D> aFirstPointAttr, aLastPointAttr;
1603   /// find the points in feature attributes
1604   FeaturePtr aBaseFeature = ModelAPI_Feature::feature(aBaseObject);
1605   std::list<AttributePtr> a2DPointAttributes = aBaseFeature->data()->attributes(
1606                                                     GeomDataAPI_Point2D::typeId());
1607   std::list<AttributePtr>::const_iterator anIt = a2DPointAttributes.begin(),
1608                                           aLast = a2DPointAttributes.end();
1609   for (; anIt != aLast; anIt++) {
1610     std::shared_ptr<GeomDataAPI_Point2D> anAttributePoint =
1611                                   std::dynamic_pointer_cast<GeomDataAPI_Point2D>(*anIt);
1612     std::shared_ptr<GeomAPI_Pnt2d> aPoint2D = anAttributePoint->pnt();
1613     std::shared_ptr<GeomAPI_Pnt> aPoint3D = sketch()->to3D(aPoint2D->x(), aPoint2D->y());
1614     if (aFirstPnt->isEqual(aPoint3D))
1615       aFirstPointAttr = anAttributePoint;
1616     else if (aLastPnt->isEqual(aPoint3D))
1617       aLastPointAttr = anAttributePoint;
1618   }
1619
1620   /// find the points in coincident features
1621   PntToAttributesMap aRefAttributes = myCashedReferences[aBaseObject];
1622   PntToAttributesMap::const_iterator
1623     aRIt = aRefAttributes.begin(), aRLast = aRefAttributes.end();
1624   for (; aRIt != aRLast; aRIt++) {
1625     std::shared_ptr<GeomDataAPI_Point2D> anAttribute = aRIt->first;
1626     std::shared_ptr<GeomAPI_Pnt> aPoint = aRIt->second;
1627     if (!aFirstPointAttr.get() && aFirstPnt->isEqual(aPoint))
1628       aFirstPointAttr = anAttribute;
1629     if (!aLastPointAttr.get() && aLastPnt->isEqual(aPoint))
1630       aLastPointAttr = anAttribute;
1631     if (aFirstPointAttr.get() && aLastPointAttr.get())
1632       break;
1633   }
1634   if (!aFirstPointAttr.get() || !aLastPointAttr)
1635     return anAttribute;
1636
1637   return isFirstAttribute ? aFirstPointAttr : aLastPointAttr;
1638 }
1639
1640 #ifdef _DEBUG
1641 std::string SketchPlugin_Split::getFeatureInfo(const std::shared_ptr<ModelAPI_Feature>& theFeature,
1642                                                const bool isUseAttributesInfo)
1643 {
1644   std::string anInfo;
1645   if (!theFeature.get()) {
1646     return "none";
1647   }
1648
1649   if (theFeature->data()->isValid())
1650     anInfo.append(theFeature->data()->name().c_str());
1651
1652   if (isUseAttributesInfo) {
1653     std::string aPointsInfo = ModelGeomAlgo_Point2D::getPontAttributesInfo(theFeature,
1654                                                              getEdgeAttributes(theFeature));
1655     /// processing of feature with point 2d attributes, like line, arc, circle
1656     if (!aPointsInfo.empty()) {
1657       anInfo += ": ";
1658       anInfo += "\n";
1659       anInfo += aPointsInfo;
1660     }
1661     else { /// process constraint coincidence, find points in ref attr attributes
1662       std::list<AttributePtr> anAttrs = theFeature->data()->attributes(
1663                                                        ModelAPI_AttributeRefAttr::typeId());
1664       std::list<AttributePtr>::const_iterator anIt = anAttrs.begin(), aLast = anAttrs.end();
1665       std::string anAttributesInfo;
1666       for(; anIt != aLast; anIt++) {
1667         if (!anAttributesInfo.empty()) {
1668           anAttributesInfo.append(", ");
1669           anAttributesInfo += "\n";
1670         }
1671         AttributePtr anAttr = *anIt;
1672         std::string aValue = "not defined";
1673         std::string aType = anAttr->attributeType();
1674         if (aType == ModelAPI_AttributeRefAttr::typeId()) {
1675           std::shared_ptr<ModelAPI_AttributeRefAttr> aRefAttr =
1676                              std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(anAttr);
1677           if (aRefAttr.get()) {
1678             if (aRefAttr->isObject()) {
1679               FeaturePtr aFeature = ModelAPI_Feature::feature(aRefAttr->object());
1680               aValue = "<object:>" + getFeatureInfo(aFeature, false);
1681             }
1682             else {
1683               AttributePtr anAttribute = aRefAttr->attr();
1684               if (anAttribute.get()) {
1685                 FeaturePtr aFeature = ModelAPI_Feature::feature(anAttribute->owner());
1686                 aValue = "<attr:>" + ModelGeomAlgo_Point2D::getPointAttributeInfo(anAttribute) +
1687                          " [" + getFeatureInfo(aFeature, false) + "]";
1688               }
1689             }
1690           }
1691         }
1692         anAttributesInfo.append("    " + anAttr->id() + ": " + aValue);
1693       }
1694       if (!anAttributesInfo.empty())
1695         anInfo = anInfo + "\n" + anAttributesInfo;
1696     }
1697   }
1698   return anInfo;
1699 }
1700 #endif