src/SketchPlugin/SketchPlugin_Arc.cpp

   1 // Copyright (C) 2014-20xx CEA/DEN, EDF R&D -->
   2
   3 // File:        SketchPlugin_Arc.cpp
   4 // Created:     26 Apr 2014
   5 // Author:      Artem ZHIDKOV
   6
   7 #include "SketchPlugin_Arc.h"
   8 #include "SketchPlugin_Sketch.h"
   9 #include <SketchPlugin_ConstraintCoincidence.h>
  10 #include <SketchPlugin_ConstraintTangent.h>
  11
  12 #include <Events_Loop.h>
  13 #include <ModelAPI_Data.h>
  14 #include <ModelAPI_ResultConstruction.h>
  15 #include <ModelAPI_AttributeDouble.h>
  16 #include <ModelAPI_AttributeRefAttr.h>
  17 #include <ModelAPI_AttributeSelection.h>
  18 #include <ModelAPI_AttributeString.h>
  19 #include <ModelAPI_Events.h>
  20 #include <ModelAPI_Validator.h>
  21 #include <ModelAPI_Session.h>
  22 #include <ModelAPI_Tools.h>
  23
  24 #include <GeomAPI_Ax2.h>
  25 #include <GeomAPI_Circ2d.h>
  26 #include <GeomAPI_Circ.h>
  27 #include <GeomAPI_Dir2d.h>
  28 #include <GeomAPI_Dir.h>
  29 #include <GeomAPI_Lin2d.h>
  30 #include <GeomAPI_Lin.h>
  31 #include <GeomAPI_Pnt2d.h>
  32 #include <GeomAPI_XY.h>
  33 #include <GeomDataAPI_Point2D.h>
  34 #include <GeomDataAPI_Dir.h>
  35 #include <GeomAlgoAPI_PointBuilder.h>
  36 #include <GeomAlgoAPI_EdgeBuilder.h>
  37 #include <GeomAlgoAPI_CompoundBuilder.h>
  38 // for sqrt on Linux
  39 #include <math.h>
  40
  41 const double tolerance = 1e-7;
  42 const double paramTolerance = 1.e-4;
  43 const double PI =3.141592653589793238463;
  44
  45 namespace {
  46   static const std::string& ARC_TYPE_CENTER_START_END()
  47   {
  48     static const std::string TYPE("CenterStartEnd");
  49     return TYPE;
  50   }
  51   static const std::string& ARC_TYPE_THREE_POINTS()
  52   {
  53     static const std::string TYPE("ThreePoints");
  54     return TYPE;
  55   }
  56
  57   static const std::string& PASSED_POINT_ID()
  58   {
  59     static const std::string PASSED_PNT("ArcPassedPoint");
  60     return PASSED_PNT;
  61   }
  62   static const std::string& RADIUS_ID()
  63   {
  64     static const std::string RADIUS("ArcRadius");
  65     return RADIUS;
  66   }
  67   static const std::string& ANGLE_ID()
  68   {
  69     static const std::string ANGLE("ArcAngle");
  70     return ANGLE;
  71   }
  72
  73   static const std::string& POINT_ID(int theIndex)
  74   {
  75     switch (theIndex) {
  76     case 1: return SketchPlugin_Arc::START_ID();
  77     case 2: return SketchPlugin_Arc::END_ID();
  78     case 3: return PASSED_POINT_ID();
  79     }
  80
  81     static const std::string DUMMY;
  82     return DUMMY;
  83   }
  84 }
  85
  86
  87
  88 SketchPlugin_Arc::SketchPlugin_Arc()
  89     : SketchPlugin_SketchEntity()
  90 {
  91   myStartUpdate = false;
  92   myEndUpdate = false;
  93   // default values
  94   myXEndBefore = 0;
  95   myYEndBefore = 0;
  96
  97   myParamBefore = 0;
  98 }
  99
 100 void SketchPlugin_Arc::initDerivedClassAttributes()
 101 {
 102   data()->addAttribute(CENTER_ID(), GeomDataAPI_Point2D::typeId());
 103   data()->addAttribute(START_ID(), GeomDataAPI_Point2D::typeId());
 104   std::shared_ptr<GeomDataAPI_Point2D> anEndAttr = std::dynamic_pointer_cast<
 105     GeomDataAPI_Point2D>(data()->addAttribute(END_ID(), GeomDataAPI_Point2D::typeId()));
 106   data()->addAttribute(EXTERNAL_ID(), ModelAPI_AttributeSelection::typeId());
 107   ModelAPI_Session::get()->validators()->registerNotObligatory(getKind(), EXTERNAL_ID());
 108
 109   AttributeBooleanPtr isInversed = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
 110     data()->addAttribute(INVERSED_ID(), ModelAPI_AttributeBoolean::typeId()));
 111
 112   // get the initial values
 113   if (anEndAttr->isInitialized()) {
 114     myXEndBefore = anEndAttr->x();
 115     myYEndBefore = anEndAttr->y();
 116   }
 117
 118   AttributeStringPtr anArcType = std::dynamic_pointer_cast<ModelAPI_AttributeString>(
 119     data()->addAttribute(ARC_TYPE(), ModelAPI_AttributeString::typeId()));
 120
 121   data()->addAttribute(PASSED_POINT_ID(), GeomDataAPI_Point2D::typeId());
 122   data()->addAttribute(TANGENT_POINT_ID(), ModelAPI_AttributeRefAttr::typeId());
 123   data()->addAttribute(RADIUS_ID(), ModelAPI_AttributeDouble::typeId());
 124   data()->addAttribute(ANGLE_ID(), ModelAPI_AttributeDouble::typeId());
 125
 126   // set after all to avoid in attributeChanged reference to not existing attributes
 127   if (!isInversed->isInitialized())
 128     isInversed->setValue(false);
 129   anArcType->setValue(ARC_TYPE_CENTER_START_END());
 130 }
 131
 132 void SketchPlugin_Arc::execute()
 133 {
 134   SketchPlugin_Sketch* aSketch = sketch();
 135   // result for the arc is set only when all obligatory attributes are initialized,
 136   // otherwise AIS object is used to visualize the arc's preview
 137   if (aSketch && isFeatureValid()) {
 138     bool hasResult = lastResult().get() != NULL;
 139
 140     // compute a circle point in 3D view
 141     std::shared_ptr<GeomDataAPI_Point2D> aCenterAttr = std::dynamic_pointer_cast<
 142         GeomDataAPI_Point2D>(data()->attribute(CENTER_ID()));
 143     // compute the arc start point
 144     std::shared_ptr<GeomDataAPI_Point2D> aStartAttr = std::dynamic_pointer_cast<
 145         GeomDataAPI_Point2D>(data()->attribute(START_ID()));
 146
 147     std::shared_ptr<GeomAPI_Pnt> aCenter(aSketch->to3D(aCenterAttr->x(), aCenterAttr->y()));
 148     // make a visible point
 149     std::shared_ptr<GeomAPI_Shape> aCenterPointShape = GeomAlgoAPI_PointBuilder::point(aCenter);
 150     std::shared_ptr<ModelAPI_ResultConstruction> aConstr1 = document()->createConstruction(
 151         data(), 0);
 152     aConstr1->setShape(aCenterPointShape);
 153     aConstr1->setIsInHistory(false);
 154     setResult(aConstr1, 0);
 155
 156     // make a visible circle
 157     std::shared_ptr<GeomDataAPI_Dir> aNDir = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 158         aSketch->data()->attribute(SketchPlugin_Sketch::NORM_ID()));
 159     std::shared_ptr<GeomAPI_Dir> aNormal = aNDir->dir();
 160     std::shared_ptr<GeomAPI_Pnt> aStartPoint(aSketch->to3D(aStartAttr->x(), aStartAttr->y()));
 161
 162     // compute and change the arc end point
 163     std::shared_ptr<GeomDataAPI_Point2D> anEndAttr = std::dynamic_pointer_cast<
 164         GeomDataAPI_Point2D>(data()->attribute(END_ID()));
 165     std::shared_ptr<GeomAPI_Circ2d> aCircleForArc(
 166         new GeomAPI_Circ2d(aCenterAttr->pnt(), aStartAttr->pnt()));
 167     std::shared_ptr<GeomAPI_Pnt2d> aProjection = aCircleForArc->project(anEndAttr->pnt());
 168     if (aProjection && anEndAttr->pnt()->distance(aProjection) > tolerance)
 169       anEndAttr->setValue(aProjection);
 170     std::shared_ptr<GeomAPI_Pnt> aEndPoint(aSketch->to3D(anEndAttr->x(), anEndAttr->y()));
 171     AttributeBooleanPtr isInversed =
 172         std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(attribute(INVERSED_ID()));
 173
 174     std::shared_ptr<GeomAPI_Dir> anXDir(new GeomAPI_Dir(aStartPoint->xyz()->decreased(aCenter->xyz())));
 175     std::shared_ptr<GeomAPI_Ax2> anAx2(new GeomAPI_Ax2(aCenter, aNormal, anXDir));
 176     std::shared_ptr<GeomAPI_Circ> aCirc(new GeomAPI_Circ(anAx2, aCenter->distance(aStartPoint)));
 177     double aParameterNew = 0.0;
 178     if(aCirc->parameter(aEndPoint, paramTolerance, aParameterNew)) {
 179       if(0 <= myParamBefore && myParamBefore <= PI / 2.0
 180         && PI * 1.5 <= aParameterNew && aParameterNew <= PI * 2.0) {
 181           isInversed->setValue(true);
 182       } else if(PI * 1.5 <= myParamBefore && myParamBefore <= PI * 2.0
 183         && 0 <= aParameterNew && aParameterNew <= PI / 2.0) {
 184           isInversed->setValue(false);
 185       }
 186     }
 187     myParamBefore = aParameterNew;
 188
 189     std::shared_ptr<GeomAPI_Shape> aCircleShape;
 190     if(!isInversed->value()) {
 191       aCircleShape = GeomAlgoAPI_EdgeBuilder::lineCircleArc(aCenter, aStartPoint, aEndPoint, aNormal);
 192     } else {
 193       aCircleShape = GeomAlgoAPI_EdgeBuilder::lineCircleArc(aCenter, aEndPoint, aStartPoint, aNormal);
 194     }
 195
 196     if (aCircleShape) {
 197       std::shared_ptr<ModelAPI_ResultConstruction> aConstr2 = document()->createConstruction(
 198           data(), 1);
 199       aConstr2->setShape(aCircleShape);
 200       aConstr2->setIsInHistory(false);
 201       setResult(aConstr2, 1);
 202     }
 203   }
 204 }
 205
 206 AISObjectPtr SketchPlugin_Arc::getAISObject(AISObjectPtr thePrevious)
 207 {
 208   SketchPlugin_Sketch* aSketch = sketch();
 209   if (aSketch) {
 210     // if the feature is valid, the execute() method should be performed, AIS object is empty
 211     if (!isFeatureValid()) {
 212       // compute a circle point in 3D view
 213       std::shared_ptr<GeomDataAPI_Point2D> aCenterAttr = std::dynamic_pointer_cast<
 214           GeomDataAPI_Point2D>(data()->attribute(CENTER_ID()));
 215
 216       std::list<std::shared_ptr<GeomAPI_Shape> > aShapes;
 217       if (aCenterAttr->isInitialized()) {
 218         std::shared_ptr<GeomAPI_Pnt> aCenter(aSketch->to3D(aCenterAttr->x(), aCenterAttr->y()));
 219
 220         std::shared_ptr<GeomDataAPI_Point2D> aStartAttr = std::dynamic_pointer_cast<
 221             GeomDataAPI_Point2D>(data()->attribute(SketchPlugin_Arc::START_ID()));
 222         std::shared_ptr<GeomDataAPI_Point2D> aEndAttr = std::dynamic_pointer_cast<
 223             GeomDataAPI_Point2D>(data()->attribute(SketchPlugin_Arc::END_ID()));
 224         AttributeStringPtr aTypeAttr = std::dynamic_pointer_cast<ModelAPI_AttributeString>(
 225             data()->attribute(ARC_TYPE()));
 226
 227         if (aStartAttr->isInitialized()) {
 228           // make a visible circle
 229           std::shared_ptr<GeomDataAPI_Dir> aNDir = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 230               aSketch->data()->attribute(SketchPlugin_Sketch::NORM_ID()));
 231           bool aHasPlane = aNDir && !(aNDir->x() == 0 && aNDir->y() == 0 && aNDir->z() == 0);
 232           if (aHasPlane) {
 233             std::shared_ptr<GeomAPI_Dir> aNormal = aNDir->dir();
 234             std::shared_ptr<GeomAPI_Pnt> aStartPoint(aSketch->to3D(aStartAttr->x(), aStartAttr->y()));
 235             std::shared_ptr<GeomAPI_Pnt> aEndPoint = aStartPoint;
 236             if (aTypeAttr && aTypeAttr->isInitialized() &&
 237                 aTypeAttr->value() == ARC_TYPE_THREE_POINTS() && aEndAttr->isInitialized() &&
 238                 aEndAttr->pnt()->distance(aStartAttr->pnt()) > tolerance) {
 239               aEndPoint = aSketch->to3D(aEndAttr->x(), aEndAttr->y());
 240               std::shared_ptr<GeomDataAPI_Point2D> aPassedAttr =
 241                 std::dynamic_pointer_cast<GeomDataAPI_Point2D>(data()->attribute(PASSED_POINT_ID()));
 242               if (!aPassedAttr->isInitialized()) { // calculate the appropriate center for the presentation
 243                 std::shared_ptr<GeomAPI_XY> aDir =
 244                   aEndAttr->pnt()->xy()->decreased(aStartAttr->pnt()->xy())->multiplied(0.5);
 245                 double x = aDir->x();
 246                 double y = aDir->y();
 247                 aDir->setX(x - y);
 248                 aDir->setY(y + x);
 249                 std::shared_ptr<GeomAPI_XY> aCenterXY = aStartAttr->pnt()->xy()->added(aDir);
 250                 aCenter = aSketch->to3D(aCenterXY->x(), aCenterXY->y());
 251               }
 252             }
 253
 254             std::shared_ptr<GeomAPI_Shape> aCircleShape = GeomAlgoAPI_EdgeBuilder::lineCircleArc(
 255                                                             aCenter, aStartPoint, aEndPoint, aNormal);
 256             if (aCircleShape)
 257               aShapes.push_back(aCircleShape);
 258           }
 259         }
 260         // make a visible point
 261         std::shared_ptr<GeomAPI_Shape> aCenterPointShape = GeomAlgoAPI_PointBuilder::point(aCenter);
 262         aShapes.push_back(aCenterPointShape);
 263       }
 264       if (!aShapes.empty()) {
 265         std::shared_ptr<GeomAPI_Shape> aCompound = GeomAlgoAPI_CompoundBuilder::compound(aShapes);
 266         AISObjectPtr anAIS = thePrevious;
 267         if (!anAIS)
 268           anAIS = AISObjectPtr(new GeomAPI_AISObject);
 269         anAIS->createShape(aCompound);
 270         anAIS->setWidth(3);
 271         return anAIS;
 272       }
 273     }
 274   }
 275   return AISObjectPtr();
 276 }
 277
 278 void SketchPlugin_Arc::move(double theDeltaX, double theDeltaY)
 279 {
 280   std::shared_ptr<ModelAPI_Data> aData = data();
 281   if (!aData->isValid())
 282     return;
 283
 284   aData->blockSendAttributeUpdated(true);
 285
 286   myStartUpdate = true;
 287   myEndUpdate = true;
 288   std::shared_ptr<GeomDataAPI_Point2D> aPoint2 = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 289       aData->attribute(SketchPlugin_Arc::START_ID()));
 290   aPoint2->move(theDeltaX, theDeltaY);
 291
 292   std::shared_ptr<GeomDataAPI_Point2D> aPoint3 = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 293       aData->attribute(SketchPlugin_Arc::END_ID()));
 294   aPoint3->move(theDeltaX, theDeltaY);
 295   myStartUpdate = false;
 296   myEndUpdate = false;
 297
 298   std::shared_ptr<GeomDataAPI_Point2D> aPoint1 = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 299       aData->attribute(SketchPlugin_Arc::CENTER_ID()));
 300   aPoint1->move(theDeltaX, theDeltaY);
 301
 302   std::shared_ptr<GeomDataAPI_Point2D> aPassedPoint =
 303       std::dynamic_pointer_cast<GeomDataAPI_Point2D>(aData->attribute(PASSED_POINT_ID()));
 304   aPassedPoint->move(theDeltaX, theDeltaY);
 305   aData->blockSendAttributeUpdated(false);
 306 }
 307
 308 bool SketchPlugin_Arc::isFixed() {
 309   return data()->selection(EXTERNAL_ID())->context().get() != NULL;
 310 }
 311
 312 bool SketchPlugin_Arc::isFeatureValid()
 313 {
 314   AttributeStringPtr anArcTypeAttr =
 315       std::dynamic_pointer_cast<ModelAPI_AttributeString>(data()->attribute(ARC_TYPE()));
 316   if (!anArcTypeAttr)
 317     return false;
 318   std::string anArcType = anArcTypeAttr->value();
 319
 320   std::shared_ptr<GeomDataAPI_Point2D> aCenterAttr = std::dynamic_pointer_cast<
 321       GeomDataAPI_Point2D>(data()->attribute(SketchPlugin_Arc::CENTER_ID()));
 322   std::shared_ptr<GeomDataAPI_Point2D> aStartAttr = std::dynamic_pointer_cast<
 323       GeomDataAPI_Point2D>(data()->attribute(SketchPlugin_Arc::START_ID()));
 324   std::shared_ptr<GeomDataAPI_Point2D> anEndAttr = std::dynamic_pointer_cast<
 325       GeomDataAPI_Point2D>(data()->attribute(SketchPlugin_Arc::END_ID()));
 326   std::shared_ptr<GeomDataAPI_Point2D> aPassedAttr = std::dynamic_pointer_cast<
 327       GeomDataAPI_Point2D>(data()->attribute(PASSED_POINT_ID()));
 328
 329   bool isValid = false;
 330   if (anArcType == ARC_TYPE_THREE_POINTS())
 331     isValid = aStartAttr->isInitialized() && anEndAttr->isInitialized() && aPassedAttr->isInitialized();
 332   else
 333     isValid = aCenterAttr->isInitialized() && aStartAttr->isInitialized() && anEndAttr->isInitialized();
 334
 335   return isValid;
 336 }
 337
 338 static inline void adjustPeriod(double& theParam)
 339 {
 340   static const double PERIOD = 2.0 * PI;
 341   while (theParam < 0.0) theParam += PERIOD;
 342   while (theParam >= PERIOD) theParam -= PERIOD;
 343 }
 344
 345 void SketchPlugin_Arc::attributeChanged(const std::string& theID)
 346 {
 347   std::shared_ptr<GeomDataAPI_Point2D> aCenterAttr = std::dynamic_pointer_cast<
 348       GeomDataAPI_Point2D>(data()->attribute(CENTER_ID()));
 349   std::shared_ptr<GeomDataAPI_Point2D> aStartAttr = std::dynamic_pointer_cast<
 350       GeomDataAPI_Point2D>(data()->attribute(START_ID()));
 351   std::shared_ptr<GeomDataAPI_Point2D> anEndAttr = std::dynamic_pointer_cast<
 352       GeomDataAPI_Point2D>(data()->attribute(END_ID()));
 353   // the second condition for unability to move external segments anywhere
 354   if (theID == EXTERNAL_ID() || isFixed()) {
 355     std::shared_ptr<GeomAPI_Shape> aSelection = data()->selection(EXTERNAL_ID())->value();
 356     // update arguments due to the selection value
 357     if (aSelection && !aSelection->isNull() && aSelection->isEdge()) {
 358       std::shared_ptr<GeomAPI_Edge> anEdge( new GeomAPI_Edge(aSelection));
 359       std::shared_ptr<GeomAPI_Circ> aCirc = anEdge->circle();
 360       if (aCirc.get()) {
 361         aStartAttr->setValue(sketch()->to2D(anEdge->firstPoint()));
 362         anEndAttr->setValue(sketch()->to2D(anEdge->lastPoint()));
 363         aCenterAttr->setValue(sketch()->to2D(aCirc->center()));
 364       }
 365     }
 366     return;
 367   }
 368   AttributeStringPtr aTypeAttr =
 369       std::dynamic_pointer_cast<ModelAPI_AttributeString>(attribute(ARC_TYPE()));
 370
 371   // this is before others since here end attribute may be changed, but with the special behavior
 372   if (aTypeAttr->value() == ARC_TYPE_TANGENT() && (theID == TANGENT_POINT_ID() || theID == END_ID())) {
 373     SketchPlugin_Sketch* aSketch = sketch();
 374     AttributeRefAttrPtr aTangPtAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 375         data()->attribute(TANGENT_POINT_ID()));
 376
 377     if (aTangPtAttr->isInitialized() && anEndAttr->isInitialized()) {
 378       // compute orthogonal direction
 379       std::shared_ptr<GeomAPI_Dir2d> anOrthoDir;
 380       std::shared_ptr<GeomDataAPI_Point2D> aTangentPoint =
 381           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(aTangPtAttr->attr());
 382       std::shared_ptr<GeomAPI_Pnt2d> aTangPnt2d = aTangentPoint->pnt();
 383       FeaturePtr aTangFeature = ModelAPI_Feature::feature(aTangentPoint->owner());
 384       std::shared_ptr<GeomAPI_Edge> aTangEdge = std::dynamic_pointer_cast<GeomAPI_Edge>(
 385           aTangFeature->lastResult()->shape());
 386       if (aTangEdge->isLine()) {
 387         std::shared_ptr<GeomAPI_Dir> aDir = aTangEdge->line()->direction();
 388         std::shared_ptr<GeomAPI_Pnt> aPnt(new GeomAPI_Pnt(aDir->x(), aDir->y(), aDir->z()));
 389         std::shared_ptr<GeomAPI_Pnt2d> aPnt2d = aSketch->to2D(aPnt);
 390         anOrthoDir = std::shared_ptr<GeomAPI_Dir2d>(new GeomAPI_Dir2d(-aPnt2d->y(), aPnt2d->x()));
 391       }
 392       else if (aTangEdge->isArc()) {
 393         std::shared_ptr<GeomAPI_Pnt> aCenter = aTangEdge->circle()->center();
 394         std::shared_ptr<GeomAPI_Pnt2d> aCenter2d = aSketch->to2D(aCenter);
 395         anOrthoDir = std::shared_ptr<GeomAPI_Dir2d>(
 396             new GeomAPI_Dir2d(aTangPnt2d->xy()->decreased(aCenter2d->xy())));
 397       }
 398
 399       // compute parameters of the middle perpendicular
 400       std::shared_ptr<GeomAPI_XY> aEndPntCoord = anEndAttr->pnt()->xy();
 401       std::shared_ptr<GeomAPI_XY> aTempDir = aEndPntCoord->decreased(aTangPnt2d->xy());
 402       std::shared_ptr<GeomAPI_Dir2d> aMidDir(new GeomAPI_Dir2d(-aTempDir->y(), aTempDir->x()));
 403       std::shared_ptr<GeomAPI_Pnt2d> aMidPnt(
 404           new GeomAPI_Pnt2d(aEndPntCoord->added(aTangPnt2d->xy())->multiplied(0.5)));
 405
 406       // compute center of arc by calculating intersection of orthogonal line and middle perpendicular
 407       std::shared_ptr<GeomAPI_Lin2d> anOrthoLine(new GeomAPI_Lin2d(aTangPnt2d, anOrthoDir));
 408       std::shared_ptr<GeomAPI_Lin2d> aMiddleLine(new GeomAPI_Lin2d(aMidPnt, aMidDir));
 409       std::shared_ptr<GeomAPI_Pnt2d> aCenter = anOrthoLine->intersect(aMiddleLine);
 410       if (aCenter) {
 411         data()->blockSendAttributeUpdated(true);
 412         aCenterAttr->setValue(aCenter);
 413         aStartAttr->setValue(aTangPnt2d);
 414         data()->blockSendAttributeUpdated(false);
 415       }
 416
 417       tangencyArcConstraints();
 418     }
 419     return;
 420   }
 421
 422   // if changed the base attributes, update all other (is necessary) without recursion
 423   if (theID == CENTER_ID() || theID == START_ID() || theID == END_ID() || theID == ARC_TYPE()) {
 424     if (!isFeatureValid())
 425       return;
 426     std::shared_ptr<GeomAPI_Pnt2d> aCenter = aCenterAttr->pnt();
 427     std::shared_ptr<GeomAPI_Pnt2d> aStart = aStartAttr->pnt();
 428     std::shared_ptr<GeomAPI_Pnt2d> anEnd = anEndAttr->pnt();
 429     double aRadius = aCenter->distance(aStart);
 430     if (aRadius < tolerance)
 431       return;
 432     std::shared_ptr<GeomAPI_Circ2d> aCircleForArc(new GeomAPI_Circ2d(aCenter, aStart));
 433
 434     data()->blockSendAttributeUpdated(true);
 435     if (theID == END_ID() && isStable()) {
 436       // The arc is under construction, so its end point projected
 437       // on the circle formed by center and start points
 438       std::shared_ptr<GeomAPI_Pnt2d> aProjection = aCircleForArc->project(anEnd);
 439       if (aProjection && anEnd->distance(aProjection) > tolerance) {
 440         anEndAttr->setValue(aProjection);
 441         anEnd = aProjection;
 442       }
 443     }
 444     // update all other attributes due to the base attributes values
 445     if (aTypeAttr->value() == ARC_TYPE_THREE_POINTS()) { // update passed point due to start, end and center
 446       if (aCenter->distance(aStart) > tolerance && aCenter->distance(anEnd) > tolerance) {
 447         // project passed point t othe circle
 448         std::shared_ptr<GeomDataAPI_Point2D> aPassedAttr =
 449           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(data()->attribute(PASSED_POINT_ID()));
 450         if (aPassedAttr->isInitialized()) {
 451           std::shared_ptr<GeomAPI_Pnt2d> aProjection = aCircleForArc->project(aPassedAttr->pnt());
 452           if (aProjection && aPassedAttr->pnt()->distance(aProjection) > tolerance) {
 453             aPassedAttr->setValue(aProjection);
 454           }
 455         } else { // initialize it by some middle - value
 456           std::shared_ptr<GeomAPI_Dir2d> aStartDir(new GeomAPI_Dir2d(
 457             aStart->xy()->decreased(aCenter->xy())));
 458           std::shared_ptr<GeomAPI_Dir2d> aEndDir(new GeomAPI_Dir2d(
 459             anEnd->xy()->decreased(aCenter->xy())));
 460           std::shared_ptr<GeomAPI_XY> aMidDirXY = aStartDir->xy()->added(aEndDir->xy());
 461           if (aMidDirXY->dot(aMidDirXY) < tolerance * tolerance) {
 462             // start and end directions are opposite, so middle direction will be orthogonal
 463             aMidDirXY->setX(-aStartDir->y());
 464             aMidDirXY->setY(aStartDir->x());
 465           }
 466           std::shared_ptr<GeomAPI_Dir2d> aMidDir(new GeomAPI_Dir2d(aMidDirXY));
 467           if ((aStartDir->cross(aMidDir) > 0) ^ !isReversed())
 468             aMidDir->reverse();
 469           std::shared_ptr<GeomAPI_XY> aPassedPnt =
 470             aCenter->xy()->added(aMidDir->xy()->multiplied(aCenter->distance(aStart)));
 471           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(data()->attribute(PASSED_POINT_ID()))->
 472             setValue(aPassedPnt->x(), aPassedPnt->y());
 473         }
 474       }
 475     }
 476     // update radius and angle
 477     AttributeDoublePtr aRadiusAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
 478         data()->attribute(RADIUS_ID()));
 479     aRadiusAttr->setValue(aRadius);
 480     AttributeDoublePtr anAngleAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
 481         data()->attribute(ANGLE_ID()));
 482     std::shared_ptr<GeomAPI_Circ2d> aCircle(new GeomAPI_Circ2d(aCenter, aStart));
 483     double aStartParam, aEndParam;
 484     aCircle->parameter(aStart, paramTolerance, aStartParam);
 485     aCircle->parameter(anEnd, paramTolerance, aEndParam);
 486     adjustPeriod(aStartParam);
 487     adjustPeriod(aEndParam);
 488     if (aTypeAttr->value() == ARC_TYPE_THREE_POINTS()) { // use the passed point for the angle calculation
 489       std::shared_ptr<GeomDataAPI_Point2D> aPassedAttr =
 490         std::dynamic_pointer_cast<GeomDataAPI_Point2D>(data()->attribute(PASSED_POINT_ID()));
 491       double aPassedParam;
 492       aCircle->parameter(aPassedAttr->pnt(), paramTolerance, aPassedParam);
 493       adjustPeriod(aPassedParam);
 494       double aNewAngle = aPassedParam >= aStartParam && aPassedParam <= aEndParam ?
 495         ((aEndParam - aStartParam) * 180.0 / PI) :
 496         ((aEndParam - aStartParam - 2.0 * PI) * 180.0 / PI);
 497       if (fabs(aNewAngle - anAngleAttr->value()) > tolerance)
 498         anAngleAttr->setValue(aNewAngle);
 499     } else {
 500       double aNewAngle = (aEndParam - aStartParam) * 180.0 / PI;
 501       if (fabs(aNewAngle - anAngleAttr->value()) > tolerance)
 502         anAngleAttr->setValue(aNewAngle);
 503     }
 504     // do not need to inform that other parameters were changed in this basis mode: these arguments
 505     // change is enough
 506     data()->blockSendAttributeUpdated(false, false);
 507     return;
 508   }
 509
 510   if (theID == PASSED_POINT_ID()) {
 511     data()->blockSendAttributeUpdated(true);
 512
 513     std::shared_ptr<GeomAPI_Pnt2d> aPoints[3];
 514     int aNbInitialized = 0;
 515     for (int i = 1; i <= 3; ++i) {
 516       std::shared_ptr<GeomDataAPI_Point2D> aCurPnt =
 517           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(attribute(POINT_ID(i)));
 518       if (aCurPnt->isInitialized())
 519         aPoints[aNbInitialized++] = aCurPnt->pnt();
 520     }
 521
 522     if (aNbInitialized == 3) {
 523       std::shared_ptr<GeomAPI_Circ2d> aCircle(
 524           new GeomAPI_Circ2d(aPoints[0], aPoints[1], aPoints[2]));
 525
 526       std::shared_ptr<GeomAPI_Pnt2d> aCenter = aCircle->center();
 527       if (aCenter) {
 528         aCenterAttr->setValue(aCenter);
 529       }
 530     }
 531     data()->blockSendAttributeUpdated(false);
 532     return;
 533   }
 534
 535   if (theID == RADIUS_ID()) {
 536     if (!aStartAttr->isInitialized() || !anEndAttr->isInitialized() || !aCenterAttr->isInitialized())
 537       return;
 538     // move center and passed point
 539     std::shared_ptr<GeomAPI_Pnt2d> aStart = aStartAttr->pnt();
 540     std::shared_ptr<GeomAPI_Pnt2d> anEnd = anEndAttr->pnt();
 541     std::shared_ptr<GeomAPI_Pnt2d> aCenter = aCenterAttr->pnt();
 542     if (aStart->distance(aCenter) < tolerance || anEnd->distance(aCenter) < tolerance)
 543       return;
 544     AttributeDoublePtr aRadiusAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
 545         data()->attribute(RADIUS_ID()));
 546     double aRadius = aRadiusAttr->value();
 547
 548     data()->blockSendAttributeUpdated(true);
 549     std::shared_ptr<GeomAPI_Dir2d> aStartDir(new GeomAPI_Dir2d(aStart->xy()->decreased(aCenter->xy())));
 550     std::shared_ptr<GeomAPI_XY> aNewStart = aStartDir->xy()->multiplied(aRadius)->added(aCenter->xy());
 551     aStartAttr->setValue(aNewStart->x(), aNewStart->y());
 552     std::shared_ptr<GeomAPI_Dir2d> anEndDir(new GeomAPI_Dir2d(anEnd->xy()->decreased(aCenter->xy())));
 553     std::shared_ptr<GeomAPI_XY> aNewEnd = anEndDir->xy()->multiplied(aRadius)->added(aCenter->xy());
 554     anEndAttr->setValue(aNewEnd->x(), aNewEnd->y());
 555     data()->blockSendAttributeUpdated(false);
 556     return;
 557   }
 558   if (theID == ANGLE_ID()) {
 559     if (!aStartAttr->isInitialized() || !aCenterAttr->isInitialized())
 560       return;
 561     AttributeDoublePtr anAngleAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
 562         data()->attribute(ANGLE_ID()));
 563     data()->blockSendAttributeUpdated(true);
 564     // move end point and passed point
 565     std::shared_ptr<GeomAPI_XY> aCenter = aCenterAttr->pnt()->xy();
 566     double anAngle = anAngleAttr->value() * PI / 180.0;
 567     double sinA = sin(anAngle);
 568     double cosA = cos(anAngle);
 569     std::shared_ptr<GeomAPI_XY> aStartDir = aStartAttr->pnt()->xy()->decreased(aCenter);
 570     std::shared_ptr<GeomAPI_XY> aDir(new GeomAPI_XY(
 571         aStartDir->x() * cosA - aStartDir->y() * sinA,
 572         aStartDir->x() * sinA + aStartDir->y() * cosA));
 573     anEndAttr->setValue(aCenter->x() + aDir->x(), aCenter->y() + aDir->y());
 574     data()->blockSendAttributeUpdated(false);
 575     return;
 576   }
 577 }
 578
 579 void SketchPlugin_Arc::setReversed(bool isReversed)
 580 {
 581   std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(attribute(INVERSED_ID()))->setValue(isReversed);
 582   myParamBefore = 0.0;
 583 }
 584
 585 bool SketchPlugin_Arc::isReversed()
 586 {
 587   return std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(attribute(INVERSED_ID()))->value();
 588 }
 589
 590 void SketchPlugin_Arc::tangencyArcConstraints()
 591 {
 592   if (!lastResult())
 593     return;
 594
 595   std::shared_ptr<GeomDataAPI_Point2D> aStartAttr =
 596       std::dynamic_pointer_cast<GeomDataAPI_Point2D>(attribute(START_ID()));
 597   AttributeRefAttrPtr aTangPtAttr = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 598       attribute(TANGENT_POINT_ID()));
 599   if (!aTangPtAttr->attr())
 600     return;
 601
 602   FeaturePtr aFeature = ModelAPI_Feature::feature(aStartAttr->owner());
 603   ObjectPtr aThisArc = aFeature->lastResult();
 604   aFeature = ModelAPI_Feature::feature(aTangPtAttr->attr()->owner());
 605   ObjectPtr aTangFeature = aFeature->lastResult();
 606
 607   // trying to find constraints to fix the tangency of the arc
 608   std::set<FeaturePtr> aCoincidence;
 609   std::set<FeaturePtr> aTangency;
 610
 611   AttributeRefAttrPtr aRefAttrA, aRefAttrB;
 612   std::set<AttributePtr> aRefs = data()->refsToMe();
 613   const std::set<AttributePtr>& aRefsToResult = lastResult()->data()->refsToMe();
 614   aRefs.insert(aRefsToResult.begin(), aRefsToResult.end());
 615   std::set<AttributePtr>::const_iterator aRefIt = aRefs.begin();
 616   for (; aRefIt != aRefs.end(); ++aRefIt) {
 617     FeaturePtr aConstrFeature = ModelAPI_Feature::feature((*aRefIt)->owner());
 618     if (aConstrFeature->getKind() == SketchPlugin_ConstraintCoincidence::ID()) {
 619       aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 620           aConstrFeature->attribute(SketchPlugin_Constraint::ENTITY_A()));
 621       aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 622           aConstrFeature->attribute(SketchPlugin_Constraint::ENTITY_B()));
 623       if ((aRefAttrA && aRefAttrA->attr() == aStartAttr) ||
 624           (aRefAttrB && aRefAttrB->attr() == aStartAttr))
 625         aCoincidence.insert(aConstrFeature);
 626     }
 627     else if (aConstrFeature->getKind() == SketchPlugin_ConstraintTangent::ID()) {
 628       aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 629           aConstrFeature->attribute(SketchPlugin_Constraint::ENTITY_A()));
 630       aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 631           aConstrFeature->attribute(SketchPlugin_Constraint::ENTITY_B()));
 632       if ((aRefAttrA && aRefAttrA->object() == aThisArc) ||
 633           (aRefAttrB && aRefAttrB->object() == aThisArc))
 634         aTangency.insert(aConstrFeature);
 635     }
 636   }
 637   // search applicable pair of constraints
 638   bool isFound = false;
 639   FeaturePtr aPrevCoincidence, aPrevTangency;
 640   std::set<FeaturePtr>::const_iterator aCIt, aTIt;
 641   for (aCIt = aCoincidence.begin(); aCIt != aCoincidence.end() && !isFound; ++aCIt) {
 642     aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 643         (*aCIt)->attribute(SketchPlugin_Constraint::ENTITY_A()));
 644     aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 645         (*aCIt)->attribute(SketchPlugin_Constraint::ENTITY_B()));
 646     AttributePtr anOtherPoint =
 647         aRefAttrA->attr() == aStartAttr ? aRefAttrB->attr() : aRefAttrA->attr();
 648     for (aTIt = aTangency.begin(); aTIt != aTangency.end() && !isFound; ++aTIt) {
 649       aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 650           (*aTIt)->attribute(SketchPlugin_Constraint::ENTITY_A()));
 651       aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 652           (*aTIt)->attribute(SketchPlugin_Constraint::ENTITY_B()));
 653       ObjectPtr anOtherObject = aRefAttrA->object() == aThisArc ?
 654           aRefAttrB->object() : aRefAttrA->object();
 655       FeaturePtr anOtherFeature = ModelAPI_Feature::feature(anOtherObject);
 656       if (anOtherPoint->owner() == anOtherFeature) {
 657         isFound = true;
 658         aPrevCoincidence = *aCIt;
 659         aPrevTangency = *aTIt;
 660       }
 661     }
 662   }
 663
 664   if (isFound) {
 665     // update previous constraints
 666     aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 667         aPrevCoincidence->attribute(SketchPlugin_Constraint::ENTITY_A()));
 668     aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 669         aPrevCoincidence->attribute(SketchPlugin_Constraint::ENTITY_B()));
 670     if (aRefAttrA->attr() == aStartAttr)
 671       aRefAttrB->setAttr(aTangPtAttr->attr());
 672     else
 673       aRefAttrA->setAttr(aTangPtAttr->attr());
 674
 675     aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 676         aPrevTangency->attribute(SketchPlugin_Constraint::ENTITY_A()));
 677     aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 678         aPrevTangency->attribute(SketchPlugin_Constraint::ENTITY_B()));
 679     if (aRefAttrA->object() == aThisArc)
 680       aRefAttrB->setObject(aTangFeature);
 681     else
 682       aRefAttrA->setObject(aTangFeature);
 683   } else {
 684     // Wait all constraints being removed, then send update events
 685     static Events_ID aDeleteEvent = Events_Loop::eventByName(EVENT_OBJECT_DELETED);
 686     bool isDeleteFlushed = Events_Loop::loop()->isFlushed(aDeleteEvent);
 687     if (isDeleteFlushed)
 688       Events_Loop::loop()->setFlushed(aDeleteEvent, false);
 689     // Remove all obtained constraints which use current arc, because
 690     // there is no information which of them were used to build tangency arc.
 691     DocumentPtr aDoc = sketch()->document();
 692     std::set<FeaturePtr> aFeaturesToBeRemoved;
 693     for (aCIt = aCoincidence.begin(); aCIt != aCoincidence.end(); ++aCIt)
 694       aFeaturesToBeRemoved.insert(*aCIt);
 695     for (aTIt = aTangency.begin(); aTIt != aTangency.end(); ++aTIt)
 696       aFeaturesToBeRemoved.insert(*aTIt);
 697     ModelAPI_Tools::removeFeaturesAndReferences(aFeaturesToBeRemoved);
 698     // Send events to update the sub-features by the solver.
 699     if (isDeleteFlushed)
 700       Events_Loop::loop()->setFlushed(aDeleteEvent, true);
 701     else
 702       Events_Loop::loop()->flush(aDeleteEvent);
 703
 704     // Wait all constraints being created, then send update events
 705     static Events_ID aCreateEvent = Events_Loop::eventByName(EVENT_OBJECT_CREATED);
 706     bool isCreateFlushed = Events_Loop::loop()->isFlushed(aCreateEvent);
 707     if (isCreateFlushed)
 708       Events_Loop::loop()->setFlushed(aCreateEvent, false);
 709
 710     // Create new constraints
 711     FeaturePtr aConstraint = sketch()->addFeature(SketchPlugin_ConstraintCoincidence::ID());
 712     aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 713         aConstraint->attribute(SketchPlugin_Constraint::ENTITY_A()));
 714     aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 715         aConstraint->attribute(SketchPlugin_Constraint::ENTITY_B()));
 716     aRefAttrA->setAttr(aStartAttr);
 717     aRefAttrB->setAttr(aTangPtAttr->attr());
 718     aConstraint->execute();
 719     ModelAPI_EventCreator::get()->sendUpdated(aConstraint, aCreateEvent);
 720
 721     aConstraint = sketch()->addFeature(SketchPlugin_ConstraintTangent::ID());
 722     aRefAttrA = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 723         aConstraint->attribute(SketchPlugin_Constraint::ENTITY_A()));
 724     aRefAttrB = std::dynamic_pointer_cast<ModelAPI_AttributeRefAttr>(
 725         aConstraint->attribute(SketchPlugin_Constraint::ENTITY_B()));
 726     aRefAttrA->setObject(aThisArc);
 727     aRefAttrB->setObject(aTangFeature);
 728     aConstraint->execute();
 729     ModelAPI_EventCreator::get()->sendUpdated(aConstraint, aCreateEvent);
 730
 731     // Send events to update the sub-features by the solver.
 732     if(isCreateFlushed)
 733       Events_Loop::loop()->setFlushed(aCreateEvent, true);
 734     else
 735       Events_Loop::loop()->flush(aCreateEvent);
 736   }
 737 }