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