From 9cc7b1cc2b55eed25553fa8a8c277e17613147d0 Mon Sep 17 00:00:00 2001 From: artemyev Date: Fri, 4 Jul 2025 10:17:52 +0700 Subject: [PATCH] feat(edge_flaps): EMAP now has sizes of F --- mesh/edge_flaps.m | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/mesh/edge_flaps.m b/mesh/edge_flaps.m index 4184f9f1..9c217584 100644 --- a/mesh/edge_flaps.m +++ b/mesh/edge_flaps.m @@ -9,12 +9,13 @@ % e=(j->i) % EI #E by 2 list of edge flap corners (see above). % uE #uE by 2 list of edge indices into V. - % EMAP #F*3 list of indices into uE, mapping each directed edge to unique + % EMAP #F by 3 list of indices into uE, mapping each directed edge to unique % unique edge in uE % - E = [F(:,2:3);F(:,[3 1]);,F(:,1:2)]; + E = [F(:,2:3);F(:,[3 1]);F(:,1:2)]; sE = sort(E,2); [uE,~,EMAP] = unique(sE,'rows'); + EMAP = reshape(EMAP,size(F)); I = (1:size(uE,1))'; [B1,E1] = ismember(uE,E,'rows'); [B2,E2] = ismember(uE,fliplr(E),'rows');