BuildBE.m

Go to the documentation of this file.

function [berefT,belrefT]=BuildBE(e2q,e2be,nq)
% function [berefT,belrefT]=BuildBE(e2q,e2be,nq)
%   Build boundaries edges with labels (non sorted)
  I=find(e2be(:,2)==1);
  be=e2q(e2be(:,1),:);
  be(I,:)=fliplr(be(I,:));
  bel=e2be(:,3);
  midpoint=e2be(:,1)+nq;
  beref=[ be(:,1) midpoint; midpoint be(:,2)];
  belref=[bel;bel];
  belrefT=belref;
  berefT=beref;
end