update

Author: BotellaA

src/geode/geometry/position.cpp

@@ -30,6 +30,7 @@
 #include <geode/geometry/basic_objects/segment.hpp>
 #include <geode/geometry/basic_objects/tetrahedron.hpp>
 #include <geode/geometry/basic_objects/triangle.hpp>
+#include <geode/geometry/bounding_box.hpp>
 #include <geode/geometry/distance.hpp>
 #include <geode/geometry/internal/position_from_sides.hpp>
 #include <geode/geometry/internal/predicates.hpp>
@@ -67,6 +68,25 @@ namespace geode
             internal::side( dot0 ), internal::opposite_side( dot1 ) );
     }
 
+    POSITION point_segment_position_exact(
+        const Point1D& point, const Segment1D& segment )
+    {
+        const auto& vertices = segment.vertices();
+        if( point == vertices[0] )
+        {
+            return POSITION::vertex0;
+        }
+        if( point == vertices[1] )
+        {
+            return POSITION::vertex1;
+        }
+        if( segment.bounding_box().contains( point ) )
+        {
+            return POSITION::inside;
+        }
+        return POSITION::outside;
+    }
+
     template < index_t dimension >
     POSITION point_triangle_position_all_zero(
         const Point< dimension >& point, const Triangle< dimension >& triangle )
@@ -305,9 +325,10 @@ namespace geode
         return GEO::PCK::aligned_3d( point0, point1, point2 );
     }
 
+    template POSITION opengeode_geometry_api point_segment_position(
+        const Point1D&, const Segment1D& );
     template POSITION opengeode_geometry_api point_segment_position(
         const Point2D&, const Segment2D& );
-
     template POSITION opengeode_geometry_api point_segment_position(
         const Point3D&, const Segment3D& );
 } // namespace geode