working on algos

src/polymesh/Mesh.hh

@@ -1120,24 +1120,48 @@ vertex_property<PropT> vertex_collection::make_property(const PropT &def_value)
     return vertex_property<PropT>(mesh, def_value);
 }
 
+template <typename PropT>
+vertex_property<PropT> const_vertex_collection::make_property(const PropT &def_value)
+{
+    return vertex_property<PropT>(mesh, def_value);
+}
+
 template <typename PropT>
 face_property<PropT> face_collection::make_property(const PropT &def_value)
 {
     return face_property<PropT>(mesh, def_value);
 }
 
+template <typename PropT>
+face_property<PropT> const_face_collection::make_property(const PropT &def_value)
+{
+    return face_property<PropT>(mesh, def_value);
+}
+
 template <typename PropT>
 edge_property<PropT> edge_collection::make_property(const PropT &def_value)
 {
     return edge_property<PropT>(mesh, def_value);
 }
 
+template <typename PropT>
+edge_property<PropT> const_edge_collection::make_property(const PropT &def_value)
+{
+    return edge_property<PropT>(mesh, def_value);
+}
+
 template <typename PropT>
 halfedge_property<PropT> halfedge_collection::make_property(const PropT &def_value)
 {
     return halfedge_property<PropT>(mesh, def_value);
 }
 
+template <typename PropT>
+halfedge_property<PropT> const_halfedge_collection::make_property(const PropT &def_value)
+{
+    return halfedge_property<PropT>(mesh, def_value);
+}
+
 inline void Mesh::register_prop(vertex_property_base *prop) const
 {
     // insert in front

src/polymesh/algorithms/geodesic_fast_marching.hh

+#pragma once
+
+#include <glm/glm.hpp>
+
+#include <polymesh/Mesh.hh>
+
+namespace polymesh
+{
+vertex_property<float> geodesic_fast_marching_linear(Mesh const& m, vertex_property<glm::vec3> const& position)
+{
+    auto dis = m.vertices().make_property<float>(-1);
+
+    // TODO
+
+    return dis;
+}
+}