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');