Merge pull request #924 from Geode-solutions/feat/cut_along_solid_facets

panquez · web-flow · commit 4f8e54de42d7 · 2024-05-20T12:19:46.000+02:00
Feat/cut along solid facets
diff --git a/include/geode/mesh/core/detail/facet_storage.h b/include/geode/mesh/core/detail/facet_storage.h
@@ -115,7 +115,7 @@ namespace geode
                 for( const auto e :
                     Range{ facet_attribute_manager_.nb_elements() } )
                 {
-                    to_delete[e] = !counter_->value( e );
+                    to_delete[e] = counter_->value( e ) == 0;
                 }
                 return delete_facets( to_delete );
             }
diff --git a/include/geode/mesh/core/grid.h b/include/geode/mesh/core/grid.h
@@ -170,4 +170,27 @@ namespace geode
         IMPLEMENTATION_MEMBER( impl_ );
     };
     ALIAS_2D_AND_3D( Grid );
+
+    namespace detail
+    {
+        static constexpr std::array< std::array< geode::local_index_t, 2 >, 12 >
+            CELL_EDGE_VERTICES_3D{ { { 0, 1 }, { 0, 2 }, { 2, 3 }, { 1, 3 },
+                { 4, 5 }, { 4, 6 }, { 6, 7 }, { 5, 7 }, { 0, 4 }, { 1, 5 },
+                { 2, 6 }, { 3, 7 } } };
+
+        static constexpr std::array< std::array< geode::local_index_t, 4 >, 6 >
+            CELL_FACET_VERTICES_3D{ { { 0, 4, 6, 2 }, { 1, 5, 7, 3 },
+                { 0, 4, 5, 1 }, { 2, 3, 7, 6 }, { 0, 1, 3, 2 },
+                { 4, 6, 7, 5 } } };
+
+        static constexpr std::array< std::array< geode::local_index_t, 3 >, 8 >
+            CELL_EDGES_AROUND_VERTEX_3D{ { { 0, 1, 8 }, { 0, 3, 9 },
+                { 1, 2, 10 }, { 2, 3, 11 }, { 4, 5, 8 }, { 4, 7, 9 },
+                { 5, 6, 10 }, { 6, 7, 11 } } };
+
+        static constexpr std::array< std::array< geode::local_index_t, 3 >, 8 >
+            CELL_FACETS_AROUND_VERTEX_3D{ { { 0, 2, 4 }, { 1, 2, 4 },
+                { 0, 3, 4 }, { 1, 3, 4 }, { 0, 2, 5 }, { 1, 2, 5 }, { 0, 3, 5 },
+                { 1, 3, 5 } } };
+    } // namespace detail
 } // namespace geode
