From 862c385c7548fb69d1674a409b3ccaeae13ddf1c Mon Sep 17 00:00:00 2001
From: azv <azv@opencascade.com>
Date: Thu, 9 Jun 2016 17:15:36 +0300
Subject: [PATCH] PlaneGCS: fix incorrect range of fillet arc on the XZ plane

---
 .../PlaneGCSSolver/PlaneGCSSolver_Storage.cpp         | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/src/SketchSolver/PlaneGCSSolver/PlaneGCSSolver_Storage.cpp b/src/SketchSolver/PlaneGCSSolver/PlaneGCSSolver_Storage.cpp
index a3f832e6e..66c47d162 100644
--- a/src/SketchSolver/PlaneGCSSolver/PlaneGCSSolver_Storage.cpp
+++ b/src/SketchSolver/PlaneGCSSolver/PlaneGCSSolver_Storage.cpp
@@ -418,6 +418,17 @@ void PlaneGCSSolver_Storage::processArc(const EntityWrapperPtr& theArc)
   if (!anArcEdge)
     return;
   anArcEdge->getRange(*aStartAngle, *aEndAngle);
+  // verify the range is correct and not shifted to an angle
+  std::shared_ptr<GeomAPI_Dir2d> aDir(new GeomAPI_Dir2d(cos(*aStartAngle), sin(*aStartAngle)));
+  std::shared_ptr<GeomAPI_Pnt2d> aCalcStartPnt(
+      new GeomAPI_Pnt2d(aCenterPnt->xy()->added(aDir->xy()->multiplied(*aRadius))));
+  if (aCalcStartPnt->distance(aStartPnt) > tolerance) {
+    std::shared_ptr<GeomAPI_Dir2d> aDirToStart(
+        new GeomAPI_Dir2d(aStartPnt->xy()->decreased(aCenterPnt->xy())));
+    double anAngle = aDir->angle(aDirToStart);
+    *aStartAngle += anAngle;
+    *aEndAngle += anAngle;
+  }
 
   // do not constraint copied arc
   if (anArcFeature->isCopy())
-- 
2.39.2