feat(AnisotropicDistanceToSurface): adding DistanceToSurface for anistropic cases.

Author: benoit2776

include/geode/mesh/helpers/aabb_surface_helpers.hpp

@@ -37,6 +37,8 @@ namespace geode
     FORWARD_DECLARATION_DIMENSION_CLASS( AABBTree );
     FORWARD_DECLARATION_DIMENSION_CLASS( SurfaceMesh );
     FORWARD_DECLARATION_DIMENSION_CLASS( TriangulatedSurface );
+    FORWARD_DECLARATION_DIMENSION_CLASS( CoordinateSystem );
+    ALIAS_2D_AND_3D( CoordinateSystem );
     template < typename T >
     class GenericMeshAABB;
 } // namespace geode
@@ -65,6 +67,26 @@ namespace geode
     };
     ALIAS_2D_AND_3D( DistanceToTriangle );
 
+    template < index_t dimension >
+    class AnisotropicDistanceToTriangle
+    {
+    public:
+        explicit AnisotropicDistanceToTriangle(
+            const TriangulatedSurface< 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 TriangulatedSurface< dimension >& mesh_;
+        const CoordinateSystem< dimension >& coordinate_system_;
+    };
+    ALIAS_2D_AND_3D( AnisotropicDistanceToTriangle );
+
     template < index_t dimension >
     class GenericMeshAABB< SurfaceMesh< dimension > >
     {

src/geode/mesh/helpers/aabb_surface_helpers.cpp

@@ -33,6 +33,7 @@
 
 #include <geode/geometry/aabb.hpp>
 #include <geode/geometry/basic_objects/triangle.hpp>
+#include <geode/geometry/coordinate_system.hpp>
 #include <geode/geometry/distance.hpp>
 #include <geode/geometry/point.hpp>
 
@@ -67,12 +68,31 @@ namespace geode
             point_triangle_distance( query, mesh_.triangle( cur_box ) ) );
     }
 
+    template < index_t dimension >
+    double AnisotropicDistanceToTriangle< dimension >::operator()(
+        const Point< dimension >& query, index_t cur_box ) const
+    {
+        const auto triangle = mesh_.triangle( cur_box );
+        const auto vertices = triangle.vertices();
+        const auto point0 = coordinate_system_.coordinates( vertices[0].get() );
+        const auto point1 = coordinate_system_.coordinates( vertices[1].get() );
+        const auto point2 = coordinate_system_.coordinates( vertices[2].get() );
+        const Triangle< dimension > triangle_in_metric_space{ point0, point1,
+            point2 };
+        const auto query_in_metric_space =
+            coordinate_system_.coordinates( query );
+        return std::get< 0 >( point_triangle_distance(
+            query_in_metric_space, triangle_in_metric_space ) );
+    }
+
     template opengeode_mesh_api AABBTree2D create_aabb_tree< 2 >(
         const SurfaceMesh2D& );
     template opengeode_mesh_api AABBTree3D create_aabb_tree< 3 >(
         const SurfaceMesh3D& );
 
     template class opengeode_mesh_api DistanceToTriangle< 2 >;
     template class opengeode_mesh_api DistanceToTriangle< 3 >;
+    template class opengeode_mesh_api AnisotropicDistanceToTriangle< 2 >;
+    template class opengeode_mesh_api AnisotropicDistanceToTriangle< 3 >;
 
 } // namespace geode