Fix coordinate position on vertical boundary

Inspired by @wwaltersdavis in #1465

Author: michaelkirk

geo/CHANGES.md

@@ -1,5 +1,9 @@
 # Changes
 
+## Unreleased
+
+- Fix `CoordinatePosition` for `Triangle` to correctly return `CoordPos::OnBoundary` for coordinate within vertical segment.
+
 ## 0.32.0 - 2025-12-05
 
 - Move `PreparedGeometry` into a new `indexed` module intended to provide index-backed geometries. `relate::PreparedGeometry` has been deprecated.

geo/src/algorithm/coordinate_position.rs

@@ -192,7 +192,12 @@ where
                 let orientation = T::Ker::orient2d(l.start, l.end, *coord);
                 if orientation == Orientation::Collinear
                     && point_in_rect(*coord, l.start, l.end)
-                    && coord.x != l.end.x
+                    // If a coordinate falls *on* one of the triangle's vertices,
+                    // it will appear both as the `end` of one segment and the `start` of the next.
+                    // We only want to count one boundary crossing in that case,
+                    // so we ignore it if it falls on `l.end`,
+                    // knowing that we'll count it on some other segments `l.start`
+                    && *coord != l.end
                 {
                     *boundary_count += 1;
                 }
@@ -445,10 +450,8 @@ where
 
 #[cfg(test)]
 mod test {
-    use geo_types::coord;
-
     use super::*;
-    use crate::{line_string, point, polygon};
+    use crate::{coord, line_string, point, polygon, wkt};
 
     #[test]
     fn test_empty_poly() {
@@ -708,12 +711,29 @@ mod test {
             triangle.coordinate_position(&coord! { x: 2.5, y: 5.0 }),
             CoordPos::OnBoundary
         );
+        assert_eq!(
+            triangle.coordinate_position(&coord! { x: 5.0, y: 10.0 }),
+            CoordPos::OnBoundary
+        );
         assert_eq!(
             triangle.coordinate_position(&coord! { x: 2.49, y: 5.0 }),
             CoordPos::Outside
         );
     }
 
+    #[test]
+    fn test_right_triangle() {
+        let triangle: Triangle = wkt!(TRIANGLE(0. 0.,10. 0.,10. 10.));
+        assert_eq!(
+            triangle.coordinate_position(&coord!(x: 0., y: 0.)),
+            CoordPos::OnBoundary
+        );
+        assert_eq!(
+            triangle.coordinate_position(&coord!(x: 10., y: 5.)),
+            CoordPos::OnBoundary
+        );
+    }
+
     #[test]
     fn test_collection() {
         let triangle = Triangle::new((0.0, 0.0).into(), (5.0, 10.0).into(), (10.0, 0.0).into());