1 // Copyright (C) 2014-20xx CEA/DEN, EDF R&D -->
3 // File: SketchPlugin_Arc.cpp
4 // Created: 26 Apr 2014
5 // Author: Artem ZHIDKOV
7 #include "SketchPlugin_Arc.h"
8 #include "SketchPlugin_Sketch.h"
9 #include <SketchPlugin_ConstraintCoincidence.h>
10 #include <SketchPlugin_ConstraintTangent.h>
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>
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>
42 const double tolerance = 1e-7;
43 const double paramTolerance = 1.e-4;
44 const double PI = 3.141592653589793238463;
47 static const std::string& POINT_ID(int 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();
55 static const std::string DUMMY;
62 SketchPlugin_Arc::SketchPlugin_Arc()
63 : SketchPlugin_SketchEntity()
65 myStartUpdate = false;
71 myParamBefore = PI * 2.0;
74 void SketchPlugin_Arc::initDerivedClassAttributes()
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());
83 AttributeBooleanPtr isInversed = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
84 data()->addAttribute(INVERSED_ID(), ModelAPI_AttributeBoolean::typeId()));
86 // get the initial values
87 if (anEndAttr->isInitialized()) {
88 myXEndBefore = anEndAttr->x();
89 myYEndBefore = anEndAttr->y();
92 AttributeStringPtr anArcType = std::dynamic_pointer_cast<ModelAPI_AttributeString>(
93 data()->addAttribute(ARC_TYPE(), ModelAPI_AttributeString::typeId()));
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());
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());
105 data()->addAttribute(CENTER_REF_ID(), ModelAPI_AttributeRefAttr::typeId());
106 ModelAPI_Session::get()->validators()->registerNotObligatory(getKind(), CENTER_REF_ID());
109 void SketchPlugin_Arc::execute()
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;
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()));
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(
129 aConstr1->setShape(aCenterPointShape);
130 aConstr1->setIsInHistory(false);
131 setResult(aConstr1, 0);
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()));
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()));
151 // compute end parameter
152 aCircleForArc->parameter(anEndAttr->pnt(), paramTolerance, myParamBefore);
154 std::shared_ptr<GeomAPI_Shape> aCircleShape;
155 if(!isInversed->value()) {
157 GeomAlgoAPI_EdgeBuilder::lineCircleArc(aCenter, aStartPoint, aEndPoint, aNormal);
160 GeomAlgoAPI_EdgeBuilder::lineCircleArc(aCenter, aEndPoint, aStartPoint, aNormal);
164 std::shared_ptr<ModelAPI_ResultConstruction> aConstr2 = document()->createConstruction(
166 aConstr2->setShape(aCircleShape);
167 aConstr2->setIsInHistory(false);
168 setResult(aConstr2, 1);
171 if (!hasResult && string(ARC_TYPE())->value() == ARC_TYPE_TANGENT()) {
172 // constraints for tangent arc have not been prepared yet, do them
173 tangencyArcConstraints();
178 AISObjectPtr SketchPlugin_Arc::getAISObject(AISObjectPtr thePrevious)
180 SketchPlugin_Sketch* aSketch = sketch();
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()));
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()));
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()));
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);
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();
231 std::shared_ptr<GeomAPI_XY> aCenterXY = aStartAttr->pnt()->xy()->added(aDir);
232 aCenter = aSketch->to3D(aCenterXY->x(), aCenterXY->y());
235 } else { // issue #1695: don't display circle if initialized only start point
236 return AISObjectPtr();
239 AttributeBooleanPtr isInversed =
240 std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(attribute(INVERSED_ID()));
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);
248 aShapes.push_back(aCircleShape);
251 // make a visible point
252 std::shared_ptr<GeomAPI_Shape> aCenterPointShape =
253 GeomAlgoAPI_PointBuilder::vertex(aCenter);
254 aShapes.push_back(aCenterPointShape);
256 if (!aShapes.empty()) {
257 std::shared_ptr<GeomAPI_Shape> aCompound = GeomAlgoAPI_CompoundBuilder::compound(aShapes);
258 AISObjectPtr anAIS = thePrevious;
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);
270 return AISObjectPtr();
273 void SketchPlugin_Arc::move(double theDeltaX, double theDeltaY)
275 std::shared_ptr<ModelAPI_Data> aData = data();
276 if (!aData->isValid())
279 bool aWasBlocked = aData->blockSendAttributeUpdated(true);
281 myStartUpdate = 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);
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;
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);
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);
307 bool SketchPlugin_Arc::isFixed() {
308 return data()->selection(EXTERNAL_ID())->context().get() != NULL;
311 bool SketchPlugin_Arc::isFeatureValid()
313 AttributeStringPtr anArcTypeAttr =
314 std::dynamic_pointer_cast<ModelAPI_AttributeString>(data()->attribute(ARC_TYPE()));
317 std::string anArcType = anArcTypeAttr->value();
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()));
328 bool isValid = false;
329 if (anArcType == ARC_TYPE_THREE_POINTS())
330 isValid = aStartAttr->isInitialized() &&
331 anEndAttr->isInitialized() && aPassedAttr->isInitialized();
333 isValid = aCenterAttr->isInitialized() &&
334 aStartAttr->isInitialized() && anEndAttr->isInitialized();
339 static inline void adjustPeriod(double& theParam)
341 static const double PERIOD = 2.0 * PI;
342 while (theParam < 0.0) theParam += PERIOD;
343 while (theParam >= PERIOD) theParam -= PERIOD;
346 void SketchPlugin_Arc::attributeChanged(const std::string& theID)
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();
358 // empty shape in selection shows that the shape is equal to context
359 ResultPtr anExtRes = selection(EXTERNAL_ID())->context();
361 aSelection = anExtRes->shape();
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();
368 aStartAttr->setValue(sketch()->to2D(anEdge->firstPoint()));
369 anEndAttr->setValue(sketch()->to2D(anEdge->lastPoint()));
370 aCenterAttr->setValue(sketch()->to2D(aCirc->center()));
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);
382 AttributeStringPtr aTypeAttr =
383 std::dynamic_pointer_cast<ModelAPI_AttributeString>(attribute(ARC_TYPE()));
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()));
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()))
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()));
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())));
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)));
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);
429 bool aWasBlocked = data()->blockSendAttributeUpdated(true);
430 aCenterAttr->setValue(aCenter);
431 aStartAttr->setValue(aTangPnt2d);
432 data()->blockSendAttributeUpdated(aWasBlocked);
435 tangencyArcConstraints();
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())
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)
450 std::shared_ptr<GeomAPI_Circ2d> aCircleForArc(new GeomAPI_Circ2d(aCenter, aStart));
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);
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);
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());
483 std::shared_ptr<GeomAPI_Dir2d> aMidDir(new GeomAPI_Dir2d(aMidDirXY));
484 if ((aStartDir->cross(aMidDir) > 0) ^ !isReversed())
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());
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()));
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);
518 double aNewAngle = (aEndParam - aStartParam) * 180.0 / PI;
519 if (!anAngleAttr->isInitialized() || fabs(aNewAngle - anAngleAttr->value()) > tolerance)
520 anAngleAttr->setValue(aNewAngle);
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;
535 // do not need to inform that other parameters were changed in this basis mode: these arguments
537 data()->blockSendAttributeUpdated(aWasBlocked, false);
541 if (theID == PASSED_POINT_ID()) {
542 bool aWasBlocked = data()->blockSendAttributeUpdated(true);
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();
553 if (aNbInitialized == 3) {
554 std::shared_ptr<GeomAPI_Circ2d> aCircle(
555 new GeomAPI_Circ2d(aPoints[0], aPoints[1], aPoints[2]));
557 std::shared_ptr<GeomAPI_Pnt2d> aCenter = aCircle->center();
559 aCenterAttr->setValue(aCenter);
562 data()->blockSendAttributeUpdated(aWasBlocked);
566 if (theID == RADIUS_ID()) {
567 if (!aStartAttr->isInitialized() || !anEndAttr->isInitialized() ||
568 !aCenterAttr->isInitialized())
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)
576 AttributeDoublePtr aRadiusAttr = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
577 data()->attribute(RADIUS_ID()));
578 double aRadius = aRadiusAttr->value();
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);
594 if (theID == ANGLE_ID()) {
595 if (!aStartAttr->isInitialized() || !aCenterAttr->isInitialized())
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);
615 void SketchPlugin_Arc::setReversed(bool isReversed)
617 std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
618 attribute(INVERSED_ID()))->setValue(isReversed);
622 bool SketchPlugin_Arc::isReversed()
624 return std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(attribute(INVERSED_ID()))->value();
627 void SketchPlugin_Arc::tangencyArcConstraints()
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())
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();
644 // trying to find constraints to fix the tangency of the arc
645 std::set<FeaturePtr> aCoincidence;
646 std::set<FeaturePtr> aTangency;
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);
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);
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) {
695 aPrevCoincidence = *aCIt;
696 aPrevTangency = *aTIt;
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());
710 aRefAttrA->setAttr(aTangPtAttr->attr());
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);
719 aRefAttrA->setObject(aTangFeature);
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);
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.
737 Events_Loop::loop()->setFlushed(aDeleteEvent, true);
739 Events_Loop::loop()->flush(aDeleteEvent);
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);
745 Events_Loop::loop()->setFlushed(aCreateEvent, false);
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);
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);
768 // Send events to update the sub-features by the solver.
770 Events_Loop::loop()->setFlushed(aCreateEvent, true);
772 Events_Loop::loop()->flush(aCreateEvent);