Feat: Add check_coplanar kwarg to Line.intersect_line

src/skspatial/objects/line.py

@@ -475,7 +475,7 @@ def distance_line(self, other: Line) -> np.float64:
 
         return distance
 
-    def intersect_line(self, other: Line, **kwargs) -> Point:
+    def intersect_line(self, other: Line, tol: np.float64 = None, **kwargs) -> Point:
         """
         Intersect the line with another.
 
@@ -485,6 +485,8 @@ def intersect_line(self, other: Line, **kwargs) -> Point:
         ----------
         other : Line
             Other line.
+        tol : np.float64, optional
+            Threshold below which values are considered zero. This value is passed to is_coplanar
         kwargs : dict, optional
             Additional keywords passed to :meth:`Vector.is_parallel`.
 
@@ -560,7 +562,7 @@ def intersect_line(self, other: Line, **kwargs) -> Point:
         if self.direction.is_parallel(other.direction, **kwargs):
             raise ValueError("The lines must not be parallel.")
 
-        if not self.is_coplanar(other):
+        if not self.is_coplanar(other, tol=tol):
             raise ValueError("The lines must be coplanar.")
 
         # Vector from line A to line B.

src/skspatial/objects/line_segment.py

@@ -106,14 +106,17 @@ def contains_point(self, point: array_like, **kwargs) -> bool:
 
         return math.isclose(similarity, -1, **kwargs)
 
-    def intersect_line_segment(self, other: LineSegment) -> Point:
+    def intersect_line_segment(self, other: LineSegment, tol: np.float64 = None) -> Point:
         """
         Intersect the line segment with another.
 
         Parameters
         ----------
         other : LineSegment
 
+        tol : np.float64, optional
+            Threshold below which values are considered zero. This value is passed to is_coplanar as tol and to is_collinear as abs_tol.
+
         Returns
         -------
         Point
@@ -145,7 +148,7 @@ def intersect_line_segment(self, other: LineSegment) -> Point:
         line_a = Line.from_points(self.point_a, self.point_b)
         line_b = Line.from_points(other.point_a, other.point_b)
 
-        point_intersection = line_a.intersect_line(line_b)
+        point_intersection = line_a.intersect_line(line_b, tol=tol, abs_tol=tol) if tol is not None else line_a.intersect_line(line_b)
 
         point_on_segment_a = self.contains_point(point_intersection)
         point_on_segment_b = other.contains_point(point_intersection)

src/skspatial/objects/triangle.py

@@ -498,12 +498,18 @@ def altitude(self, vertex: str) -> Line:
 
         raise ValueError("The vertex must be 'A', 'B', or 'C'.")
 
-    def orthocenter(self) -> Point:
+    def orthocenter(self, tol: np.float64 = None) -> Point:
         """
         Return the orthocenter of the triangle.
 
         The orthocenter is the intersection point of the three altitudes.
 
+        Parameters
+        ----------
+        tol : np.float64, optional
+            Threshold below which values are considered zero. This value is passed to is_coplanar as tol and to is_collinear as abs_tol.
+
+
         Returns
         -------
         Point
@@ -523,7 +529,7 @@ def orthocenter(self) -> Point:
         line_alt_a = self.altitude('A')
         line_alt_b = self.altitude('B')
 
-        return line_alt_a.intersect_line(line_alt_b)
+        return line_alt_a.intersect_line(line_alt_b, tol=tol, abs_tol=tol) if tol is not None else line_alt_a.intersect_line(line_alt_b)
 
     def classify(self, **kwargs: float) -> str:
         """