feat(AnisotropicDistanceToEdge): adding a new DistanceToEdge object in an anisotropic space.

Author: benoit2776

include/geode/mesh/helpers/aabb_edged_curve_helpers.hpp

@@ -36,6 +36,8 @@ namespace geode
     FORWARD_DECLARATION_DIMENSION_CLASS( Point );
     FORWARD_DECLARATION_DIMENSION_CLASS( AABBTree );
     FORWARD_DECLARATION_DIMENSION_CLASS( EdgedCurve );
+    FORWARD_DECLARATION_DIMENSION_CLASS( CoordinateSystem );
+    ALIAS_2D_AND_3D( CoordinateSystem );
     template < typename T >
     class GenericMeshAABB;
 } // namespace geode
@@ -63,6 +65,25 @@ namespace geode
     };
     ALIAS_1D_AND_2D_AND_3D( DistanceToEdge );
 
+    template < index_t dimension >
+    class AnisotropicDistanceToEdge
+    {
+    public:
+        explicit AnisotropicDistanceToEdge( const EdgedCurve< dimension >& mesh,
+            const CoordinateSystem< dimension >& coordinate_system )
+            : mesh_( mesh ), coordinate_system_( coordinate_system )
+        {
+        }
+
+        [[nodiscard]] double operator()(
+            const Point< dimension >& query, index_t cur_box ) const;
+
+    private:
+        const EdgedCurve< dimension >& mesh_;
+        const CoordinateSystem< dimension >& coordinate_system_;
+    };
+    ALIAS_2D_AND_3D( AnisotropicDistanceToEdge );
+
     template < index_t dimension >
     class GenericMeshAABB< EdgedCurve< dimension > >
     {

src/geode/mesh/helpers/aabb_edged_curve_helpers.cpp

@@ -33,6 +33,7 @@
 
 #include <geode/geometry/aabb.hpp>
 #include <geode/geometry/basic_objects/segment.hpp>
+#include <geode/geometry/coordinate_system.hpp>
 #include <geode/geometry/distance.hpp>
 #include <geode/geometry/point.hpp>
 #include <geode/geometry/projection.hpp>
@@ -65,6 +66,20 @@ namespace geode
         return point_segment_distance( query, segment );
     }
 
+    template < index_t dimension >
+    double AnisotropicDistanceToEdge< dimension >::operator()(
+        const Point< dimension >& query, index_t cur_box ) const
+    {
+        const auto segment = mesh_.segment( cur_box );
+        const auto vertices = segment.vertices();
+        const auto point0 = coordinate_system_.coordinates( vertices[0].get() );
+        const auto point1 = coordinate_system_.coordinates( vertices[1].get() );
+        const auto query_in_metric_space =
+            coordinate_system_.coordinates( query );
+        geode::Segment< dimension > segment_in_metric_space{ point0, point1 };
+        return point_segment_distance( query_in_metric_space, segment );
+    }
+
     template opengeode_mesh_api AABBTree1D create_aabb_tree< 1 >(
         const EdgedCurve1D& );
     template opengeode_mesh_api AABBTree2D create_aabb_tree< 2 >(
@@ -75,4 +90,7 @@ namespace geode
     template class opengeode_mesh_api DistanceToEdge< 1 >;
     template class opengeode_mesh_api DistanceToEdge< 2 >;
     template class opengeode_mesh_api DistanceToEdge< 3 >;
+
+    template class opengeode_mesh_api AnisotropicDistanceToEdge< 2 >;
+    template class opengeode_mesh_api AnisotropicDistanceToEdge< 3 >;
 } // namespace geode