replace f by ff!

and use facet_edges() helper function
This commit is contained in:
Jane Tournois 2024-05-06 17:50:23 +02:00
parent f8ff5f6f1c
commit 30a66089b8
1 changed files with 3 additions and 11 deletions

View File

@ -365,19 +365,11 @@ private:
const Facet ff = facets.front(); const Facet ff = facets.front();
facets.pop_front(); facets.pop_front();
const typename C3t3::Cell_handle ch = f.first; const Vector_3& ref = fnormals[ff];
const std::array<std::array<int, 2>, 3> edges for (const Edge& ei : facet_edges(ff.first, ff.second, tr))
= {{ {{(ff.second + 1) % 4, (ff.second + 2) % 4}}, //edge 1-2
{{(ff.second + 2) % 4, (ff.second + 3) % 4}}, //edge 2-3
{{(ff.second + 3) % 4, (ff.second + 1) % 4}} //edge 3-1
}}; //vertex indices in cells
const Vector_3& ref = fnormals[f];
for (const std::array<int, 2>& ei : edges)
{ {
Edge edge(ch, ei[0], ei[1]);
if (std::optional<Facet> neighbor if (std::optional<Facet> neighbor
= find_adjacent_facet_on_surface(f, edge, c3t3)) = find_adjacent_facet_on_surface(ff, ei, c3t3))
{ {
const Facet neigh = *neighbor; //already a canonical_facet const Facet neigh = *neighbor; //already a canonical_facet
if (fnormals[neigh] == CGAL::NULL_VECTOR) //check it's not already computed if (fnormals[neigh] == CGAL::NULL_VECTOR) //check it's not already computed