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