diff --git a/include/geode/mesh/core/meshes_mapping.h b/include/geode/mesh/core/meshes_mapping.h
@@ -0,0 +1,41 @@
+/*
+ * Copyright (c) 2019 - 2024 Geode-solutions
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a copy
+ * of this software and associated documentation files (the "Software"), to deal
+ * in the Software without restriction, including without limitation the rights
+ * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+ * copies of the Software, and to permit persons to whom the Software is
+ * furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ *
+ */
+
+#pragma once
+
+#include <geode/basic/mapping.h>
+
+#include <geode/mesh/common.h>
+
+namespace geode
+{
+    using ElementsMapping = GenericMapping< index_t, index_t >;
+
+    struct MeshesElementsMapping
+    {
+        ElementsMapping vertices;
+        ElementsMapping edges;
+        ElementsMapping polygons;
+        ElementsMapping polyhedra;
+    };
+} // namespace geode
diff --git a/include/geode/mesh/helpers/detail/cut_along_solid_facets.h b/include/geode/mesh/helpers/detail/cut_along_solid_facets.h
@@ -0,0 +1,63 @@
+/*
+ * Copyright (c) 2019 - 2024 Geode-solutions
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a copy
+ * of this software and associated documentation files (the "Software"), to deal
+ * in the Software without restriction, including without limitation the rights
+ * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+ * copies of the Software, and to permit persons to whom the Software is
+ * furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ *
+ */
+
+#pragma once
+
+#include <absl/types/span.h>
+
+#include <geode/basic/pimpl.h>
+
+#include <geode/mesh/common.h>
+#include <geode/mesh/core/meshes_mapping.h>
+
+namespace geode
+{
+    FORWARD_DECLARATION_DIMENSION_CLASS( SolidMesh );
+    FORWARD_DECLARATION_DIMENSION_CLASS( SolidMeshBuilder );
+    ALIAS_3D( SolidMesh );
+    ALIAS_3D( SolidMeshBuilder );
+} // namespace geode
+
+namespace geode
+{
+    namespace detail
+    {
+        class opengeode_mesh_api CutAlongSolidFacets
+        {
+        public:
+            CutAlongSolidFacets(
+                const SolidMesh3D& mesh, SolidMeshBuilder3D& builder );
+            ~CutAlongSolidFacets();
+
+            /*
+             * Cuts the solid along given facets, and returns the mapping on
+             * vertices and facets.
+             */
+            MeshesElementsMapping cut_solid_along_facets(
+                absl::Span< const PolyhedronFacet > facets_list );
+
+        private:
+            IMPLEMENTATION_MEMBER( impl_ );
+        };
+    } // namespace detail
+} // namespace geode
diff --git a/include/geode/model/helpers/detail/cut_along_internal_surfaces.h b/include/geode/model/helpers/detail/cut_along_internal_surfaces.h
@@ -43,7 +43,6 @@ namespace geode
         class opengeode_model_api CutAlongInternalSurfaces
         {
         public:
-            CutAlongInternalSurfaces( BRep& model );
             CutAlongInternalSurfaces( const BRep& model, BRepBuilder& builder );
             ~CutAlongInternalSurfaces();
 
diff --git a/src/geode/mesh/CMakeLists.txt b/src/geode/mesh/CMakeLists.txt
@@ -117,6 +117,7 @@ add_geode_library(
         "helpers/detail/component_identifier.cpp"
         "helpers/detail/create_mesh.cpp"
         "helpers/detail/curve_merger.cpp"
+        "helpers/detail/cut_along_solid_facets.cpp"
         "helpers/detail/debug.cpp"
         "helpers/detail/element_identifier.cpp"
         "helpers/detail/solid_merger.cpp"
@@ -298,6 +299,7 @@ add_geode_library(
         "core/detail/geode_elements.h"
         "core/detail/vertex_cycle.h"
         "helpers/detail/curve_merger.h"
+        "helpers/detail/cut_along_solid_facets.h"
         "helpers/detail/component_identifier.h"
         "helpers/detail/create_mesh.h"
         "helpers/detail/debug.h"
diff --git a/src/geode/mesh/builder/solid_mesh_builder.cpp b/src/geode/mesh/builder/solid_mesh_builder.cpp
@@ -293,49 +293,49 @@ namespace
             }
             const auto position_it =
                 absl::c_find( facet_vertices, polyhedron_vertex );
-            if( position_it != facet_vertices.end() )
+            if( position_it == facet_vertices.end() )
             {
-                geode::PolyhedronFacetVertices facet_vertices_id;
-                facet_vertices_id.reserve( nb_facet_vertices );
-                for( const auto& v : facet_vertices )
-                {
-                    facet_vertices_id.emplace_back(
-                        solid.polyhedron_vertex( v ) );
-                }
+                continue;
+            }
+            geode::PolyhedronFacetVertices facet_vertices_id;
+            facet_vertices_id.reserve( nb_facet_vertices );
+            for( const auto& v : facet_vertices )
+            {
+                facet_vertices_id.emplace_back( solid.polyhedron_vertex( v ) );
+            }
 
-                const auto position = static_cast< geode::index_t >(
-                    std::distance( facet_vertices.begin(), position_it ) );
+            const auto position = static_cast< geode::index_t >(
+                std::distance( facet_vertices.begin(), position_it ) );
 
-                if( solid.are_facets_enabled() )
+            if( solid.are_facets_enabled() )
+            {
+                auto facets = builder.facets_builder();
+                facets.update_facet_vertex(
+                    facet_vertices_id, position, new_vertex_id );
+            }
+
+            if( solid.are_edges_enabled() )
+            {
+                auto edges = builder.edges_builder();
+                const auto next =
+                    position + 1 == nb_facet_vertices ? 0 : position + 1;
+                std::array< geode::index_t, 2 > next_edge_vertices{
+                    facet_vertices_id[position], facet_vertices_id[next]
+                };
+                if( next_edge_vertices[0] < next_edge_vertices[1] )
                 {
-                    auto facets = builder.facets_builder();
-                    facets.update_facet_vertex(
-                        facet_vertices_id, position, new_vertex_id );
+                    edges.update_edge_vertex(
+                        next_edge_vertices, 0, new_vertex_id );
                 }
-
-                if( solid.are_edges_enabled() )
+                const auto prev =
+                    position == 0 ? nb_facet_vertices - 1 : position - 1;
+                std::array< geode::index_t, 2 > previous_edge_vertices{
+                    facet_vertices_id[prev], facet_vertices_id[position]
+                };
+                if( previous_edge_vertices[0] < previous_edge_vertices[1] )
                 {
-                    auto edges = builder.edges_builder();
-                    const auto next =
-                        position + 1 == nb_facet_vertices ? 0 : position + 1;
-                    std::array< geode::index_t, 2 > next_edge_vertices{
-                        facet_vertices_id[position], facet_vertices_id[next]
-                    };
-                    if( next_edge_vertices[0] < next_edge_vertices[1] )
-                    {
-                        edges.update_edge_vertex(
-                            next_edge_vertices, 0, new_vertex_id );
-                    }
-                    const auto prev =
-                        position == 0 ? nb_facet_vertices - 1 : position - 1;
-                    std::array< geode::index_t, 2 > previous_edge_vertices{
-                        facet_vertices_id[prev], facet_vertices_id[position]
-                    };
-                    if( previous_edge_vertices[0] < previous_edge_vertices[1] )
-                    {
-                        edges.update_edge_vertex(
-                            previous_edge_vertices, 1, new_vertex_id );
-                    }
+                    edges.update_edge_vertex(
+                        previous_edge_vertices, 1, new_vertex_id );
                 }
             }
         }
diff --git a/src/geode/mesh/core/solid_mesh.cpp b/src/geode/mesh/core/solid_mesh.cpp
@@ -1146,6 +1146,7 @@ namespace geode
         }
         return result;
     }
