#3278 fixed

BHoM · Feb 12, 2024 · 19144dc · 19144dc
1 parent 4e5ef3c
commit 19144dc
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/Geometry_Engine/Compute/FitLine.cs b/Geometry_Engine/Compute/FitLine.cs
@@ -64,8 +64,16 @@ public static Line FitLine(this IEnumerable<Point> points, double tolerance = To
                 double sumXY = points.Sum(x => x.X * x.Y);
                 double p1 = sumX * sumX - sumY * sumY - n * (sumXsq - sumYsq);
                 double p2 = n * sumXY - sumX * sumY;
-                double b = (p1 + Math.Sqrt(p1 * p1 + 4 * p2 * p2)) / (2 * p2);
-                Vector dir = new Vector { X = 1, Y = b };
+
+                Vector dir;
+                if (Math.Abs(p2) > tolerance)
+                {
+                    double b = (p1 + Math.Sqrt(p1 * p1 + 4 * p2 * p2)) / (2 * p2);
+                    dir = new Vector { X = 1, Y = b };
+                }
+                else
+                    dir = p1 < 0 ? Vector.XAxis : Vector.YAxis;
+
                 return new Line { Start = C, End = C + dir };
             }
             else