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