+
     template < index_t dimension >
     absl::optional< index_t > SolidMesh< dimension >::polyhedron_around_edge(
         const std::array< index_t, 2 >& vertices ) const
diff --git a/src/geode/mesh/helpers/detail/cut_along_solid_facets.cpp b/src/geode/mesh/helpers/detail/cut_along_solid_facets.cpp
diff --git a/src/geode/model/helpers/detail/cut_along_internal_surfaces.cpp b/src/geode/model/helpers/detail/cut_along_internal_surfaces.cpp
diff --git a/tests/mesh/test-convert-solid.cpp b/tests/mesh/test-convert-solid.cpp

Original file line number	Diff line number	Diff line change
`@@ -115,7 +115,7 @@ namespace geode`
`115`	`115`	`for( const auto e :`
`116`	`116`	`Range{ facet_attribute_manager_.nb_elements() } )`
`117`	`117`	`{`
`118`		`- to_delete[e] = !counter_->value( e );`
	`118`	`+ to_delete[e] = counter_->value( e ) == 0;`
`119`	`119`	`}`
`120`	`120`	`return delete_facets( to_delete );`
`121`	`121`	`}`
Original file line number	Diff line number	Diff line change
`@@ -43,7 +43,6 @@ namespace geode`
`43`	`43`	`class opengeode_model_api CutAlongInternalSurfaces`
`44`	`44`	`{`
`45`	`45`	`public:`
`46`		`- CutAlongInternalSurfaces( BRep& model );`
`47`	`46`	`CutAlongInternalSurfaces( const BRep& model, BRepBuilder& builder );`
`48`	`47`	`~CutAlongInternalSurfaces();`
`49`	`48`
Original file line number	Diff line number	Diff line change
`@@ -1146,6 +1146,7 @@ namespace geode`
`1146`	`1146`	`}`
`1147`	`1147`	`return result;`
`1148`	`1148`	`}`
	`1149`	`+`
`1149`	`1150`	`template < index_t dimension >`
`1150`	`1151`	`absl::optional< index_t > SolidMesh< dimension >::polyhedron_around_edge(`
`1151`	`1152`	`const std::array< index_t, 2 >& vertices ) const`